grads处理多个ctl文件和nc文件解析Word格式.docx

1、 两位天 h1 1或者2位时 h2 2位时例如： 原文件其中之一的文件名为gdas2006050812f00，且所有文件只有天和时的变化 那么新描述文件的文件名为：gdas200605%d2%h2f00另外如果源文件里有index项的话，需要修改其idx的文件名，假设改成fnl.idx。并用在dos下用gribmap函数生成一个新的idx文件。gribmap -e -i fnl.ctl（加绝对路径）open fnl.ctl就可以打开所有文件。*若想要提取从1951-2006年56年nc文件中的某些数据，一个一个处理非常麻烦，这里介绍种较为简易的方法。例如想提取6-8月的位势高度资料。reini

2、tt5=1951*作文件名循环while（t5=2006）set gxout fwriteset fwrite D:sichuanhgt1%t5%.datsdfopen e:ncep1hgthgt.nct3=t5-1950*判断是否为闰年if（t3=2|t3=6|t3=10|t3=14|t3=18|t3=22|t3=26|t3=30|t3=34|t3=38|t3=42|t3=46|t3=50|t3=54） to=153else to=152endif t4=to+91 while（to=t4） set t to t1=1 while（t1=12）set z t1set lon 80 140se

3、t lat 15 55d hgt t1=t1+1 endwhile to=to+1*这里必须先观点上述运行的文件，grads最多同时可以打开20个文件左右。t5=t5+1endwhile这样可以提取你想要的年数据，然后你大可运用fortran对数据进行随心所欲的处理。能否直接生成一个文件还正在探索中。批量读取nc数据，用你的方法成功了，谢谢！直接配个批量描述的ctl就可以了有一批nc数据，一个月一个文件，现将文件名改为：197901.nc，197902.nc，依次类推，对二进制的数据知道写ctl文件来进行批处理运算，那么nc数据应该怎么做呢？试过了写ctl文件，sdfopen *%y4%m2.

4、nc，year=1978while（year=2011）month=01while（monthhgt 17t,z,y,xmean Daily Geopotential height10.endvars复制代码光看着正常不行啊，需要和原始的图对比验证了才能确定是对的吧。所以在grads里面用sdfopen命令打开hgt.1948.nc画第一层第一个时次的图，再用xdfopen打开你编写的ctl，也画第一层第一个时次的图看看。我擦！对不上啊，看来还是有问题啊。说明上面的ctl是有问题的，还得改进才行。其实到这步已经基本成功了，仔细看看叠加起来的那两张图，几乎是一样的，只是南北的方向是反的。那就好办

5、了，在ctl里加上options yrev，告诉grads南北要反向不就行了么。于是最终的ctl出来了，如下3.options yrev4.* yrev表示y轴反向5.undef -9996.xdef lon 144 linear 0 2.5 7.ydef lat 73 linear -90 2.5 8.zdef level 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 109.tdef time 366 linear 00Z01Jan1948 1440mn10.vars 1 11.hgt=12.e

6、ndvars再用xdfopen打开画图，这回就一模一样了啊。 说明成功了！那下面的批量描述就太简单了，比如我要批量描述1948和1949两年的，算一下一个闰年一个平年，一共有时次366+365=731，那么就修改ctl吧dset F:/ncep/daily/hgt.%y4.nctitle mean daily NMC Reanalysis options yrev* yrev表示y轴反向*undef -999xdef lon 144 linear 0 2.5 ydef lat 73 linear -90 2.5 zdef level 17 levels 1000 925 850 700 600

7、 500 400 300 250 200 150 100 70 50 30 20 10tdef time 731 linear 00Z01Jan1948 1440mnvars 1 hgt=就这么简单，只要用%y4代表四位年，把总的时次改成731，再增加一句options template就可以了。然后就可以利用xdfopen命令打开画图和原来的对比一下，一样一样的吧，可以批量描述了吧！注意：有人可能注意到我把缺测那一项给注释掉了，其实这个是完全不需要的，去掉或者改成任何值，都不影响。这是问什么？不细说了，因为我没仔细研究（比较晚了，我要去睡觉了，以后再说，或者有人知道可以告诉我）其他的资料都同

