- 积分
- 10617
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-10-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
原始的nc资料是在欧洲中心下的,三个变量下在一个nc里了,查看过原始资料里面的变量名分别是u,v,t,现在想计算里面的额v变量,但是总是读不对,这个脚本以前读只含一个变量的nc时没有问题,不知道是不是多个变量的读法不同,有人可以帮忙解答吗。
脚本如下
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
in1 =addfile("/home/hp/uvt.nc","r")
time = in1->time
lev = in1->level
lon = in1->longitude
lat = in1->latitude
hgt = short2flt(in1->v(2:133,:,:,:))
tt = cd_calendar(hgt&time,0)
print(hgt)
h = hgt(0:32,:,:,:)
h = 0.0
fud = addfile ("vcha.nc","c")
;-------------------------------------------------------------------------------------------------
h = new((/33,17,181,360/),float)
do k=0,32
h(k,:,:,:)=(hgt(k,:,:,:)+hgt(k+1,:,:,:)+hgt(k+2,:,:,:)+hgt(k+3,:,:,:))/4.0
end do
eld = new((/17,181,360/),float)
elg = new((/17,181,360/),float)
elq = new((/17,181,360/),float)
eld(:,:,:)=(h(1,:,:,:)+h(9,:,:,:)+h(10,:,:,:)+h(11,:,:,:)+h(12,:,:,:)+h(14,:,:,:)+h(16,:,:,:)+h(24,:,:,:)+h(26,:,:,:)+h(28,:,:,:)+h(31,:,:,:)+h(32,:,:,:))/12.0
elg(:,:,:)=(h(0,:,:,:)+h(5,:,:,:)+h(7,:,:,:)+h(13,:,:,:)+h(15,:,:,:)+h(17,:,:,:)+h(18,:,:,:)+h(20,:,:,:)+h(21,:,:,:)+h(22,:,:,:)+h(25,:,:,:))/11.0
elq(:,:,:)=eld(:,:,:)-elg(:,:,:)
copy_VarCoords(hgt(0,:,:,:),elq)
fud->V = elq
end
|
|