ncl做散度

yang69can

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

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

风往北吹

我觉得你画的两张图都是一样的，但是显示的范围有差异，仔细看右上角的高低值平行的“条”就知道了。

yang69can

风往北吹 发表于 2014-12-3 09:28
我觉得你画的两张图都是一样的，但是显示的范围有差异，仔细看右上角的高低值平行的“条”就知道了。

但是唯一变的是脚本中的紫色的一两句计算，其他的都没有变呀！范围为什么会变呢？

风往北吹

yang69can 发表于 2014-12-3 09:46
但是唯一变的是脚本中的紫色的一两句计算，其他的都没有变呀！范围为什么会变呢？

第一种方案应该是没错的，可能是在赋值传递过程中坐标信息丢失了吧。

yang69can

风往北吹 发表于 2014-12-3 10:02
第一种方案应该是没错的，可能是在赋值传递过程中坐标信息丢失了吧。

但是为什么第二种没有丢失，第一种丢失了呢？从数学上看是同样的计算思路呀！所以这里我不太懂，不知道您可以帮我解决下这个问题

dataiyang

yang69can 发表于 2014-12-3 10:05
但是为什么第二种没有丢失，第一种丢失了呢？从数学上看是同样的计算思路呀！所以这里我不太懂，不知道您 ...

ncl中等号的功能跟一般的还是有一些差异的，根据之前的解释，给第一种方案重新在赋一个经纬度属性，如果跟第二种是一样的结果的话应该就是坐标信息不对了！

dataiyang

yang69can 发表于 2014-12-3 10:05
但是为什么第二种没有丢失，第一种丢失了呢？从数学上看是同样的计算思路呀！所以这里我不太懂，不知道您 ...

ncl中等号的功能跟一般的还是有一些差异的，根据之前的解释，给第一种方案重新在赋一个经纬度属性，如果跟第二种是一样的结果的话应该就是坐标信息不对了！

centaurus

虽然不知道为什么会这样，看着就像镜子的两面，图沿着45°纬度倒过来了。不过呈东北-西南走向的第一个图应该是对的吧。希望有高人能解答

吴闫521

ncl作图真心漂亮，，，可惜不会啊