#!/usr/bin/env python
# -*- coding: utf8 -*-
"""Diagnostics de maree"""
from vcmq import cdms2, P, curve2, savefigs, data_sample
from vacumm.tide.filters import demerliac, godin
from vacumm.tide.filters import extrema, zeros
from vacumm.tide.marigraph import Marigraph
from vacumm.tide.station_info import StationInfo
# Stations
station = StationInfo('Brest')
print station.attributes()
print station.name, station.longitude
print 'Niveau moyen a Brest:', station.nm
# Read sea level at Brest
f = cdms2.open(data_sample("tide.sealevel.BREST.mars.nc"))
sea_level = f('sea_level')
f.close()
# Surcotes/decotes
cotes, tide = demerliac(sea_level, get_tide=True) # -> ESSAYER GODIN
kwp = dict(date_fmt='%b', date_locator='month')
curve2(sea_level, 'b', show=False, figsize=(15, 4), **kwp)
curve2(tide, 'r', **kwp)
curve2(cotes, figsize=(15, 4), **kwp)
# Extremas
slzoom1 = sea_level(('2006-10-01', '2006-10-02'))[::4] # Toutes les heures
bm, pm = extrema(slzoom1, spline=True, ref='mean') # -> SANS SPLINES
# -*- coding: utf8 -*-
# On cherche directement 'Bre' pour Brest
from vacumm.tide.station_info import StationInfo
station = StationInfo('Bre')
# ->:
#Chargement de la station suivante :
# Nom : Ile de Brehat (Port-Clos)
# Position : 3.0°W / 48.9°N
# Zone : http://www.shom.fr/fr_page/fr_act_oceano/maree/zone6_9.htm
# BM45 : 3.8
# BM95 : 1.35
# NM : 5.89
# PBM : 0.1
# PHM : 11.6
# PM45 : 8
# PM95 : 10.4
# ZERO_HYDRO : -5.4
#Definition des termes accessible avec .definitions()
# Loupé !
# On récupère la première station trouvée lors de l'initialisation
# Mais bon, on vérifie quand même
print station.attributes()
# -> ['igs', 'psmsl', 'uhslc', 'gloss', 'shom', 'legos', 'latitude',
# 'longitude', 'zone', 'phm', 'pm95', 'pm45', 'nm', 'bm45',
# 'bm95', 'pbm', 'zero_hydro']
print station.name, station.longitude
# -> Ile de Brehat (Port-Clos) -3.0
# On peut se servir de station pour continuer à chercher
# car le fichier est déjà chargé
#!/usr/bin/env python
# -*- coding: utf8 -*-
"""Diagnostics de maree"""
from vcmq import cdms2, P, curve2, savefigs, data_sample
from vacumm.tide.filters import demerliac, godin
from vacumm.tide.filters import extrema, zeros
from vacumm.tide.marigraph import Marigraph
from vacumm.tide.station_info import StationInfo
# Stations
station = StationInfo('Brest')
print station.attributes()
print station.name, station.longitude
print 'Niveau moyen a Brest:', station.nm
# Read sea level at Brest
f = cdms2.open(data_sample("tide.sealevel.BREST.mars.nc"))
sea_level = f('sea_level')
f.close()
# Surcotes/decotes
cotes, tide = demerliac(sea_level, get_tide=True) # -> ESSAYER GODIN
kwp = dict(date_fmt='%b', date_locator='month')
curve2(sea_level, 'b', show=False, figsize=(15, 4), **kwp)
curve2(tide, 'r', **kwp)
curve2(cotes, figsize=(15, 4), **kwp)
# -*- coding: utf8 -*-
# On cherche directement 'Bre' pour Brest
from vacumm.tide.station_info import StationInfo
station = StationInfo('Bre')
# ->:
#Chargement de la station suivante :
# Nom : Ile de Brehat (Port-Clos)
# Position : 3.0°W / 48.9°N
# Zone : http://www.shom.fr/fr_page/fr_act_oceano/maree/zone6_9.htm
# BM45 : 3.8
# BM95 : 1.35
# NM : 5.89
# PBM : 0.1
# PHM : 11.6
# PM45 : 8
# PM95 : 10.4
# ZERO_HYDRO : -5.4
#Definition des termes accessible avec .definitions()
# Loupé !
# On récupère la première station trouvée lors de l'initialisation
# Mais bon, on vérifie quand même
print station.attributes()
# -> ['igs', 'psmsl', 'uhslc', 'gloss', 'shom', 'legos', 'latitude',
# 'longitude', 'zone', 'phm', 'pm95', 'pm45', 'nm', 'bm45',
# 'bm95', 'pbm', 'zero_hydro']
print station.name, station.longitude
# -> Ile de Brehat (Port-Clos) -3.0
# On peut se servir de station pour continuer à chercher
# car le fichier est déjà chargé
