【精细版】用GrADS转换nc数据

artmunich

                                         用GrADS转换nc数据
                                                      ——by Artmunich from bbs.06climate.com
    kittyhare 在一篇帖子里向我们初步介绍了使用GrADS转化nc为grd格式（摸我去看看~），我想能不能再详细点，于是借鉴以前自己见到的一个教程，更详细的介绍一下这方面的知识，希望大家指正。
   先给出一个单层的二进制文件转化的gs

'reinit'
'sdfopen f:/data/nc/1980/air.1980.nc'
'set gxout fwrite'
'set fwrite f:/data//air.1980.bin'
'set lon 0 357.5'
'set lat -90 90'
'set lev 1000'
'set t 1 640'
'd air'
'reinit'

    但我们知道nc文件是按照经度、纬度、高度、变量、时次顺序排列，要转nc文件，需要靠循环，由于高度的不连续性，我们可以在时间循环里面把每层的高度写出来。
'reinit'
'sdfopen f:/data/nc/1980/air.1980.nc'
'set gxout fwrite'
'set fwrite f:/data//air.1980.bin'
'set lon 0 357.5'
'set lat -90 90'
t=1
while(t<=365)
'set t' t
'set lev 1000'
'd air'
'set lev 925'
'd air'
*You can continue to write hgt
t=t+1
endwhile
'reinit'
   当然，如果用z坐标系，这样z就是从1到17的，在时间的循环里嵌套z=1;while(z<=17)我认为也是可行的，不过自己没试过，有用过的朋友可以告诉一下结果。

尽头的尽头

辛苦了啊，不过手上没有NC资料，先收下来

传说中的谁

把高度循环的加上来吧，我懒得贴了

zsmnantong

顶一下，楼主辛苦，谢谢分享

artmunich

@传说中的谁 基情奉献
鼓掌

while(t<=8)
'set t 't''
z=1
while(z<=21)
'set z ' z''
'set lon 90 125'
'set lat 10 35'
'define t=tmpprs'
'define rh=rhprs'
'define 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-10'
'define 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 thetse'
z=z+1
endwhile
t=t+1
endwhile

传说中的谁

artmunich 发表于 2011-8-22 20:14

                               
登录/注册后可看大图

@传说中的谁 基情奉献
鼓掌

错别字啊有木有

artmunich

传说中的谁 发表于 2011-8-22 20:17

                               
登录/注册后可看大图

错别字啊有木有

嘿嘿，没有错别字
看样子咱论坛还是能召唤的么~~哇哈哈

传说中的谁

artmunich 发表于 2011-8-22 20:18

                               
登录/注册后可看大图

嘿嘿，没有错别字
看样子咱论坛还是能召唤的么~~哇哈哈

P啊，我自己屁颠屁颠跑进来的，看看你有木有好好干活

artmunich

传说中的谁 发表于 2011-8-22 20:19

                               
登录/注册后可看大图

P啊，我自己屁颠屁颠跑进来的，看看你有木有好好干活

版主你太小气了。。。我这么用心的干活发帖，都木有一毛钱。。。桑心呐。。。

传说中的谁

artmunich 发表于 2011-8-22 21:16

                               
登录/注册后可看大图

版主你太小气了。。。我这么用心的干活发帖，都木有一毛钱。。。桑心呐。。。

忘记啦