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"
