用GRADS打开NC文件

map999

ncfile = NetCDFFile('e:/data/ncep/hgt.mon.mean.nc')
hgt = ncfile.variables['hgt'][:,5,:37,:]  #北半球500hpa
lats = ncfile.variables['lat'][:37]
lons = ncfile.variables['lon'][:]
ts= ncfile.variables['time'][:]
lev = ncfile.variables['level'][:]
ncfile.close()
ny,nx = hgt.shape[1:]
nt = 2010-1948
h_summer = np.empty( ( nt, ny, nx ), 'f' )
for yr in range( 1948, 2010 ):
    t = (yr-1948)*12
    mean = hgt[ t+5:t+8,:,:].mean( axis=0 )
    h_summer[yr-1948,:,:] = mean
ncfile = NetCDFFile( 'c:/x.nc', mode='w' )
ncfile.createDimension( 'lon', nx )
ncfile.createDimension( 'lat', ny )
ncfile.createDimension( 'time',nt )
x = ncfile.createVariable( 'lon', 'f', ('lon',) )
y = ncfile.createVariable( 'lat', 'f', ('lat',) )
t = ncfile.createVariable( 'time','i', ('time',))
v = ncfile.createVariable( 'hgt', 'f', ('time','lat','lon') )
#---units需要正确设置，否则grads sdfopen就会出现SDF file has no discernable X coordinate的错误
y.units = 'degrees_north'
x.units = 'degrees_east'
t.units = 'years since 1-1-1 0:0:0'
v.units = 'm'
x[:] = lons[:]
y[:] = lats[:]
t[:] = np.arange( 1948,2010,dtype='i' )
v[:,:,:] = h_summer[:,:,:]
ncfile.close()

mofangbao

没明白啥意思啊、、、、你确定这是GrADS的语法，还是发个帖子拿积分？

whmooncool

确定是grads吗？

MeteoInfo

典型的Python语法呀

maplesnow82

MeteoInfo 发表于 2012-3-12 20:29

                               
登录/注册后可看大图

典型的Python语法呀

对啊，不是grads吧？

xjliutao

这个可以输出吗？

wwwaaasss

{:5_214:}{:5_214:}

zdqpla

好东西啊，谢谢分享！

李振宇

mofangbao 发表于 2012-3-12 18:53
没明白啥意思啊、、、、你确定这是GrADS的语法，还是发个帖子拿积分？

完全是错误的语法

dgbest

我差点复制粘贴了