ncl做散度

yang69can · 发表于 2014-12-3 08:21:24

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 yang69can 于 2014-12-3 08:22 编辑

我想问下各位，我用同样的脚本，只是中间处理数据有一句不一样，为什么画出图差距如此之大？（输出来的变量数据值时一样的）
脚本一、

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"
begin
      a = addfile("./grib2006071600.grib1","r")
      u = a->U_GRD_3_ISBL({200},:,:)
      v = a->V_GRD_3_ISBL({200},:,:)

      duv = uv2dvG_Wrap(u,v)/4


      scale = 1.e05
      duv = duv*scale

      wks = gsn_open_wks("png","divergerce200-2")
      res = True
      res@gsnDraw = False
      res@gsnFrame = False
      res@gsnPaperOrientation = "Landscape"

      res@tmYROn = False
      res@tmXTOn = False
      res@tmBorderThicknessF = 8.0

      mpres = res
      mpres@mpFillOn = False
      mpres@mpOutlineOn = True
      mpres@mpOutlineSpecifiers = "China:states"
      mpres@mpDataBaseVersion = "MediumRes"
      mpres@mpDataSetName = "Earth..4"
      mpres@mpGeophysicalLineThicknessF = 2
      mpres@mpNationalLineThicknessF = 2
      mpres@mpUSStateLineThicknessF = 2

      mpres@mpMinLatF = 15
      mpres@mpMaxLatF = 55
      mpres@mpMinLonF = 75
      mpres@mpMaxLonF = 135

      cnres = res
      cnres@cnFillOn = True
      cnres@gsnSpreadColors = True
      cnres@lbLabelAutoStride = True
      cnres@lbOrientation = "Vertical"
      cnres@lbLabelFontHeightF = 0.01
      cnres@lbPerimThicknessF = 0.01
      cnres@lbBoxLinesOn = False
      cnres@pmLabelBarWidthF = 0.05
      cnres@cnLinesOn = False


      cnres@cnLevelSelectionMode = "ManualLevels"
      cnres@cnMinLevelValF = -8
      cnres@cnMaxLevelValF = 8
      cnres@cnLevelSpacingF = 0.5

      cnres@gsnLeftString = "divergence"+"  " +"height:10m"
      cnres@gsnRightString = "1/s"+"  "+1.0/scale
      cnres@gsnCenterString = "06-07-19-06"

      map = gsn_csm_map(wks,mpres)
      plot = gsn_csm_contour(wks,duv,cnres)

      overlay(map,plot)

      draw(map)
frame(wks)
end

脚本二：
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"
begin
      a = addfile("./grib2006071600.grib1","r")
      u = a->U_GRD_3_ISBL({200},:,:)
      v = a->V_GRD_3_ISBL({200},:,:)


   duv = uv2dvG_Wrap(u,v)
      duv = duv/4

      scale = 1.e05
      duv = duv*scale

      wks = gsn_open_wks("png","divergerce200-2")
      res = True
      res@gsnDraw = False
      res@gsnFrame = False
      res@gsnPaperOrientation = "Landscape"

      res@tmYROn = False
      res@tmXTOn = False
      res@tmBorderThicknessF = 8.0

      mpres = res
      mpres@mpFillOn = False
      mpres@mpOutlineOn = True
      mpres@mpOutlineSpecifiers = "China:states"
      mpres@mpDataBaseVersion = "MediumRes"
      mpres@mpDataSetName = "Earth..4"
      mpres@mpGeophysicalLineThicknessF = 2
      mpres@mpNationalLineThicknessF = 2
      mpres@mpUSStateLineThicknessF = 2

      mpres@mpMinLatF = 15
      mpres@mpMaxLatF = 55
      mpres@mpMinLonF = 75
      mpres@mpMaxLonF = 135

      cnres = res
      cnres@cnFillOn = True
      cnres@gsnSpreadColors = True
      cnres@lbLabelAutoStride = True
      cnres@lbOrientation = "Vertical"
      cnres@lbLabelFontHeightF = 0.01
      cnres@lbPerimThicknessF = 0.01
      cnres@lbBoxLinesOn = False
      cnres@pmLabelBarWidthF = 0.05
      cnres@cnLinesOn = False


      cnres@cnLevelSelectionMode = "ManualLevels"
      cnres@cnMinLevelValF = -8
      cnres@cnMaxLevelValF = 8
      cnres@cnLevelSpacingF = 0.5

      cnres@gsnLeftString = "divergence"+"  " +"height:10m"
      cnres@gsnRightString = "1/s"+"  "+1.0/scale
      cnres@gsnCenterString = "06-07-19-06"

      map = gsn_csm_map(wks,mpres)
      plot = gsn_csm_contour(wks,duv,cnres)

      overlay(map,plot)

      draw(map)
frame(wks)
end