#!/bin/bash

# trace 1 carte de la variable dS
#
# exemple des réponses sous forme d'un fichier a rediriger dans 
# exemple-6cartes.dat
#================================================================
cp $1 mask.grd  # celui du snapshot a tracer

#timetitle=$2    #celui sur le dessin
color_name='COLOR/mask.cpt'
iso_name='ISO/mask.iso'
gl_name='ISO/gl_LGM.iso'
iso_pres='ISO/surf.iso'

# caracteristiques generales du dessin
rm -f .gmtdefaults4
gmtset ANOT_FONT_SIZE 14p # taille des annotations sur les axes
gmtset HEADER_FONT_SIZE 16p # taille des titres sur les axes
gmtset LABEL_FONT_SIZE 16p


# map size and  position
plot_size=20
size='-JX'$plot_size # taille en cm du cadre

# ranges
range0='-R-2797500/2992500/-2797500/2797500'  # LBq
step='-I15000'
#range='-R-2790/2790/-2790/2790'  #can be smaller than range0 for zooms

# pour tracer les meridiens eet paralleles
#mapproject << END -Js0/-90/-71/1:50000000 -Fk -R/0/360/-90/-60 -C -I
#> 2790 0
#> END
#
#steropol='-R0/360/-90./-64.7233 -JS0/-90/'$plot_size' -Bg30/g5'
steropol='-R225/-54.7922/45./-54.7922r -JS0/-90/'$plot_size' -Bg30/g5'

# pos : position et tics du cadre
# scalpos : position de l échelle de couleur

# placement des cartes
titre_var='(km)'
title=' surface elevation (m)'
pos='-X0.5 -Y0.5'
Btics='-Bf1000a2000:km:/f1000a2000:km:WeSn'

# placement echelle de couleur

#scalpos='-X20 -Y8'
#scalsize='-D2/2.5/13/0.3' # taille de l'échelle

scalpos='-X8 -Y17'
scalsize='-D2/2.5/13/0.3h' # taille de l'échelle


#Lecture des fichiers nc + correction des axes
#===============================================

# grdedit S.grd $range0  # ne semble pas marcher
grd2xyz mask.grd -ZBL | xyz2grd -ZBL -Gmask.grd $range0 $step

grdmath mask.grd mgrounded.grd ADD = mask.grd
# -------------------------------------------- LEGENDE GENERALE
#pstext -K   -JX21/29.6 -R0/21/0/29.6 -X0 -Y0 <<END > POST
#11 5   15  0. 0  10 title
#END


grdimage mask.grd $range $size  $pos  -K -C$color_name  > POST

grdcontour mask.grd  $range $size  -O -K  -C$iso_name -W4/70/70/70 >> POST
grdcontour S_gl.grd  $range $size  -O -K  -C$gl_name  -W6/0/0/255 >> POST


# passe en stereo pour tracer mediens et paralleles
gmtset BASEMAP_FRAME_RGB 150/150/150

#psbasemap -R0/360/-90./-64.7233 -JS0/-90/$plot_size -Bg30/g5 -O >> POST
psbasemap $steropol -O  >> POST
gmtset BASEMAP_FRAME_RGB 0/0/0

rm -f .gmtdefaults4
rm -f .gmtcommand

#pstext -O -K  $range $size <<EOF >> POST
#-2100 -2400  30  0.  0  10 $timetitle yr
#EOF

#psscale  -C$color_name -L $scalsize -E $scalpos   -O  >> POST






gmtset ANOT_FONT_SIZE 14p # remet la taille des annotations à la valeur par défaut

#=================================================================================