请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2599|回复: 1

[图形美化] 求助:绘制整层水汽通量图像和文献中图像不一样

[复制链接]

新浪微博达人勋

发表于 2023-1-14 11:12:48 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 WOT 于 2023-1-14 11:16 编辑

本人想绘制1982-2020年5~8月的整层水汽通量图像,用的是necp月平均数据,但是图像和文献中不相符,箭头长短也不一样,部分阴影区域也不一样。文献中图像描述为” Climatology of the vertically integrated water vapor flux (vector) and its magnitude (shading) in the wet season from 1982 to 2020.“是不是文献中绘制的阴影区域是水汽通量的量级而我画的是水汽通量结果,所以导致我这个图像有偏差?或者是数据不对呢?如果要用逐日是数据,气压是否也要用逐日数据呢?
'reinit'
'sdfopen E:\GrADS\testwww\uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\vwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\shum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\pres.mon.mean.nc'


'set parea 1.5 10.0 1.5 8.0'
'set map 1 1 3'
'set mpdset cnworld cnriver'


'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 10'
'set ylint 10'


'set t 413 869'
'set z 1 8 '
'define qu5=uwnd*shum.3/1000'
'define qv5=vwnd.2*shum.3/1000'
'set t 413'
'set z 1 8 '
'define qud5=ave(qu5,t+0,t=869,12)'
'define qvd5=ave(qv5,t+0,t=869,12)'
'define qdiv5=hdivg(qud5,qvd5)'


'set t 414 870'
'set z 1 8 '
'define qu6=uwnd*shum.3/1000'
'define qv6=vwnd.2*shum.3/1000'
'set t 414'
'set z 1 8 '
'define qud6=ave(qu6,t+0,t=870,12)'
'define qvd6=ave(qv6,t+0,t=870,12)'
'define qdiv6=hdivg(qud6,qvd6)'


'set t 415 871'
'set z 1 8 '
'define qu7=uwnd*shum.3/1000'
'define qv7=vwnd.2*shum.3/1000'
'set t 415'
'set z 1 8 '
'define qud7=ave(qu7,t+0,t=871,12)'
'define qvd7=ave(qv7,t+0,t=871,12)'
'define qdiv7=hdivg(qud7,qvd7)'


'set t 416 872'
'set z 1 8 '
'define qu8=uwnd*shum.3/1000'
'define qv8=vwnd.2*shum.3/1000'
'set t 416'
'set z 1 8 '
'define qud8=ave(qu8,t+0,t=872,12)'
'define qvd8=ave(qv8,t+0,t=872,12)'
'define qdiv8=hdivg(qud8,qvd8)'

'set z 1 8 '
'define qud=(qud5+qud6+qud7+qud8)/4'
'define qvd=(qvd5+qvd6+qvd7+qvd8)/4'
'define qdiv=(qdiv5+qdiv6+qdiv7+qdiv8)/4'


'set z 1'
'define qus=vint(pres.4(z=1),qud,300)'
'define qvs=vint(pres.4(z=1),qvd,300)'


'define shuiqi=mag(qus,qvs)'


*'define sandu=hdivg(qus,qvs)'




'set lon 40 120'
'set lat 20 60'


'set gxout shaded'
'd shuiqi'
'set gxout barb'
'set ccolor 1'
'set gxout vector'
'd qus;qvs'
'cbarn 1.2 0'


'draw title 1000-300hpa vapor flux and divergence'


'printim E:\GrADS\test\44.png  white x2000 y1800'
'disable print'


文献中图像

文献中图像

gradsz绘制的图像

gradsz绘制的图像
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-1-17 08:33:59 | 显示全部楼层
色标都不一致,肯定不一致
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表