wiki:ValidateOrchidee

Version 14 (modified by mmaipsl, 13 years ago) (diff)

--

Simulations

Voir IGCMGContinentIPSLCM5

  • Modifications dans 1.9.4_AR5 versus 1.9.4 : humrel dans la formule vbeta3, phenologie corrigée et harvest ajouté (récoltes). Ces modifs ont été introduites dans la branche "temporaire" : orchidee_1_9_4_AR5 et sont dans mod.def.
  • vbeta2 : modif retenue pour le couplé et commitée sur orchidee_1_9_4_AR5. 20 ans Off-line avec CHECK WATER BALANCE à TRUE dans l'hydrologie ont été effectués la et on montre la comparaison avec la version orchidee_1_9_4 ici et plus précisément
  • 51 ans avec land use ont tournés ici.
  • La modif 1264 de LMDZ doit permettre de changer de fichier en fin d'année.
  • Faire le même test en couplé (sans aerosols et en calendrier vrai)
  • A noter pour démarrer en 1995, il faut mettre VEGET_YEAR à 151 dans le fichier orchidee.def. Voir : [797]
  • On doit prendre le nouveau fichier de LAND USE pour l'AR5 : pour l'instant (à faire : déplacemenent/renommage dans IGCM) /dmnfs/cont003/p86manci/IGCM/BC/SRF/OL2/ORCHIDEE_PFTmap_1850to2005_AR5_LUHa.rc2.nc
  • Monitoring : attention à la gestion des points mixtes continents/océan dans les calculs de bilan (il faut rajouter CONTFRAC dans les MONITORINGs SRF). Pour l'instant, il vaut mieux regarder les cartes ATM (justes) plutot que les cartes SRF.

DGVM + LAND USE

Test "SpF" 1950 CRU 1° with or not LAND USE and with or not DGVM after a spinup equilibrium in a zoom in the south west of France.

See DODS LAND-USE+DGVM.

  • SpFORC_1 : spinup run. spinup algorithm is here :
    1. sechiba and stomate : duree_inistomate=10
    2. Loop configuration for spin-up : n_iter=1
    3. orchidee (ORC_1) : duree_sechiba=150
    4. teststomate : duree_stomate=0
    5. forcesoil : duree_carbonsol=10000
    6. Finalization : duree_final=50
  • SpFHIST : 10 Years historical run after the spinup without LAND USE nor DGVM.
  • SpFHISTLU : 10 Years historical run with LAND USE.
  • SpFHISTLUDGVM : same with LAND USE and DGVM.
  • SpFHISTLUDGVM2 : same with DGVM with last correction of N. Viovy on tmin_crit and tcm_crit.

ATLAS des différences entre les versions

Etude couplé versus forcé LMDzOR versus Off-line : Dans le DODS ssipsl, on pourra trouver à la page
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/ETUDE_COUPLE_FORCE
les ATLAS de différences normalisés suivant.

ORCHIDEE OFF-LINE

Dans ce paragraphe, on trouvera des ATLAS de comparaison sur la seconde décennie d'une simulation off-line ORCHIDEE seul entre la version orchidee_1_9_4 et la version orchidee_1_9_4_AR5 avec démarrage de la phénologie corrigé.

LMDzOR

Dans ce paragraphe, on trouvera des ATLAS de comparaison entre différentes versions de LMDzOR et des forcés ORCHIDEE off-line.

  • Comparaison (seconde décénie) Tag 1.6 :
    LMDzOR (versions LMDZ4_V3_4 et orchidee_1_9_4) versus LMDzOR (versions LMDZ4_V3_4 et orchidee_1_6) (STOMATE activé, RVEG à 1) dans
    ATLAS-DIFF-INITSS-LUCIDpdINIT
  • Comparaison (seconde décénie) on-line - off-line avec STOMATE :
    LMDzOR (versions LMDZ4_V3_4 et orchidee_1_9_4) versus ORCHIDEE off-line (STOMATE activé, RVEG à 1) dans
    ATLAS-DIFF-INITSS-OOINITSS_NEW
  • Comparaison (seconde décénie) on-line - off-line avec STOMATE et RVEG modifié :
    LMDzOR (versions LMDZ4_V3_4 et orchidee_1_9_4) versus ORCHIDEE off-line (STOMATE activé, RVEG à 0.5) dans
    ATLAS-DIFF-INITSSRV-OOINITSSRV_NEW

