- 积分
- 1564
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-3-23
- 最后登录
- 1970-1-1
|
发表于 2022-3-7 20:02:01
|
显示全部楼层
begin
diri = "/mnt/e/data/hgt_daily/"
files = systemfunc("ls " + diri + "hgt.*.nc" ) ;"文件名相同的部分,不一样的地方用*代替")
f = addfiles(files, "r")
ListSetType(f, "cat")
;var = new((/10,14610,37,73/), float)
;hh_100 = f[:]->hgt(:,{100},{0:90},{0:180});要合并的数据变量
;hh_150 = f[:]->hgt(:,{150},{0:90},{0:180});要合并的数据变量
hh_200 = f[:]->hgt(:,{250},:,:);要合并的数据变量
;hh_250 = f[:]->hgt(:,{250},:,:);要合并的数据变量
;hh_300 = f[:]->hgt(:,{300},{0:90},{0:180});要合并的数据变量
;hh_400 = f[:]->hgt(:,{400},{0:90},{0:180});要合并的数据变量
;hh_500 = f[:]->hgt(:,{500},{0:90},{0:180});要合并的数据变量
;hh_600 = f[:]->hgt(:,{600},{0:90},{0:180});要合并的数据变量
;hh_700 = f[:]->hgt(:,{700},{0:90},{0:180});要合并的数据变量
;hh_850 = f[:]->hgt(:,{850},{0:90},{0:180});要合并的数据变量
; hh_850 = f[:]->hgt(:,{850},:,:);要合并的数据变量
;hh = f[:]->$var(i)$(:,:,:,:)
;print(files)
;================================================================数据筛选
;================================================================数据输出
fout = addfile("/mnt/e/h250.nc", "c") ; create 文件用以输出
fout@title = "250hpa_hgt" ; 给文件加上一个全域属性
fout->hgt = hh_200 ; 将 hh 存储到 fout 里 ,变量名命名为 h
;var = new((/14610,10,37,73/), float)
;var(:,0,:,:) = hh_100
;var(:,1,:,:) = hh_150
;var(:,2,:,:) = hh_200
;var(:,3,:,:) = hh_250
;var(:,4,:,:) = hh_300
;var(:,5,:,:) = hh_400
;var(:,6,:,:) = hh_500
;var(:,7,:,:) = hh_600
;var(:,8,:,:) = hh_700
;var(:,9,:,:) = hh_850
;fname = "/mnt/d/hgt_250.nc"
;system("rm -f "+fname)
;fout := addfile(fname,"c")
;fout->hgt=hh_250
end
|
|