关于oacres函数下define变量名的问题

linda霖霖

我想要实现两个或多个时次的降水量分布画在同一张图上，第一个时次的可以顺利出图。但从第二个时次，问题就出现了define这里了，define的变量名只能是rgrid吗？oacres括号内的，g.3、rain.4 应该没有问题吧？求大神们解答

river

你明白后面的点几 点几 是什么意思吗？不知道的话请先去了解清楚

linda霖霖

river 发表于 2018-4-27 21:40
你明白后面的点几 点几 是什么意思吗？不知道的话请先去了解清楚

应该是所在open的数据文件吧，第14行rain.2就是表示打开的第二个文件20.ctl里的rain。
然后后面的g.3就是23grid.ctl里面的g，rain.4就是23.ctl里面的rain。

它没有error提示，说明gs应该没问题。提示是少于2个站点没法插值，有点疑惑，第一部分（20grid.ctl，20.ctl）是能顺利出图的，但第二部分（23grid.ctl，23.ctl）就不能，尝试过把第二部分单独拿出来，是能够正常出图的。另外，尝试过在每个文件打开下一行分别加set dfile 1，...，set dfile 4语句，则Define error:  Define not yet valid for station data.Default file is a station data file。

我在想会不会是因为第一部分和第二部分的站点分布不一致（站点资料本身该时次有降水则有显示，没降水则没有这个站点），第二部分的降水站点分布与第一部分不一样，所以插值结果也就不一样，所以导致了第二部分不能正常出图？

zhi726

似乎应该是维度匹配的问题，你查一查两者的时间维度，还有空间定义是不是与前面的一致？

linda霖霖

zhi726 发表于 2018-4-28 12:55
似乎应该是维度匹配的问题，你查一查两者的时间维度，还有空间定义是不是与前面的一致？

时间维度应该没问题，都是该时次的1h降水量。空间定义范围也一样。分别单独做20时、23时的降水量都没有问题，但想要同时出在一张图就不行，只出了20时的，23时的则如上述提示那样。不知道原因在哪咯

zhi726

linda霖霖 发表于 2018-4-28 13:03
时间维度应该没问题，都是该时次的1h降水量。空间定义范围也一样。分别单独做20时、23时的降水量都没有问 ...

我是说第三个文件和第4个文件里的时间是否一致，就是g.3和grain.4的时间是不是一致？

linda霖霖

zhi726 发表于 2018-4-29 04:07
我是说第三个文件和第4个文件里的时间是否一致，就是g.3和grain.4的时间是不是一致？

是一致的，把第三、第四个文件单独拿出来作图是没有问题的。

问题已经解决了，在第一部分后面加了close 2和close 1语句，就顺利出图了。