from timeserie.io.ascii.sirocco.SiroccoWriter import SiroccoWriter
from timeserie.io.ascii.puertos.PuertosTideReader import PuertosTideReader
import logging
import numpy as np
if __name__ == "__main__":
print("Transform tide data to SIROCCO format")
logging.basicConfig(format='[%(levelname)s] %(message)s',level=logging.INFO)
outDir="/tmp/"
## AJACCIO_ASPRETTO Données REFMAR
reader = RefmarReader('/home/retf/work/fieldsites/med-cruesim/observations/sea-level/raw-data/REFMAR/Ajaccio_aspretto.txt')
serie = TimeSerie(reader,'H','2010-01-01','2016-01-01');
serie.name_station = "AJACCIO_ASPRETTO"
serie.data_source = "SHOM / OCA / Marine nationale"
serie.vertical_datum = "zero_hydrographique"
serie.x_coord="8.76284981"
serie.y_coord="41.92279816"
logging.info(str(serie.name_station))
writer = SiroccoWriter(outDir+str(serie.name_station)+'_tide_'+str(serie.time_range[0].strftime("%Y"))+'_to_'+str(serie.time_range[serie.time_range.size-1].strftime("%Y"))+'.dat');
writer.write_tide(serie)
print 'End of programm'
for file in modelFiles:
filename = basename(file)
m = re.search('^([A-Z]+[\-\_]*[A-Z]*)\-[a-z]+.*', filename)
if m:
name_station = m.group(1)
logging.info("Station : "+str(name_station))
obsFiles = glob.glob(obsDir+name_station+"_*.dat")
if len(obsFiles) == 1:
# modèle
modelReader = WW3Reader(file)
model = TimeSerie(modelReader,start='2013-03-03',end='2013-04-01',freq='30min');
# observations
obsReader = SiroccoReader(obsFiles[0])
obs = TimeSerie(obsReader,start='2013-03-03',end='2013-04-01',freq='H');
# join sur les données
modelData = model.read_data()
obsData = obs.read_data()
result = obsData.join(modelData, how='inner',lsuffix='_obs', rsuffix='_model')
var_cor = correlation(result.sea_surface_wave_significant_height_obs,result.sea_surface_wave_significant_height_model)
var_bias = bias(result.sea_surface_wave_significant_height_obs,result.sea_surface_wave_significant_height_model)
var_rmse = rmse(result.sea_surface_wave_significant_height_obs,result.sea_surface_wave_significant_height_model)
var_si = si(result.sea_surface_wave_significant_height_obs,result.sea_surface_wave_significant_height_model)
if __name__ == "__main__":
print ("Transform waves data to SIROCCO format")
logging.basicConfig(format="[%(levelname)s] %(message)s", level=logging.INFO)
outDir = "/tmp/"
#
# Data from Puertos
#
sourceDirPuertos = "/home/retf/work/fieldsites/med-cruesim/observations/waves/raw-data/PUERTOS/"
## CaboBegur
reader = PuertosOceanReader(sourceDirPuertos + "REDEXT_T_HIS_CaboBegur.dat")
serie = TimeSerie(reader, "H", "2010-01-01", "2016-01-01")
serie.name_station = "CABO-BEGUR"
serie.data_source = "Puertos del Estado"
serie.meta_data = "Burst = 26 min"
serie.x_coord = "3.65"
serie.y_coord = "41.92"
logging.info(str(serie.name_station))
writer = SiroccoWriter(
outDir
+ str(serie.name_station)
+ "_wave_"
+ str(serie.time_range[0].strftime("%Y"))
+ "_to_"
+ str(serie.time_range[serie.time_range.size - 1].strftime("%Y"))
+ ".dat"
)
