source: branches/GRISLIv3/DOC/LeBrocq-to-Grisli.txt @ 383

Donnees d'entree pour GRISLI a partir des cartes de Le Brocq


La construction des fichiers se fait dans MAPS/Antarctique-LeBrocq

prendre la carte lsrf2 (base de la glace) et topg2 (bed) qui utilisent le creusement dans recovery.
   faire lsrf2.grd et topg2.grd avec 
           grdmath ALBMAPv1.nc?lsfr2 1 MUL = lsrf2.grd

* prendre usrf (surface) et soustraire firn pour travailler en glace equivalent
           grdmath ALBMAPv1.nc?usrf ALBMAPv1.nc?firn SUB = usrf-ice.grd 

         et lui changer le header
           grdedit usrf-ice-equ.grd -D=/=/"upper ice surface elevation in ice equivalent "/=/=/=/=

* pour l'épaisseur surface - base de la glace   : thick2-ice-equ.grd
            grdmath usrf-ice-equ.grd lsrf2.grd SUB 0 MAX  = thick2-ice-equ.grd
    et changer le header

* recalculer la surface a partir de l'epaisseur pour ne pas avoir de correction firn dans les zones sans glace
           grdmath thick2-ice-equ.grd lsrf2.grd ADD = usrf-ice-equ-final.grd

* prendre l'accumulation d'Arthern (m equ. glace)   acca. mais mettre +2 m dans l'ocean
           grdmath ALBMAPv1.nc?acca -9999 NAN 1.999 XOR = acca.grd

* le flux geothermique de Shapiro (parce qu'il est sur toute la zone)
           grdmath ALBMAPv1.nc?ghfsr 1 MUL = ghfsr.grd

* prendre la temperature en surface de Comiso (temp)
           grdmath ALBMAPv1.nc?temp 1 MUL = temp.grd
        
* Pour le masque : choisir mask_plus : 
-  0 : ocean libre
-  1 : pose avec glace 
-  2 : ice shelves
-  3 : outcrops
-  4  : ice streams 

exemple de grdinfo

usrf-ice-equ-final.grd: Title: ice thickness in ice equiv (incl. curvebed in Recovery)
usrf-ice-equ-final.grd: Command: grdmath thick2-ice-equ.grd lsrf2.grd ADD = usrf-ice-equ-final.grd
usrf-ice-equ-final.grd: Remark: 
usrf-ice-equ-final.grd: Gridline node registration used
usrf-ice-equ-final.grd: Grid file format: nf (# 18)
usrf-ice-equ-final.grd: x_min: -2.7975e+06  x_max: 2.9975e+06 x_inc: 5000 name: x-coordinate [meter] nx: 1160
usrf-ice-equ-final.grd: y_min: -2.7975e+06 y_max: 2.7975e+06 y_inc: 5000 name: y-coordinate [meter] ny: 1120
usrf-ice-equ-final.grd: z_min: 0 z_max: 4043 name: ice thickness in ice equivalent  [meter]

-R-2.7975e+06/2.9975e+06/-2.7975e+06/2.7975e+06 -I5000

Finalement : liste des grd a utiliser

acca.grd               : accum
temp.grd               : temperature en surface
ghfsr.grd              : flux geothermique
mask_plus.grd          : masque
lsrf2.grd              : base de la glace
thick2-ice-equ.grd     : epaisseur
topg2.grd              : socle et bathymetrie
usrf-ice-equ-final.grd : surface


Passage en 15 km

Bornes 

La carte de Le Brocq n'est pas centrée sur le pole et je garde cet aspect

 xmin_grisli = -2797500             
 xmax_grisli =  2992500            
 ymin_grisli = -2797500           
 ymax_grisli =  2797500  

makegrid-15.sh qui appelle average-5-to-15km-LeBrocq.f90  
nx = 387   ny = 374

Pour le masque, je fais sample-5-15

Carte Le Brocq 15 km -> LBq-15


-I15000 -R-2797500/2992500/-2797500/2797500


Masque pour eviter les iles : 
Avoir un fichier i,j des points a eviter.
Les noter sous ncview (sans oublier d'ajouter 1 a i et j)
puis utiliser ij_grisli2mask.sh qui produit un masque 1 pour ces points, 0 ailleurs

ij_grisli2mask.sh points-iles-a-enlever.dat
grdmath masque-glace_15km.grd mask-from-ij.grd SUB = mask-glace-ss-iles.grd

inverser ce masque : 
grdmath mask-glace-ss-iles.grd 1 SUB ABS = mask-glace-ss-iles.grd

en faire un ZBL
grd2xyz mask-glace-ss-iles.grd -ZBL > mask-glace-ss-iles_ZBL_15km.dat


Finalement faire un masque : mask_plus_ss_iles_ZBL_15km.dat  
colonne 1 : mask_plus
colonne 2 mask glace sans iles (0 pour glace , 1 ailleurs)
colonne 3 et 4 i j
(comme ca on peut intervenir directement dans le fichier pour corriger une ile qui aurait des pb)


longitude et latitude 
mapproject xy.dat  -Js0/-90/-71/1:50000000 -F -R/0/360/-90/-60 -C -I > lonlat.dat 
awk '{print $2}' lonlat.dat |xyz2grd -ZBL -I15000 -R-2797500/2992500/-2797500/2797500 -Glat_LBq_15km.grd
awk '{print $1}' lonlat.dat |xyz2grd -ZBL -I15000 -R-2797500/2992500/-2797500/2797500 -Glong_LBq_15km.grd