加载中…NCL--批量读取FNL日资料后reshape
(2020-06-28 19:46:53) | 分类： NCL的小tips |
FNL的grib2资料无时间维，批量读取后会对气压维或纬度维进行整合
可以利用reshape函数，根据文件数重新改写变量维度，增加时间维
后来发现，上面是一种很麻烦的方法。。。
在addfiles后面加一句
ListSetType(file,"join")
即可增加一个数目为文件个数的新维度
code:
DATADir=(/"/public/home/ncep_daily/"+iyear+"/"/)
FILE_1= systemfunc (" ls -1 " + DATADir +
"fnl_*"+itime+"_00")
f_1 =
addfiles(FILE_1+".grib1","r")
% ListSetType(file,"join")
t_a=addfiles_GetVar (f_1,FILE_1+".grib1","TMP_3_ISBL")
;;;reshape data
size_t=dimsizes(t_a)
size_p_3=dimsizes(p_3)
file_num=size_t(0)/size_p_3
t=reshape(t_a,
(/file_num,size_p_3,size_t(1),size_t(2)/))
喜欢
0
赠金笔