level=logging.INFO)
# Read file
reader = SymphonieReader2015(
'/work/sciences/projects/WWB-2017/Manicouagan/configuration_V2015/TStra_N/OFFLINE/grid.nc',
'/work/sciences/projects/WWB-2017/Manicouagan/configuration_V2015/TStra_N/GRAPHIQUES/20090301_071217.nc'
)
coverageOrig = TimeLevelCoverage(reader)
# InterpolatorCore.INTERPOLATION_METHOD = "linear";
InterpolatorCore.INTERPOLATION_METHOD = "nearest"
depths = np.array([0.0, 10.0, 50.0, 100.0])
coverage = TimeLevelCoverage(
CoverageInterpolator(coverageOrig, 0.001, 0.001, depths))
#Renommer toutes les noms des fonctions
writer = DefaultWriter(coverage, '/tmp/symphonie_regular.nc')
writer.write_variable_baroclinic_sea_water_velocity()
writer.write_variable_barotropic_sea_water_velocity()
writer.write_variable_sea_surface_height_above_mean_sea_level()
writer.write_variable_wind_stress()
writer.write_variable_sea_water_temperature()
writer.write_variable_sea_water_salinity()
# writer.write_variable_wet_binary_mask()
writer.write_variable_2D_sea_binary_mask()
writer.write_variable_wind_10m()
writer.write_variable_mesh_size()
#writer.write_variable_sea_surface_wave_significant_height()
logging.basicConfig(format='[%(levelname)s] %(message)s',level=logging.INFO)
# Read file
reader = MercatorReader('/home/retf/work/fieldsites/med-cruesim/modelling/mercator/grid/mercator_grid.nc',
'/home/retf/work/fieldsites/med-cruesim/modelling/mercator/netcdf/pool/ext-PSY2V4R4_1dAV_20120706_20120707_grid2D_R20120711.nc',
'/home/retf/work/fieldsites/med-cruesim/modelling/mercator/netcdf/pool/ext-PSY2V4R4_1dAV_20120706_20120707_gridT_R20120711.nc',
'/home/retf/work/fieldsites/med-cruesim/modelling/mercator/netcdf/pool/ext-PSY2V4R4_1dAV_20120706_20120707_gridU_R20120711.nc',
'/home/retf/work/fieldsites/med-cruesim/modelling/mercator/netcdf/pool/ext-PSY2V4R4_1dAV_20120706_20120707_gridV_R20120711.nc')
coverage = TimeLevelCoverage(reader);
LevelCoverage.LEVEL_DELTA = 50.0 # On change le delta par défaut pour trouver une couche
interpolator = CoverageInterpolator(coverage,0.01,0.01,'/home/retf/work/fieldsites/med-cruesim/modelling/mercator/netcdf/pool/test_resample.nc') # résolution voulue en degrès
interpolator.resample_variable_current_at_level(155.0) # on peut donner la profondeur en mètres positifs (plus proche) ou en indice de couche verticale.
interpolator.resample_variable_ssh()
interpolator.close()
print 'End of programm'
if __name__ == "__main__":
print("Transform/Interpole Symphonie to GMT")
logging.basicConfig(format='[%(levelname)s] %(message)s',level=logging.INFO)
# Read file
reader = WW3Reader('/work/sciences/projects/ArcEmeraude_2018/waves/Re_ Bon___/ww3.20150105_T405_2s.nc')
coverageOrig = TimeCoverage(reader);
# InterpolatorCore.INTERPOLATION_METHOD = "linear";
InterpolatorCore.INTERPOLATION_METHOD = "nearest";
depths = np.array([])
coverage = TimeCoverage(CoverageInterpolator(coverageOrig, 0.1, 0.1, depths))
#Renommer toutes les noms des fonctions
writer = DefaultWriter(coverage, '/tmp/ww3_regular.nc')
writer.write_variable_barotropic_sea_water_velocity()
writer.write_variable_sea_surface_height_above_mean_sea_level()
writer.write_variable_2D_sea_binary_mask()
writer.write_variable_sea_surface_wave_significant_height()
writer.write_variable_sea_surface_wave_mean_period()
# writer.write_variable_sea_surface_wave_peak_period()
writer.write_variable_sea_surface_wave_from_direction()
writer.write_variable_sea_surface_wave_to_direction()
writer.write_variable_atmosphere_momentum_flux_to_waves()
writer.write_variable_waves_momentum_flux_to_ocean()
writer.write_variable_sea_surface_wave_stokes_drift_velocity()
from coverage.operator.interpolator.CoverageInterpolator import CoverageInterpolator
from coverage.io.netcdf.symphonie.SymphonieReader import SymphonieReader
import logging
if __name__ == "__main__":
print("Transform/Interpole Symphonie to GMT")
logging.basicConfig(format='[%(levelname)s] %(message)s',level=logging.INFO)
# Read file
reader = SymphonieReader('GRAPHIQUES/grid.nc',
'GRAPHIQUES/20110320_000002.nc')
coverage = TimeLevelCoverage(reader);
interpolator = CoverageInterpolator(coverage,0.02,0.02,'GRAPHIQUES/20110220_regular.nc',
[]) # résolution voulue en degrès
#interpolator = CoverageInterpolator(coverage,0.02,0.02,'GRAPHIQUES/20110220_regular.nc',
# [0.0,10.0,50.0,150.0,200.0,250.0,300.0,350.0,400.0,450.0,500.0,550.0,600.0,650.0,700.0,750.0,850.0,900.0,950.0,1000.0,1100.0,1200.0,1300.0,1400.0,1500.0,1600.0,1700.0,1800.0,1900.0,2000.0,2300.0,2400.0,2500.0,2600.0,2700.0,2850.0,2870.0]) # résolution voulue en degrès
#interpolator.resample_variable_current_at_depths()
interpolator.resample_variable_salinity_at_depths()
#interpolator.resample_variable_ssh()
#interpolator.resample_variable_current()
#interpolator.resample_variable_2D_mask()
#interpolator.resample_variable_waves_dir()
#interpolator.resample_variable_hs()
#interpolator.resample_variable_wind()
#interpolator.resample_variable_mesh_size()
#interpolator.resample_variable_wetmask()
interpolator.close()
from coverage.TimeCoverage import TimeCoverage
from coverage.operator.interpolator.CoverageInterpolator import CoverageInterpolator
from coverage.io.netcdf.ww3.WW3Reader import WW3Reader
import logging
if __name__ == "__main__":
"""
Cette routine permet d'interpoler des résultats de WaveWatchIII sur une grille régulière.
"""
print ("Transform/Interpole WW3 to GMT ")
logging.basicConfig(format="[%(levelname)s] %(message)s", level=logging.INFO)
# Read file
reader = WW3Reader("/home/retf/work/fieldsites/med-cruesim/modelling/waves/med/outputs/netcdf/ww3.201103.nc")
coverage = TimeCoverage(reader)
interpolator = CoverageInterpolator(
coverage,
0.01,
0.01,
"/home/retf/work/fieldsites/med-cruesim/modelling/waves/med/outputs/netcdf/regular/ww3.201103.nc",
) # résolution voulue en degrès
interpolator.resample_variable_current()
interpolator.resample_variable_ssh()
interpolator.close()
print "End of programm"