8、理了，思路都差不多，变通一下就可以了。1.首先编写可以正确描述一个资料的ctl，并能正确出图。2.修改该ctl里dest后的文件名，使用合适的替代格式替代（具体的格式看实用手册）3.增加options template，这个是批量描述必须加上的一句4.计算出所要批量描述的文件的总时次，修改ctl里的总时次5.画图进行验证，如果不对，再根据具体情况做出相应的修改6.大功告成，可以用了7.需注意的一点就是options yrev，这是可选项，要根据资料实际情况来使用。今天这个帖子有点儿罗嗦了，大家捡有用的看就行。其实说那么多，无非是想和一些人说我们这个年代的人已经是站在前人的肩膀上了，很多东西别人

9、都已经有了经验，分享给你你就要好好利用，不要把时间都浪费在发帖和等回复上面。很多东西，只要自己有一定的基础，在一定合理的范围内就会找到那个答案的GrADS grads批量处理nc文件求区域平均 复制链接 |关注本帖取消最新回复取消置顶回复取消最新编辑stefan_scofield stefan_scofield 当前离线 积分4181贡献精华在线时间小时注册时间2013-7-18最后登录1970-1-1窥视卡 雷达卡 电梯直达 楼主 海温资料是1965-2013年的逐月资料（每年的4,5,6月），我需要求每年4,5,6月nino3.4区域的平均海温，然后三个月再平均，就是每年得一个值。win7

10、 批量处理nc文件，合成一个dat文件时遇到问题，请大家帮指点一下！set fwrite g:mls1.20.32.datt5=1=20）sdfopen g:mlsh2o.t5set x 1 144set y 1 91set t 1z=1while（z=25）zd h2o.t5 z=z+1disable fwritegrads提取某个范围内的资料时注意的一个问题热度 4已有 524 次阅读2011-11-19 16:10 | 资料 这两天被一个资料的提取程序弄的晕乎乎的，今天终于发现问题所在了。提取资料的时候，经度的范围即使原来就是全球的，也不能写成 set lon 0 360，得写成set

11、lon 0 357.5，就这个2.5（格点间距）把我折腾了这么久，下次一定要注意了！报错：Data request warning: request is completely outside file limitsgs如下：set fwrite d:/sst/4-6nino.grdyy=1965while（yypcp 0 t,z,y,x precipitationENDVARSgs文件如下：open d:/trmmdata/work/2000.ctlset lat 0 60set lon 0 30set lev 1000d ave（pcp,t=1,t=12）用GrADS转换nc数据by Ar

12、tmunich from kittyhare 在一篇帖子里向我们初步介绍了使用GrADS转化nc为grd格式（摸我去看看），我想能不能再详细点，于是借鉴以前自己见到的一个教程，更详细的介绍一下这方面的知识，希望大家指正。先给出一个单层的二进制文件转化的gssdfopen f:/data/nc/1980/air.1980.ncset fwrite f:/data/air.1980.binset lon 0 357.5set lat -90 90set t 1 640d air 但我们知道nc文件是按照经度、纬度、高度、变量、时次顺序排列，要转nc文件，需要靠循环，由于高度的不连续性，我们可以在时

13、间循环里面把每层的高度写出来。t=1while（t=365）set t tset lev 925*You can continue to write hgtt=t+1当然，如果用z坐标系，这样z就是从1到17的，在时间的循环里嵌套z=1;=17）我认为也是可行的，不过自己没试过，有用过的朋友可以告诉一下结果。高度循环的：=8）t=21） zset lon 90 125set lat 10 35define t=tmpprsdefine rh=rhprsdefine es=（6.112*exp（17.67*（t-273.15）/（t-29.65）define q=rh*（0.62197*es/（lev-es）/100.define e=lev*q/（0.62197+q）+1e-10define tlcl=55.0+2840.0/（3.5*log（t）-log（e）-4.805）define theta=t*pow（1000./lev）,（0.2854*（1.0-0.28*q）define thetse=theta*exp（3376./tlcl）-2.54）*q*（1.0+0.81*q）d thetsez=z+1