Python isdepthup Examples

Example #1

File: courses_grids.py Project: jfleroux/vacumm

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

Example #2

File: showvar.py Project: jfleroux/vacumm

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

Example #3

File: courses_grids.py Project: acoat/vacumm

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

Example #4

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