利用FNL资料画高空图

珮瑶瑶

利用FNL资料画高空图，绘制了高度场 温度场  风场 还有湿度分布。高度场 温度场  风场 都比较正常，但是出现了大范围的高湿区，不知道问题出在哪里？gs文件如下
'open d:\NCDownloads\fnl_20180227_00_00.ctl'

'set mpdset cnworld cnriver'
'set map 8 1 2'
'set grads off'
'set grid off'
'set lon 60 140'
'set lat 10 70'
'set xlint 10'
'set ylint 10'
"set xlpos 0 b"
"set ylpos 0 l"
"set xlopts 1 1 0.15"
"set ylopts 1 1 0.15"


"set lev 850"'set gxout contour'
'set clopts 0'
"set gxout shaded"
'set ccols 4'
"set cint 10"
'set cmin 50'
'set cmax 100'
*'set clab on'
*'set clopts 1 2 0.1'
'd ave(RHprs,t=1,t=1)'
"cbarn 1 1"

"set lev 500"
'set gxout contour'
"set cint 4"
"set ccols 8"
'set cthick 6'
"set clab on"
'set clopts 1 2 0.1'
'd ave(HGTprs/10,t=1,t=1)'


"set lev 500"
'set gxout contour'
"set cint 3"
"set ccols 9"
'set cthick 6'
"set clab on"
'set clopts 2 2 0.1'
'd ave(Tmpprs-273.15,t=1,t=1)'


"set lev 500"
"set gxout barb"
"set digsize 0.05"
"set ccolor 1"
"set cthick 4"
'define u=ave(UGRDprs*2.5,t=1,t=1)'
'define v=ave(VGRDprs*2.5,t=1,t=1)'
'd skip(u,3);skip(v,3)'


'draw title 500Hpa 271200z FEB 2018'
"gxprint D:/500850hpa.png x1600 y1200 white"
"reinit"

hbl1021

直接display高度温度风就可以了吧，应该没必要用ave函数求平均

珮瑶瑶

hbl1021 发表于 2021-2-15 12:54
直接display高度温度风就可以了吧，应该没必要用ave函数求平均

这里用的是单时次的资料，直接set t 1应该就可以了  我这里是参考了前人的一个帖子。但是湿度大和这个好像没有关系

黑马骑士

除了求平均没太懂，脚本没问题。从图上可以看出，高湿位于两个气旋内，而且是850hpa低空的，可能本身数据就是如此分布的呢

珮瑶瑶

黑马骑士 发表于 2021-2-15 22:01
除了求平均没太懂，脚本没问题。从图上可以看出，高湿位于两个气旋内，而且是850hpa低空的，可能本身数据就 ...

我画了500百帕的呢，也是大片的高湿区，和这个850的差不多。南方的可以理解，北方那么大范围，似乎有点说不过去吧

珮瑶瑶

这个是500的相对湿度，gs中相对湿度的地方 850改成了500，两张图几乎一模一样

黑马骑士

珮瑶瑶 发表于 2021-2-15 22:03
我画了500百帕的呢，也是大片的高湿区，和这个850的差不多。南方的可以理解，北方那么大范围，似乎有点说 ...

北方高湿也可能啊，只要天气形势配合。而且北方温度低啊，相对湿度大，又不表示绝对水汽多

珮瑶瑶

黑马骑士 发表于 2021-2-15 22:10
北方高湿也可能啊，只要天气形势配合。而且北方温度低啊，相对湿度大，又不表示绝对水汽多

你说的有道理。我再换一天的资料看看