IPSLCM4_v1

Dans ce paragraphe, on trouvera des ATLAS de différences entre des configuration de l'ancien et du nouveau couplé.

Le couplé récent ST06 utilisé pour ces comparaisons est décris à la page ContinentIPSLCM5.

IPSLCM5

Dans ce paragraphe, on compare le couplé ST06 avec des runs couplés et forcés récents :

LMDzOR

En LMDZ4OR_V2 : (juste SECHIBA)
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSEC/MONITORING/
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSECKO/MONITORING/

Vous avez :

  • Les "spinups" INIT d'initialisation (sans ATLAS) : 20 ans avec l'année

PFT=1992, HADISST=1972

INITSS avec STOMATE et INITSEC juste SECHIBA. (sans ATLAS ...)

  • Les runs PD005 : on repart de la dernière année INIT et on fait 30 ans

avec HADISST=73-

SEC et SECR => juste avec SECHIBA sans et avec les nouveaux RESTARTS SS et SSR => la même chose avec STOMATE en +

Les Simulations

  • LMDZ4OR_v2 sans routage :

http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSR/MONITORING/
à comparer à avec routage :
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSS/MONITORING/

  • RVEG=0.5 (au lieu de 1 en COUPLE LMDZ4OR_v2) :

http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSRV/MONITORING/
à comparer avec
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSS/MONITORING/

  • OK_CO2=n SECHIBA

http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSECKO/MONITORING/
à comparer avec
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSEC/MONITORING/

Simulations with 0.5 for RVEG_PFT for class of vegetations (and other at 1) : http://dods.idris.fr/rces533/LMDZ4OR_v2/?C=M;O=Dp

  • INITSSTr

Tr : 2 & 3 (Tropiquaux) à 0.5
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSTr/

  • INITSSBo

Bo : 8 & 9 (Boreaux) à 0.5
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSBo/

  • INITSSTe

Te : 4,5 & 6 (Temperes) à 0.5
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSTe/

  • INITSSHe

He : 10,11,12 & 13 (Herbaces) à 0.5
http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSHe/

  • La nouvelle transpiration de Nicolas avec LMDZ4OR_v2 peut être regardée là :

http://dods.idris.fr/rces533/LMDZ4OR_v2/INITSSV/

Les ATLAS de diff

entre ces simulations sont là :
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/

On trouve par exemple les différences suivantes :

On peut aussi les télécharger en local :

L'atlas comparant ORCHIDEE 1.9.4 On-line et Off-line est celui-ci :
ORCHIDEE on-line - off-line

Spinups

To do so, we first run the model at a resolution of 2°×2° until the C pools reach equilibrium based on transient climate data during the period of 1901-1910 and the 1860 atmospheric CO2 concentration of 286.05 ppm. The model was then run to 1900 with the variable climate of the period 1901-1910 and corresponding atmospheric CO2 concentration data during 1860-1900. This state was used as the initial condition

The meteorological data used to run ORCHIDEE include air temperature, precipitation, wet day frequency, diurnal temperature range, cloud cover, relative humidity of the air, and wind speed. Monthly data sets, with a spatial resolution of 0.5°×0.5° for 1901-2002, were supplied by the Climatic Research Unit (CRU), School of Environmental Sciences, University of East Anglia, U.K. (Mitchell and Jones, 2005). Data on annual atmospheric CO2 concentration were obtained from Rayner et al. (2005).

