请问cdo如何输出nc文件中的一个变量到新nc文件中

chongzika

張顧.煒Zane

ncl 可以
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

        system ("rm ud.nc vd.nc")

        in1 =addfile("/scratch/student/bozk/ex/vwnd.mon.mean.nc","r")
        in2 =addfile("/scratch/student/bozk/ex/uwnd.mon.mean.nc","r")
               time  = in1->time
               lev   = in1->level
               lon   = in1->lon
               lat   = in1->lat

        v  = short2flt(in1->vwnd(:,:,:,:))
        u  = short2flt(in2->uwnd(:,:,:,:))

          fud   = addfile ("ud.nc","c")
          fvd   = addfile ("vd.nc","c")

        div = uv2dvG(u,v)

          ud    = new ( dimsizes(u), typeof(u), u@_FillValue )
        vd    = new ( dimsizes(v), typeof(v), v@_FillValue )

          dv2uvg(div,ud,vd) ; div  ==> divergent  wind components

          copy_VarCoords(u, ud )
          copy_VarCoords(v, vd )
          ud@long_name  = "Zonal Divergent Wind"
          ud@units      = u@units
          vd@long_name  = "Meridional Divergent Wind"
          vd@units      = v@units

          fud->u = ud
          fvd->v = vd
end

chongzika

張顧.煒Zane 发表于 2015-1-23 22:18
ncl 可以
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/ncls ...

ok非常感谢