#-----------------------------------------------------------
# A) on coupe sur la zone Groenland, et on "affine" a 30 sec
grdcut ../IN_OUTPUT/IN_INIT/ETOPO1_Ice_g.grd -G../IN_OUTPUT/TEMP/etopo1min_surface_groen.grd -R/-110/20/50/85
grdsample ../IN_OUTPUT/TEMP/etopo1min_surface_groen.grd -G../IN_OUTPUT/TEMP/etopo30s_surface.grd -I30c
rm ../IN_OUTPUT/TEMP/etopo1min_surface_groen.grd
#-----------------------------------------------------------


#-----------------------------------------------------------
# B) on passe du grd en lon/lat/z en un ASCII x/y/z
# on utilise la projection de Bamber
grd2xyz ../IN_OUTPUT/TEMP/etopo30s_surface.grd -R-110/20/50/85  > ../IN_OUTPUT/TEMP/etopo30s_surface.lonlatz
mapproject ../IN_OUTPUT/TEMP/etopo30s_surface.lonlatz  -Js-39/90/71/1:50000000 -Fk -R/-110/20/50/85 -C   > ../IN_OUTPUT/TEMP/etopo30s_surface.xyz.txt
#rm ../IN_OUTPUT/TEMP/etopo30s_surface.lonlatz
rm ../IN_OUTPUT/TEMP/etopo30s_surface.grd
#-----------------------------------------------------------


#-----------------------------------------------------------
# C) on recree des fichiers grd et xyz (ordonne)
# Grille finale plus grande que celle de Bamber
xyz2grd ../IN_OUTPUT/TEMP/etopo30s_surface.xyz.txt -G../IN_OUTPUT/TEMP/surface_xyz_5km.grd -I5 -R-990/895/-3515/-550
grd2xyz ../IN_OUTPUT/TEMP/surface_xyz_5km.grd > ../IN_OUTPUT/TEMP/surface_xyz_5km.txt
rm ../IN_OUTPUT/TEMP/etopo30s_surface.xyz.txt
#-----------------------------------------------------------