Version 65 (modified by mafoipsl, 17 years ago) (diff) |
---|
Page to follow climatology of IPSL model at different resolution
- Voir Test Orchidée là : Climatology
Resolutions tests
Les tests de résolution se font sur des simulations de contrôle, aérosols constants sur valeurs 1980.
ORCA2xLMD9671
2L20
- Atlas are here : http://mc2.ipsl.jussieu.fr/PHP/experience.php?exp=2L20
ORCA2xLMD14496
Sur mercure, toute est avec les nouveaux scripts
HH20
- Atlas are here : http://mc2.ipsl.jussieu.fr/PHP/testing.php?exp=HH20
- See wiki page (ipsl/vivaipsl) : http://wiki.ipsl.jussieu.fr/wiki_ipsl/CoupléHauteResolution or attached file : attachment:"HH20.pdf"
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)
- readsulfate : LMDZ4/libf/phylmd/readsulfate.F ligne 106-107 pour faire un contrôle actuel
- 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'
- readsulfate.F
- 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.
- physiq.def
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
- Démarrage :
- 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 |
- Atlas de comparaison HH20/VV20 : voir par décennie
- Atlas de comparaison 2L20/VV20 : voir par décennie
- Atlas de comparaison 2L20/HH20 : 7 décennies
ORCHIDEE tests
Tests couplé orchidee_1_3_beton / orchidee1.5
- Atlas de comparaison 2L20/FH01 : http://dods.idris.fr/rces452/ATLAS_COMPARAISONS/2L20_FH01
- Atlas de comparaison 2L24/FH02 : http://dods.idris.fr/rces452/ATLAS_COMPARAISONS/2L24_FH02
- Atlas de comparaison 2L24/FH04 : http://dods.idris.fr/rces452/ATLAS_COMPARAISONS/2L24_FH04
- 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)
- HH20.pdf (127.2 KB) - added by mafoipsl 17 years ago.
-
config.card
(4.1 KB) -
added by mafoipsl 17 years ago.
config.card used for HH20B
Download all attachments as: .zip