底图区域不对 经纬度不显示

追风的阳光

longlivehj 发表于 2014-4-21 21:36
用wrf_user_unstagger来做插值

谢谢！。。。。。。。。。

追风的阳光

longlivehj 发表于 2014-4-21 21:36
用wrf_user_unstagger来做插值

你好，我用官网上的脚本试着做出了涡度的图像，上图左边是母网格的图像，边界上有问题，但中心区域还是好的，右边是子网格（位置大概处于母网格中心），子网格的图像明显错误，但我又不知道是什么问题，你能帮我再看看吗？谢谢！下面附上我的脚本：
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"

begin
  in  = addfile("/home/zssapr/WRF3.5.1/20110710/data/wrf_data/wrfout_d02_2011-07-10_00:00:00","r")                             ; open netcdf file

  U = in->U
  ua = wrf_user_unstagger(U,U@stagger)
  V = in->V
  va = wrf_user_unstagger(V,V@stagger)



  scale = 1.e05
  vrt   = ua                                      ; retain coordinates
  vrt   = uv2vrG(ua,va)*scale
  vrt@long_name = "vorticity"
  vrt@units     = "scaled"


  wks  = gsn_open_wks("ps","vort4")               ; specifies a ps plot
  gsn_define_colormap(wks,"BlWhRe")              ; choose colormap

  res                   = True                   ; plot mods desired
  res@cnFillOn          = True                   ; turn on color
  res@gsnSpreadColors   = True                   ; use full range of color map
  res@lbLabelAutoStride = True                   ; nice label bar labels
  res@cnLinesOn         = True                  ; no contour lines

  res@cnLevelSelectionMode = "ManualLevels"       ; set manual contour levels
  res@cnMinLaevelValF       = -4                   ; set min contour level
  res@cnMaxLevelValF       =  4                   ; set max contour level
  res@cnLevelSpacingF      =  0.5                 ; set contour spacing



  plot = gsn_csm_contour(wks,vrt(20,9,:,:),res) ; create plot
*******************************************
end

追风的阳光

还是调下区域，增大scale的值。。。

longlivehj

追风的阳光 发表于 2014-4-23 11:17
你好，我用官网上的脚本试着做出了涡度的图像，上图左边是母网格的图像，边界上有问题，但中心区域还是 ...

用下面两个函数直接求
wrf_avo
wrf_pvo