ETAPE 0 (rien a changer ici)
=======
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
OK

=======================================================================
ETAPE 1 (a refaire pour toute nouvelle grille modele)
========
Cette etape consiste a interpoler sur la  grille ORCA2 (aux  points T)
les champs B, SIGMAsed, qui ont ete crees dans un meme fichier netcdf
lors de l'etape 0  (fichier condmag.nc).
Pour ce faire, on suit pas a pas la doc de la routine SAXO 
Voir toute la DOC de cette etape dans:
http://www.lodyc.jussieu.fr/opa/INTERP/index.html

On lance le programme comme ceci (cf la doc)
idl
@init_interp
interpolation

On va realiser 2 runs successifs de la routine d'interpolation: un
passage pour le champ magnetique B (un fichier Br_ORCA2.nc est genere)
et un deuxieme run  pour le champ de conductance (un  fichier
cond_sed_ORCA2.nc est cree)

Il faut donc faire quelque changements de configuration dans le
fichiers suivant:
>> editer interp.pro=> changer le nom de la variable du  fichier condmag.nc a interpoler (fich cree sous matlab), L91:
   zdata=ncdf_lec('/usr/work/sur/fvi/OPA/geomag/condmag.nc',var='Br')
>>editer naminterp.pro
L67: ;   data_name = 'cond_sed'

>>init_path.pro
L44   outputfile = 'cond_sed_ORCA2.nc'

Le meshmask est different pour chaque grille model, penser a pointer
vers le bon fichier (idem pour l'etape 3 plus bas)

==========================================================================
ETAPE 2: (a refaire pour toute nouvelle grille modele)
========
C'est l'etape de calcul des forcages geomagnetiques a partir des
sorties ORCA et des fichiers Br et conductance. Le calcul est fait sur
la grille modele pour que le calcul de divergence soit aussi precis
que possible. Ce point est vraiment critique.

Techniquement, tout se passe dans le fichier  forcagequimarche.pro
On lance  idl puis:

@init2
.r forcagequimarche.pro
forcagequimarche,'1993','01'

et idem pour les autres annees: forcagequimarche,'1994','01'


Pour info, voila ce que fait le programme:
-----------------------------------------
1-  calcul d'une matrice 3D de conductivite sigma=f(T,S) en fonction
de T,S de ORCA2, sur la grille T.  Formule assez  empirique, en
siemens/m que j'ai pondu apres regression: 
sigma(T,S)=0.02047780622061 + 0.00273147624197*T + 0.00035133182334*T^2 + 0.09139808809909*S + 0.00241425798890*S*T -0.00023998958774*S^2

2- Calcul du champ 2D de conductance de l'ocean par integration verticale 
        /
SIGMAoc=| sigma dz
        /

3- definir SIGMA=SIGMAsed+SIGMAoc


4-  calculer 
           /
Bustar= B* | sigma * u /SIGMA dz
           /
           /
Bvstar= B* | sigma * v /SIGMA dz
           /

5- calculer div(Bustar,Bvstar)  ---> utilise les  routines SAXO
(tres legerement modifiees ici pour certaines d'entre elles (divfred,
gradfred)



========================================================================
ETAPE 3 (a refaire pour chaque nouvelle grille)
========
On fait une Reverse  interpolation de la grille modele vers la grille
geographique des fichiers crees a l'etape precedente (Divergence,
Sigma, annee par annee)
Pour cette etape, suivre pas a pas l'aide en ligne  sur:
http://www.lodyc.jussieu.fr/opa/REVTERP/

Il faut penser a:
1- editer /usr/work/sur/fvi/OPA/geomag/INTERP2/naminterp2.pro
et verifier le nom de la variable que l'on veut interpoler: p. ex
name_input_field='Sigma'
2- editer  /usr/work/sur/fvi/OPA/geomag/INTERP2/init_path2.pro  
changer les  datafile (sortie de l'etape 2), outputfile et le cas
echeant   input_dir, output_dir inputmesh weightfile

Pour lancer, aller dans le repertoire INTERP2:
cd /usr/work/sur/fvi/OPA/geomag/INTERP2

Puis:
> idl
> @init
> interp2