Python ZonalMeansの例

コード例 #1

bg= vcs.test.support.bg
x=vcs.init()
# Open data file
f=cdms2.open('test_data.nc')
O2=f("O2",level=slice(0,5))
print O2.shape,O2.mask
print 'Level:',O2.getLevel()
f.close()

# Open grid file
f=cdms2.open('test_grid.nc')

# 1st with just bounds
blon = f('bounds_lon')
blat = f('bounds_lat')
mask = f('mask')
g=mask.getGrid()
O2.setAxis(-1,mask.getAxis(-1))
O2.setAxis(-2,mask.getAxis(-2))
O2.setGrid(g)

#sys.exit()
area = f('area')

res=ZonalMeans.compute(O2,area=area,delta_band=5)

print res[1].shape
x.plot(res[1][:,0],bg=bg)

vcs.test.support.check_plot(x)

コード例 #2

ファイル: test_zonal.py プロジェクト: NCPP/uvcdat-devel

# Adapted for numpy/ma/cdms2 by convertcdms.py
import cdms2,sys,ZonalMeans
cdms2.setAutoBounds("on")
cdms2.axis.level_aliases.append('depth')
import vcs
import vcs.test.support
bg= vcs.test.support.bg
x=vcs.init()
# Open data file
f=cdms2.open('test_data.nc')
O2=f("O2",level=slice(0,5))
print O2.shape,O2.mask
print 'Level:',O2.getLevel()
f.close()
# Open grid file
f=cdms2.open('test_grid.nc')

# 1st with just bounds
blon = f('bounds_lon')
blat = f('bounds_lat')
#mask = 1 - f('mask')
area = f('area')

res=ZonalMeans.compute(O2,bounds_lon=blon,bounds_lat=blat,area=area,delta_band=5)

print res[1].shape
x.plot(res[1][:,0],bg=bg)
vcs.test.support.check_plot(x)

コード例 #3

ファイル: test_zonal_4.py プロジェクト: ptbremer/uvcdat-devel

# Adapted for numpy/ma/cdms2 by convertcdms.py
import cdms2
import sys, MV2
import ZonalMeans

delta_lon = .5
delta_lat = .5
nlat = int(180. / delta_lat)
nlon = int(360 / delta_lon)
s = MV2.ones((nlat, nlon), typecode='f')
print delta_lon, nlon, delta_lat, nlat
g = cdms2.createUniformGrid(-89.5, nlat, delta_lat, 0., nlon, delta_lon)
s.setAxis(-1, g.getLongitude())
s.setAxis(-2, g.getLatitude())
s.setGrid(g)
print s.shape
print s.getGrid()
print s.getLongitude()
print s.getLatitude()
zm = ZonalMeans.compute(s, delta_band=5)
print zm

コード例 #4

ファイル: test_zonal.py プロジェクト: ptbremer/uvcdat-devel

cdms2.setAutoBounds("on")
cdms2.axis.level_aliases.append('depth')
import vcs
import vcs.test.support
bg = vcs.test.support.bg
x = vcs.init()
# Open data file
f = cdms2.open('test_data.nc')
O2 = f("O2", level=slice(0, 5))
print O2.shape, O2.mask
print 'Level:', O2.getLevel()
f.close()
# Open grid file
f = cdms2.open('test_grid.nc')

# 1st with just bounds
blon = f('bounds_lon')
blat = f('bounds_lat')
#mask = 1 - f('mask')
area = f('area')

res = ZonalMeans.compute(O2,
                         bounds_lon=blon,
                         bounds_lat=blat,
                         area=area,
                         delta_band=5)

print res[1].shape
x.plot(res[1][:, 0], bg=bg)
vcs.test.support.check_plot(x)

コード例 #5

ファイル: test_zonal_2.py プロジェクト: ptbremer/uvcdat-devel

# Adapted for numpy/ma/cdms2 by convertcdms.py
import cdms2,ZonalMeans,sys,os


## 
f=cdms2.open(os.path.join(cdms2.__path__[0],'..','..','..','..','sample_data','clt.nc'))
s=f('u')
print s.shape
res = ZonalMeans.compute(s)
import vcs
import vcs.test.support
bg= vcs.test.support.bg
x=vcs.init()
x.plot(res[1],bg=bg)
vcs.test.support.check_plot(x)
lats=[-90,-80,-70,-60,-50,-40,-30,-20,-10,-5,0,5,10,20,30,40,50,60,70,80,90]
res = ZonalMeans.compute(s,delta_band=lats)
## print res

コード例 #6

ファイル: test_zonal_4.py プロジェクト: NCPP/uvcdat-devel

# Adapted for numpy/ma/cdms2 by convertcdms.py
import cdms2
import sys,MV2
import ZonalMeans

delta_lon=.5
delta_lat=.5
nlat=int(180./delta_lat)
nlon=int(360/delta_lon)
s=MV2.ones((nlat,nlon),typecode='f')
print delta_lon,nlon,delta_lat,nlat
g=cdms2.createUniformGrid(-89.5, nlat, delta_lat, 0., nlon, delta_lon)
s.setAxis(-1,g.getLongitude())
s.setAxis(-2,g.getLatitude())
s.setGrid(g)
print s.shape
print s.getGrid()
print s.getLongitude()
print s.getLatitude()
zm=ZonalMeans.compute(s,delta_band=5)
print zm