一些常用的grads命令与ncl命令对应关系

majianyong

lz,请教个问题：
我有一个nc文件，变量维数是[time,z,lat,lon]，分别是366天，1个高度，360个纬度，720个经度，现在我想用write_matrix语句输出所有时次的数据(我知道这个命令只能输出2维，所以想写个循环输出)，循环写成如下格式：

a=addfile("2008_prcp.nc","r")
prcp=a->prcp
time=a->time
nlon=dimsizes(longitude)
opt=True
fmtf= nlon+"f8.2"
opt@fout="2008_prcp_all.dat"

do n=0,dimsizes(time)-1
write_matrix (prcp(n,0,:,:), fmtf, opt)
end do
-----------------------------为什么最后只输出第366天的值，要怎么改这个循环才能把所有的天数数据都输出呢？？求解答，很感谢！

游子

majianyong 发表于 2013-11-8 13:54

                               
登录/注册后可看大图

lz,请教个问题：
我有一个nc文件，变量维数是[time,z,lat,lon]，分别是366天，1个高度，360个纬度，720个经 ...

这是因为write-matrix只支持矩阵输出，只能输出二维场，你是用循环时，后一个将覆盖前一个。循环几次就覆盖几次。所以最后只能输出一个时次、

游子

游子 发表于 2013-11-8 14:54

                               
登录/注册后可看大图

这是因为write-matrix只支持矩阵输出，只能输出二维场，你是用循环时，后一个将覆盖前一个。循环几次就覆 ...

解决方法如下：
a=addfile("2008_prcp.nc","r")
prcp=a->prcp
time=a->time
nlon=dimsizes(longitude)
opt=True
fmtf= nlon+"f8.2"
opt@fout="2008_prcp_all.dat"

dim = dimsizes(prcp)

prcp_2D = onedtond(ndtooned(prcp(:,0,:,:)), (/dim(0)*dim(1), dim(2)/)
write_matrix (prcp_2D, fmtf, opt)


这就是“降维”的思想。。。

游子

majianyong 发表于 2013-11-8 13:54

                               
登录/注册后可看大图

lz,请教个问题：
我有一个nc文件，变量维数是[time,z,lat,lon]，分别是366天，1个高度，360个纬度，720个经 ...

dim = dimsize(prcp(:,0,:,:))

majianyong

游子 发表于 2013-11-8 15:01

                               
登录/注册后可看大图

dim = dimsize(prcp(:,0,:,:))

多谢lz提点！dim = dimsize(prcp(:,0,:,:))应该是dim = dimsizes(prcp(:,0,:,:))吧！
还有你写的这句(/dim(0)*dim(1), dim(2)/)相当于/366天×360lat,720lon/是吧？！这样利用marix输出的是366天的360×720个格点么？我感觉是不是要写成(/dim(0),  dim(1)*dim(2)/)相当于/366天 , 360lat×720lon/比较合理？因为在NCL不是列优先么？

游子

majianyong 发表于 2013-11-9 18:29

                               
登录/注册后可看大图

多谢lz提点！dim = dimsize(prcp(:,0,:,:))应该是dim = dimsizes(prcp(:,0,:,:))吧！
还有你写的这句(/d ...

OK，你理解了就好。挺简单的。

majianyong

游子 发表于 2013-11-9 18:39

                               
登录/注册后可看大图

OK，你理解了就好。挺简单的。

抱头痛哭。。。刚学用起来还是很生疏，困难重重！

游子

万事开头难，熟能生巧。

abc258886

好啊，顶啊！

hd4908

简单的程序  不过顶一下