立即注册 登录
气象家园 返回首页

山色in的个人空间 http://bbs.06climate.com/?134279 [收藏] [复制] [分享] [RSS]

日志

将多个逐年数据编写成一个nc数据 ncl

已有 71 次阅读2023-5-23 23:35 | NCL, 数据处理

begin
Files=systemfunc("ls evspsbl_Amon_EC-Earth3-Veg-LR_historical_r1i1p1f1_gr_*01-*12.nc")
f=addfiles(Files,"r")
;printVarSummary(f)
evp0=f[0]->evspsbl(:,{24:43},{73:107})
evp1=f[1]->evspsbl(:,{24:43},{73:107})
evp2=f[2]->evspsbl(:,{24:43},{73:107})
evp3=f[3]->evspsbl(:,{24:43},{73:107})
evp4=f[4]->evspsbl(:,{24:43},{73:107})
evp5=f[5]->evspsbl(:,{24:43},{73:107})
evp6=f[6]->evspsbl(:,{24:43},{73:107})
evp7=f[7]->evspsbl(:,{24:43},{73:107})
evp8=f[8]->evspsbl(:,{24:43},{73:107})
evp9=f[9]->evspsbl(:,{24:43},{73:107})
evp10=f[10]->evspsbl(:,{24:43},{73:107})
evp11=f[11]->evspsbl(:,{24:43},{73:107})
evp12=f[12]->evspsbl(:,{24:43},{73:107})
evp13=f[13]->evspsbl(:,{24:43},{73:107})
evp14=f[14]->evspsbl(:,{24:43},{73:107})
evp15=f[15]->evspsbl(:,{24:43},{73:107})
evp16=f[16]->evspsbl(:,{24:43},{73:107})
evp17=f[17]->evspsbl(:,{24:43},{73:107})
evp18=f[18]->evspsbl(:,{24:43},{73:107})
evp19=f[19]->evspsbl(:,{24:43},{73:107})
evp20=f[20]->evspsbl(:,{24:43},{73:107})
evp21=f[21]->evspsbl(:,{24:43},{73:107})
evp22=f[22]->evspsbl(:,{24:43},{73:107})
evp23=f[23]->evspsbl(:,{24:43},{73:107})
evp24=f[24]->evspsbl(:,{24:43},{73:107})
evp25=f[25]->evspsbl(:,{24:43},{73:107})
evp26=f[26]->evspsbl(:,{24:43},{73:107})
evp27=f[27]->evspsbl(:,{24:43},{73:107})
evp28=f[28]->evspsbl(:,{24:43},{73:107})
evp29=f[29]->evspsbl(:,{24:43},{73:107})
evp30=f[30]->evspsbl(:,{24:43},{73:107})
evp31=f[31]->evspsbl(:,{24:43},{73:107})
evp32=f[32]->evspsbl(:,{24:43},{73:107})
evp33=f[33]->evspsbl(:,{24:43},{73:107})
evp34=f[34]->evspsbl(:,{24:43},{73:107})
evp35=f[35]->evspsbl(:,{24:43},{73:107})
fout=addfile("evspsbl_EC-Earth3-Veg-LR_197901-201412.nc","c")
filedimdef (fout, "time", -1, True)

;array_append_record将逐年提取出的数据的时间维组合在一起
e1=array_append_record(evp0,evp1,0)
e2=array_append_record(e1,evp2,0)
e3=array_append_record(e2,evp3,0)
e4=array_append_record(e3,evp4,0)
e5=array_append_record(e4,evp5,0)
e6=array_append_record(e5,evp6,0)
e7=array_append_record(e6,evp7,0)
e8=array_append_record(e7,evp8,0)
e9=array_append_record(e8,evp9,0)
e10=array_append_record(e9,evp10,0)
e11=array_append_record(e10,evp11,0)
e12=array_append_record(e11,evp12,0)
e13=array_append_record(e12,evp13,0)
e14=array_append_record(e13,evp14,0)
e15=array_append_record(e14,evp15,0)
e16=array_append_record(e15,evp16,0)
e17=array_append_record(e16,evp17,0)
e18=array_append_record(e17,evp18,0)
e19=array_append_record(e18,evp19,0)
e20=array_append_record(e19,evp20,0)
e21=array_append_record(e20,evp21,0)
e22=array_append_record(e21,evp22,0)
e23=array_append_record(e22,evp23,0)
e24=array_append_record(e23,evp24,0)
e25=array_append_record(e24,evp25,0)
e26=array_append_record(e25,evp26,0)
e27=array_append_record(e26,evp27,0)
e28=array_append_record(e27,evp28,0)
e29=array_append_record(e28,evp29,0)
e30=array_append_record(e29,evp30,0)
e31=array_append_record(e30,evp31,0)
e32=array_append_record(e31,evp32,0)
e33=array_append_record(e32,evp33,0)
e34=array_append_record(e33,evp34,0)
evp=array_append_record(e34,evp35,0)
fout->evspsbl=evp

end

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 立即注册

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部