- 积分
- 2228
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-26
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2012-7-29 08:04:20
|
显示全部楼层
hzwjy 发表于 2012-7-28 15:02
dims = dimsizes(A)
ntimes = 30
因为我是增加高度那个维,因此我是这样写的:
dims = dimsizes(x1)
ntim = dims(0)
nlat = dims(1)
nlon = dims(2)
scale = new ((/ntim,4,nlat,nlon/) , typeof(x1))
scale!0 = "time"
scale!1 = "lev"
scale!2 = "lat"
scale!3 = "lon"
scale&lev = (/ 1000, 100, 10, 1 /)
但是这个生产的数据用grads打不开,知道是不符合grads标准的写法,但是在这个文件中赋了一些属性,总是有错误,因此我又根据兰溪的帖子写了个新的ncl脚本,如下:
begin
f = addfile("./scale.nc","w")
f->lev@units = "generic"
f->lev@units = "level"
f->lev@axis = "z"
f->lat@units = "degrees_north"
f->lon@units = "degrees_east"
end
这样两个先后运行,就可以了,现在请教如果想生成的第一个nc文件就符合标准,应该怎么改呢? |
|