Tags/196/ExternalisationParameters: externalisation

Nathalie
1) Reperer tous les parametres potentiellement modifiables du code et les en extraire pour pouvoir
les modifier dans le run.def sans avoir a recompiler le modele
- have a set of PFTs that are associated to potentially different physics/equations but as restricted as possible
==> exemples Ciais
2) pour les PFTs existants, spatialiser les parametres de facon a pouvoir commencer a prendre en compte une notion de 'varietes'. (Voir comment on articule ca avec la lecture de la carte de vegetation)
3) pour les PFTs existants, definir de nouveaux parametres 'descriptifs' permettant de les definir en externe. Par exemple 'tree', 'C3', 'deciduous', ..... (Voir comment on articule ca avec la lecture de la carte de vegetation)
4) Toujours pour les PFTs existants, autoriser le choix du nombre de PFTs (par exemple 11 si pas d'agriculture, 13 si agriculture) a etre modifie en fonction de la simulation.
5) Et enfin pouvoir augmenter le nombre de PFTs mais cela releve tout d'abord des travaux des uns et des autres.

Nicolas
Il serait peut être sympa de pouvoir rajouter la possibilité de lire une carte (spatialisée) de
correspondance entre les PFTs du modÚles et ceux de la carte de végétation lue. Comme cela on
pourrait par exemple continuer à utiliser une carte 13 PFT existante mais définir que le PFT culture
C4 (13 de la carte de veg) en Europe c'est du maïs (par ex pft 15 du modÚle) et du mil en Afrique
(par ex pft 16 du modÚle).

les fichiers textes existent déjà : les orchidee.def
=> on peut (doit !) les découper. 2 propositions :
1) orchidee_ol.def, sechiba.def, stomate.def
2) forcing_orchidee.def, physic_orchidee.def

Pour externalisé les paramÚtres :
1) Nombre de PFT flottant => quid de la carte de végÚt ?? => plusieurs jeux de carte de végÚt ?
2) Les paramÚtres choisit sont transformés 1D (/pft) en 2D (/pts/pft) => spacialisation
3) Mise en oeuvre info :
On refait presque tout getincom !! ==> transformation des listes de paramÚtre en xml ??
===> On lit le paramÚtre toujours comme une chaîne de caractÚre 

exemples de getin :
i) RVEG_PFT = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.

ii) SOILTYPE_FILE = ../surfmap/soils_param.nc

iii) 
SECHIBA_VEG__01 = 0.2
SECHIBA_VEG__02 = 0.0
SECHIBA_VEG__03 = 0.0
SECHIBA_VEG__04 = 0.0
SECHIBA_VEG__05 = 0.0
SECHIBA_VEG__06 = 0.0
SECHIBA_VEG__07 = 0.0
SECHIBA_VEG__08 = 0.0
SECHIBA_VEG__09 = 0.0
SECHIBA_VEG__10 = 0.8
SECHIBA_VEG__11 = 0.0
SECHIBA_VEG__12 = 0.0
SECHIBA_VEG__13 = 0.0


SECHIBA_VEG__06 = 0.3
SECHIBA_VEG = 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13, 1./13

SECHIBA_VEG = ../fichier_SECHIBA_VEG.nc
ou
SECHIBA_VEG__06 = ../fichier_SECHIBA_VEG__06.nc

spacialisation => 
DIMENSION(kjpindex,nvm) :: SECHIBA_VEG
sinon 
DIMENSION(1,nvm) :: SECHIBA_VEG
DIM1_SECHIBA_VEG=kjpindex/1

toto(1:kjpindex,1:nvm)=tata(1:kjpindex,1:nvm)*SECHIBA_VEG(1:nvm) 
=>
toto(1:kjpindex,1:nvm)=tata(1:kjpindex,1:nvm)*SECHIBA_VEG(1:DIM1_SECHIBA_VEG,1:nvm) 

Q10 ??
DPU

 Fonction englobant 
a) un getin + une lecture de carte normalisée ( projection + valeur en cas de missing ==> flag qui
bétonne les missings -> plantages)
b) valeur par défaut / valeur par variétés==méta-classe (nom à trouver) => tableau de correspondance à rajouter
==> On lit le paramÚtre toujours comme une chaîne de caractÚre 
On initialise avec les valeurs par défaut (/variétés ?)
Si le getin est une valeur simple, on écrit tout le tableau avec cette valeur
Si le getin est une liste de valeur => on écrit par pft donnée
Si le getin est un chemin vers un fichier : lecture de cette carte + utilisation de la projection
   => quid des PFTs manquantes dans la carte (variété ?)
   => quid des points de grille manquants ?? 
vcmax_op__18="MonfichierPFT18"
vcmax_op="MonfichierPFT18"

Flant Functionnal Type  => PFT
Plant Functionnal class => PFC   
Meta Class of Vegetation
BioGechemical Class

====> restarter ces paramÚtres ??!!


implications :
1) changer le nombre de PFT => supprimer tous les paramÚtres tableaux par PFT pour les relire
   pendant le run
2) les méta-classes (les PFTs actuelles ??) => définir un tableau de correspondance systématique et
initialiser tous les tableaux de sous-classe avec une fonction utilisant systématiquement ce
tableau.

exemple :
=========
  vcmax_opt_m(2:nvm_m) =     &
 & (/      65.,     65.,     35.,     45.,     55.,     35.,      &
 &         45.,     35.,     70.,     70.,    70.,    70.      /)

appel de la fct de dim pour récupérer dim1_vcmax_opt
allocation 
ALLOCATE( vcmax_opt(1:dim1_vcmax_opt,2:nvm), ERR )
initialisation :
DO j=1:dim1_vcmax_opt
  vcmax_opt(j,2:nvm) = Transforma_meta_class(vcmax_opt_m(2:nvm_m))
ENDDO
CALL lecture de carte :
  CALL getin_p("vcmax_opt", vcmax_opt)