1 | #!/bin/bash |
---|
2 | |
---|
3 | #Petit script super simple pour tracer le Groenland |
---|
4 | |
---|
5 | # On part du fichier hz, S, H, B enroule j a l'exterieur |
---|
6 | # nx=41 ny= 71 |
---|
7 | |
---|
8 | |
---|
9 | file=$1 |
---|
10 | |
---|
11 | RANGE=-R0/1600/0/2800 |
---|
12 | PROJEC=-JX12/21 |
---|
13 | |
---|
14 | # creation des fichiers grd |
---|
15 | tail +5 $file | awk '{print $1}' |xyz2grd -Gsurf.grd -I40 $RANGE \ |
---|
16 | -ZBL |
---|
17 | |
---|
18 | grdsample surf.grd -Gsurf-fin.grd $RANGE -I10 |
---|
19 | |
---|
20 | #grd2xyz surf-fin.grd -ZBL > surf-fin.z |
---|
21 | |
---|
22 | tail +5 $file | awk '{print $3}' |xyz2grd -Gbed.grd -I40 $RANGE \ |
---|
23 | -ZBL |
---|
24 | |
---|
25 | grdsample surf.grd -Gbed-fin.grd $RANGE -I10 |
---|
26 | |
---|
27 | tail +5 $file | awk '{print $2}' |xyz2grd -GH.grd -I40 $RANGE \ |
---|
28 | -ZBL |
---|
29 | |
---|
30 | grdsample H.grd -GH-fin.grd $RANGE -I5 |
---|
31 | |
---|
32 | #grd2xyz H-fin.grd -ZBL > Hfin.z |
---|
33 | |
---|
34 | #paste surf-fin.z Hfin.z > SH.z |
---|
35 | |
---|
36 | #awk '{ if ($2 <=10 ) print $1/10000 ; else print $1 }' SH.z \ |
---|
37 | # | xyz2grd -Gcombine-fin.grd -I10 $RANGE -ZBL |
---|
38 | |
---|
39 | # tail +5 $file | awk '{ if ($2 <=1 ) print $1/10000 ; else print $1 }' \ |
---|
40 | # |xyz2grd -Gcombine.grd -I40 $RANGE -ZBL |
---|
41 | # grdsample combine.grd -Gcombine-fin.grd $RANGE -I10 -Q |
---|
42 | |
---|
43 | rm contour*xyz |
---|
44 | # la mer en bleu |
---|
45 | psxy <<EOF -P $RANGE $PROJEC -G140/102/255 -K > POST |
---|
46 | 0 0 |
---|
47 | 1600 0 |
---|
48 | 1600 2800 |
---|
49 | 0 2800 |
---|
50 | 0 0 |
---|
51 | EOF |
---|
52 | |
---|
53 | grdcontour bed.grd -Ccote.iso $RANGE $PROJEC -D -S20 > BIDON |
---|
54 | psclip contour*xyz $RANGE $PROJEC -O -K >> POST |
---|
55 | |
---|
56 | # la Terre en vert |
---|
57 | psxy <<EOF -P $RANGE $PROJEC -G000/255/000 -K -O >> POST |
---|
58 | 0 0 |
---|
59 | 1600 0 |
---|
60 | 1600 2800 |
---|
61 | 0 2800 |
---|
62 | 0 0 |
---|
63 | EOF |
---|
64 | |
---|
65 | # enleve le clip terre |
---|
66 | psclip -C -K -O >> POST |
---|
67 | rm contour*xyz |
---|
68 | |
---|
69 | |
---|
70 | #grdimage bed-fin.grd -Ba500f100 -Cmer-grass.cpt -P $RANGE $PROJEC -K -X2 -Y1 > POST |
---|
71 | |
---|
72 | grdcontour H-fin.grd -Ccontour.iso -W4/255/0/0 $RANGE $PROJEC -D -S20 > bidon |
---|
73 | |
---|
74 | psclip contour*xyz $RANGE $PROJEC -O -K >> POST |
---|
75 | |
---|
76 | #grdimage surf-fin.grd -Ba500f100 -Cglace-bleue.cpt -P $RANGE $PROJEC -K -O >> POST |
---|
77 | |
---|
78 | # grdview combine.grd -Ba500f100 -Cglace-sur-land.cpt -Qi -P $RANGE $PROJEC -K -X2 -Y1 > POST |
---|
79 | |
---|
80 | #grdcontour H-fin.grd -Csurf.iso $RANGE $PROJEC -O -K >> POST |
---|
81 | |
---|
82 | |
---|
83 | #grdview surf-fin.grd -Cglace-bleue.cpt -Qs -P $RANGE $PROJEC -K -O >> POST |
---|
84 | grdimage surf-fin.grd -Cglace-bleue.cpt -P $RANGE $PROJEC -K -O >> POST |
---|
85 | grdcontour surf-fin.grd -C500 -A1000 $RANGE $PROJEC -O -K >> POST |
---|
86 | |
---|
87 | psclip -C -K -O >> POST |
---|
88 | # grdcontour surf.grd -Ccontour.iso $RANGE $PROJEC -O -K >> POST |
---|
89 | |
---|
90 | |
---|
91 | |
---|
92 | #grdcontour H-fin.grd -Ccontour.iso -W4/255/0/0 $RANGE $PROJEC -O -S20 >> POST |
---|
93 | psxy contour*xyz $RANGE $PROJEC -O -K -W4/255/0/0 >> POST |
---|
94 | |
---|
95 | # gates |
---|
96 | psxy gate*.txt $RANGE $PROJEC -O -K -W15 >> POST |
---|
97 | |
---|
98 | # trace les meridiens et paralleles |
---|
99 | # lon lat des coins bas gauche -55.82902 58.77798 haut droit : 12.45875 80.15318 |
---|
100 | |
---|
101 | #psbasemap -JS-44/90/12 -R-55.82902/58.77798/12.45875/80.15318r -Ba10g10/a5g5 -O >> POST |
---|
102 | pscoast -JS-44/90/12 -R-55.82902/58.77798/12.45875/80.15318r -Ba10g10/a5g5 -O -W2 >> POST |
---|