Rq +3°C par doublement de CO2 (attention a coherence SST et CO2 si on veut avoir tt le deltaT)

Rq Sveltlana a fait 0 ; 15 Ma ; 30Ma ; 42 Ma ; 55Ma avec des CO2 variables (x0 x2 x3 x4 et cold or warm orbits)
/ccc/work/cont003/dsm/p25botsy/LMDZ/modipsl/config/LMDZOR_v5/CR55Ma4x/

Eemien de Pascale runs couplés avec veget landuse de present , sans doute fichiers de lancements ont été perdus dans migration
dans ces runs c le changement d'insolation qui créé le chgmt C (sans doute pas de changement de CO2 pris en compte). Ces chgmts d'insolation créent un gradient interhémisphérique et des changements de saisonnalité tres forts.

Pliocene (mails de Camille CONTOUX 27/02/2018)
“PlioMIP1” date de 2012 (voir Contoux et al. 2012).
simu LMDZOR atmosphere-only Pliocène :
/ccc/store/cont003/dsm/p25camil/IGCM_OUT/LMDZOR/PROD/CLIM/Pcomp4
Normalement les restarts sont encore là, j'ai vérifié seulement dans ATM. Il s'agit d'une simu LMDZ5 en 96x95x39, avec LAI interactif. Si c'est juste pour des tests, tu peux prendre cela, c'est la seule simu “atmosphere-only” qui ait fait partie d'un MIP.
Les SSTs et la végétation de la simul LMDZOR Pcomp4 ont été générées à partir de bases de données.

Si tu veux vraiment regarder la chimie au Pliocène, il vaut mieux prendre les restart d'une simu couplée, les températures de l'océan seront plus réalistes.
Dans ce cas tu peux aller voir là :
/ccc/store/cont003/dsm/p25camil/IGCM_OUT/IPSLCM5A/DEVT/pliocene/Pcalt4/RESTART
lancé depuis /ccc/work/cont003/dsm/p25camil/work_titane/work/cont003/p25camil/modipsl/config/IPSLCM5A/Pcalt4REDO/

En gros la simul LMDZOR est plus chaude que la simu CM5A car les SSTS imposées viennent d'une carte de SST construite à partir de données. La végét imposée est la même pour les deux simus, et est également une carte construite à partir de données.
Si tu veux faire un truc “propre”, je te conseille de prendre plutôt les sorties du couplé, car les sorties de LMDZOR n'ont pas vraiment été exploitées depuis, et des problèmes subsistent concernant la base de données avec laquelle a été construite la carte des SSTs imposées. Du coup si un jour tu veux le publier, tu risquerais de te faire embêter la-dessus. Avec la sortie du couplé, c'est béton, elle a été exploitée et décrite dans de nombreux papiers issus de l'intercomparaison.
Par contre, dans LMDZOR tu as un réchauffement pour important à cause des SSTs imposées, donc peut être que tu verras un meilleur signal pour tes tests.

Maintenant on est à “PlioMIP2”, avec que des simus couplées. Pour cela, voir avec Ning, les résultats sont en cours d'analyse.
dans la simu PlioMIP2 de Ning, il y a des changements de la configuration océanique qui n'existaient pas dans PlioMIP1. Dans mes simus il n'y a pas de changements du land/sea mask, alors que dans celles de Ning il y a fermeture de Bering Stratit + Canadian Archipelago.

- fichiers de restart doivent contenir les champs chimiques de nos fichiers de restarts habituels (restartnitrate) et les champs physiques (de Svetlana par exemple) voir ajout_field.nco dans BRICOLAGE_RESTART*
- Climato d'ozone pour la strato dans INCA, correction sur le premier niveau du modele (mettre une pression bidon mais suffisamment élevée), cf mm repertoire que ci-dessus. PAS ABSOLUMENT CERTAINE QUE CE SOIT UN PROBLEME car a été diagnostiqué dans un eperiode de bugs multiples
- préciser pour lmdz qu'on est en bucket en commentant dans config.card les lignes débutant par SBG et SRF (pas besoin de recompiler)
- préciser dans COMP/lmdz.card un config.def adequat ( par ex config.def_preind pour 1850) pour les LL-GHG et paramètres orbitaux
- les SST et autres paramètres de surface sont définis dans le fichier limit.nc (dans le runs de Svetlana les SSt albedo, rugosité de surface, etc ont été précalculés par FOAM+LPJ puis remappés avec LMDz pr créer un fichier au limite)

basé sur discussion ave JB Ladant et testé dans un run présenté ici
Quand on est sur périodes sans glace sur pole sud, ca pose problème a ORCHIDEE. dans ce cas il faut modifier les fichier start.nc et limit.nc (cf BRICOLAGE_RESTART*) et modifier 3 routines de LMDz pour qu'il sache quece n'est pa un vrai point de glace (hydrol.F90 surf_lan_ice_mod.F90 et fonte_neige_mod.F90)

*BRICOLAGE_RESTART= $WORKDIR/LMDZORINCA_v6.2/modipsl/config/LMDZORINCA_v6/Eocene0/BRICOLAGE_RESTART/

isoprene dans fichier eclai 40.7
dans fichier eocene 71.5

Pour avoir les emissions sur la bonne grille, il faut indiquer le masque a INCAsflx. Pour cela, récuperer le fichier /ccc/store/cont003/dsm/p25botsy/IGCM_OUT/LMDZ/B55Ma4x/ATM/Output/Grid/B55Ma4x_grilles_gcm.nc qui contient la variable mask.
Puis “ncks -A -v lonu,lonv lmdz_grid_96x95.nc B55Ma4x_grilles_gcm.nc” si lonu et lonv ne debutent par par -180.

Les concentrations de d'ODS et GHG pour REPROBUS sont a changer dans les fichiers rangés à /$STORE/STRATO_INPUT_PALEO/RCP2.6CCMI_REFC2V2-PALEO/GHG_ODS/
RQ: Gauss 2006 CH3Cl and CH3Br were set to 550 and 6 pptv, respectively, in both the 2000 and 1850 simulations.

NMHC emission ratio Min Max Eocene 1.8 2.4 Miocene Pliocene 1.3 LGM 0.45 1.14 Preind 0.88 1.12 Present (2010) 1 1.54 RCP 0.74 1.27

CO emission ratio Eocene 1.8 2.5 Miocene Pliocene 2.1 LGM 0.36 0.65 Preind 0.72 1.276 Present (2010) 1.75 2.77 RCP 0.46 1.62

NOx emission ratio Min Max Eocene 1.5 1.7 Miocene Pliocene 1.3 LGM 0.41 0.62 Preind 0.77 1.45 Present (2010) 2.4 4.19 RCP 0.46 1.49