1 | Les codes sont dans ~MAPS/Calcul-vit-bil |
---|
2 | |
---|
3 | vitbil_LeBrocq_to_Grisli-15km.f90 |
---|
4 | |
---|
5 | - Le principe est de partir d'un champ estime et de le corriger en ne |
---|
6 | touchant que les flux sortants. Le flux de départ est sur les noeuds majeurs. |
---|
7 | |
---|
8 | - A partir de ce flux, on calcule quelles sont les faces sortantes. |
---|
9 | - quand toutes les vitesses sont nulles, on devrait supposer que c'est |
---|
10 | un point sortant avec 1/4 sur chaque face. |
---|
11 | XpartW(i,j)=0.25 ... |
---|
12 | |
---|
13 | - Ce qui reste conserve sur chaque face est : |
---|
14 | (flux sortant par cette face)/(flux total sortant) |
---|
15 | |
---|
16 | - Ca correspond a fixer l'azimut du flux sortant |
---|
17 | |
---|
18 | Comme champ de départ |
---|
19 | ------------------------ |
---|
20 | on peut partir des champs "pente" ou "driving |
---|
21 | stress" tels qu'ils sont calcules par La Brocq. |
---|
22 | |
---|
23 | |
---|
24 | On peut aussi utiliser les vitesses observees (E. Rignot). |
---|
25 | Avantage, on peut faire les shelves |
---|
26 | |
---|
27 | La difficulte est de melanger les deux. |
---|
28 | |
---|
29 | Approche travailler sur les azimuths : |
---|
30 | -------------------------------------- |
---|
31 | - comparer les cartes d'azimuth dans les deux cas. |
---|
32 | - prendre ceux de Rignot où ils existent |
---|
33 | - prendre ceux de Le Brocq sinon |
---|
34 | - interpoler ceux de Rignot sur les ice shelves |
---|
35 | - supposer perpendiculaire au front pour les bords d'ice shelves |
---|
36 | |
---|
37 | Raffinement "outcrop". |
---|
38 | ----------------------- |
---|
39 | |
---|
40 | Comment définir un outcrop ? Attention, il faudrait une définition qui soit compatible avec |
---|
41 | la méthode qui sera utilisée dans Grisli pour les traiter. |
---|
42 | |
---|
43 | Prendre tous les noeuds H < 10. |
---|
44 | -> si tous les noeuds autour sont plus bas -> outcrop |
---|
45 | -> si les noeuds autour qui sont plus haut |
---|
46 | |
---|
47 | les noeuds au dessus sont aussi des outcrop. |
---|
48 | -> noeuds qui ne font que donner sont d'abord tagges outcrop. |
---|
49 | |
---|