linint_Wrap缺测值

zhangkf_7

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/csm/shea_util.ncl"
begin
        f = addfile("/home/zkf/paper/Gt/data/sst_new.nc","r")
        sst  = f->sst
        time = f->time
        lat  = f->latitude({-90:90})
        lon  = f->longitude({-180:180})
        
        sst!0 = "time"
        sst!1 = "lat"
        sst!2 = "lon"
        sst&time = time
        sst&lat  = lat
        sst&lon  = lon
        ;print(sst&lon)
        ;printVarSummary(sst)     
;************************************************
; interpolate to new grid
;***********************************************
    newlat = fspan(-90.,90.,180)
    newlon = fspan(0.,360.,360)
    newlat@units = "degrees_north"
    newlon@units = "degrees_east"
         
    newsst = linint2_Wrap(lon,lat,sst,True,newlon,newlat,0)
   
    newsst!1   ="lat"   
    newsst!2   ="lon"
    newsst&lat = newlat
    newsst&lon = newlon
以上是我的脚本代码
然后图是：

                               
登录/注册后可看大图

只有东经才出现图，且也有缺测值，但西经就全是。
麻烦各路大神请教

zhangkf_7

file:///home/zkf/paper/Gt/plot/ave.ps

packard

我看不到图。但按你的描述，你的原始数据是-180到180的，但 newlon 是0到360的。可以试试都弄成 -180到180的吧。
做完插值后，假如你真的想要 0 到360 那样的。可以看看 函数lonFlip （https://www.ncl.ucar.edu/Document/Functions/Contributed/lonFlip.shtml）能不能帮你。注意这个函数0和360不能同时重复出现，比如 fspan(-180.,179.,360)

zhangkf_7

packard 发表于 2017-7-27 16:41
我看不到图。但按你的描述，你的原始数据是-180到180的，但 newlon 是0到360的。可以试试都弄成 -180到180 ...

嗯嗯，十分感谢，昨天还一直想着有没有这样的函数，结果我是用较笨的方法对它进行转换;
脚本如下：
;数据的转换，传统方式（复杂：将元数据保留）
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/csm/shea_util.ncl"
begin
        f = addfile("/home/zkf/paper/Gt/data/sst_new.nc","r")
        sst  = f->sst
        time = f->time
        lat  = f->latitude
        lon  = f->longitude
        
        sst!0 = "time"
        sst!1 = "lat"
        sst!2 = "lon"
        sst&time = time
        sst&lat  = lat
        sst&lon  = lon
        
;************************************************
; interpolate to new grid
;***********************************************
        mtime = dimsizes(time)
        mlat = dimsizes(lat)
        mlon = dimsizes(lon)
        
        array = new((/mtime,mlat,mlon/),float)
        xlon = fspan(0.5,359.5,360)
        copy_VarMeta(lon,xlon)
        array(:,:,0:(mlon/2-1)) = sst(:,:,(mlon/2):(mlon-1))
        array(:,:,(mlon/2):(mlon-1)) = sst(:,:,0:(mlon/2-1))
        
        array!0 = "TIME"
        array!1 = "LAT"
        array!2 = "LON"
        array&TIME = time
        array&LAT = lat
        array&LON = xlon


        fout = addfile("/home/zkf/paper/Gt/data/sst.nc","c")
        fout->SST = array
               
end
不过这个脚本也让我对数组更加熟悉一点。