Exemple #1
0
 def gridL(self,dlalo=0.5,nij=10):
     """provide grid and rec values for other tests"""
     grtyp='L'
     grref=grtyp
     la0 = 0.-dlalo*(nij/2.)
     lo0 = 180.-dlalo*(nij/2.)
     ig14 = (ig1,ig2,ig3,ig4) =  rpnstd.cxgaig(grtyp,la0,lo0,dlalo,dlalo)
     axes = (None,None)
     hasAxes = 0
     ij0 = (1,1)
     doCorners = 0
     (la,lo) = Fstdc.ezgetlalo((nij,nij),grtyp,(grref,ig1,ig2,ig3,ig4),axes,hasAxes,ij0,doCorners)
     grid = rpnstd.RPNGrid(grtyp=grtyp,ig14=ig14,shape=(nij,nij))
     return (grid,la,lo)
Exemple #2
0
 def gridZL(self,dlalo=0.5,nij=10):
     """provide grid and rec values for other tests"""
     (g1,la1,lo1) = self.gridL(dlalo,nij)
     x_axis_d = lo1[:,0].reshape((lo1.shape[0],1)).copy('FORTRAN')
     y_axis_d = la1[0,:].reshape((1,la1.shape[1])).copy('FORTRAN')
     grtyp='L'
     la0 = 0.
     lo0 = 180.
     ig14 = (ig1,ig2,ig3,ig4) =  rpnstd.cxgaig(grtyp,0.,0.,1.,1.)
     ip134 = (1,2,1,1)
     g1.ig14 = ig14
     x_axis = rpnstd.RPNRec(x_axis_d,rpnstd.RPNMeta())
     y_axis = rpnstd.RPNRec(y_axis_d,rpnstd.RPNMeta())
     grid = rpnstd.RPNGrid(grtyp='Z',ig14=ip134,shape=(nij,nij),g_ref=g1,xyaxis=(x_axis,y_axis))
     return (grid,la1,lo1)