请教有关tlogp图绘制的问题

cloud_alive

楼主现在需要用grads，利用模式输出结果画各高度层温度的tlogp图。本人是初学grads，中间遇到各种问题，请各位大神指教~之前看有位大神的帖子，所以写了如下代码：
'reinit'
'open c:/Users/root/Desktop/grapesdata/post1.ctl_2013060700'
'set lon 123'
'set lat 23'
'set t 9'
'set lev 1000 100'
'define t=tmpprs-273.16'
'define rH=rhprs'
'td=t-((14.55+0.114*t)*(1-0.01*RH) + pow((2.5+0.007*t)*(1-0.01*RH),3) + (15.9+0.37*T)*pow((1-0.01*RH),14))' *计算露点温度
'define uv=mag(vgrdprs,ugrdprs)' *计算风速
'define pi=3.14159'
if (ugrdprs>=0.)
if (vgrdprs>=0.)
'define dir=180+atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (vgrdprs<=0.)
if (ugrdprs>=0.)
'define dir=270+atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (ugrdprs<=0.)
if (vgrdprs>=0.)
'define dir=180-atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (vgrdprs<0.)
if (ugrdprs<=0.)
'define dir=270-atan2(ugrdprs,vgrdprs)*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(t,td,uv,dir)


在这之后将plotskew.gs中的复制到了上面代码的后面
运行了这个文件后，显示如下图所示：
根本没有办法出图，请各位帮帮忙，谢谢大家啦

运行结果

传说中的谁

首先你得弄清需要哪些要素，这些要素在你得数据文件里是什么字段，在示例里又是什么字段

river

提示很明显了，你的ctl里面的变量名称是tmpprs和rhprs吗

cloud_alive

传说中的谁 发表于 2015-4-14 19:25
首先你得弄清需要哪些要素，这些要素在你得数据文件里是什么字段，在示例里又是什么字段

好的，请问要是我的数据文件里没有相对湿度和露点温度，是不是不太好画tiogp图呢

youyou

从提示来看，你换成GRADS2.0版本的试试，因为以前我也遇到过这样的提示，好像就是版本的问题。

cloud_alive

river 发表于 2015-4-14 20:58
提示很明显了，你的ctl里面的变量名称是tmpprs和rhprs吗

发现变量名是不对，但是改了之后还是运行不了

cloud_alive

youyou 发表于 2015-4-15 09:18
从提示来看，你换成GRADS2.0版本的试试，因为以前我也遇到过这样的提示，好像就是版本的问题。

好的，请问你有安装包么

river

cloud_alive 发表于 2015-4-15 09:56
发现变量名是不对，但是改了之后还是运行不了

变量名不对，改了以后还是运行不了？那改完以后运行的错误提示呢？
你觉得挤牙膏一点儿一点儿的，是不是挺累的呢

cloud_alive

river 发表于 2015-4-15 12:20
变量名不对，改了以后还是运行不了？那改完以后运行的错误提示呢？
你觉得挤牙膏一点儿一点儿的，是不是 ...

艾玛，是挺累的呢，sorry~
嘿嘿，改了变量名之后，运行程序，结果变成了file:///C:\Users\root\AppData\Roaming\Tencent\Users\1976816205\QQ\WinTemp\RichOle\3CACA%F7{]~RUT]WR((T1~9.png

运行结果

1. 'open c:/Users/root/Desktop/grapesdata/post1.ctl_2013060700''set lon 123''set lat 23''set t 1''set lev 1000 100''define tem=t-273.16''define rH=qv'*'td=tem-((14.55+0.114*tem)*(1-0.01*RH) + pow((2.5+0.007*tem)*(1-0.01*RH),3) + (15.9+0.37*tem)*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'endifendifif (v<=0.)if (u>=0.)'define dir=270+atan2(u,v)*180/pi'endifendifif (u<=0.)if (v>=0.)'define dir=180-atan2(u,v)*180/pi'endifendifif (v<0.)if (u<=0.)'define dir=270-atan2(u,v)*180/pi'endifendif*if(u > 0.);'define dir=90-180/pi*atan2(v,u)';endif*if(u < 0.);'define dir=270-180/pi*atan2(v,u)';endifrc=plotskew(t,uv,dir)

'open c:/Users/root/Desktop/grapesdata/post1.ctl_2013060700''set lon 123''set lat 23''set t 1''set lev 1000 100''define tem=t-273.16''define rH=qv'*'td=tem-((14.55+0.114*tem)*(1-0.01*RH) + pow((2.5+0.007*tem)*(1-0.01*RH),3) + (15.9+0.37*tem)*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'endifendifif (v<=0.)if (u>=0.)'define dir=270+atan2(u,v)*180/pi'endifendifif (u<=0.)if (v>=0.)'define dir=180-atan2(u,v)*180/pi'endifendifif (v<0.)if (u<=0.)'define dir=270-atan2(u,v)*180/pi'endifendif*if(u > 0.);'define dir=90-180/pi*atan2(v,u)';endif*if(u < 0.);'define dir=270-180/pi*atan2(v,u)';endifrc=plotskew(t,uv,dir)

这是修改后的gs文件，提示第十行有错，但不知道错出在哪里，麻烦给看看呗

river

cloud_alive 发表于 2015-4-15 14:29
艾玛，是挺累的呢，sorry~
嘿嘿，改了变量名之后，运行程序，结果变成了

没法看，这样谁知道第十行是哪一行