- 积分
- 1831
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-6
- 最后登录
- 1970-1-1
|
NCL
系统平台: |
linux |
问题截图: |
|
问题概况: |
我做出了某一层散度区域平均的时间序列图,发现平均散度量级达到1.e-3之多,而先前做的剖面图量级正常(1.e-5) |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
2 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
俩张图的散度都乘了1.e05,剖面图上散涡的量级基本在1.e-05上,而前图(折线)的散度量级却非常大,请问这是怎么回事呢?
下面是折线图的脚本:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
begin
in = addfile("/home/zssapr/WRF3.5.1/20110710/WRFV3/test/em_real/wrfout_d02_2011-07-05_00:00:00","r")
FNM =in->FNM
U =in->U
ua = wrf_user_unstagger(U,U@stagger)
V = in->V
va = wrf_user_unstagger(V,V@stagger)
lat = in->XLAT(0,:,0)
lon = in->XLONG(0,0,:)
T =in->Times
scale=1.e05
dv = new ( dimsizes(ua), typeof(ua), ua@_FillValue )
dv@long_name = "divengence"
dv@units = "scaled"
dv = uv2dv_cfd (ua,va,lat,lon, 2)
ud = new ( dimsizes(ua), typeof(ua), ua@_FillValue )
vd = new ( dimsizes(va), typeof(va), va@_FillValue )
copy_VarCoords(ua, ud )
copy_VarCoords(ua, vd )
dv2uvg(dv,ud,vd)
ud@long_name = "Zonal Divergent Wind"
ud@units = ua@units
vd@long_name = "Meridional Divergent Wind"
vd@units = va@units
dv2 = uv2dv_cfd (ud,vd,lat,lon, 2)*scale
vr = new ( dimsizes(ua), typeof(ua), ua@_FillValue )
vr@long_name = "vorticity"
vr@units = "scaled"
vr = uv2vr_cfd (ua,va,lat,lon, 2)
uort = new ( dimsizes(ua), typeof(ua), ua@_FillValue )
vort = new ( dimsizes(va), typeof(va), va@_FillValue )
copy_VarCoords(ua, uort )
copy_VarCoords(ua, vort )
vr2uvg(vr,uort,vort) ; div ==> rotational wind components
uort@long_name = "Zonal Divergent Wind"
uort@units = ua@units
vort@long_name = "Meridional Divergent Wind"
vort@units = va@units
vr2 = uv2vr_cfd (uort,vort,lat,lon, 2)*scale
Dv2 =dv2(:,25,:,:)
Vr2 =vr2(:,25,:,:)
Dv2abs =abs(Dv2)
Vr2abs =abs(Vr2)
avedv3 = dim_avg_n_Wrap(Dv2abs(:,:,: ), (/1,2/))
avevr3 = dim_avg_n_Wrap(Vr2abs(:,:,: ), (/1,2/))
wks = gsn_open_wks("ps","vrdv_lev")
res = True
res@tiMainString = "Vort&Div Profile Plot" ; add title
res@xyDashPatterns = 12 ; choose dash patterns
res@vpWidthF = 0.5 ; change aspect ratio of plot
res@vpHeightF = 0.3
res2 =True
res2@xyDashPatterns = 0
res2@vpWidthF = 0.5 ; change aspect ratio of plot
res2@vpHeightF = 0.3
plot2 =gsn_csm_y (wks,avedv3(::24),res)
end
|
|