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

river

   最近浏览很多帖子，学到不少东西。咱们论坛高手很多啊!在这儿先赞一个！但是也看到很多人问剖面图的画法，论坛里相关的帖子啊 gs文件也都很多。今天我也来凑个数，主要是为了抛砖引玉，希望论坛的大大们可以出个专题总结一下这类图的画法。

  进入正题：画剖面图的话，常常需要的有时间--高度剖面图、固定纬度(经度)的经度(纬度)--高度剖面度、延任意直线的剖面图，还有····额，想不出来了，智商有限啊。思路比较简单，根据不同的需要，固定其中的一个值，让其他两个量变化，然后出图就行了。  但是有些时候往往总是画不出来，各种出错。
  我感觉最大的一个原因就是我们在学习的时候往往都是比较简单的例子，有一个物理量的多时次多层次的资料，然后在gs里固定经度（纬度），让时间或者高度变化，然后display变量，出图。这就造成一个错觉就是一旦要画某个剖面图，第一步就是要固定某个变量。但是我们出图需要的物理量往往是需要计算的。而计算过程中是需要这个量变化的，所以一旦先固定了某个变量，那计算过程一定会出错。

所以大家在画剖面图的时候，一定不要急于固定某一个量，要等计算过程完成了，在出图的时候再固定这个量，这样就不会出错啦。

下面给出一个我自己画时间高度剖面图的例子：

'reinit'
'sdfopen f:\2012\uwnd.2012.nc'
'sdfopen f:\2012\vwnd.2012.nc'
'sdfopen f:\2012\rhum.2012.nc'
'sdfopen f:\2012\shum.2012.nc'
'sdfopen f:\2012\omega.2012.nc'
'sdfopen f:\2012\air.2012.nc'
'set t 709 721'
*-----水平散度------
'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)'
*--------GRADS设置------------
'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'

*---------输出图形1------------
'enable print f:\2012\1\qdiv.gmf'
'set lon 125'
'set lat 35'
'set zlog on'
*'set xlint 4'
'set ylevs 1000 925 850 700 600 500 400 300'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'set gxout shaded'
'set cint 0.3'
'set cmax 0'
'd qdiv'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.3'
'set cthick 11'
*'set clskip 1'

'd qdiv'
'writehz 0.4 4 气压（hPa） 1 2 1 1.3 0.6 90 11'
'writehz 1 8.1 水汽通量散度（经度125纬度35） 1 2.5 1 1.3 0.6 0 11'
'writehz 0.7 1.0 6月26日 1 2 1 1.3 0.6 0 7'
'writehz 2.8 1.0 6月27日 1 2 1 1.3 0.6 0 7'
'writehz 5.9 1.0 6月28日 1 2 1 1.3 0.6 0 7'
'writehz 9.1 1.0 6月29日 1 2 1 1.3 0.6 0 7'

'print'
'disable print'
'printim f:/2012/1/qdiv.gif white'
'c'
*---------输出图形2------------
'enable print f:\2012\1\jia.gmf'
'set grads off'
'set lon 125'
'set lat 35'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'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'
*'set clskip 1'
'd eqt'

'writehz 0.4 4 气压（hPa） 1 2 1 1.3 0.6 90 7'
'writehz 1 8.1 假相当位温（经度125纬度35） 1 2.5 1 1.3 0.6 0 7'
'writehz 0.7 1.0 6月26日 1 2 1 1.3 0.6 0 7'
'writehz 2.8 1.0 6月27日 1 2 1 1.3 0.6 0 7'
'writehz 5.9 1.0 6月28日 1 2 1 1.3 0.6 0 7'
'writehz 9.1 1.0 6月29日 1 2 1 1.3 0.6 0 7'

'print'
'disable print'
'printim f:/2012/1/jia.gif white'
'c'
*---------输出图形3------------
'enable print f:\2012\omega.gmf'
'set lev 1000 100'
'set grads off'
'set lon 125'
'set lat 35'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'set gxout shaded'
'set cint 0.05'
'set cmax 0'
'd omega.5'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.05'
'set cthick 11'
'd omega.5'
'writehz 0.4 4 气压（hPa） 1 2 1 1.3 0.6 90 7'
'writehz 1 8.1 垂直速度（经度125纬度35） 1 2.5 1 1.3 0.6 0 7'
'writehz 0.7 1.0 6月26日 1 2 1 1.3 0.6 0 7'
'writehz 2.8 1.0 6月27日 1 2 1 1.3 0.6 0 7'
'writehz 5.9 1.0 6月28日 1 2 1 1.3 0.6 0 7'
'writehz 9.1 1.0 6月29日 1 2 1 1.3 0.6 0 7'
'print'
'disable print'
'printim f:/2012/1/omega.gif white'
'c'
*---------输出图形4------------
'enable print f:\2012\1\div.gmf'
'set grads off'
'set lon 125'
'set lat 35'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'set gxout shaded'
'set cint 0.5'
'set cmax 0'
'd div'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.5'
'set cthick 11'
*'set clskip 1'
'd div'
'writehz 0.4 4 气压（hPa） 1 2 1 1.3 0.6 90 7'
'writehz 1 8.1 散度（经度125纬度35） 1 2.5 1 1.3 0.6 0 7'
'writehz 0.7 1.0 6月26日 1 2 1 1.3 0.6 0 7'
'writehz 2.8 1.0 6月27日 1 2 1 1.3 0.6 0 7'
'writehz 5.9 1.0 6月28日 1 2 1 1.3 0.6 0 7'
'writehz 9.1 1.0 6月29日 1 2 1 1.3 0.6 0 7'
'print'
'disable print'
'printim f:/2012/1/div.gif white'
;



我出的图


如果要画多时次的固定纬度(经度)的经度(纬度)--高度剖面图的话，还要再加上时间循环，这个给个例子放在附件里吧。 poumian.gs (2.05 KB, 下载次数: 797)

2012-10-18 21:20 上传

点击文件名下载附件
阅读权限: 10

如果不是NCEP资料的话，道理是一样的，只要转成grads能识别的资料，配好ctl文件，就都能画。
至于任意直线的那个还是等高人吧，我就会这点儿啦。

如果有不对的地方或者更好的方法还请大家赐教啊{:soso__8961432591078930798_3:}

另外有关沿任意直线画剖面的做法请看我的另一个帖子啊 地址：http://bbs.06climate.com/forum.php?mod=viewthread&tid=10966

hustlh

楼主，真是细心了，辛苦了！

最近正好在做剖面图！

Aires

没画过，楼主辛苦了~

言深深

剖面图一般是难以掌握的，很多时候绘制的也不好看，楼主的整的很不错

flyingbird2008

,楼主辛苦了！

mofangbao

整理的挺全的，几个经典的、问的多的要素都写出来了，值得借鉴~

whmooncool

收起来以后用~~

风荷

收起来备用

river

言深深 发表于 2012-10-19 08:15

                               
登录/注册后可看大图

剖面图一般是难以掌握的，很多时候绘制的也不好看，楼主的整的很不错

谢谢，以后还得多交流，很多东西还不会弄

竹岭

正需要学习呢,楼主辛苦了