1 | #!/bin/bash |
---|
2 | |
---|
3 | # trace 1 carte de la variable dS |
---|
4 | # |
---|
5 | # exemple des réponses sous forme d'un fichier a rediriger dans |
---|
6 | # exemple-6cartes.dat |
---|
7 | #================================================================ |
---|
8 | cp $1 mask.grd # celui du snapshot a tracer |
---|
9 | |
---|
10 | #timetitle=$2 #celui sur le dessin |
---|
11 | color_name='COLOR/mask.cpt' |
---|
12 | iso_name='ISO/mask.iso' |
---|
13 | gl_name='ISO/gl_LGM.iso' |
---|
14 | iso_pres='ISO/surf.iso' |
---|
15 | |
---|
16 | # caracteristiques generales du dessin |
---|
17 | rm -f .gmtdefaults4 |
---|
18 | gmtset ANOT_FONT_SIZE 14p # taille des annotations sur les axes |
---|
19 | gmtset HEADER_FONT_SIZE 16p # taille des titres sur les axes |
---|
20 | gmtset LABEL_FONT_SIZE 16p |
---|
21 | |
---|
22 | |
---|
23 | # map size and position |
---|
24 | plot_size=20 |
---|
25 | size='-JX'$plot_size # taille en cm du cadre |
---|
26 | |
---|
27 | # ranges |
---|
28 | range0='-R-2797500/2992500/-2797500/2797500' # LBq |
---|
29 | step='-I15000' |
---|
30 | #range='-R-2790/2790/-2790/2790' #can be smaller than range0 for zooms |
---|
31 | |
---|
32 | # pour tracer les meridiens eet paralleles |
---|
33 | #mapproject << END -Js0/-90/-71/1:50000000 -Fk -R/0/360/-90/-60 -C -I |
---|
34 | #> 2790 0 |
---|
35 | #> END |
---|
36 | # |
---|
37 | #steropol='-R0/360/-90./-64.7233 -JS0/-90/'$plot_size' -Bg30/g5' |
---|
38 | steropol='-R225/-54.7922/45./-54.7922r -JS0/-90/'$plot_size' -Bg30/g5' |
---|
39 | |
---|
40 | # pos : position et tics du cadre |
---|
41 | # scalpos : position de l échelle de couleur |
---|
42 | |
---|
43 | # placement des cartes |
---|
44 | titre_var='(km)' |
---|
45 | title=' surface elevation (m)' |
---|
46 | pos='-X0.5 -Y0.5' |
---|
47 | Btics='-Bf1000a2000:km:/f1000a2000:km:WeSn' |
---|
48 | |
---|
49 | # placement echelle de couleur |
---|
50 | |
---|
51 | #scalpos='-X20 -Y8' |
---|
52 | #scalsize='-D2/2.5/13/0.3' # taille de l'échelle |
---|
53 | |
---|
54 | scalpos='-X8 -Y17' |
---|
55 | scalsize='-D2/2.5/13/0.3h' # taille de l'échelle |
---|
56 | |
---|
57 | |
---|
58 | #Lecture des fichiers nc + correction des axes |
---|
59 | #=============================================== |
---|
60 | |
---|
61 | # grdedit S.grd $range0 # ne semble pas marcher |
---|
62 | grd2xyz mask.grd -ZBL | xyz2grd -ZBL -Gmask.grd $range0 $step |
---|
63 | |
---|
64 | grdmath mask.grd mgrounded.grd ADD = mask.grd |
---|
65 | # -------------------------------------------- LEGENDE GENERALE |
---|
66 | #pstext -K -JX21/29.6 -R0/21/0/29.6 -X0 -Y0 <<END > POST |
---|
67 | #11 5 15 0. 0 10 title |
---|
68 | #END |
---|
69 | |
---|
70 | |
---|
71 | grdimage mask.grd $range $size $pos -K -C$color_name > POST |
---|
72 | |
---|
73 | grdcontour mask.grd $range $size -O -K -C$iso_name -W4/70/70/70 >> POST |
---|
74 | grdcontour S_gl.grd $range $size -O -K -C$gl_name -W6/0/0/255 >> POST |
---|
75 | |
---|
76 | |
---|
77 | # passe en stereo pour tracer mediens et paralleles |
---|
78 | gmtset BASEMAP_FRAME_RGB 150/150/150 |
---|
79 | |
---|
80 | #psbasemap -R0/360/-90./-64.7233 -JS0/-90/$plot_size -Bg30/g5 -O >> POST |
---|
81 | psbasemap $steropol -O >> POST |
---|
82 | gmtset BASEMAP_FRAME_RGB 0/0/0 |
---|
83 | |
---|
84 | rm -f .gmtdefaults4 |
---|
85 | rm -f .gmtcommand |
---|
86 | |
---|
87 | #pstext -O -K $range $size <<EOF >> POST |
---|
88 | #-2100 -2400 30 0. 0 10 $timetitle yr |
---|
89 | #EOF |
---|
90 | |
---|
91 | #psscale -C$color_name -L $scalsize -E $scalpos -O >> POST |
---|
92 | |
---|
93 | |
---|
94 | |
---|
95 | |
---|
96 | |
---|
97 | |
---|
98 | gmtset ANOT_FONT_SIZE 14p # remet la taille des annotations à la valeur par défaut |
---|
99 | |
---|
100 | #================================================================================= |
---|