test leafagecrit_tab[pft=12-13] =

  • 90 : Nouveau driver : (terminé)

http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb1ORC_1/MONITORING/

Ancien driver : (terminé)

http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb2ORC_1/MONITORING/

  • 120 : Nouveau driver : (en cours)

http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb3ORC_1/MONITORING/

Ancien driver : (terminé)

http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarbORC_1/MONITORING/

après FORCESOIL, 50 ans ORCHIDEE final :
leafagecrit_tab[pft=12-13] =

  • 90 : Nouveau driver : (terminé)

http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb1FIN/MONITORING/

  • 120 : Ancien driver : (terminé)

http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarbFIN/MONITORING/

Appel slowproc à la fin du pas de temps : (terminé)
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarbiORC_1/MONITORING/
à comparer avec
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb1ORC_1/MONITORING/

Enfin j'ai lancé un SPINUP SANS PARTIR des RESTART de SHILONG : (en cours)
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb4ORC_1/MONITORING/
à comparer avec
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb1ORC_1/MONITORING/

J'ai lancé un spinup 20, 2°, en séparant les composantes SECHIBA et STOMATE pour bien voir les MONITORING carbone. C'est là :
/work/cont003/p86manci/AGRI90/modipsl/modeles/ORCHIDEE_OL/SPINUP_CARBON_VALIDE_deux
http://dods.extra.cea.fr/data/p86manci/OL2/V90deuxORC_1*-9/MONITORING/

La simulation V90deux est quasi à l'équilibre : 0.09106 GtC/an.
Voici la figure
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/V90deux/NEP-SPINUP-1_66.gif

