Python resol Examples

Example #1

File: courses_grids.py Project: jfleroux/vacumm

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(...)
print coord2slice(grid, lon=(4, 8), lat=(44, 46)) # -> TESTEZ SUR GRIDC
# -> TESTEZ VARSEL


# Transformer
gridcr = curv2rect(gridc)       # -> TESTEZ AVEC VOTRE GRILLE NON RECT 
ax = create_lon([350, 0, 10.])
print monotonic(ax)
print xshift(var, 2)[0]         # -> COMPAREZ A VAR ET VERIFIEZ AXE
# -> TESTEZ XEXTEND
grido = rotate_grid(grid, 30)   # -> TRACEZ LES LONGITUDE (PCOLOR)
grids = shiftgrid(grid, 1, -1)  # -> VERIFIEZ PUIS TESTEZ SUR GRIDC + PASSER DE T À U
gride = extendgrid(grid, iext=(2, 3))   # -> VERIFIER PUIS TESTEZ LE MODE 


# Exploiter
print resol(grid)               # -> TESTEZ EN METRES ET SUR AXE
print depth2dz(dep)
print get_closest(lon2d, lat2d, 2.3, 1.2)


# Utilitaires sut les coordonnées
xx, yy = meshgrid(xx, yy[:, 0])
xxb, yyb = meshcells(xx, yy)      # -> EN 1D?

Example #2

File: courses_regrid.py Project: acoat/vacumm

#!/usr/bin/env python
# -*- coding: utf8 -*-
"""Le regrillage"""


from vcmq import N, cdms2, curve2, map2, regrid1d, regrid2d, data_sample, \
    create_grid, regrid_method, CDATRegridder, resol, add_grid, P


# Lecture
f = cdms2.open(data_sample('mfs.nc'))
select = dict(lat=(42.5, 43.5), lon=(3, 4))
ssti = f('votemper', level=slice(0, 1), squeeze=1, **select)
f.close()
gridi = ssti.getGrid()
dxi, dyi = resol(gridi)


# Nouvelle grille
factor = 0.634                                                # -> CHANGEZ LE FACTEUR
dxo = dxi*factor
dyo = dyi*factor
grido = create_grid((3.13, 3.8, dxo), (42.62, 43.4, dyo))


# Quelle méthode ?
method = regrid_method(gridi, grido)
print method


# Regrillage

Example #3

File: courses_regrid.py Project: jfleroux/vacumm

    create_grid,
    regrid_method,
    CDATRegridder,
    resol,
    add_grid,
    P,
)


# Lecture
f = cdms2.open(data_sample("mfs.nc"))
select = dict(lat=(42.5, 43.5), lon=(3, 4))
ssti = f("votemper", level=slice(0, 1), squeeze=1, **select)
f.close()
gridi = ssti.getGrid()
dxi, dyi = resol(gridi)


# Nouvelle grille
factor = 0.634  # -> CHANGEZ LE FACTEUR
dxo = dxi * factor
dyo = dyi * factor
grido = create_grid((3.13, 3.8, dxo), (42.62, 43.4, dyo))


# Quelle méthode ?
method = regrid_method(gridi, grido)
print method


# Regrillage

Example #4

File: courses_grids.py Project: acoat/vacumm

# 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(...)
print coord2slice(grid, lon=(4, 8), lat=(44, 46))  # -> TESTEZ SUR GRIDC
# -> TESTEZ VARSEL

# Transformer
gridcr = curv2rect(gridc)  # -> TESTEZ AVEC VOTRE GRILLE NON RECT
ax = create_lon([350, 0, 10.])
print monotonic(ax)
print xshift(var, 2)[0]  # -> COMPAREZ A VAR ET VERIFIEZ AXE
# -> TESTEZ XEXTEND
grido = rotate_grid(grid, 30)  # -> TRACEZ LES LONGITUDE (PCOLOR)
grids = shiftgrid(grid, 1,
                  -1)  # -> VERIFIEZ PUIS TESTEZ SUR GRIDC + PASSER DE T À U
gride = extendgrid(grid, iext=(2, 3))  # -> VERIFIER PUIS TESTEZ LE MODE

# Exploiter
print resol(grid)  # -> TESTEZ EN METRES ET SUR AXE
print depth2dz(dep)
print get_closest(lon2d, lat2d, 2.3, 1.2)

# Utilitaires sut les coordonnées
xx, yy = meshgrid(xx, yy[:, 0])
xxb, yyb = meshcells(xx, yy)  # -> EN 1D?