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

kust55

剖面面图很难画，谢谢1

ziqiangbuxi

对这个方面不知道，开始学习，楼主总结的太好了。

liy72427

好帖子，学习了

liy72427

之前老是出错，现在豁然明了了哈哈哈

青衣

感谢楼主，很不错

拖延症

撒花点赞

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
而我想得到的是类似这样的效果图：

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