grads处理多个ctl文件和nc文件.docx

《grads处理多个ctl文件和nc文件.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、grads处理多个ctl文件和nc文件  2011-10-1021:03:59

2、  分类：grads学习

3、  标签：

4、举报

5、字号大中小 订阅       下载LOFTER我的照片书  

6、用grads处理多个相同格式的数据时若单个单个处理非常麻烦，当文件非常多的时候是单个处理是不实际的。下面介绍一种方法;  第一步，在这种情况下可以重新写一个ctl描述文件，其文件变量都和已知的ctl相同，若原来的n文件只是时间不同，那么新描述文件的时间维数是所有原文件的时间的和。同样，若其他维数不同时也用同样的方法处理。   第二

7、步，在第一行之后添加一行：optionstemplate表示多个时间序列原始数据文件想用一个描述文件统一地描述。这些原数据的原文件名由dset定义的形势命名文件名。  第三步，修改dset的文件名。原路径不变，把文件名用％表示。其中：  ％y2 代表两位数年  ％y4 代表四位数年  ％m1 代表一位或者两位数的月  ％m2 代表两位数月（用0补齐1位数）  ％mc 3个字符月份的缩写  ％d1 1或2位天  ％d2 两位天  ％h1 1或者2位时  ％h2 2位时例如：  原文件其中之一的文件名为gdas200

8、6050812f00，且所有文件只有天和时的变化  那么新描述文件的文件名为：gdas200605%d2%h2f00     另外如果源文件里有index项的话，需要修改其idx的文件名，假设改成fnl.idx。并用在dos下用gribmap函数生成一个新的idx文件。gribmap-e-ifnl.ctl（加绝对路径）openfnl.ctl就可以打开所有文件。 *****************************************************************************

9、*****************************************************若想要提取从1951-2006年56年nc文件中的某些数据，一个一个处理非常麻烦，这里介绍种较为简易的方法。例如想提取6-8月的位势高度资料。'reinit't5=1951*作文件名循环while(t5<=2006) 'setgxoutfwrite' 'setfwriteD:sichuanhgt1'%t5%'.dat' 'sdfopene:cep1hgthgt.'%t5%'.nc' t3=t5-

10、1950*判断是否为闰年 if（t3=2

11、t3=6

12、t3=10

13、t3=14

14、t3=18

15、t3=22

16、t3=26

17、t3=30

18、t3=34

19、t3=38

20、t3=42

21、t3=46

22、t3=50

23、t3=54) to=153 else to=152 endif   t4=to+91  while(to<=t4)    'sett'to     t1=1     while(t1<=12)       'setz't1       'setlon80140'       'setlat1555'       'dhgt'    

24、    t1=t1+1        endwhile     to=to+1  endwhile*这里必须先观点上述运行的文件，grads最多同时可以打开20个文件左右。  'reinit' t5=t5+1endwhile'reinit'这样可以提取你想要的年数据，然后你大可运用fortran对数据进行随心所欲的处理。 能否直接生成一个文件还正在探索中。bbs.06climate./forum.php?mod=viewthread&tid=7310&extra=&page=1批量读取nc数据，用你的方法成功了，

25、！！！直接配个批量描述的ctl就可以了有一批nc数据，一个月一个文件，现将文件名改为：197901.nc，197902.nc，依次类推，对二进制的数据知道写ctl文件来进行批处理运算，那么nc数据应该怎么做呢？试过了写ctl文件，sdfopen***%y4%m2.nc，year=1978while(year<=2011)month=01  while(month<=12)'sdfopen***'year''month'.nc'...  month=month+1  endwhileyear=year+1end

26、while实我也是糊里糊涂的解决了。。。ctl文件如下：dset  ^%y4%m2.ncundef1e+15optionstemplatetitleMERRAdatadtypeNetCDFydef  144linear-901.25xdef  288linear-1801.25zdef  21levels1000975950925900875850825800775750