1 | Methode effectivement utilisee pour les initial guess des vitesses: |
---|
2 | |
---|
3 | Ca parait plus simple de prendre le run dans la partie flottante et la vitesse de Le Brocq dans la partie posee |
---|
4 | |
---|
5 | Ne garder que la partie shelves des runs |
---|
6 | grdmath Uy-LBq15km_A15spn61.grd mask_shelves_15km-NAN.grd MUL = Uy-LBq15km_A15spn61-shelves.grd |
---|
7 | |
---|
8 | Merger les deux |
---|
9 | grdmath Uy-LBq15km_A15spn61-shelves.grd Uy-LeBrocq.grd AND = Uy-merge-LeBrocq-A15spn61.grd |
---|
10 | |
---|
11 | Verifier l'apparence de l'angle et de la norme |
---|
12 | |
---|
13 | angle |
---|
14 | grdmath Ux-merge-LeBrocq-A15spn61.grd Uy-merge-LeBrocq-A15spn61.grd ATAN2 R2D = angle-merge-LeBrocq-A15spn61.grd |
---|
15 | |
---|
16 | norme en changeant les 9999 en -9999 en passant par le Nan |
---|
17 | |
---|
18 | grdmath Ux-merge-LeBrocq-A15spn61.grd Uy-merge-LeBrocq-A15spn61.grd HYPOT 9999 NAN -9999 XOR = norme-merge-LeBrocq-A15spn61.grd |
---|
19 | |
---|
20 | faire le masque glace |
---|
21 | |
---|
22 | |
---|
23 | Construction du fichier input_SHbmMkULBq-ux-uy-15km.dat |
---|
24 | |
---|
25 | pour nourir vitbil_LeBrocq_to_Grisli-15km-avr2010.f90 |
---|
26 | |
---|
27 | Il tourne 900000 iterations et produit des vitbil_iter_900000 |
---|
28 | |
---|
29 | Pour visialiser une iteration : |
---|
30 | |
---|
31 | tail -n +5 vitbil_iter_100000 |awk '{print $2*$4}' |xyz2grd -ZBL -R-2797500/2992500/-2797500/2797500 -I15000 -Ghdot_100000.grd |
---|
32 | |
---|
33 | 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 |
---|
34 | |
---|
35 | H,masque,bm,hdot,uxmaj,uymaj,ux,uy |
---|
36 | |
---|
37 | Il ne semble pas y avoir un gros gain passé 300000 iterations |
---|
38 | |
---|
39 | Remarque importante |
---|
40 | |
---|
41 | Il y a des problemes dans les zones ou il y a beaucoup de Nunataks (typiquement le glacier Lambert) : |
---|
42 | |
---|
43 | * 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. |
---|
44 | |
---|
45 | |
---|
46 | Nouvelle Procedure : |
---|
47 | |
---|
48 | * On passe en flux en multipliant les vitesses merge (Ux-merge-LeBrocq-A15spn61.grd) par l'epaisseur |
---|
49 | - fluX-merge-LeBrocq-A15spn61.grd et fluY-merge-LeBrocq-A15spn61.grd |
---|
50 | |
---|
51 | * On met en Nan tous les 0 |
---|
52 | - grdmath fluX-merge-LeBrocq-A15spn61.grd 0 NAN = fluX-merge-LeBrocq-A15spn61-NAN.grd |
---|
53 | |
---|
54 | * On filtre avec une moyenne mobile sur 30 km |
---|
55 | - grdfilter fluY-merge-LeBrocq-A15spn61-NAN.grd -D0 -Fb30000 -GfluY-merge-LeBrocq-A15spn61-Fb30.grd |
---|
56 | |
---|
57 | * On calcule la direction |
---|
58 | - grdmath fluX-merge-LeBrocq-A15spn61-Fb30.grd fluY-merge-LeBrocq-A15spn61-Fb30.grd ATAN2 R2D = angle-flux-merge-Fb30.grd |
---|
59 | |
---|
60 | * On recalcule Ux et Uy en utilisant la direction lissee |
---|
61 | - grdmath angle-flux-merge-Fb30.grd 90 SUB COSD norme-merge-LeBrocq-A15spn61.grd MUL -9999 AND = Ux-merge-dirFb30.grd |
---|
62 | - grdmath angle-flux-merge-Fb30.grd -90 SUB SIND norme-merge-LeBrocq-A15spn61.grd MUL -9999 AND = Uy-merge-dirFb30.grd |
---|
63 | |
---|
64 | * On fait les ZBL |
---|
65 | - grd2xyz -ZBL Ux-merge-dirFb30.grd >Ux-merge-dirFb30_ZBL_15km.dat |
---|
66 | |
---|
67 | * et on re-paste le fichier input_SHbmMkULBq-ux-uy-15km.dat |
---|
68 | |
---|
69 | |
---|
70 | |
---|
71 | En utilisant les resultats d'un run GRISLI |
---|
72 | |
---|