# Lecture du niveau de la mer
import cdms2
from vacumm.config import data_sample
f = cdms2.open(data_sample('tide.sealevel.BREST.mars.nc'))
sea_level = f('sea_level', time=('2006-10', '2006-10-07'))
f.close()
# Initialisation de l'objet marégraphique
from vacumm.tide.marigraph import Marigraph
mg = Marigraph(sea_level, verbose=True)
# Extraction du signal de maree
tide = mg.tide(tide_filter='demerliac')
# On peut aussi specifier le filtre avec :
mg.set_tide_filter('demerliac') # On peut aussi choisir 'godin'
# Calcul surcotes/decotes
cotes = mg.cotes()
# Calcul des pleines et basses mers
ref = 'mean'
highs = mg.highs(ref=ref)
lows = mg.lows(ref=ref)
zeros = mg.zeros(ref=ref)
# On plot
import pylab as P
P.rc('font', size=10)
P.figure(figsize=(5, 5))
kwplot = dict(date_fmt='%d/%m', date_locator='day', show=False, hspace=.2, left=.15)
# Lecture du niveau de la mer
import cdms2
from vacumm.config import data_sample
f = cdms2.open(data_sample('tide.sealevel.BREST.mars.nc'))
sea_level = f('sea_level', time=('2006-10', '2006-10-07'))
f.close()
# Initialisation de l'objet marégraphique
from vacumm.tide.marigraph import Marigraph
mg = Marigraph(sea_level, verbose=True)
# Extraction du signal de maree
tide = mg.tide(tide_filter='demerliac')
# On peut aussi specifier le filtre avec :
mg.set_tide_filter('demerliac') # On peut aussi choisir 'godin'
# Calcul surcotes/decotes
cotes = mg.cotes()
# Calcul des pleines et basses mers
ref = 'mean'
highs = mg.highs(ref=ref)
lows = mg.lows(ref=ref)
zeros = mg.zeros(ref=ref)
# On plot
import pylab as P
P.rc('font', size=10)
P.figure(figsize=(5, 5))
kwplot = dict(date_fmt='%d/%m',
