grads批量合并有规律的二进制文件（如NC逐日资料等）和批量提取特定时间的特定要素

白兔糖0217

非常有用，多谢楼主！

白兔糖0217

我也参考着用，请问楼主，
'reinit'
yy=2014
mm=07
dd=01
while(dd<32)
'open D:/Fortran/ARWpost/test/'yy''mm''dd'.ctl'   
为什么老是在这open一行出错？

river

白兔糖0217 发表于 2014-10-23 15:26
我也参考着用，请问楼主，
'reinit'
yy=2014

yy和mm要循环吗，不循环的话直接写在文件明上就行了。还有就是缺乏grads基础，单引号双引号不是乱用的
'open D:/Fortran/ARWpost/test/201407'dd'.ctl'   

白兔糖0217

river 发表于 2014-10-23 17:02
yy和mm要循环吗，不循环的话直接写在文件明上就行了。还有就是缺乏grads基础，单引号双引号不是乱用的
' ...

这样试过，不行唉

白兔糖0217

白兔糖0217 发表于 2014-10-23 18:17
这样试过，不行唉

logina:/global/research/iumhouyy/ARWpost> vim mywrf.gs
      1 'reinit'
      2 dd=10
      3 while(d<12)
      4 'open /global/research/iumhouyy/ARWpost/201406'dd'.ctl'
      5 'set gxout fwrite'
      6 *'set fwrite /global/research/iumhouyy/grads/201406'dd'.grd'
      7 'set lat 32.0'
      8 'set lon 118.8'
      9 it=1
     10 while(it<5)
     11 'set t 'it
     12 'q time'
     13 m=substr(result,13,3)
     14 d=substr(result,11,2)
     15 h=substr(result,8,2)
     16 'set fwrite /global/research/iumhouyy/grads/201406'%dd%''h'.grd'
     17 i=1
     18 while(i<48)
     19 'set z 'i
     20 'd pressure'
     21 'd tk'
     22 'd tc'
     23 'd td'
     24 'd rh'
     25 i=i+1
     26 endwhile
     27 'disable fwrite'
     28 ‘reinit’
     29 ;
这是我修改后的。无错误，无.grd生成

river

白兔糖0217 发表于 2014-10-23 18:26
logina:/global/research/iumhouyy/ARWpost> vim mywrf.gs
      1 'reinit'
      2 dd=10

你确定没错？三个while循环，只有一个endwhile，而且还有两个根本没有增量，你说没错怎么可能，是grads崩溃了吧······

白兔糖0217

river 发表于 2014-10-23 21:32
你确定没错？三个while循环，只有一个endwhile，而且还有两个根本没有增量，你说没错怎么可能，是grads崩 ...

多谢指正！已经知道所有错误了。谢谢

天之边

我就喜欢这种把程序摆出来的好楼主

__明月__

总觉得这次用grads很奇怪。。第一，reinit之后，yy=1979就执行不了，不用说后面的了。。囧
第二，我不能同时打开多个necp的nc额。。。。。要去投奔Fortran了。。囧

budapestw

感谢楼主 学习一下