- 积分
- 557
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-11-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
ncl是不是对写入的数据大小有限制,当删掉几个变量就可以写入,或者写进两个nc文件里就不报错,可是现在要写进一个文件,求大神帮助。
错误提示为:
ncendef: ncid 131072: NetCDF: One or more variable sizes violate format constraints
fatal:NetCDF: Operation not allowed in define mode: error attempting to write variable (tc) to file
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 196 in file tiqu3.ncl
ncclose: ncid 131072: NetCDF: One or more variable sizes violate format constraints
脚本如下:就是简单的将wrfout数据读出来再写入另外一个文件。
t=0
; do nt=3,ntimes-1,3
do nt=3,4
qc(t,:,:,:) = wrf_user_getvar(aa,"QCLOUD",nt)
qr(t,:,:,:) = wrf_user_getvar(aa,"QRAIN",nt)
qi(t,:,:,:) = wrf_user_getvar(aa,"QICE",nt)
qs(t,:,:,:) = wrf_user_getvar(aa,"QSNOW",nt)
xlat(t,:,:)= wrf_user_getvar(aa,"XLAT",nt)
xlong(t,:,:)= wrf_user_getvar(aa,"XLONG",nt)
psfc(t,:,:) = wrf_user_getvar(aa,"PSFC",nt)*0.01
tg(t,:,:)= wrf_user_getvar(aa,"TSK",nt)
tc(t,:,:,:) = wrf_user_getvar(aa,"tc",nt)
u(t,:,:,:) = wrf_user_getvar(aa,"ua",nt)
v(t,:,:,:) = wrf_user_getvar(aa,"va",nt)
ter(t,:,:) = wrf_user_getvar(aa,"ter",nt)
qv(t,:,:,:) = wrf_user_getvar(aa,"QVAPOR",nt)
h(t,:,:,:) = wrf_user_getvar(aa,"height",nt)
if(i .ne. 4 )then
qg(t,:,:,:) = wrf_user_getvar(aa,"QGRAUP",nt)
qt(t,:,:,:) = qc(t,:,:,:)+qr(t,:,:,:)+qi(t,:,:,:)+qs(t,:,:,:)+qg(t,:,:,:)
else
qt(t,:,:,:) = qc(t,:,:,:)+qr(t,:,:,:)+qi(t,:,:,:)+qs(t,:,:,:)
end if
t= t+1
end do
fout->qc=qc
fout->qr=qr
fout->qi=qi
fout->qs=qs
fout->xlat=xlat
fout->xlong=xlong
if(i .ne. 4 )then
fout->qg=qg
end if
fout->tc=tc
fout->ter=ter
fout->u=u
fout->v=v
fout->qv=qv
fout->h=h
fout->psfc=psfc
fout->tg=tg
fout->qt=qt
print("success: test"+a(i))
end do
|
|