du SPINUP de la NEP (ou NEE ou CO2FLUX, c'est selon ;-))...

Enfin nous avons comparé le CO2FLUX (HF) au CO2FLUX_MONTHLY qui donne le même résultat sur la figure
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/V90deux/NEP-CO2FLUX-CO2FLUX_MONTHLY-66.gif .
Vous trouvez aussi les stocks de C dans

  • la végétation (519 GtC)

http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/V90deux/TOTAL_M-cveg-SPINUP.gif

  • et le sol (1295 GtC)

http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/V90deux/TOTAL_SOIL_CARB-csoil-SPINUP.gif .

C'est mieux par rapport au SPINUP de LOOP ! (stocks C dans la végèt = 280GtC et dans le sol = 1070 GtC).
Et donc nous estimons qu'il est préférable de lancer un premier HISTORIQUE maintenant avec le calcul du FAPAR et de voir si un forcesoil sur 365 jours et non 12 mois améliore ces stocks de carbone.

SECHIBA :
300 ans OK_CO2=n
http://dods.extra.cea.fr/data/p86manci/OL2/SpinSECKOORC_1/MONITORING/
800 ans OK_CO2=y
http://dods.extra.cea.fr/data/p86manci/OL2/SpinSECORC_1/MONITORING/


Changement de driver : deux runs avec les mêmes paramètre (le 2 est l'ANCIEN driver)
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb1ORC_1/MONITORING/
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarb2ORC_1/MONITORING/
Ca change beaucoup !!
http://igcmg.lsce.ipsl.fr/monitoring/tmp/fegg_plot01_Sd9Bev_prod/


slowproc inversé : (même paramètres, même nouveau driver que SpinCarb1ORC_1
http://dods.extra.cea.fr/data/p86manci/OL2/SpinCarbiORC_1/MONITORING/
c'est vraiment pareil en forcé, avec STOMATE :
http://igcmg.lsce.ipsl.fr/monitoring/prod.php?name=tmp/fegg_plot01_9IoWz7


En préparation (qui ne tourne pas encore :-( ) :
avec l'ancien driver /OL2/SSCARB
/work/cont003/p86manci/PB_WEATHER/modipsl/modeles/ORCHIDEE_OL/OOL_SEC_STO_CARBON
avec le nouveau et le paramètre à 90 : /OL2/SSCARB1
/work/cont003/p86manci/ORCHIDEE_1_9_4/modipsl/modeles/ORCHIDEE_OL/OOL_SEC_STO_CARBON

Couplé

IPSLCM4_v2

IPSLCM5

SEE ALSO : IGCMGContinentIPSLCM5

Arnaud Caubel a testé le tag ORCHIDEE_1_9_3 dans le couple IPSLCM5 en resolution 96x71 sur une simu de 50 ans.

1) La simu s'appelle : VAORC193. Les resultats sont la : http://dods.extra.cea.fr/data/p86caub/IPSLCM5/VAORC193 2) Cette simu est a comparer avec la simu de reference IPSLCM5 BACKUP.

Les seules differences entre BACKUP et VAORC193 sont

  • les tag orchidee : 1_9_2 pour BACKUP et 1_9_3 pour VAORC193
  • les orchidee.def

Des atlas de comparaisons entre les 2 simus sont la :
http://dods.extra.cea.fr/data/p86mart/ANALYSES/IPSLCM5/BACKUP_VAORC193
Les comparaisons de monitoring entre les 2 simus sont la :
http://dods.extra.cea.fr/data/p86caub/MONITORING/fegg_plot01_zO3DlJ_prod/

Old comparaisons

For comparaison of OLD ORCHIDEE version 1.4, 1.5, 1.6 1.8 and 1.9 (SEQUENTIAL an PARALLEL), please see :
http://dods.ipsl.jussieu.fr/orchidee/COMPARAISONS_ORCHIDEE_VERSIONS/

  • ATLAS_DIFF_1.4-1.5_WG_4
  • ATLAS_DIFF_1.4-1.6
  • ATLAS_DIFF_1.6_OLD_SA_corr_1.8_OLD_SA_corr
  • ATLAS_DIFF_1.8_AV_corr_1.9_AV
  • ATLAS_DIFF_1.8_OLD_SA_corr_NEW_SA_1.9
  • ATLAS_DIFF_1.8_corr_AV_1.6_corr_SA
  • ATLAS_DIFF_1.9_SEQ_PARA_SA
  • ATLAS_DIFF_NEW_1.9_AV-SA
  • ATLAS_DIFF_vbeta1

Most of all, see README file in this directory.

List of Nathalie De Noblet simulations :
http://dods.ipsl.jussieu.fr/orchidee/Liste_Simuls.doc

FLUXNET sites

Old runs on FLUXNET sites :
http://dods.ipsl.jussieu.fr/orchidee/ATLAS/

  • atlas_AB_SECHIBA
  • atlas_CP_SECHIBA
  • atlas_Cabauw_SECHIBA
  • atlas_FL_SECHIBA
  • atlas_LW_SECHIBA
  • atlas_MA_SECHIBA
  • atlas_PO_SECHIBA
  • atlas_SH_SECHIBA
  • atlas_mc2_3
  • atlas_t2m.alb.spread
  • atlas_t2m.normal
  • atlas_t2m.swnet
  • atlas_t2m.swnet.spread
  • atlas_t30m.swnet
  • atlas_t30m.swnet.spread

Comparaisons of ORCHIDEE versions over well known Harvard (HV) FLUXNET site :
http://dods.ipsl.jussieu.fr/orchidee/SANORCHIDEE/SORTIE_FLUXNET/

  • HV13s
  • HV13sFIN
  • HV16hs
  • HV16hsFIN
  • HV193s
  • HV193sFIN
  • HVMPLAI
  • HVMPLAIFIN
  • HVOKCO2
  • HVOKCO2FIN
  • HVOLD
  • HVOLDFIN
  • HVSECH
  • HVSECHFIN
  • HVSECH_LAI
  • HVSECH_LAIFIN
  • HVSECH_LAI_CW
  • HVSECH_LAI_CWFIN
  • HVSECH_LAI_CW_T
  • HVSECH_LAI_CW_TFIN
  • HVSECH_LAI_T
  • HVSECH_LAI_TFIN

See protocol of runs in EXPLICATION file in this directory.

Differencies of SPRED_SPREC parameter

This parameter give (in Offline ORCHIDEE runs) the frequency of the precipitation over the forcing time step.
http://dods.ipsl.jussieu.fr/orchidee/COMPARAISONS_TESTS_SPRED_PREC/

  • ATLAS_DIFF_SPRED_1_3
  • ATLAS_DIFF_SPRED_1_7
  • ATLAS_DIFF_SPRED_1_9
  • ATLAS_DIFF_SPRED_1_12

See protocol of runs in EXPLICATION.txt file in this directory.

Programm of simulations for AR5 validations

1) globaux / OFF-LINE (force par ??? sur 20 ans par exemple)
1-1 les 2 versions IPCC-AR4 vs IPCC-AR5 (sequentiel pour le moment). LAI lu, ancienne carte de vegetation.
1-2 meme chose que dans 1-1 mais avec ok_stomate (toujours ancienne carte de vegetation).
1-3 comme 1-2 mais avec version IPCC-AR5 seulement et nouvelle carte de vegetation (pour montrer que l'effet de cette carte n'est pas dramatique).

2) globaux / ON-LINE / LMDZOR (version LMDZ ?) / CDRAG calcule dans LMDZ
2-1 les 2 versions avec LAI lu et ancienne carte
2-2 les 2 versions avec ok-STOMATE=yes
2-3 juste AR5 avec nouvelle carte de vegetation
Dans les 2 cas, la version OFF-LINE sera celle qui sera utilisee comme 'donnee' pour la validation et la comparaison dans les atlas quand on ne dispose pas de donnees.

3) global ... comme dans 2-3, mais avec CDRAG calcule dans ORCHIDEE
donc la version off-line est la meme que dans le point 2-3. Cela nous permettra de voir la sensibilite au cdrag. Car en fait on ne compare pas vraiment des choses comparables puisque certains 'processus' ne sont pas calcules de la meme facon on-line et off-line et c'est un point que j'aimerai soulever au cspole.

