利用NCEP资料画剖面图的一点儿经验

☆梧桐花开☆

'reinit'
'open e:\ncep\fnl.ctl'
'set t 5 13'
'set lon 115 117'
'set lat 20 40'
'define prs=lev'
'set ccolor 1'
'define  l1=TMPprs'
if(l1>263.0);
'define a1=0.622*6.11*EXP(17.26*(TMPprs-273.16)/(TMPprs-35.86))'
'define b1=prs-0.278*EXP(17.26*(TMPprs-273.16)/(TMPprs-35.86))'
'define  qs1=a1/b1'
endif
if(l1<=263.0);
'define a1=0.622*6.11*EXP(21.87*(TMPprs-273.16)/(TMPprs-7.66))'
'define b1=prs-0.278*EXP(21.87*(TMPprs-273.16)/(TMPprs-7.66))'
'define qs1=a1/b1'
endif
'define q1=qs1*RHprs'
'define e1=lev*q1/100./(0.62197+q1/100.0)'
'define tk1=55.0+2840.0/(3.5*log(TMPprs)-log(e1)-4.805)'
'define pot1=TMPprs*pow((1000/lev),(0.2854*(1.0-0.28*q1/100.0)))'
'define ept1=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0))'
'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 116'
'set lat 35'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set xlabs 0720|0802|0808|0814|0820|0902|0908|0914|0920'
'set gxout contour'
'set cint 3'
'set cthick 6'
'set clskip 1'
'd ept1'
'printim e:\ncep\fnl\'0902'.gif gif white'


请大家帮我看看假相当位温的时间高度剖面，为什么画不出来？

river

☆梧桐花开☆ 发表于 2013-10-23 00:30

                               
登录/注册后可看大图

'reinit'
'open e:\ncep\fnl.ctl'
'set t 5 13'

你的高度变化呢，高度不变化还能画剖面呢·····
还有有问题要提到点子上，画不出来一般就是有错误，错误提示呢？直接贴一个脚本程序出来，短的别人还能看一下，长的谁没事儿一句一句的去看啊

Figo

楼主辛苦了！

cat1106123

这个厉害了！

德意志

真是好东西，这几天就急等这个东西，楼主好人一生平安~~

376760020

顶之~~~~~

clare

请教LZ你的中文是怎么显示出来的呢？用的家园的集成grads但是使用writehz中字是乱码

river

clare 发表于 2013-11-9 07:31

                               
登录/注册后可看大图

请教LZ你的中文是怎么显示出来的呢？用的家园的集成grads但是使用writehz中字是乱码

乱码？好奇怪啊，你上一张图片看看

1079652110

好牛

1079652110

努力学习中，谢谢了