爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 5468|回复: 1

时间序列图中物理量量级不对

[复制链接]
发表于 2014-6-4 15:17:40 | 显示全部楼层 |阅读模式
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

密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-6-5 15:59:16 | 显示全部楼层
就是说为何将风场分解为散度风和涡旋风后再求散度,散度会有这么大的变化呢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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