1 | ETAPE 0 (rien a changer ici) |
---|
2 | ======= |
---|
3 | 1) lire B (champ magnetique), et les conductances: cond_sed, cond_oc, cond_tot, dans le fichier condmag.nc- ce fichier a ete genere sous matlab M_lib4/geomag/MOED2D/ORCA_writenetcdf.m |
---|
4 | OK |
---|
5 | |
---|
6 | ======================================================================= |
---|
7 | ETAPE 1 (a refaire pour toute nouvelle grille modele) |
---|
8 | ======== |
---|
9 | Cette etape consiste a interpoler sur la grille ORCA2 (aux points T) |
---|
10 | les champs B, SIGMAsed, qui ont ete crees dans un meme fichier netcdf |
---|
11 | lors de l'etape 0 (fichier condmag.nc). |
---|
12 | Pour ce faire, on suit pas a pas la doc de la routine SAXO |
---|
13 | Voir toute la DOC de cette etape dans: |
---|
14 | http://www.lodyc.jussieu.fr/opa/INTERP/index.html |
---|
15 | |
---|
16 | On lance le programme comme ceci (cf la doc) |
---|
17 | idl |
---|
18 | @init_interp |
---|
19 | interpolation |
---|
20 | |
---|
21 | On va realiser 2 runs successifs de la routine d'interpolation: un |
---|
22 | passage pour le champ magnetique B (un fichier Br_ORCA2.nc est genere) |
---|
23 | et un deuxieme run pour le champ de conductance (un fichier |
---|
24 | cond_sed_ORCA2.nc est cree) |
---|
25 | |
---|
26 | Il faut donc faire quelque changements de configuration dans le |
---|
27 | fichiers suivant: |
---|
28 | >> editer interp.pro=> changer le nom de la variable du fichier condmag.nc a interpoler (fich cree sous matlab), L91: |
---|
29 | zdata=ncdf_lec('/usr/work/sur/fvi/OPA/geomag/condmag.nc',var='Br') |
---|
30 | >>editer naminterp.pro |
---|
31 | L67: ; data_name = 'cond_sed' |
---|
32 | |
---|
33 | >>init_path.pro |
---|
34 | L44 outputfile = 'cond_sed_ORCA2.nc' |
---|
35 | |
---|
36 | Le meshmask est different pour chaque grille model, penser a pointer |
---|
37 | vers le bon fichier (idem pour l'etape 3 plus bas) |
---|
38 | |
---|
39 | ========================================================================== |
---|
40 | ETAPE 2: (a refaire pour toute nouvelle grille modele) |
---|
41 | ======== |
---|
42 | C'est l'etape de calcul des forcages geomagnetiques a partir des |
---|
43 | sorties ORCA et des fichiers Br et conductance. Le calcul est fait sur |
---|
44 | la grille modele pour que le calcul de divergence soit aussi precis |
---|
45 | que possible. Ce point est vraiment critique. |
---|
46 | |
---|
47 | Techniquement, tout se passe dans le fichier forcagequimarche.pro |
---|
48 | On lance idl puis: |
---|
49 | |
---|
50 | @init2 |
---|
51 | .r forcagequimarche.pro |
---|
52 | forcagequimarche,'1993','01' |
---|
53 | |
---|
54 | et idem pour les autres annees: forcagequimarche,'1994','01' |
---|
55 | |
---|
56 | |
---|
57 | Pour info, voila ce que fait le programme: |
---|
58 | ----------------------------------------- |
---|
59 | 1- calcul d'une matrice 3D de conductivite sigma=f(T,S) en fonction |
---|
60 | de T,S de ORCA2, sur la grille T. Formule assez empirique, en |
---|
61 | siemens/m que j'ai pondu apres regression: |
---|
62 | sigma(T,S)=0.02047780622061 + 0.00273147624197*T + 0.00035133182334*T^2 + 0.09139808809909*S + 0.00241425798890*S*T -0.00023998958774*S^2 |
---|
63 | |
---|
64 | 2- Calcul du champ 2D de conductance de l'ocean par integration verticale |
---|
65 | / |
---|
66 | SIGMAoc=| sigma dz |
---|
67 | / |
---|
68 | |
---|
69 | 3- definir SIGMA=SIGMAsed+SIGMAoc |
---|
70 | |
---|
71 | |
---|
72 | 4- calculer |
---|
73 | / |
---|
74 | Bustar= B* | sigma * u /SIGMA dz |
---|
75 | / |
---|
76 | / |
---|
77 | Bvstar= B* | sigma * v /SIGMA dz |
---|
78 | / |
---|
79 | |
---|
80 | 5- calculer div(Bustar,Bvstar) ---> utilise les routines SAXO |
---|
81 | (tres legerement modifiees ici pour certaines d'entre elles (divfred, |
---|
82 | gradfred) |
---|
83 | |
---|
84 | |
---|
85 | |
---|
86 | ======================================================================== |
---|
87 | ETAPE 3 (a refaire pour chaque nouvelle grille) |
---|
88 | ======== |
---|
89 | On fait une Reverse interpolation de la grille modele vers la grille |
---|
90 | geographique des fichiers crees a l'etape precedente (Divergence, |
---|
91 | Sigma, annee par annee) |
---|
92 | Pour cette etape, suivre pas a pas l'aide en ligne sur: |
---|
93 | http://www.lodyc.jussieu.fr/opa/REVTERP/ |
---|
94 | |
---|
95 | Il faut penser a: |
---|
96 | 1- editer /usr/work/sur/fvi/OPA/geomag/INTERP2/naminterp2.pro |
---|
97 | et verifier le nom de la variable que l'on veut interpoler: p. ex |
---|
98 | name_input_field='Sigma' |
---|
99 | 2- editer /usr/work/sur/fvi/OPA/geomag/INTERP2/init_path2.pro |
---|
100 | changer les datafile (sortie de l'etape 2), outputfile et le cas |
---|
101 | echeant input_dir, output_dir inputmesh weightfile |
---|
102 | |
---|
103 | Pour lancer, aller dans le repertoire INTERP2: |
---|
104 | cd /usr/work/sur/fvi/OPA/geomag/INTERP2 |
---|
105 | |
---|
106 | Puis: |
---|
107 | > idl |
---|
108 | > @init |
---|
109 | > interp2 |
---|
110 | |
---|
111 | |
---|
112 | |
---|