#!/bin/bash
# script gmt pour convertir la carte de sediment sur la grille hemin40 de GRISLI Nord



awk '{print $5,$6}' /homel/dumas/CODE/version-6-Hudson-Hemin40-oct2008/INPUT/HEMIN40/coord-nord-40km.dat > lon_lat_hemin40.dat
awk '{print $1,$2}' /homel/dumas/CODE/version-6-Hudson-Hemin40-oct2008/INPUT/HEMIN40/coord-nord-40km.dat > ij_hemin40.dat

grdtrack lon_lat_hemin40.dat -Gsedmap.grd > sediment_hemin40.dat

# passage de l'epaisseur de sediment en metre
paste ij_hemin40.dat sediment_hemin40.dat | awk '{print $1,$2,$5*1000}' > sediment_ij_hemin40.dat


xyz2grd sediment_ij_hemin40.dat -I1 -Gsediment_hemin40.grd -R1/241/1/241

#makecpt -Csealand -T0/100/10 -I -Z > sediment2.cpt

psscale -P -Csediment2.cpt -D8/22/15/0.5h -B5:"Epaisseur de Sédiment":/:m: -L -V -K > sediment_hemin40.ps

grdimage sediment_hemin40.grd -Csediment2.cpt -P -R1/241/1/241 -JX15/15 -O -K -Y3 >> sediment_hemin40.ps

pscoast -R-47.12/16.366/137.72/34.48r -Bg45/g10/a10g10 -JA0/90/15 -Dc -P -W2 -O >> sediment_hemin40.ps

ps2raster -Tf sediment_hemin40.ps

gv sediment_hemin40.pdf