谁有合并nc文件的ncl？

zongqing123

{:eb502:}{:eb502:}

zongqing123

cdo 不太方便

FrankieLJY

likailing2008 发表于 2013-10-13 19:02
用cdo吧，这个东西处理数据什么的最方便了

你好我用cdo的mergetime合并文件，但是一直提醒
Warning (cdf_set_var) : Inconsistent variable definition for XLAT_V!
Warning (cdf_set_var) : Inconsistent variable definition for XLONG_V!
不知道层主遇到过这个问题没有

likailing2008

FrankieLJY 发表于 2017-3-29 10:02
你好我用cdo的mergetime合并文件，但是一直提醒
Warning (cdf_set_var) : Inconsistent variable defini ...

没遇到过这种情况，会不会是你两个数据的经纬度坐标使用的不一样？如果是建议先插值统一一下坐标

FrankieLJY

likailing2008 发表于 2017-3-30 08:36
没遇到过这种情况，会不会是你两个数据的经纬度坐标使用的不一样？如果是建议先插值统一一下坐标

你好，我用的应该都是lambert投影，如果要转换投影的话，请问有什么方法呢

likailing2008

FrankieLJY 发表于 2017-4-4 20:05
你好，我用的应该都是lambert投影，如果要转换投影的话，请问有什么方法呢

太久没研究了，有点忘了。但是不知你用NCL不？要是用我这里有一个插值程序，在CDO不好用的时候的别的方法。

FrankieLJY

likailing2008 发表于 2017-4-5 08:18
太久没研究了，有点忘了。但是不知你用NCL不？要是用我这里有一个插值程序，在CDO不好用的时候的别的方法 ...

你好 请问你的程序可以分享一下吗 谢谢

likailing2008

FrankieLJY 发表于 2017-4-5 10:15
你好 请问你的程序可以分享一下吗 谢谢

; ***********************************************
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
    FILES = systemfunc (" ls -1 "+"psl*.nc ")        ;从这里开始到下面都是读取文件的命令
    numFILES = dimsizes(FILES)

        do k = 0,numFILES-1  

  
         b = addfile("bcc-csm1-1.nc","r")         ;把所有读入的文件写入一个叫做a的文件中       
         lon2 := b->lon
         lat2 := b->lat
         f = addfile(FILES(k),"r")
                 aa=strlen(FILES(k))
         bb= str_get_cols(FILES(k),0,aa-4)
         lon := f->lon
         lat := f->lat
         psl := f->psl
         rcmpsl := linint2 (lon,lat,psl(:,:,:), True, lon2,lat2, 0)

  
        rcmpsl!0="time"
        rcmpsl!1="lat"
        rcmpsl!2="lon"

      rcmpsl&lat=lat2
      rcmpsl&lon=lon2
          rcmpsl&time=psl&time

  fout   = addfile(bb+"_remapbil.nc","c")
  fout->psl=rcmpsl

  end do

end

不懂的请微博私信我