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" )