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)
# 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)
# 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
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)
# 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
# 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