[[PageOutline]] = Tests IPSLCM5 avec 39 niveaux verticaux = == Conteste des tests de réglage en 39 niveaux verticaux == - On part d'un '''LMDZOR''' 96x95x19 à peu près équilibré ('''X8L1900''', '''X8''' la machine et '''L19''' le nombre de niveaux). - On en décline un LMDZOR '''96x95x39''' avec les mêmes '''sources''' et les mêmes '''.def''' ('''X8L3900''', '''X8''' la machine et '''L39''' le nombre de niveaux), ('''P6L3900''', '''P6''' la machine Power6 et '''L39''' le nombre de niveaux) - Les mêmes def à ceci près que X8L39* a '''ok_strato=y''' et une '''nouvelle discrétisation verticale'''. - On n'utilise pas les aérosols. ('''ok_ade=n ; ok_aie=n''') - Début de simulation en 1980-01-01, start sechiba issue de X8L1900 pour les X8L39* - Les fichiers .def et les start* dans ce [http://dods.extra.cea.fr/data/p86denv/LMDZOR/START_DEF_39/ Kit de démarrage] L3900. '''Bientôt sur svn'''. - '''Compilé''' pour le couplé, '''tourné''' en forcé. - '''On fait varier un certain nombre de paramètres dans le but de savoir comment permettre un équilibrage (X8L3900 --> X8L3913)''' - - Révision LMDZ : 1165 de lmdz-dev + ORCHIDEE_1_9_2 - '''On fait varier un certain nombre de paramètres dans le but de comprendre comment éviter le plantage du modèle dans les hautes couches atmosphériques (X8L3914 --> X8L39??)''' - Révision LMDZ : 1179 de lmdz-dev + ORCHIDEE_1_9_2 : dépendance verticale de l'intensité de la diffusion, activation de la paramétrisation de Hines dans la "sponge layer" stratosphérique. - '''On utilise une première combinaison de paramètres pour équilibrer le modèle et on refait varier un certain nombre de paramètres pour affiner l'équilibrage. X8L3925 --> X8L39??), entre autre la diffusion''' - Révision LMDZ : 1185 de lmdz-dev + ORCHIDEE_1_9_2 : dépendance verticale de l'intensité de la diffusion, activation de la paramétrisation de Hines dans la "sponge layer" stratosphérique. || '''PREMIER ROUND D'ÉQUILIBRAGE et Accès dods''' || Résolutions || day_step || '''Strato/Hines/NewDiffusion/TopBound''' || tetagdiv || tetagrot || tetatemp || pmagic || '''Autres Diffs // référence L1900''' || Flux net au sommet || Durée || || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L1900/ X8L1900] || 96x95x'''19''' || 480 ||'''N'''/N/N/N|| 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || || -0.7W/m² || 721 mois && STOP || || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3900/ X8L3900] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || || -7.3W/m² ||82 mois && HGARDFOU|| || [http://dodsp.idris.fr/rces452/LMDZOR/P6L3900/ P6L3900] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || || -7.6W/m² || '''348 mois''' && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3901/ X8L3901] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || cld_lc_[lsc][con]=1.3 10-4 || -4.1W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3902/ X8L3902] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || cld_tau_[lsc][con]=1800. || -2.0W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3903/ X8L3903] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || ffalv_[lsc][con]=0.5 || -3.9W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3904/ X8L3904] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || coef_eva=1e-5 || -5.9W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3905/ X8L3905] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || ratqshaut=0.2 || -5.9W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3906/ X8L3906] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || rad_froid=50 || -7.3W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3907/ X8L3907] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || rad_chau[\1][\2]=14 || -5.9W/m² || 57 mois && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3908/ X8L3908] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || epmax=0.995 || -6.0W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3909/ X8L3909] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || nbapp_rad=48 || -5.8W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3910/ X8L3910] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || '''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.5''' || +1.7W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3911/ X8L3911] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || '''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.75''' || -0.9W/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3912/ X8L3912] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || '''0.00''' || '''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.5''' || +3.7/m² || 73 mois && STOP|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3913/ X8L3913] || 96x95x'''39''' || 480 || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || '''0.00''' || '''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.75''' || +1.1W/m² || 73 mois && STOP|| || '''TEST DE STABILITÉ,EXPÉRIENCES,DÉVELOPPEMENTS''' || || || || || || || || || || || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3914/ X8L3914] || 96x95x'''39''' || '''720''' || Y/N/N/N || 5400(it=2)|| 7200(it=2) || 7200(it=2) || 0.02 || || -7.0W/m² || 139 mois && HGARDFOU|| || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3915/ X8L3915] || 96x95x'''39''' || 480 || Y/N/N/N || '''3600'''(it=2)|| '''5400'''(it=2)|| '''5400'''(it=2)|| 0.02 || || -7.7W/m² || '''320 mois''' && HGARDFOU|| || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3916/ X8L3916] || 96x95x'''39''' || 480 || Y/N/N/N || 5400'''(it=1)'''|| 7200(it=2) || 7200(it=2) || 0.02 || || -8.4W/m² || 157 mois && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3917/ X8L3917] || 96x95x'''39''' || 480 || Y/N/N/N || '''3600(it=1)'''|| 7200(it=2) || 7200(it=2) || 0.02 || || -8.4W/m² || 128 mois && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3918/ X8L3918] || 96x95x'''39''' || 480 || Y/N/N/N || 5400'''(it=1)'''|| '''5400'''(it=2) || '''5400'''(it=2) || 0.02 || || -8.4W/m² || '''477 mois''' && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3919/ X8L3919] || 96x95x'''39''' || 480 || Y/N/N/N || 5400'''(it=1)'''|| '''5400'''(it=2) || 7200(it=2) || 0.02 || || -8W/m² || 19 mois && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X8L3920/ X8L3920] || 96x95x'''39''' || 480 || Y/N/N/N || 5400'''(it=1)'''|| 7200(it=2) || '''5400'''(it=2) || 0.02 || || -8.3W/m² || 283 mois && HGARDFOU|| || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3921/ X8L3921] || 96x95x'''39''' || 480 ||'''Y/Y/Y'''/N|| '''3600'''(it=2)|| '''5400'''(it=2)|| '''5400'''(it=2)|| 0.02 || || -7.5W/m² || '''1649 mois & Running''' || || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3922/ X8L3922] || 96x95x'''39''' || 480 ||'''Y/N/Y'''/N|| '''3600'''(it=2)|| '''5400'''(it=2)|| '''5400'''(it=2)|| 0.02 || || -7.5W/m² || '''799 mois & HGARDFOU''' || || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3923/ X8L3923] || 96x95x'''39''' || 480 ||'''Y/Y/Y'''/N|| '''5400(it=1)'''|| '''5400'''(it=2)|| '''5400'''(it=2)|| 0.02 || || -8.1W/m² || '''1478 mois & Running''' || || [http://dodsp.idris.fr/rces452/LMDZOR/X8L3924/ X8L3924] || 96x95x'''39''' || 480 ||'''Y/N/Y'''/N|| '''5400(it=1)'''|| '''5400'''(it=2)|| '''5400'''(it=2)|| 0.02 || || -8.1W/m² || '''493 mois''' & HGARDFOU || || '''DEUXIÈME ROUND D'ÉQUILIBRAGE''' || || || || || || || || || || || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X9L3925/ X9L3925] || 96x95x'''39''' || 480 ||'''Y/Y/Y'''/N|| '''5400(it=1)'''|| '''5400'''(it=2)|| '''5400'''(it=2)|| '''0.0'''||'''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83''' || 0.05W/m² || '''Running''' || || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X9L3926/ X9L3926] || 96x95x'''39''' || 480 ||'''Y/N/Y'''/N|| '''5400(it=1)'''|| '''5400'''(it=2)|| '''5400'''(it=2)|| '''0.0'''||'''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83''' || 0.17W/m² || '''127 mois''' && HGARDFOU|| || [http://dods.extra.cea.fr/data/p86denv/LMDZOR/X9L3927/ X9L3927] || 96x95x'''39''' || 480 ||'''Y/Y/Y/Y'''|| '''5400(it=1)'''|| '''5400'''(it=2)|| '''5400'''(it=2)|| '''0.0'''||'''cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83''' || 0.06W/m² || '''Running''' || == Commentaires / Mémos / '''pour la partie équilibrage''' == - Des simulations d'équilibrage, '''X8L3900''' est la plus froide avec '''-7,3 W/m²'''. - ''' X8L3900 plante peu après les 6 années de simulation.''' --> X8L3914 et X8L3915 ... pour se rassurer et comprendre. - ''' P6L3900 plante peu après 29 années de simulation. Elle a tourné sur Power6'''. - À part P6L3900, toutes les simulations ont tourné sur '''SX8 IDRIS ou SX8-R CCRT'''. Les simulations LMDZOR sont '''identiques entre SX8 et SX8-R'''. - Les mieux du point de vu de l'équilibrage étant la '''X8L3911''' '''(-1 W/m²)''', et la '''X8L3913''' (+1 W/m²)'''. - '''En réduisant le temps de vie des nuages bas''' (cld_tau_lsc/con=1800, X8L3902) on regagne presque tout et les CRFsw ne sont pas trop bête. En plus, ce sont les constantes qu'on avait dans le modèle dans des réglages plus anciens. cf : [http://igcmg.lsce.ipsl.fr/monitoring/tmp/fegg_plot02_lmXHvW_prod/ X8L1900 vs X8L3900 vs X8L3901 vs X8L3902 vs X8L3903 vs X8L3904] - On a pu gagner les 2W/m² manquant avec plusieurs autres réglages. On a choisi les '''ffalv_[lsc][con] à 0.5 et 0.75 en combinaison des cld_tau_[lsc][con]=1800'''. --> X8L3910 et X8L3911 cf : [http://igcmg.lsce.ipsl.fr/monitoring/tmp/fegg_plot02_0JrsnE_prod/ X8L1900 vs X8L3900 vs X8L3902 vs X8L3903 vs X8L3910 vs X8L3911] - On cherche par ailleurs à '''estimer''' la valeur de pmagic permettant de corriger le "coté transitoire" du contrôle ---> X8L3912 et X8L3913 - Le passage 0.02 à 0.0 en pmagic fait gagner 2W/m² lors des deux tests X8L3912 et X8L3913 - On gagne donc '''1W/m²''' par pas de '''0.01''' en pmagic en forcé. cf : [http://igcmg.lsce.ipsl.fr/monitoring/tmp/fegg_plot02_Ylwl6y_prod/ X8L1900 vs X8L3900 vs X8L3910 vs X8L3911 vs X8L3912 vs X8L3913] - Le passage à nitergdiv=1 (la préférence de Fred) fait perdre 0.6 W/m² - On peut imaginer que la combinaison '''nitergdiv=2+cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.66+ *pmgagic=0.02*''' nous donne un flux net au sommet de '''0W/m²''' - On peut imaginer que la combinaison '''nitergdiv=2+cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.75+ *pmgagic=0.01*''' nous donne un flux net au sommet de '''0W/m²''' - On peut imaginer que la combinaison '''nitergdiv=2+cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=1.0+ *pmgagic=0.0*''' nous donne un flux net au sommet de '''0W/m²''' - On peut imaginer que la combinaison '''nitergdiv=1+cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83+ *pmgagic=0.0*''' nous donne un flux net au sommet de '''0W/m²''' == Ce qu'il reste à faire pour ces réglages en forcé avant passage couplé == - Estimation de l'impact sur l'équilibrage LMDZOR de '''l'ajout des aérosols''' période actuelle (sulfatés principalement). En attente de Yves. - '''Reréglage en conséquence'''. Jean-Louis valide l'ajout des aérosols en regardant leurs forçages radiatifs. - Estimation du déséquilibrage que devra compenser pmagic pour refléter l'aspect transitoire d'un '''réglage sur la période actuel'''. - Réadapter les réglages (à pmagic constant, avec tous les aérosols) une fois '''inclus les nouveaux coefficients cdrag''' (Pascale, Frédéric) == Commentaires / Mémos / '''pour la partie hgardfou dans les hautes couches''' == - L'arrêt hgardfou se produit '''systématiquement''' aux hautes latitudes, pendant l'hiver (boréal ou austral), dans la '''couche 35''' du modèle. - L'arrêt hgardfou se produit '''systématiquement''' juste '''sous la "sponge layer"''' qui est active dans les '''4 dernières couches''' du modèle par défaut. - Ajout d'une '''dépendance verticale''' de l'intensité de la '''diffusion'''. (Plus de diffusion en haut). Depuis la révision '''1176''' lmdz. - On a activé la '''paramétrisation de Hines''' dans la "sponge layer" des hautes couches stratosphériques (+7 à 8 % en calcul) ; '''ok_hines=y''' == Voir aussi la table de Fred des simulations de juillet 2009 == [http://dodsp.idris.fr/rlmd282/LMDZAR5/R96x95x39A2/HTML/tab_simulations.html Table des simulations de juillet 2009] == Tests aerosols == * Sources : * LMDZ 1196 avec modifs lecture aerosols (minval/maxval) * Orchidee 1 9 4 * 96x95x39 * C'est là : /scratch/cont003/p86maf/CMIP5_LMDZOR/PARA_SX9_CM5_R1196_39L/modipsl/config/LMDZOR === flag_aerosol=1 === * Lecture SO4, fichiers interpollés par Fred et mis au CCRT là : * /dmnfs/cont003/p86maf/IGCM/BC/ATM/LMDZOR/LMD9695/AR5/HISTORIQUE/SO4.run1980.nc -> SO4.2000.2000.96x95.nc * SO4.run.nat.nc -> SO4.1750.1750.96x95.nc * 5 tests (1 an, puis 10 ans), lev_histday=3, lev_histmth=4 pour avoir les diagnostics aerosols * CLD10 idem CLD10 Sébastien * LMDACT : ok_ade=y, period=actuel * LMDIACT : idem + ok_aie=y * LMDPRE : ok_ade=y, period=preind * LMDIPRE : idem + ok_aie=y * Temps calcul : * CLD10 : 700-900 s avec diagnostics (600 s sans) (elapsed time) * LMDACT et LMDPRE : 850-1000 s avec diagnostics (800 s sans) (elapsed time) * LMDIACT et LMDIPRE : 1000-1200 s avec diagnostics (830 s sans) (elapsed time) * Accès aux fichiers || ok_ade=n ok_aie=n || actuel |||| preind |||| || || ok_ade=y ok_aie=n || ok_ade=y ok_aie=y || ok_ade=y ok_aie=n || ok_ade=y ok_aie=y || || CLD10 || LMDACT || LMDIACT || LMDPRE || LMDIPRE || || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/CLD10/MONITORING monitoring CLD10] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDACT/MONITORING monitoring LMDACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDIACT/MONITORING monitoring LMDIACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDPRE/MONITORING monitoring LMDPRE] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDIPRE/MONITORING monitoring LMDIPRE] || || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/CLD10/ATM/Output sorties MO et DA CLD10] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDACT/ATM/Output sorties MO et DA LMDACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDIACT/ATM/Output sorties MO et DA LMDIACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDPRE/ATM/Output sorties MO et DA LMDPRE] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LMDIPRE/ATM/Output sorties MO et DA LMDIPRE] || === flag_aerosol=6 === * Lecture AIBCM, AIPOMM, ASBCM, ASPOMM, ASSSM, ASSO4M, CIDUSTM, CSSO4M, CSSSM, SO4, SSSSM : fichiers interpollés par Fred et mis au CCRT là : * /dmnfs/cont003/p86maf/IGCM/BC/ATM/LMDZOR/LMD9695/AR5/HISTORIQUE/XXX.run1980.nc -> XXX.2000.2000.96x95.nc * XXX.run.nat.nc -> XXX.1750.1750.96x95.nc * 2 tests (1 an), lev_histday=3, lev_histmth=4 pour avoir les diagnostics aerosols * LM6DACT : ok_ade=y, period=actuel * LM6DIACT : idem + ok_aie=y * Temps calcul : * LM6DACT : 1150 s avec diagnostics (elapsed time) * LMD6IACT : 1180 s avec diagnostics (elapsed time) * Accès aux fichiers || actuel |||| preind (a venir) |||| || ok_ade=y ok_aie=n || ok_ade=y ok_aie=y || ok_ade=y ok_aie=n || ok_ade=y ok_aie=y || || LM6DACT || LM6DIACT || LM6DPRE || LM6DIPRE || || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DACT/MONITORING monitoring LM6DACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DIACT/MONITORING monitoring LM6DIACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DPRE/MONITORING monitoring LM6DPRE] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DIPRE/MONITORING monitoring LM6DIPRE] || || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DACT/ATM/Output sorties MO et DA LM6DACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DIACT/ATM/Output sorties MO et DA LM6DIACT] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DPRE/ATM/Output sorties MO et DA LM6DPRE] || [http://dods.extra.cea.fr/data/p86maf/LMDZOR/LM6DIPRE/ATM/Output sorties MO et DA LM6DIPRE] || == Tests en couplé == * Sources : * IPSLCM5_v1 * LMDZ 1196 * Orchidee 1 9 4 * NEMO * Avec les mêmes réglages LMDZ que CLD10. === 96x95x39 === * SX9 * C5993A * 10 ans finis * 100 ans en route (8 juillet 2009) * 4 cpus , 1000 s par mois environ (elapsed time), 9 Go de mémoire * C'est là : /work/cont003/p86maf/CMIP5/modipsl/config/IPSLCM5/C5993A * Voir suivi là : [http://dods.extra.cea.fr/data/p86maf/IPSLCM5/C5993A/MONITORING/ Monitoring C5993A] * Arrêté après 40 ans car c'est le même que : http://dods.extra.cea.fr/data/p86fair/IPSLCM5/LR04/MONITORING/ === 144x142x39 === * SX9 * C514143A * 10 ans finis * 100 ans en route (8 juillet 2009) * 8 cpus, 1500 s par mois environ (elapsed time), 22 Go de mémoire * C'est là : /work/cont003/p86maf/CMIP5/C514143/modipsl/config/IPSLCM5/C514143A * Voir suivi là : [http://dods.extra.cea.fr/data/p86maf/IPSLCM5/C514143A/MONITORING/ Monitoring C14143A]