如何将多个nc文件中包含的相同维度的二维数据按行或按列进行合并

DCC

现在我想通过利用ncl实现9个nc文件的数据拼接，其中absc532在文件中是一个600*48（高度*时间）的数组，我想做个时序分布的填充图，也就是想将9个文件中的absc532拼接成600*（48*9）的数组；但是我用ListSetType(nc_file,"cat")后发现，数据是在时间维度上进行了拼接，也就是形成（600*9）*48的数组，出来的图也就成了附件里这个样子。所以想请教大家该如何解决这个问题   我的程序文本如下


load "/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

files = systemfunc("ls /mnt/f/Lidar/data/Japan/*")+".nc" ;
nc_file = addfiles(files,"r");
ListSetType(nc_file,"cat");

time     = nc_file[0]->time                   ; 时间
height   = nc_file[0]->height_a               ; 高度
absc532  = nc_file[:]->absc532                ; 532nm后散

absc532 = absc532 * 1000              ;
printVarSummary(absc532);
absc532@units = "(1/km/sr)"           ;
ext532@units = "(1/km)"               ;


;*********************************
;absc532nm  contour
;*********************************
wks = gsn_open_wks("PNG","/mnt/f/Lidar/pictures/absc532_fukuoka_Nov")      ; send graphics to PNG file.
res                = True                 ; plot模式要求
res@tiMainString   = "absc532nm,Fukuoka,2014,11,01"    ;设置标题
res@cnFillOn       = True                  ; turn on color fill
res@cnFillPalette  = "wh-bl-gr-ye-re"      ; choose colormap
res@cnLinesOn      = False                 ;
res@lbOrientation  = "Vertical"

res@cnLevelSelectionMode = "ManualLevels"      ;set manual contour levels
res@cnMinLevelValF       =  0                  ; set min contour
res@cnMaxLevelValF       =  0.03               ; set max contour
res@cnLevelSpacingF      =  0.001              ; set contour space
plot = gsn_csm_contour(wks, absc532, res)     ; Draw a contour plot.
end

戏子丶

丁总？？？？

cang

请问，问题解决了吗？

DCC

cang 发表于 2019-1-8 10:36
请问，问题解决了吗？

最后是用cdo实现的nc文件拼接，然后用ncl画的图 僵硬的操作

凝慕Oo

可以先改变纬度顺序，cat应该是合并最左侧的一维。如原始数据a是600*48，高度*时间，先改成b=a(time|:,height|:)..再合并b。