请教前辈TLnp绘图中的问题

wygohappy

请教各位大神帮忙，使用资料ERA5 nc格式数据画tlnp，报错

附上代码：
reinit'
'sdfopen F:\data\202008\20200807_0815.nc'
'set grid off'
'set grads off'
'set lev 1000 10'
'set lat 45.1'
'set lon 115.3'
tt=1
while(tt<=2)
'c'
'set grads off'
'set t 'tt' '
'define tmp=t-273.16'
'define RH=r'
'td=tmp-((14.55+0.114*tmp)*(1-0.01*RH) + pow((2.5+0.007*tmp)*(1-0.01*RH),3) + (15.9+0.37*tmp)*pow((1-0.01*RH),14))'
'define uv=mag(v,u)'
'define pi=3.14159'
if (u>=0.)
if (v>=0.)
'define dir=180+atan2(u,v)*180/pi'
endif
endif
if (v<=0.)
if (u>=0.)
'define dir=270+atan2(u,v)*180/pi'
endif
endif
if (u<=0.)
if (v>=0.)
'define dir=180-atan2(u,v)*180/pi'
endif
endif
if (v<0.)
if (u<=0.)
'define dir=270-atan2(u,v)*180/pi'
endif
endif
*if(ugrdprs > 0.);'define dir=90-180/pi*atan2(vgrdprs,ugrdprs)';endif
*if(ugrdprs < 0.);'define dir=270-180/pi*atan2(vgrdprs,ugrdprs)';endif
rc=plotskew(tmp,td,uv,dir)
'printim D:\data\TlnP\'%tt%'.gif  white'
tt=tt+1
endwhile

后续代码为自带探空gs文件内容粘贴，行数太多省略
function plotskew(sndtemp,snddewp,sndspd,snddir)
........

出图：

提示报错如下：

其中报错行代码：
40:   rc=plotskew(tmp,td,uv,dir)
878:     PTemp=LiftDry(tloop,1000,_pmin,1,_pmin,_pmax)     
1836:    temp=starttk*_powpres.subscr*powstart-273.15

wygohappy

已解决，set lev 的问题，层次改为1000 100就好了。

西南反气旋

楼主您好，请问“后续代码为自带探空gs文件内容粘贴，行数太多省略
function plotskew(sndtemp,snddewp,sndspd,snddir)”，这个自带的探空GS在哪里可以找到

djg

西南反气旋 发表于 2022-6-14 09:29
楼主您好，请问“后续代码为自带探空gs文件内容粘贴，行数太多省略
function plotskew(sndtemp,snddewp,sn ...

OpenGrADS 文件夹下的 \Contents\Resources\Scripts

sk378179128

请问下出了图以后怎么看时间啊？