- 积分
- 28
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-11-25
- 最后登录
- 1970-1-1
|
发表于 2017-4-18 12:08:11
|
显示全部楼层
'reinit'
'sdfopen E:\daxuewenjian\air.2015.nc'
'sdfopen E:\daxuewenjian\rhum.2015.nc'
'sdfopen E:\daxuewenjian\uwnd.2015.nc'
'sdfopen E:\daxuewenjian\vwnd.2015.nc'
'enable print E:\daxuewenjian\shuiqi.gmf'
'set mpdset cnworld'
'set lon 105 120'
'set lat 18 25'
'set lev 925'
t=275
while(t<=365)
'set t 't''
'set grads off'
'define p=lev'
'define e=exp(17.67*(air-273.15)/(air-29.65))*(rhum.2/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.3*q/9.8'
'define q2=vwnd.4*q/9.8'
'set gxout shaded'
'set cint 20'
'd mag(q1,q2)'
'set gxout barb'
'set ccolor 4'
'd uwnd.3*2.5;vwnd.4*2.5'
'cbarn 1 0'
t=t+1
endwhile
'print'
'disable print'
'c'。
这是我的gs,希望大神帮忙改一下。
我只要2015年10月4号前后的图就可以了,一层的水汽,不需要平均。 |
|