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

青衣

感谢楼主，很不错

拖延症

撒花点赞

dzn

太棒了，非常感谢！

CDH

版主您好，我现在就是像绘制一个简单的固定维度，经度和高度变化的剖面图，我下载了你的poumian.gs，仿照这个写了一个处理自己的数据的gs.我分别贴上我的ctl 和gs
ctl:
dset F:/FNL/sidr/profilemap/Temperature/fnl_200711%d2_%h2_00.grib1.spasub.hui139497
options template
index F:/FNL/sidr/profilemap/Temperature/fnl_20071110_00_00.grib1.spasub.hui139497.idx
undef 9.999E+20
title F:/FNL/sidr/profilemap/Temperature/fnl_20071110_00_00.grib1.spasub.hui139497
*  produced by grib2ctl v0.9.12.5p33k
dtype grib 255
options yrev
ydef 36 linear 0.000000 1
xdef 61 linear 60.000000 1.000000
tdef 28 linear 00Z10nov2007 6hr
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 1
TMPprs 26 11,100,0 ** Temp. [K]
ENDVARS


gs:
'reinit'
'open F:\FNL\sidr\profilemap\Temperature\fnl_200711.ctl'
'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'

i=1
while(i<=28)
'set t 'i''
'set lev 1000 300'
'set lon 60 120'
'set lat 10'
'set zlog on'
'set gxout contour'
'd TMPprs'
'printim F:\FNL\sidr\profilemap\Temperature\'i'.gif white'
'c'
i=i+1
endwhile
;
但是我得到的图是这样的，不知道是哪里出的问题呢？期待您的回复。万分感谢。

CDH

而我想得到的是类似这样的效果图：

river

CDH 发表于 2015-11-17 22:15
而我想得到的是类似这样的效果图：

那你得先搞清楚别人画的是什么变量啊，难道就是普通的温度？

CDH

river 发表于 2015-11-17 22:23
那你得先搞清楚别人画的是什么变量啊，难道就是普通的温度？

您的意思是说，我这样的FNL的数据只能画出来这样的结果？

CDH

river 发表于 2015-11-17 22:23
那你得先搞清楚别人画的是什么变量啊，难道就是普通的温度？

您的意思是说，我这样的FNL的数据只能画出来这样的结果？我现在的问题是因为我选择的数据的问题？谢谢您的回复。

river

CDH 发表于 2015-11-17 22:26
您的意思是说，我这样的FNL的数据只能画出来这样的结果？我现在的问题是因为我选择的数据的问题？谢谢您 ...

我没那么说···你要画温度随高度变化那么剧烈的，那好歹得有个什么天气过程之类的吧。没什么特别的不就是那样温度随高度降低么···

CDH

river 发表于 2015-11-17 22:50
我没那么说···你要画温度随高度变化那么剧烈的，那好歹得有个什么天气过程之类的吧。没什么特别的不就 ...

嗯，就是的，我后来也想明白您的意思了。下载的这段资料是南亚的一次台风sidr，我想分析一下台风的垂直结构，还下载了垂直速度等资料。
不知道您对分析台风垂直结构这方面有没有相关的了解，我也再查一下，是不是要用其他的资料才行，这个1°*1°的应该不能显示出台风的结构。
谢谢您这么晚还回复我，感谢。