- 积分
- 25437
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-12-13
- 最后登录
- 1970-1-1
![[wet510] 粉丝数:55 微博数:3 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
绘制多年夏季平均场的水汽通量,结果出来的等值线图不连贯,断断续续的,不知是脚本问题还是数据问题,向大家请教,期待高手指导!!- 'reinit'
- 'set grads off'
- 'set grid off'
- 'set timelab off'
- 'sdfopen E:\nc_data_monthly\air.mon.mean.nc'
- 'sdfopen E:\nc_data_monthly\rhum.mon.mean.nc'
- 'sdfopen E:\nc_data_monthly\uwnd.mon.mean.nc'
- 'sdfopen E:\nc_data_monthly\vwnd.mon.mean.nc'
- 'sdfopen E:\nc_data_monthly\pres.sfc.mon.mean.nc'
- 'set lon 115 145'
- 'set lat 35 60'
- 'set z 1 8'
- 'set t 157 756'
- *file='D:\NEdraw\vapor_flux\cv_ave\ave.dat'
- *'set gxout fwrite'
- *'set fwrite 'file
- 'define prs=lev'
- 'define es=6.112*exp(17.2693822*(air.1-273.16)/(air.1-29.65))'
- 'define qs=0.622*es/(prs-0.378*es)'
- 'define qv=rhum.2*qs/100'
- 'define uq=qv*uwnd.3*1000'
- 'define vq=qv*vwnd.4*1000'
- 'define uqs=vint(pres.5(lev=0),uq,300)'
- 'define vqs=vint(pres.5(lev=0),vq,300)'
- 'set z 1'
- 'define sq=mag(uqs,vqs)'
- 'set t 6 8'
- 'define sq1=ave(sq,t+156,t=756,12)'
- 'modify sq1 seasonal'
- 'define sq2=sum(sq1,t=162,t=164)'
- 'd sq2'
复制代码 下面是出的图:
丑了点,见谅
|
|