用ncep温度和湿度资料算相对位温时出现空白区域

echo要读书

GS文件如下：

'sdfopen f:\data\air.mon.mean.nc'
'sdfopen f:\data\rhum.mon.mean.nc'
'set display color white'
'c'

'set xlopts 1 8 0.25'
'set ylopts 1 8 0.25'
'set clopts 1 8 0.25'
'set lon 90 130'
'set lat 20 50'
'set lev 850'
'set xlevs 90 95 100 105 110 115 120 125 130'
'set ylevs 20 25 30 35 40 45 50'


'set t 390'
'set lev 850'
'define prs=850'
'define es=(6.112*exp(17.67*(air.1(t=390)-273.15)/(air.1(t=390)-29.65)))'
'define q=rhum.2(t=390)*(0.62197*es/(prs-es))/100'
'define theta=air.1(t=390)*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'set cint 0.5'
'set ccolor 1'
'set cthick 16'
'set grid off'
'set grads off'
'd theta'                   <----单独画湿度和温度的时候只改了这一句。
'enable print g:\data\theta.gmf'
'print'
'disable print'

玉儿

我用你的程序试了下，图是正常的。
F:\theta.gmf

river

程序看着没问题

seaou

学到很多，经常来学习了

seaou

大家在这里讨论和学习，真是难得！

echo要读书

玉儿 发表于 2014-4-14 09:02
我用你的程序试了下，图是正常的。

我重复试了很多次都有问题..能否私联一下？我给你发站内信咯.

玉儿

echo要读书 发表于 2014-4-15 01:33
我重复试了很多次都有问题..能否私联一下？我给你发站内信咯.

好的，我去看下。