Methode effectivement utilisee pour les initial guess des vitesses: 
 
Ca parait plus simple de prendre le run dans la partie flottante et la vitesse de Le Brocq dans la partie posee

Ne garder que la partie shelves des runs
grdmath Uy-LBq15km_A15spn61.grd mask_shelves_15km-NAN.grd MUL = Uy-LBq15km_A15spn61-shelves.grd

Merger les deux
grdmath Uy-LBq15km_A15spn61-shelves.grd Uy-LeBrocq.grd AND = Uy-merge-LeBrocq-A15spn61.grd

Verifier l'apparence de l'angle et de la norme

angle
grdmath  Ux-merge-LeBrocq-A15spn61.grd Uy-merge-LeBrocq-A15spn61.grd  ATAN2 R2D = angle-merge-LeBrocq-A15spn61.grd

norme en changeant les 9999 en -9999 en passant par le Nan

grdmath Ux-merge-LeBrocq-A15spn61.grd Uy-merge-LeBrocq-A15spn61.grd HYPOT 9999 NAN -9999 XOR = norme-merge-LeBrocq-A15spn61.grd

faire le masque glace 


Construction du fichier input_SHbmMkULBq-ux-uy-15km.dat 

pour  nourir vitbil_LeBrocq_to_Grisli-15km-avr2010.f90

Il tourne 900000 iterations et produit des vitbil_iter_900000

Pour visialiser une iteration :

tail -n +5 vitbil_iter_100000 |awk '{print $2*$4}' |xyz2grd -ZBL -R-2797500/2992500/-2797500/2797500 -I15000 -Ghdot_100000.grd 

tail -n +5 vitbil_iter_100000 |awk '{print $2*($5*$5+$6*$6)**0.5}' |xyz2grd -ZBL -R-2797500/2992500/-2797500/2797500 -I15000 -GUmag_100000.grd

H,masque,bm,hdot,uxmaj,uymaj,ux,uy

Il ne semble pas y avoir un gros gain passé 300000 iterations

Remarque importante

Il y a des problemes dans les zones ou il y a beaucoup de Nunataks (typiquement le glacier Lambert) :

* ce serait mieux de faire une moyenne sur les flux plutôt que sur les vitesses. Comme un noeud de faible epaisseur tend a avoir une vitesse forte, il est sur-represente.


Nouvelle Procedure : 

* On passe en flux en multipliant les vitesses merge (Ux-merge-LeBrocq-A15spn61.grd) par l'epaisseur  
- fluX-merge-LeBrocq-A15spn61.grd et fluY-merge-LeBrocq-A15spn61.grd

* On met en Nan tous les 0 
- grdmath fluX-merge-LeBrocq-A15spn61.grd 0 NAN = fluX-merge-LeBrocq-A15spn61-NAN.grd

* On filtre avec une moyenne mobile sur 30 km
- grdfilter fluY-merge-LeBrocq-A15spn61-NAN.grd -D0 -Fb30000 -GfluY-merge-LeBrocq-A15spn61-Fb30.grd

* On calcule la direction
- grdmath fluX-merge-LeBrocq-A15spn61-Fb30.grd fluY-merge-LeBrocq-A15spn61-Fb30.grd ATAN2 R2D = angle-flux-merge-Fb30.grd
	
* On recalcule Ux et Uy en utilisant la direction lissee 
- grdmath angle-flux-merge-Fb30.grd 90 SUB COSD norme-merge-LeBrocq-A15spn61.grd MUL -9999 AND  = Ux-merge-dirFb30.grd
- grdmath angle-flux-merge-Fb30.grd -90 SUB SIND norme-merge-LeBrocq-A15spn61.grd MUL -9999 AND = Uy-merge-dirFb30.grd

* On fait les ZBL
- grd2xyz -ZBL Ux-merge-dirFb30.grd >Ux-merge-dirFb30_ZBL_15km.dat

* et on re-paste le fichier input_SHbmMkULBq-ux-uy-15km.dat 



En utilisant les resultats d'un run GRISLI