wiki:IPSLCMIP5/Centennal/Table3

Version 51 (modified by aclsce, 13 years ago) (diff)

--

Table 3 : Baseline simulations for model evaluation and for understanding historical and paleoclimates.

3,1 : piControl

piControl2 :

  • p86caub : Voir MONITORING piControl2
  • Moyenne saisonnière 100 ans disponible. Par exemple pour atmosphère : ATM/Analyse/SE_100Y
  • Atlas des différences (sauf coupe avec 39 niveaux) avec CM5MRPIRC1 disponible : CM5MRPIRC1_piControl2/SE_1950_2049_2000_2099/
  • Réalisée avec :
    • 1800 à 1850 : IPSLCM5A figé le 16 juillet 2010 (rev 1417 pour LMDZ)
    • Après 1850 : IPSLCM5A figé le 26 juillet 2010 (rev 1420 pour LMDZ) (Vérification résultats identiques)
      • post traitements sur cesium
      • Aga dans histmth
      • version COSP v1.3 dans LMDZ
      • correction prise en compte raz_date
  • Fréquence Outputs COSP :
    • monthly de 1800 à (2005?) 2020 puis COSP desactivé
    • monthly + daily : années 1 à 50 (1850 à 1900) et 121 à 170 (1970 à 2020)
  • En machine sur 8 CPUs depuis le 19 juillet 2010. Démarre en 1800-01-01 à partir de piControl1 à la date de 2049-12-31.
    • Repassé sur 4 CPUs lors de l'activation de COSP daily (1970)
  • Frequence Outputs :
    • lev_histmth=2 de 1800 à 2369 puis lev_histmth=5 à partir de 2370.
  • Comment reproduire exactement la simulation piControl2 ?
    • Récupérer un IPSLCM5A à jour (./model IPSLCM5A) et utiliser la configuration config/IPSLCM5A/piControl.
    • Utiliser le fichier climoz_LMDZ_1855.nc qui a servi à la simulation piControl2 en modifiant le COMP/lmdz.card (plus d'infos ) :
       ListNonDel= (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/AR5/HISTORIQUE/aerosols_11YearsClim_1855_v3.nc, aerosols.nat.nc),\
      -           (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/AR5/HISTORIQUE/climoz_LMDZ_1855_v2.nc, climoz_LMDZ.nc)
      +           (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/AR5/HISTORIQUE/climoz_LMDZ_1855.nc, climoz_LMDZ.nc)
      
    • Utiliser les mêmes états de démarrage que la simulation piControl2. Voir config.card

piControl1 :

  • p86maf
  • Réalisée avec la version IPSLCM5A du modèle figée le 17 juin 2010. Voir IPSLCM5A.
  • Fichiers de configuration enregistrés là : browser:CMIP5/piControl/piControl1
  • Expérience piControl de CMIP5. Démarre en 1800-01-01 à partir de CM5PIRC11 (3049-12-31). Voir : MONITORING piControl1
  • Fini en 2049.
  • La carte de végétation a changé lors du démarrage de CM5PIRC11. C'est l'année 1860 qui a été prise pour démarrer cet état initial.
  • Constante solaire = 1365.6537
  • Attention!
    • volume océan calculé sans partial steps
  • Souci dans les axes de temps de certains fichiers océan. Voir piControl1B pour récupérer les fichiers.

piControl1B :

  • p86maf
  • duplicata de piControl1
  • 1890-2029 pour remplacer les fichiers Output océan, glace et pisces.
  • Les solver.stat sont identiques
  • Les fichiers Output OCE, ICE et MBG ont été transférés dans la simu piControl1. Les TS OCE, ICE et MBG relancées sur cesium. Pas toutes OK. Les anciennes TS gardées sous Analyse_KO.

3,2 : historical

v3.historical4

v3.historical3

v3.historical2

v3.historical1

  • Vérifications carte végétation :
    • Evolution des part des 13 PFT dans v3.historical1
    • Evolution des parts des 13 PFT dans le fichier dtentrée
  • Comment reproduire exactement la simulation v3.historical1 ?
    • Récupérer un IPSLCM5A à jour (./model IPSLCM5A) et utiliser la configuration config/IPSLCM5A/historical.
    • Utiliser les mêmes états de démarrage que la simulation v3.historical1. Voir config.card

historical5

  • p86caub : Voir MONITORING historical5
  • Réalisée avec IPSLCM5A extrait et compilé le 27 août 2010
  • copie historical4 + démarrage décalé de 10 ans : départ de piControl2 à la date 1879-12-31
  • Sera arrêté en 19951231 (pour l'instant). En attente fichiers forcage 1995-2005.

historical4

  • p86caub : Voir MONITORING historical4
  • Réalisée avec IPSLCM5A extrait et compilé le 27 août 2010. Arrêté en 19951231 (en attente fichiers forçage 1995-2005, plus d'infos ) puis prolongé le 22/10/2007.
  • copie historical2 + démarrage décalé de 10 ans : départ de piControl2 à la date 1869-12-31
  • Sera arrêté en 19951231 (pour l'instant). En attente fichiers forcage 1995-2005.

historical3

  • p86manci : Voir MONITORING historical3
  • Réalisée avec IPSLCM5A extrait et compilé le 17 août 2010
  • copie historical2 + démarrage décalé de 10 ans :
  • partir du piControl2 à la date 1859-12-31
  • Sera arrêté en 19951231 (pour l'instant). En attente fichiers forcage 1995-2005.

historical2

  • p86caub : Voir MONITORING historical2
  • Réalisé avec IPSLCM5A figé le 26 juillet 2010 (rev 1420 pour LMDZ). Arrêté en 19951231 (en attente fichiers forcage 1995-2005, plus d'infos ) puis prolongé le 22/10/2007.
    • Aga
    • post-traitements sur cesium
    • version COSP v1.3 dans LMDZ
    • correction prise en compte raz_date
    • Mise à jour du fichier SOLARANDVOLCANOES_1610_2008.txt. Vous le trouverez disponible dans le répertoire qui suit (sur mercure): /dmnfs/cont003/p86ipsl/IGCM/BC/ATM/IPSLCM5A/IPCC_AR5/HISTORIQUE/. La version précédente de ce fichier est nommée SOLARANDVOLCANOES_1610_2008_V0.txt. Vous pouvez également le télécharger txt. Ou, Téléchargez ce fichier au format excel xls
    • Ce qui a changé : nous avons calculé la moyenne des valeurs « Sato volcanic radiative forcing (annual) » entre 1850 et 2003. La moyenne vaut -0,250 W/m2. On a simplement souhaité compenser ce forçage. D'où, une modification de la formule utilisée dans le fichier excel (en attaché également). Donc, suivant ce rainsonnement, les moyennes des valeurs SOLAR et résultat sur la période 1850 à 2003 devraient être les mêmes (après vérification, c'est sensiblement le cas).
    • Forcage solaire incluant les volcans : SOLARANDVOLCANOES_1610_2008.txt
      • Annual solar constant evolution, les figures correspondantes pdf
  • Fréquence Outputs COSP :
    • monthly tout le temps
    • monthly + daily : années 1979 à 2005
  • En machine sur 8 CPUs depuis le 24 juillet 2010. Démarre en 1850-01-01 à partir de piControl2 à la date de 1849-12-31.
  • note annexe de Martial : j'ai fait un test de relance de historical2 avec la carte de végètation de 1850 (et pas 1860 comme c'est récupéré par le piControl2 - voir les infos de cette simu). Voir le mutli-monitoring suivant : http://dods.ipsl.jussieu.fr/orchidee/CMIP5/VALID_histo2-histo3-histo5 différence relance carte de végétation 1850. Après analyse avec JLD, ces figures montrent que changer la carte entraîne un choc dans le couplé historique et que la correction des flux de land cover change apparaît minime.

historical1 :

  • p86maf
  • Mêmes binaires que piControl1.
  • Répertoire : historical
  • Fichier de forcage solaire utilisé : SOLAR_1610_2008.txt
  • Pas de volcans
  • En machine depuis le 21 juin. Voir: MONITORING historical1
  • Démarrage (config.card) RestartJobName=piControl1 RestartDate=1849-12-31 RestartPath=/dmnfs/cont003/IGCM_OUT/p86maf/IGCM_OUT/IPSLCM5A/PROD/piControl
  • Attention!
    • soumission par 24 mois et fréquence des rebuild 24 mois
    • TS 2D, monitoring, SE et atlas soumis sur la frontale
    • Séries temporelles (TS) soumises sur cesium. Voir détail dans log.
  • Souci dans les axes de temps de certains fichiers océan

Infos

Inutile depuis libIGCM v1_8 : tous les post-traitements sont faits sur cesium

  • 6/7/2010 : Explications pour tourner les TS complètes sur cesium :
    • A faire la première fois :
      • Préparer le terrain:
        cd $WORKDIR ; mkdir -p POST_CMIP5 ; 
        cd POST_CMIP5 ; svn_ano ; cd modipsl/util ; ./model libIGCM ; ./ins_job ; 
        
    • A faire pour chaque nouvelle simulation :
      cd $WORKDIR/POST_CMIP5 ; mkdir -p historical/MASIMU ; 
      cd POST_CMIP5/historical/MASIMU
      
      • Récupérer la description des post-traitements de la config sur mercure. xxxx est le répertoire où modipsl a été installé :
        scp -r mercure:xxx/modipsl/config/IPSLCM5A/historical/COMP .
        scp -r mercure:xxx/modipsl/config/IPSLCM5A/historical/POST .
        scp mercure:xxx/modipsl/config/IPSLCM5A/historical/config.card .
        
      • Programmer les post-traitements TS complets décrits dans les fichiers xxx.card.FULLTS :
        cd COMP
        cp lmdz.card lmdz.card_old ; cp lmdz.card.FULLTS lmdz.card
        cp opa9.card opa9.card_old ; cp opa9.card.FULLTS opa9.card
        cp orchidee.card orchidee.card_old ; cp orchidee.card.FULLTS orchidee.card
        cp stomate.card stomate.card_old ; cp stomate.card.FULLTS stomate.card
        cd ..
        
      • Installer les jobs à soumettre :
        cp ../../modipsl/libIGCM/create_ts.job .
        cp ../../modipsl/libIGCM/create_se.job .
        cp ../../modipsl/libIGCM/rebuild_fromWorkdir.job .
        cp ../../modipsl/libIGCM/TimeSeries_Checker.job  .
        
      • Vérifier/modifier dans !TimeSeries_Checker.job les variables :
        libIGCM=${libIGCM:=/work/cont003/p86maf/POST_CMIP5/modipsl/libIGCM}  # OK a priori
        POST_DIR=${POST_DIR:=$WORKDIR/OutScript/histobig4} # A choisir
        EXPERIMENT=${EXPERIMENT:=historical} # A modifier
        JobName=${JobName:=historical1} # A modifier
        CARD_DIR=${CURRENT_DIR} # A modifier
        
    • A faire à chaque fois que l'on veut relancer le calcul des TS complets :
      • Mettre à jour DateEnd dans config.card avec la date de fin de décennie à traiter
        vi config.card
        DateEnd=1859-12-31
        
      • Lancer en interactif !TimeSeries_Checker.job, en gardant la trace pour consultation ultérieure. xxx peut être la date traitée:
        ./TimeSeries_Checker.job |tee OUT_TSC_xxxx
        
      • L'ensemble des vérifications des TS est affiché. Pour celles qui manquent, un job est lancé sur cesium. Cela fait autour de 30 jobs soumis. Voir le nombre par exemple ainsi :
        grep Submitted OUT_TSC_xxx |wc -l
        
      • Les scripts de sortie sont disponibles dans le répertoire POST_DIR : $WORKDIR/OutScript/historical1
  • 5/7/2010 :
    • pour pouvoir tourner sans souci de quotas, il faut lancer les TS sur cesium (à l'aide de ./TimeSeriesChecker?), garder sur mercure, les rebuild, les TS monthly, SE et les atlas et monitoring.
    • ncwa : de temps en temps les rebuild s'arrêtent sur ces lignes là, il faut relancer le rebuild en question (et les TS et SE) à partir du fichier (1D_grid_T la plupart du temps) :
      IGCM_sys_ncwa : -Oh -a record -B mask3D==1 file1.nc file2.nc
      %NQSII(INFO): Batch job received signal SIGKILL. (Killed)
      
    • axe des temps des fichiers icemod avec valeurs nulles
    • les fichiers icemod peuvent avoir des champs nuls mais pas systématiquement lorsque le pas de temps est nul. (Vérification systématique en cours)

Attachments (5)

Download all attachments as: .zip