Context Navigation

sediment_hemin40.gmt @ 359

Last change on this file since 359 was 4, checked in by dumas, 10 years ago
initial import GRISLI trunk
File size: 1.0 KB

Line
1	#!/bin/bash
2	# script gmt pour convertir la carte de sediment sur la grille hemin40 de GRISLI Nord
3
4
5
6	awk '{print $5,$6}' /homel/dumas/CODE/version-6-Hudson-Hemin40-oct2008/INPUT/HEMIN40/coord-nord-40km.dat > lon_lat_hemin40.dat
7	awk '{print $1,$2}' /homel/dumas/CODE/version-6-Hudson-Hemin40-oct2008/INPUT/HEMIN40/coord-nord-40km.dat > ij_hemin40.dat
8
9	grdtrack lon_lat_hemin40.dat -Gsedmap.grd > sediment_hemin40.dat
10
11	# passage de l'epaisseur de sediment en metre
12	paste ij_hemin40.dat sediment_hemin40.dat \| awk '{print $1,$2,$5*1000}' > sediment_ij_hemin40.dat
13
14
15	xyz2grd sediment_ij_hemin40.dat -I1 -Gsediment_hemin40.grd -R1/241/1/241
16
17	#makecpt -Csealand -T0/100/10 -I -Z > sediment2.cpt
18
19	psscale -P -Csediment2.cpt -D8/22/15/0.5h -B5:"Epaisseur de Sédiment":/:m: -L -V -K > sediment_hemin40.ps
20
21	grdimage sediment_hemin40.grd -Csediment2.cpt -P -R1/241/1/241 -JX15/15 -O -K -Y3 >> sediment_hemin40.ps
22
23	pscoast -R-47.12/16.366/137.72/34.48r -Bg45/g10/a10g10 -JA0/90/15 -Dc -P -W2 -O >> sediment_hemin40.ps
24
25	ps2raster -Tf sediment_hemin40.ps
26
27	gv sediment_hemin40.pdf

Note: See TracBrowser for help on using the repository browser.