wiki:Climatology

Version 65 (modified by mafoipsl, 17 years ago) (diff)

--

Page to follow climatology of IPSL model at different resolution

Resolutions tests

Les tests de résolution se font sur des simulations de contrôle, aérosols constants sur valeurs 1980.

ORCA2xLMD9671

2L20

ORCA2xLMD14496

Sur mercure, toute est avec les nouveaux scripts

HH20

HH20B

  • mercure :
  • HH20B - comparable à HH20 - arrêt à 10 ans - voir HH20C plus loin :
    • mercure
    • 3 corrections : fonte calotte, calving, convection à 1.02 (cvl_corr dans physiq.def)
    • orchidee beton
    • nouveaux scripts :
      • voir mode d'emploi là : http://wiki.ipsl.jussieu.fr/wiki_ipsl/IGCMG/libIGCM/DocUtilisateur
      • EXP00/COMP/lmdz.card :
        • fichiers aerosols : garder 1980 et .nat seulement
        • suppression fichier flux_iceberg
      • config.card pris modèle sur config.card.OldName. Voir "attachment:config.card"
      • cartes PBS dans EXP00/Job_HH20B
         #PBS -l cputim_job=24:00:00   # limite en temps total
         #PBS -l restim_prc=20:00:00   # limite en temps par process
         #PBS -l memsz_job=10gb      # limite memoire
         ...
         PeriodNb=24
        
      • restart depuis HHSTART recopié de rhodes vers mercure (pas de restart Orchidee, restart depuis fichiers ancienne nomenclature) :
        /dmnfs/p86maf/SORTIES_CPL_IPSL/HHSTART/
        
        ./Atm/Restart:
        total 23812
        -rw-r--r--  1 p86maf dsm  8826164 2005-07-01 20:17 HHSTART_20901230_restart.nc
        -rw-r--r--  1 p86maf dsm 15552736 2005-07-01 20:17 HHSTART_20901230_restartphy.nc
        
        ./Cpl/Restart:
        total 3560
        lrwxrwxrwx  1 p86maf dsm      28 2007-06-01 18:46 HHSTART_20901230_flxat.nc -> HHSTART_flxatmos_20901230.nc
        lrwxrwxrwx  1 p86maf dsm      28 2007-06-01 18:46 HHSTART_20901230_sstoc.nc -> HHSTART_sstocean_20901230.nc
        -rw-r--r--  1 p86maf dsm 2124400 2005-06-30 11:49 HHSTART_flxatmos_20901230.nc
        -rw-r--r--  1 p86maf dsm 1519344 2005-06-30 11:45 HHSTART_sstocean_20901230.nc
        
        ./Lim/Restart:
        total 12504
        -r--r--r--  1 p86maf dsm 12800355 2004-11-22 15:55 HHSTART_20901230_restart_icemod
        
        ./Oce/Restart:
        total 88988
        -r--r--r--  1 p86maf dsm       38 2004-11-22 15:54 HHSTART_20901230_EMPave.dat
        -r--r--r--  1 p86maf dsm 91119608 2004-11-22 15:54 HHSTART_20901230_restart.nc
        
        
    • modifs sources nécessaires (comme HH20) :
      • readsulfate : LMDZ4/libf/phylmd/readsulfate.F ligne 106-107 pour faire un contrôle actuel
         IF (iyr .lt. 4850) THEN
             cyear='1980'
        
      • filtrez : LMDZ4/libf/filtrez/parafilt.h (à automatiser?)
         c 144 96 19 non-zoom:
               PARAMETER (nfilun=16, nfilus=16, nfilvn=16, nfilvs=16)
        
      • iceberg : Commenter dans LMDZ4/libf/phylmd/interface_surf.F90 la lecture du fichier iceberg (obsolète)
         !!    open (12,file='flux_iceberg',form='formatted',status='old')
         !!    read (12,*) coeff_iceberg
         !!    close (12)
        
    • lancement le 1er juin 2007

HH20BB

  • mercure
  • Tests sur nouveaux poids.
    • Dans le Job : R_BC_TEST=/dmnfs/p25mart/IGCM/BC
    • Dans COMP/oasis.card : {R_BC_TEST}
  • Plantage au bout de 8 ans sur un point d'océan à salinité négative : i=117, j=75 (Amazone).

HH20C

  • mercure : /home/p86maf/HH20C/modipsl/config/IPSLCM4_v1_OASIS3/EXP20C
  • En fait, le paramêtre cvl_corr n'était pas utilisé dans HH20B.

Le 11 juin 2007 : ajout dans physiq.F de :

         CALL conema3 (dtime,
 ...
          ENDIF ! ok_cvl
c Correction precip
          rain_con = rain_con * cvl_corr
          snow_con = snow_con * cvl_corr 
  • c'est la seule différence entre HH20B et HH20C
  • relance pour 10 ans
  • relance pour 100 ans, compte tenu des soucis de la simu avec les nouveaux poids HH20B3.

HH20B3

  • mercure
  • Reprise de la simulation avec nouveaux poids (HH20BB).
  • Idem HH20C
  • Plantage au bout de 2 ans.
  • Dans OPA/SRC_ORCA/flxrnf.F. Reprise modifs Olivier. Shéma upstream autour des embouchures.
diff /scratchdir/p86caub/CONTNEW/modipsl/modeles/OPA/SRC_ORCA/flxrnf.F flxrnf.F
flxrnf.F flxrnf.F                                                             <
143,148d142
<           IF (nrunoff == -2) THEN
<              IF(lwp)WRITE(numout,*) ' '
<              IF(lwp)WRITE(numout,*) 'Simulation with Upstream '
<              IF(lwp)WRITE(numout,*) '    at river mouth       '
<              IF(lwp)WRITE(numout,*) ' '
<           ENDIF
284,311c278
<       IF ((nrunoff == 2 .OR. nrunoff == -2 ) .AND. kt == nit000 ) THEN
<           IF ( nrunoff == -2 ) then
<               clname='runoff_1m_nomask' 
<               CALL flinopen(clname,nizoom,jpi,njzoom,jpj,.false.,ipi,ipj
<      $            ,ipk,zlon,zlat,zlev,jpmois,istep,zdate0,rdt,numrof) 
< C   ... Title, dimensions and tests
< #  if ! defined key_coupled
<               IF(iman /= jpmois) THEN
<                   IF(lwp)WRITE(numout,*) ' '
<                   IF(lwp)WRITE(numout,*) 'problem with time coordinates'
<                   IF(lwp)WRITE(numout,*) ' iman ',iman,' jpmois ',jpmois
<                   STOP 'flxrnf.F'
<               ENDIF
<               IF(lwp)WRITE(numout,*) iman,istep,zdate0,rdt,numrof
<               IF(lwp)WRITE(numout,*) 'numrof=', numrof
<               IF(lwp)WRITE(numout,*) 'jpmois=', jpmois
<               IF(lwp)WRITE(numout,*) 'rdt=', rdt
< #  endif
<               IF(ipi /= jpidta .and. ipj /= jpjdta .and. ipk /= 1) THEN
<                   IF(lwp)WRITE(numout,*) ' '
<                   IF(lwp)WRITE(numout,*) 'problem with dimensions'
<                   IF(lwp)WRITE(numout,*) ' ipi ',ipi,' jpidta ',jpidta
<                   IF(lwp)WRITE(numout,*) ' ipj ',ipj,' jpjdta ',jpjdta
<                   IF(lwp)WRITE(numout,*) ' ipk ',ipk,' =? 1'
<                   STOP 'flxrnf.F'
<               ENDIF
<               IF(lwp)WRITE(numout,*) 'ipi=', ipi,' ipj=',ipj,' ipk=',ipk
<           ENDIF
---
>       IF (nrunoff == 2 .AND. kt == nit000) THEN

HOL20 : ORCA2xLMD14496

  • mercure : /home/p86maf/HOL20/modipsl/config/IPSLCM4_v1_OASIS3/EXP00
  • Comparable à HH20 : en -6000 - arrêt à 10 ans - voir HOL20C plus loin :
    • lancé le 4 juin 2007
    • sur mercure : /home/p86maf/HOL20/modipsl/config/IPSLCM4_v1_OASIS3/EXP00
    • Démarrage à partir de état initial et Levitus
    • fichiers résultats : /dmnfs/p86maf/IGCM_OUT/IPSLCM4_v1_OASIS3/HOL20
    • modifications dans les sources LMDZ :
      • readsulfate.F
         IF (iyr .lt. 4850) THEN
                   cyear='.nat'
        
    • modifications dans les fichiers d'entrée :
      • physiq.def
         # parametres climatiques
         # pour 6ka
         R_ecc  = 0.018682
         R_peri = 0.87
         R_incl = 24.105
         solaire = 1365.
         co2_ppm = 280.
         CH4_ppb = 650.
         N2O_ppb = 270.
         CFC11_ppt = 0.
         CFC12_ppt = 0.
        

HOL20C

  • mercure : /home/p86maf/HOL20/modipsl/config/IPSLCM4_v1_OASIS3/EXP20C
  • En fait, le paramêtre cvl_corr n'était pas utilisé dans HOL20.

Le 11 juin 2007 : ajout dans physiq.F de :

         CALL conema3 (dtime,
 ...
          ENDIF ! ok_cvl
c Correction precip
          rain_con = rain_con * cvl_corr
          snow_con = snow_con * cvl_corr 
  • c'est la seule différence entre HOL20 et HOL20C
  • relance pour 10 ans
  • relance pour 100 ans. COmpte tenu des pbs de poids dans HH20B3.

ORCA2xLMD144142

VV20 : ORCA2xLMD144142

  • brodie : /workdir/rech/psl/rpsl003/VV20/modipsl/config/IPSLCM4_v1_OASIS3/EXP00
  • Etapes :
    • Démarrage :
      • Choix de la résolution : LMD144142
      • LMDZOR simulation (tags orchidee "beton" et LMDZ4 HEAD ), Ionela récupère o2a.nc à la bonne résolution et fait tourner en forcé les années AMIP: HAFO1
    • Olivier prépare les états initiaux du coupleur (poids et fichiers initiaux faits)
    • Sébastien fait les fichiers aérosols à la bonne résolution (fait et accessible sur CCRT et IDRIS)
    • Modifications dans les sources de LMDZ4 :
      • parafilt.h contient les valeurs : 24 23 24 24
    • Marie-Alice lance le couplé VV20
  • Etat au 18 janvier 2007 :
    • 10 ans de VV20 ont tourné. Plantage en février de la 11ème année (hgardfou). Il y a eu un plantage en février de la 4ème année sur hgardfou. On a relancé ce seul mois de février avec les paramêtres tetagdiv, tetagrot et tetatemp divisés par 2 (3600, 5400 et 5400).
    • Inquiétude sur temps de calcul : un mois de couplé VV20 (IPCC) prend environ 3500 secondes de temps CPU pour la partie LMDZ-ORCHIDEE. Pour le forcé (HEAD), c'est le double de temps. Souci de vectorisation inhibée réglé depuis.
  • Etat au 25 janvier 2007 :
    • 20 ans de VV20 ont tourné. Les paramêtres tetagdiv, tetagrot et tetatemp ont été divisés par 2 (et sont à 3600, 5400 et 5400) depuis le mois de février 1870.
  • Etat au 5 février 2007 :
    • relance car arrêt le 29 janvier, inexpliqué. Relance au 1892 12 1.
    • chaque executable oasis était gardé. rm de tous sauf les 2 premiers. Date bizarre à décembre 2007. modifiée à la main pour être comme les autres executable : cd /u/rech/psl/rpsl003/SORTIES_CPL_IPSL/VV20/Exe ; touch 01081700 VV20_18600101_18600130_oasis ; touch 01081710 VV20_18600201_18600230_oasis. OK
  • Etat au 19 février 2007 :
    • relance car VV20 s'était arrêté en time limit le 14 février à la date de février 1899
  • 21 mars 2007 : 50 ans simulés avec atlas y compris océans additionnels et synchro sur le site mc2
  • 22 mars 2007 : Relance pour atteindre 100 ans au moins (demande Eric Guilyardi et Frédéric Hourdin)
  • 4 avril 2007 : Plantage CPU limit le 1er avril, relance le 4 avril 2007
  • 18 avril 2007 : Poursuite au delà des 100 ans (voir mail Frédéric Hourdin 17 avril 2007)
  • 23 avril 2007 : relance après CPU limit exceeded le 22 avril 2007(1958 7 1)
  • 11 mai 2007 : arrêt. Bilan d'eau à vérifier. 1996 12 1. Lancement VV21, avec correctif calotte oté. Voir plus bas.
  • Resolution ORCA2xLMD144142 ajouté pour test dans les configurations IPSLCM4_v1_OASIS3 (HEAD et CHILI)
  • Configuration IPSLCM4_v1_6 (22 janvier 2007 modifié 15 février 2007 suite réunion cplipsl) :
    • But : mise en place de couplés en production avec physique revue, pour test seulement
    • avec Orchidee 1_6 : attention aux paramêtres, fichier LAI corrigé, ...
    • LMDZ4 - V3 validé pour Orchidee 1_6 (Voir Laurent Fairhead)
    • pas de changement pour le reste par rapport à IPSLCM4_v1_OASIS3 :
      • OASIS3
      • OPA - LIM
      • Résolutions LMDZ : 72x45, 96x71, 144x96 et 144x142
      • Scripts : anciens et nouveaux

VV20BB

  • mercure : /home/p86maf/VV20BB/modipsl/config/IPSLCM4_v1_OASIS3/EXPBB
  • tests sur nouveaux poids
    • Dans le Job : R_BC_TEST=/dmnfs/p25mart/IGCM/BC
    • Dans COMP/oasis.card : {R_BC_TEST}
  • plantage au bout de 3 ans au même point océan que HH20BB : i=117, j=75 (Amazone).

VV20B3

  • mercure : /home/p86maf/VV20BB/modipsl/config/IPSLCM4_v1_OASIS3/EXP20B3
  • Reprise de la simulation avec nouveaux poids (VV20BB).
  • Idem HH20C :
    • Dans OPA/SRC_ORCA/flxrnf.F. Reprise modifs Olivier. Schéma upstream autour des embouchures.
diff /scratchdir/p86caub/CONTNEW/modipsl/modeles/OPA/SRC_ORCA/flxrnf.F flxrnf.F
... Voir HH20B3 plus haut

VV20P : ORCA2xLMD144142 - LMDZ et ORCHIDEE parallèles

  • brodie : /workdir/rech/psl/rpsl003/VV20P/modipsl/config/IPSLCM4_v2_PAR/EXP00
  • anciens scripts
  • 3 mois de VV20P ont tourné. Plantage sur un point dans l'océan avec salinité négative. C'est l'embouchure du Zambèze. Même plantage sur 1 proc que sur 4 procs. Même plantage mais 6 jours après sur mercure. Recherche en cours sur la physique. Aide bienvenue. Couplé 96x71 plante aussi.
  • En attente (juin 2007)

VV21 : ORCA2xLMD144142

  • brodie : /workdir/rech/psl/rpsl003/VV20/modipsl/config/IPSLCM4_v1_OASIS3/EXP21
  • Idem VV20 - mêmes binaires
  • ok_lic_melt à .true.
  • 10 ans pour voir

Comparaisons 2L20, HF2, HH20, HH20B, HAFO1, VV20, VV20P

Tableau récapitulatif :

Infos 2L20 HF2 HH20 HH20B HAFO1 VV20 VV20P
Comment mercure - nouveaux scripts - 3 correctifs appliqués test couplé parallèle
Résolution Atm 96x71 144x96 144x96 144x96 144x142 144x142 144x142
Résolution Oce ORCA2 = = = = = =
Tag Atm LMDZ 3.3 IPSL-CM4_IPCC_branch = = Head IPSL-CM4_IPCC_branch V3_test
Tag Orchidee orchidee_1_3 orchidee_1_3_1 = beton beton = ORCHIDEE_PARA
SECHIBA_QSINT 0.02 = = 0. 0. 0. -
RSOL_CSTE - - - 33000. 33000. = -
RVEG_CSTE - - - 0.5 0.5 = -
cvl_corr - - - 1.02 - - -
traceur.def none yes = = = = =
co2_ppm 348. = = = = = =
Aérosols constant 1980 = = = sans constant 1980 constant 1980
day_step 480 960 - plantage au bout de 3 ans pour 480 = 720 720 720 720
tetagdiv 10800. 3600. = = (7200.) 3600. =
tetagrot 18000. 3600. = = (10800.) 5400. =
tetatemp 18000. 3600. = = (10800.) 5400. =
Etat initial Atm create_etat0_limit = = = = 30/12/1980 HAFO1 en inhibant la lecture de <DELTAT> (VV20_etat0 a tourné 1 an à partir état initial) 30/12/1980 HAFO1 - plante sur pression sol négative dans integrd
Etat initial Oce Levitus - décennie 2090 de 2L20 - idem HH20 idem HH20
Fichier LAI V0 = = V1 V1 V1 V1
Sorties LMDZ #define histhf, histday, histmth, histins, histISCCP #define histhf, histday, histmth, histins #define histins, histhf, histday, histmth, histmthNMC, histISCCP #define histhf, histday, histmth, histins #define histhf, histday, histmth, histins
atlas 2L20 HF2 HH20 VV20
moyennes annuelles OK - OK - OK -
status ended ended ended running running running depuis etat0 plantage au 3ème mois, depuis HAFO1 plantage integrd

ORCHIDEE tests

Tests couplé orchidee_1_3_beton / orchidee1.5

  • Tableau récapitulatif des runs. FH01 est un run actuel issu de orchidee_1_3 + béton. Les autres sont des runs pré-industriels issus de variantes de orchidee_1_5.
  • FH01, FH05 ET FH06 ont le physiq.def et le orchidee.def tunés par Frédéric.
  • FH02 et FH04 ont le physiq.def et le orchidee.def du CVS couplé.
Infos FH01 FH02 FH04 FH05 FH06 FH07 FH08 FH09 FH10 FH11 FH12
Comment correction d'eau OM+JLD de septembre 2006
Résolution Atm 96x71x19 = = = = = = = = = =
Résolution Oce ORCA2 = = = = = = = = = =
Tag Atm LMDZ 3.3 IPSL-CM4_IPCC LMDZ4 IPSL-CM4_IPCC_branch = = = = = = = = =
Tag Orchidee orchidee_1_3 + beton orchidee_1_5 = = = = = = = orchidee_1_8 orchidee_1_8
Albedo Orchidee orchidee_1_3 = = = = = = = Albedo de orchidee_1_6 = =
LMDZ-T2M-Q2M natif = Nathalie natif = = Nathalie Nathalie Nathalie Nathalie Nathalie
CDRAG_from_GCM TRUE = = = FALSE TRUE = = = = =
SECHIBA_QSINT 0. 0.02 0.02 0. 0. 0. 0.02 0.02 0.02 0.02 0.02
RSOL_CSTE 33000. 33000. dans constantes_veg.f90 = = = = = = = = =
RVEG_CSTE 0.5 1 1 0.5 = = = = = = =
epmax .995 .99 .99 .995 = = = = = = =
cld_tau_lsc 1800. 3600. 3600. 1800. = = = = = = =
cld_tau_con 1800. 3600. 3600. 1800. = = = = = = =
fact_cldcon 2. 1. 1. 2. = = = = = = =
ratqshaut .25 .33 .33 .25 = = = = = = =
day_step 480 = = = = = = = =
tetagdiv 10800. = = = = = = = = = =
tetagrot 18000. = = = = = = = = = =
tetatemp 18000. = = = = = = = = = =
traceur.def none yes = = = = = = = = =
co2_ppm 348. 280. 280. 280. 280. 348. = = = = =
Aérosols constant 1980 pré-industriel pré-industriel pré-industriel pré-industriel constant 1980 = = = = =
Etat initial Atm 30 décembre 2209 2L20 30 décembre 2659 2L24 = = = 30 décembre 2209 2L20 = create_etat0 = = =
Etat initial Oce 30 décembre 2209 2L20 30 décembre 2659 2L24 = = = 30 décembre 2209 2L20 = Levitus = = =
Fichier LAI V0 V1 = = = = = = = = =
slowproc.f90 CVS 1_3 CVS 1_5 = = = = = = = = =
Sorties LMDZ #define histday, histmth, histNMC = = = = = = = = = =
atlas FH01 FH02 FH04 FH05 FH06 FH07 FH08 FH09 FH10 FH11 FH12

Attachments (2)

Download all attachments as: .zip