dep = create_dep((0., 10)) # -> AFFICHEZ LES INFOS xx, yy = N.meshgrid(N.arange(5.), N.arange(4.)) lon2d, lat2d = create_axes2d(xx, yy) ii = lon2d.getAxis(1) # - grille grid = create_grid(lon, lat) # -> ESSAYEZ AVEC LON EXPLICITE gridc = create_grid(lon2d, lat2d) # Verifier print islon(lon) print isgrid(grid) # -> TEST PARAM CURV=... print isrect(gridc) # -> CREEZ GRILLE NON RECT ET RETESTER print isdepthup(dep) # -> TESTEZ EN CHANGEANT ATTRIBUT POSITIVE ET VALEURS print isregular(lon) # Affecter var = MV2.ones(grid.shape) set_grid(var, grid) varc = MV2.ones(gridc.shape) set_grid(varc, gridc) # -> VERIFIEZ ID DES AXES # Récupérer mygrid = get_grid(gridc) # -> TESTEZ AVEC (LON,LAT) ET PARAMS STRICT ET INTERCEPT mylon2d = get_axis(gridc, 1) # -> COMPAREZ AVEC .GETAXIS()
# Read variable vv = [f(vname, **kwread) for vname in args.var] if args.dvar: depth = f(args.dvar, **kwread) f.close() if args.verbose: for v in vv: v.info() order = vv[0].getOrder() ndim = vv[0].ndim # Reduce rank when needed if ndim>2: for i in xrange(ndim-2): if order[i]=='z': islice = -1 if isdepthup(var[0].getAxis(0)) else 0 else: islice = 0 vv = [v[islice] for v in vv] order = vv[0].getOrder() ndim = vv[0].ndim if args.dvar and depth.ndim>ndim: for i in xrange(depth.ndim-ndim): depth = depth[0] # Operators if args.add: for v in vv: v[:] += args.add if args.mul: for v in vv:
lat = create_lat(N.arange(43, 50.)) dep = create_dep((0., 10)) # -> AFFICHEZ LES INFOS xx, yy = N.meshgrid(N.arange(5.), N.arange(4.)) lon2d, lat2d = create_axes2d(xx, yy) ii = lon2d.getAxis(1) # - grille grid = create_grid(lon, lat) # -> ESSAYEZ AVEC LON EXPLICITE gridc = create_grid(lon2d, lat2d) # Verifier print islon(lon) print isgrid(grid) # -> TEST PARAM CURV=... print isrect(gridc) # -> CREEZ GRILLE NON RECT ET RETESTER print isdepthup(dep) # -> TESTEZ EN CHANGEANT ATTRIBUT POSITIVE ET VALEURS print isregular(lon) # Affecter var = MV2.ones(grid.shape) set_grid(var, grid) varc = MV2.ones(gridc.shape) set_grid(varc, gridc) # -> VERIFIEZ ID DES AXES # Récupérer mygrid = get_grid( gridc) # -> TESTEZ AVEC (LON,LAT) ET PARAMS STRICT ET INTERCEPT mylon2d = get_axis(gridc, 1) # -> COMPAREZ AVEC .GETAXIS() # Sélectionner print coord2slice(lon, lon=(2, 4.5)) # -> COMPAREZ AVEC .GETINTERVALEXT(...)
# Read variable vv = [f(vname, **kwread) for vname in args.var] if args.dvar: depth = f(args.dvar, **kwread) f.close() if args.verbose: for v in vv: v.info() order = vv[0].getOrder() ndim = vv[0].ndim # Reduce rank when needed if ndim > 2: for i in xrange(ndim - 2): if order[i] == 'z': islice = -1 if isdepthup(var[0].getAxis(0)) else 0 else: islice = 0 vv = [v[islice] for v in vv] order = vv[0].getOrder() ndim = vv[0].ndim if args.dvar and depth.ndim > ndim: for i in xrange(depth.ndim - ndim): depth = depth[0] # Operators if args.add: for v in vv: v[:] += args.add if args.mul: for v in vv: