传说版主，画出来的斜温图没有绿色湿绝热线，也没有状态曲线

海蓝蓝

就是用你的代码，还有FNL一度的资料。但是，画出来的斜温图没有绿色湿绝热线，也没有状态曲线。
用的是opengrads2.0，xp系统。
画出来的图和meehooqq画出来的这个一样：http://bbs.06climate.com/data/attachment/forum/201111/22/1304519ljn1x164l6x4b8l.jpg
脚本附上了，可否给看一眼？ tlnp.gs (62.01 KB, 下载次数: 46)

2012-1-18 09:37 上传

点击文件名下载附件

脚本如下（不用下载那个附件了）：
'reinit'
'open fnl_20110428_12_00.ctl'
'set lon 123'
'set lat 23'
'set t 1'
'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)

函数function plotskew略

marklandy

斜温图作图的地方，顺便翻译下：
* ------------------------ Output Options --------------------------------
*
* All variables below are logical .. 1=yes, 0=no, unless otherwise     ！1是画出，0是不画
* specified.
*
* DrawBarb = Draw wind barbs along right side of plot                    ！画风向杆
* DrawThet = Draw dry adiabats                                                    ！画干绝热线
* DrawThtw = Draw moist pseudo-adiabats                                    ！画湿绝热线
* DrawTemp = Draw temperature lines                                           ！画等温线（斜的）
* DrawMix  = Draw mixing ratio lines                                               ！画混合比线
* DrawTSnd = Draw temperature sounding                                    ！画温度线（层结曲线）
* DrawDSnd = Draw dewpoint sounding                                         ！画露点线
* DrawRH   = Draw relative humidity sounding                                ！画相对湿度线
* DrawPrcl = Draw parcel path from surface upward                       ！画状态曲线
* DrawPMax = Draw parcel path from most unstable level upward  ！从最不稳定层画状态曲线（要在DrawPrcl=1的条件下）
* DrawIndx = Display stability indices & CAPE                                  ！是否显示CAPE
* DrawHeli = Calculate and display absolute and storm-relative helicity   ！是否计算风暴螺旋度
* DrawHodo = Draw hodograph                                                         ！右上风速图
* DrawPLev = Draw Pressure Levels                                                   ！高度hpa
* DrawZLev = Draw height levels and lines                                          ！高度km
*            0 = no lines
*            1 = above ground level (AGL)
*            2 = above sea level (ASL)
* DrawZSTD = Draw Height levels using standard atm lapse rate           
* LblAxes  = Label the x,y axes (temperature, pressure,mixing ratio)               ！xy轴标注
*
* ThtwStop = Pressure level at which to stop drawing Theta-w lines                 ！停止画theta-w图高度
* MixStop  = Pressure level at which to stop drawing Mixratio lines                   ！停止画混合比图高度

传说中的谁

plotskew.gs里面有一些可调参数的，不知道你调过没有？

topmad

笨儿　　你好厉害啊

海蓝蓝

传说中的谁 发表于 2012-1-18 17:49

                               
登录/注册后可看大图

plotskew.gs里面有一些可调参数的，不知道你调过没有？

没有调。暂时看不懂呢。
你还记得调了那些地方？

传说中的谁

海蓝蓝 发表于 2012-1-18 22:04

                               
登录/注册后可看大图

没有调。暂时看不懂呢。
你还记得调了那些地方？

想不起来了……建议看看吧，虽然是英文，但大概应该能看懂的

传说中的谁

topmad 发表于 2012-1-18 20:39

                               
登录/注册后可看大图

笨儿　　你好厉害啊

给笨笨邮一只烤全羊过来

topmad

传说中的谁 发表于 2012-1-18 23:29

                               
登录/注册后可看大图

给笨笨邮一只烤全羊过来

刚拷出来的才好吃，过来吧，还有妹子陪你喝酒

feiji158

我也是这样的，不知如何解决！

海蓝蓝

feiji158 发表于 2012-1-20 14:50

                               
登录/注册后可看大图

我也是这样的，不知如何解决！

看来带有普遍性啊。版主来给解决一下！

风荷

改drawthtw＝0为drawthtw＝1画是绝热线