4) FLUXNET (apres avoir revu quels sites utiliser avec P. Ciais)
4-1 faire AR4 & AR5 avec LAI lu
4-2 idem mais avec ok_stomate.
Je ne sais pas s'il est interessant de mettre les 4 sur le meme graphique, sinon : AR4 vs AR5 (pour LAIlu, puis ok_stomate), AR4 : LAIlu vs ok_stomate et meme chose pour AR5 (cela va nous permettre de voir si ok_stomate ameliore le modele).

Images

quelques graphes : !VegetFrac1 = PFT 1 à la fin du SPINUP (ORCHIDEE 800 ans)

  • pour l'ancien (SpinCarbORC_1)

http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/VegetFrac1_SpinCarbORC_1.gif

  • et le nouveau (!SpinCarb1ORC_1)

http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/VegetFrac1_SpinCarb1ORC_1.gif

et le MaxVegetFrac pour comparer le sol nu
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/MaxVegetFrac1_SpinCarbORC_1.gif
et le sol avec "plus de plantes car elles sont toutes mortes" !

Le restart de SHILONG (veget_SHILONG_rest_1)
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/veget_SHILONG_rest_1.gif
est quand même différent de celui (veget_194_rest) obtenu
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/veget_194_rest_1.gif
à la fin du !SpinCarbFIN (après forcesoil et 50ans).

Enfin le fracnobio
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/fracnobio.gif
est OK mais le snownobio
http://dods.ipsl.jussieu.fr/orchidee/VALIDS_194/snow_nobio_tot_restart_SHILONG.gif
n'est pas du tout controlé dans SECHIBA (toutes les simus) et peut-être très positif ou très négatif !! en off-line, c'est embêtant !