#!/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)