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