- 积分
- 15640
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-9-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
先附脚本程序:
'reinit''sdfopen D:\Research_Room\Data\ThirdGroup\nc\multiple\201705_0.5_all_fundamental_variebles.nc'*z represents Geopotential, t represents Temperature, q represents Specific Humidity,*w represents Verticle volocity, u and v represents U and V component of wind, r represents Relative Humidity.* t=16~31(2017.05.05 ~ 2017.05.08)'set grid on 3''set grads off''set parea 1 10 1 7.6''set xlopts 1 6 0.2''set ylopts 1 6 0.2''set x 1''set y 1'*Starting to draw RH''set ylpos -0.0 l''set vrange 10 70''set yaxis 10 70 5''set gxout line''set lev 1000''set t 20 28''set ccolor 2''set cmark 3''set cthick 6''define averh=aave(r,lon=116,lon=121,lat=33,lat=37)''d averh''define aa=ave(averh,t=20,t=28)''set cmark 0''set cstyle 2''set ccolor 2''d aa'*Starting to draw SH''set ylpos -0.0 r''set vrange 0.001 0.013''set yaxis 0.001 0.013 0.001''set gxout line''set lev 1000''set t 20 28''set ccolor 1''set cmark 3''set cthick 6''define averh=aave(q,lon=116,lon=121,lat=33,lat=37)''d averh''define aa=ave(averh,t=20,t=28)''set cmark 0''set cstyle 2''set ccolor 1''d aa''draw title Average RH & SH of 5.6 2:00 ~ 5.8 2:00 in DTRZ''printim D:\Research_Room\Output\picture\meteorology\SH_RH\RH\RH-time_series.png white';
出图效果:
这里有一个问题:右坐标轴为绝对湿度,感觉数值太小,有文献表示
ECMWF数据资料:
数值在10^-1(Unit:g/g),而我用的ECMWF的资料画出来的是10^-3(Unit:kg/kg),不明白为什么,一般大气中绝对湿度的数值跨度很大吗?
|
|