- 积分
- 468
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2022-12-3
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
这是我的脚本,我想绘制1990——2020年850hpa的水汽散度图像,但是我现在只能绘制某一年一个月的图像,数据时间是1948年1月1日到2022年11月1日
求哪位大神帮我改写一下脚本呢,求求了,救救孩子吧
'reinit'
'sdfopen E:\GrADS\testwww\uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\vwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\rhum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\shum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\air.mon.mean.nc'
'set mpdset cnworld'
'set parea 1 9.5 0.5 8'
'set grid off'
'set grads off'
'set lat 30 60'
'set lon 70 100'
'set t 876'
'set lev 850'
'define p=lev'
'define e=exp(17.67*(air.5-273.15)/(air.5-29.65))*(rhum.3/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.1*q/9.8'
'define q2=vwnd.2*q/9.8'
'define p=lev'
'define e2=(6.112*exp((17.67*(air.5-273.16))/(air.5-29.65)))'
'define qs=(0.622*e2/(p-0.378*e2))'
'define q3=qs*rhum.3/100'
'define qdiv=hdivg(q3*uwnd.1/9.8,q3*vwnd.2/9.8)'
'set xlpos -10'
'set ylpos -10'
'set gxout shaded'
'd qdiv'
'run cbarn.gs'
'd mag(q1,q2)'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type b -position m -label on -interval 10 -sinterval 2 -size 0.1 -voffset -0.1 -suffix `3.'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type L -position o -start 30 -end 60 -suffix `3.'
'set gxout barb'
'set digsiz 0.1'
'set cthick 6'
'set ccolor 1'
'd uwnd.1;vwnd.2'
'printim E:\GrADS\testwww\6.png white '
;
|
|