- 积分
- 27079
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-5
- 最后登录
- 1970-1-1
|
发表于 2017-4-5 14:22:03
|
显示全部楼层
; ***********************************************
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
不懂的请微博私信我 |
|