弱弱的问一个问题，关于NCEP再分析资料的试用。。。。

小cc爱吐槽

grads新人一枚~~~
想画一张假相当位温的剖面图，下载了air.2012.nc和rhum.2012.nc，均是每6小时的数据。按照论坛上提供的一些方法，我编写了gs文件，运行时出错了。原因是air文件等压面层是17层，而rhum是8层的，不能一起用，请问这个问题如何解决？

jeydragon

我觉得你应该将出错信息贴出来，以及你的gs文件拿出来，这样大家才好有的放矢

小cc爱吐槽

jeydragon 发表于 2012-12-3 10:04

                               
登录/注册后可看大图

我觉得你应该将出错信息贴出来，以及你的gs文件拿出来，这样大家才好有的放矢

thes.gs (782 Bytes, 下载次数: 4)

2012-12-3 10:14 上传

点击文件名下载附件 好吧，真心求助啊

mofangbao

楼主的gs要这样贴
'reinit'
'sdfopen c:\2012\air.2012.nc'
'sdfopen c:\2012\rhum.2012.nc'
'set lev 1000 100'
if(air>263.0);
'define a1=0.622*6.11*exp(17.26*(air-273.16)/(air-35.86))'
'define b1=lev-0.278*exp(17.26*(air-273.16)/(air-35.86))'
'define qs1=a1/b1'
endif
if(air<=263.0);
'define a1=0.622*6.11*EXP(21.87*(air-273.16)/(air-7.66))'
'define b1=lev-0.278*EXP(21.87*(air-273.16)/(air-7.66))'
'define qs1=a1/b1'
endif
'define q1=qs1*rhum.2'
'define e1=lev*q1/100./(0.62197+q1/100.0)'
'define tk1=55.0+2840.0/(3.5*log(air)-log(e1)-4.805)'
'define pot1=air*pow((1000/lev),(0.2854*(1.0-0.28*q1/100.0)))'
'define ept1=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0))'
'set lon 116'
'set lat 32 42'
'set t 773'
'set cint 3'
'd ept1'
'printim c:\shi.gif gif white'
;

小cc爱吐槽

mofangbao 发表于 2012-12-3 13:27

                               
登录/注册后可看大图

楼主的gs要这样贴
'reinit'
'sdfopen c:\2012\air.2012.nc'

哎呀，谢谢谢谢！新人哈，好多不对的，包涵包涵下~~~

小cc爱吐槽

哎~~~哪位给个解释啊，为什么相对湿度的资料\rhum.2012.nc的层数和别的不一致呢？为什么这么特别啊？这样怎么用啊？{:eb328:}{:eb328:}

小cc爱吐槽

q file 后，相对湿度文件Zsize=8，而其他文件都是17 ，这样就不能在一起用了吗？尤其是计算假相当位温和水汽通量散度的时候，需要打开好几个nc文件啊，这样到底怎么解决呢？

小cc爱吐槽

我问的这个问题很弱智吗？还是怎么了？怎么一直没有人回答？呼叫版主啊，新人在这里就这么难混吗？

小cc爱吐槽

{:eb306:}

river

既然楼主自己知道shum没有300以上的数据，那就算不了300以上的假相当位温呗。把高度设置成1000到300不就行了么