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

liudage1983

正好学习一下，谢谢

qxtlyf

楼主辛苦了,

yifu422

请问楼主我基本用的是你的gs,只是改了路径，数据也都分别画过图，没有错，为什么grads运行报错，而且显示的界面有图，不过x轴的日期不显示，而且生成的gmf也打不开。（另外请问为什么要生成gif图形？）
另附上gs文件：
'reinit'
'sdfopen d:\ncep\uwnd.2015.nc'
'sdfopen d:\ncep\vwnd.2015.nc'
'sdfopen d:\ncep\rhum.2015.nc'
'sdfopen d:\ncep\shum.2015.nc'
'sdfopen d:\ncep\omega.2015.nc'
'sdfopen d:\ncep\air.2015.nc'
'set t 101 120'
'set lev 1000 100'
'define div=hdivg(uwnd.1,vwnd.2)*1e5'
'set lev 1000 300'
'define qdiv=hdivg(uwnd.1*shum.4*1e3/9.8,vwnd.2*shum.4*1e3/9.8)*(1e5)'
'define qu=(uwnd.1*shum.4*1e3/9.8)'
'define qv=(vwnd.2*shum.4*1e3/9.8)'
'define br=(0.1158*log(rhum.3)-1.5332787)*air.6-789.92816'                              
'define cr=789.92816*air.6'                                                         
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'                                               
'define pc2=lev*pow(tc2/air.6,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 grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'enable print d:\2015snow\out\jia.gmf'
'set grads off'
'set lon 116.5'
'set lat 31.8'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300'
'set xlabs 08|14|20|02|08|14|20|02|08|14|20|02|08|14|20|02|08|14|20|02|'
'set gxout contour'
'set cint 3'
'set cthick 6'
'set clskip 1'
'd eqt'
'set ylpos 0 r'
'set gxout contour'
'set cint 3'
'set cthick 11'
'd eqt'
'writehz 0.4 4 气压（hPa） 1 2 1 1.3 0.6 90 7'
'writehz 1 8.1 假相当位温（六安 经度116.5 纬度31.8） 1 2.5 1 1.3 0.6 0 11'
writehz 0.7 1.0 1月26日 1 2 1 1.3 0.6 0 7'
'writehz 2.7 1.0 1月27日 1 2 1 1.3 0.6 0 7'
'writehz 5.7 1.0 1月28日 1 2 1 1.3 0.6 0 7'
'writehz 9.0 1.0 1月29日 1 2 1 1.3 0.6 0 7'
'writehz 12.0 1.0 1月30日 1 2 1 1.3 0.6 0 7'
'print'
'disable print'
;

yifu422

不好意思，自己已经发现错误了，原来是有一条命令忘了个引号，运行时还是报错，不过出图倒是正常了。

river

yifu422 发表于 2015-3-28 09:16
不好意思，自己已经发现错误了，原来是有一条命令忘了个引号，运行时还是报错，不过出图倒是正常了。

数据有问题，怎么会有负数呢？

liushihua

嘻嘻，，又捡到一个有用的！！谢谢楼主

笑哈乐呵

好东东，收藏了

19921219

真的挺有用的

jr921102

赞一个~好好学习剖面图的画法

云大学子

楼主，你帮了我大忙了