（已解决）画假相当位温图出现空缺是什么原因？

tszhaoweiwen

大家好，我画假相当位温图出现空缺是什么原因？

GS脚本如下：
'reinit'
'open C:\Users\sta288\Desktop\fnl1013\fnl08.ctl'
'set t 3'
'set grads off'
'set grid off'
'set csmooth on'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set grads off'
'set lon 112.5'
'set lat 0 40'
'set lev 1000 300'
'set zlog on'
'set ylevs 1000 925 900 850 800 750 700 650 600 550 500 450 400 350 300 '

'define br=(0.1158*log(RHprs)-1.5332787)*TMPprs-789.92816'                              
'define cr=789.92816*TMPprs'                                                         
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'                                               
'define pc2=lev*pow(tc2/TMPprs,3.5005574)'                                                        
'define etc2=6.1078*pow(273.16/tc2,5.1347779)*exp(3148.0973*(tc2-273.16)/(126.06334*tc2))'   
'define wc2=0.622*etc2/(pc2-etc2)'                                                        
'define ed2=tc2*pow(1000.0/(pc2-etc2),0.28765924)'                                          
'define lc2=2500.79-2.3697*(tc2-273.16)'                                                  
'define eqt=ed2*exp(wc2*lc2/1.0048/tc2)'

'set gxout shaded'
'set cthick 6'
'set clskip 1'
'set clevs 320 330 340 350 360'
'd eqt'
'run cbarn.gs 1 0'
'set ylpos 0 r'
'set gxout contour'
'set cthick 11'
'set cint 4'
'd eqt'

'set gxout contour'
'set clevs 302 304 306 308 310 312'
'd eqt'
'writehz 1 8.1 假相当位温 1 2.5 1 1.3 0.6 0 7'
'printim D:\eqt2.png white'

小小人

我以前做风场图的时候也出现过，我出现的原因是由于首先数据缺失，然后是写fortran程序转化为.grd文件时出现了错误，我给这位的一点自己小小的拙见是检查数据和程序。

river

你的截图太小了，家园的水印把重要的东西给挡住了啊。看着应该是log函数的问题，肯定是有数据超出了它的定义域。另外还想说的就是楼主用的fnl资料的话，应该可以画到100hPa的吧

空谷幽竹

还能出现这样的空缺啊  数据缺失吧 或者就是画图软件的问题了

兰北

log函数的问题，应该是有数据超出了它的定义域

或者是gs文件格式的问题，有空格啊什么的，你再一点点复制一遍，估计就好了

wuwei2163

数据的问题吧 或者用程序处理的时候出错了

zxhhuanyue

这个gs对不对呢？从别处看到的好像不一样啊

zxhhuanyue

。。。。减我钱，呜呜

木雨

以前也遇到过类似的情况，但始终没有解决

天道酬勤

怎么解决的，请楼主解释一下吧。。。