wiki:ResolutionIPSLCM5

Tests IPSLCM5 avec 39 niveaux verticaux

Contexte 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*
  • 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
X8L1900 96x95x19 480 N/N/N/N 5400(it=2) 7200(it=2) 7200(it=2) 0.02 -0.7W/m² 721 mois && STOP
X8L3900 96x95x39 480 Y/N/N/N 5400(it=2) 7200(it=2) 7200(it=2) 0.02 -7.3W/m² 82 mois && HGARDFOU
P6L3900 96x95x39 480 Y/N/N/N 5400(it=2) 7200(it=2) 7200(it=2) 0.02 -7.6W/m² 348 mois && HGARDFOU
X8L3901 96x95x39 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
X8L3902 96x95x39 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
X8L3903 96x95x39 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
X8L3904 96x95x39 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
X8L3905 96x95x39 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
X8L3906 96x95x39 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
X8L3907 96x95x39 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
X8L3908 96x95x39 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
X8L3909 96x95x39 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
X8L3910 96x95x39 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
X8L3911 96x95x39 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
X8L3912 96x95x39 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
X8L3913 96x95x39 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
X8L3914 96x95x39 720 Y/N/N/N 5400(it=2) 7200(it=2) 7200(it=2) 0.02 -7.0W/m² 139 mois && HGARDFOU
X8L3915 96x95x39 480 Y/N/N/N 3600(it=2) 5400(it=2) 5400(it=2) 0.02 -7.7W/m² 320 mois && HGARDFOU
X8L3916 96x95x39 480 Y/N/N/N 5400(it=1) 7200(it=2) 7200(it=2) 0.02 -8.4W/m² 157 mois && HGARDFOU
X8L3917 96x95x39 480 Y/N/N/N 3600(it=1) 7200(it=2) 7200(it=2) 0.02 -8.4W/m² 128 mois && HGARDFOU
X8L3918 96x95x39 480 Y/N/N/N 5400(it=1) 5400(it=2) 5400(it=2) 0.02 -8.4W/m² 477 mois && HGARDFOU
X8L3919 96x95x39 480 Y/N/N/N 5400(it=1) 5400(it=2) 7200(it=2) 0.02 -8W/m² 19 mois && HGARDFOU
X8L3920 96x95x39 480 Y/N/N/N 5400(it=1) 7200(it=2) 5400(it=2) 0.02 -8.3W/m² 283 mois && HGARDFOU
X8L3921 96x95x39 480 Y/Y/Y/N 3600(it=2) 5400(it=2) 5400(it=2) 0.02 -7.5W/m² 2400 mois & Completed
X8L3922 96x95x39 480 Y/N/Y/N 3600(it=2) 5400(it=2) 5400(it=2) 0.02 -7.5W/m² 799 mois & HGARDFOU
X8L3923 96x95x39 480 Y/Y/Y/N 5400(it=1) 5400(it=2) 5400(it=2) 0.02 -8.1W/m² 2400 mois & Completed
X8L3924 96x95x39 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
X9L3925 96x95x39 480 Y/Y/Y/N 5400(it=1) 5400(it=2) 5400(it=2) 0.0cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83 -0.10W/m² 2400 mois && Completed
X9L3926 96x95x39 480 Y/N/Y/N 5400(it=1) 5400(it=2) 5400(it=2) 0.0cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83 0.17W/m² 127 mois && HGARDFOU
X9L3927 96x95x39 480 Y/Y/Y/Y 5400(it=1) 5400(it=2) 5400(it=2) 0.0cld_tau_[lsc][con]=1800.+ffalv_[lsc][con]=0.83 -0.05W/m² 2400 mois && Completed

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²).
  • 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 : 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
  • Le passage à nitergdiv=1 (la préférence de Fred) fait perdre 0.6 W/m²
  • On peut imaginer que les combinaisons suivantes donnent un flux net au sommet de 0W/m²
    • nitergdiv=2 + cld_tau_[lsc][con]=1800. + ffalv_[lsc][con]=0.66 + pmgagic=0.02
    • nitergdiv=2 + cld_tau_[lsc][con]=1800. + ffalv_[lsc][con]=0.75 + pmgagic=0.01
    • nitergdiv=2 + cld_tau_[lsc][con]=1800. + ffalv_[lsc][con]=1.0 + pmgagic=0.0
    • nitergdiv=1 + cld_tau_[lsc][con]=1800. + ffalv_[lsc][con]=0.83 + pmgagic=0.0
  • La combinaison retenue :
    • solaire=1361 + nbapp_rad=24 + epmax=0.995 + nitergdiv=1 + cld_tau_[lsc][con]=1800. + ffalv_[lsc][con]=0.92 + pmgagic=0.008

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
  • Extension de la couche éponge, la couche éponge est active dans les couches de pression plus faible que 100 fois la pression de la dernière couche : iflag_top_bound=2

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). Les fichiers actuels et pré-industriels sont sur les comptes communs. Cf simulations ci dessous.
  • 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 : Un pmagic de 0.008 dixit Jean Louis. Soit 0.8 W/m² en terme de flux net au sommet.
  • Réadapter les réglages (à pmagic constant, avec tous les aérosols) une fois inclus les nouveaux coefficients cdrag (Pascale, Frédéric)

Simulations de juillet 2009 : troisième round d'équilibrage/réglage

  • Sources :
    • LMDZ 1196
    • Orchidee 1_9_4
    • 96x95x39
    • C'est là : /work/cont003/p86denv/PARA_SX9_CM5_R1196_39L/modipsl/config/LMDZOR
  • Rappel de la combinaison retenue :
    • solaire=1361 + nbapp_rad=24 + epmax=0.995 + nitergdiv=1 + cld_tau_[lsc][con]=1800. + ffalv_[lsc][con]=0.92 + pmgagic=0.008

Objectifs principaux

  • On part du nouveau réglage issue des tests précédents (actuellement sur svn pour LMDZOR et IPSLCM5_v1)
  • Le job de référence (REFERENCE) et le test CLD10 sont identiques. ÉTONNANT. On a déjà vu que rad_froid n'avait pas d'effet, mais rad_chau[1][2] devaient en avoir. Bug dans le script de création d'ensembles. Rerun CLD10.
RUN atlas+ .defStatusfnettopstoplcrfswcrfnetevapr
REFERENCE X Completed 0.2523 247.0 246.7 -41.34 -17.53 2.973 2.974
DRAG1 X f_cdrag_oce=1. Completed -1.389 246.3 247.7 -42.22 -18.58 3.134 3.134
DRAG2 X f_cdrag_oce=1. cld_lc_lsc=2.e-4 cld_lc_con=2.e-4 Completed 0.5476 248.4 247.9 -40.11 -16.67 3.127 3.128
DRAG3 X f_rugoro=1. Completed 0.08636 246.6 246.6 -41.54 -17.61 2.980 2.980
DISS1 X tetagdiv=7200. Hgardfou 0.4293 247.0 246.5 -41.43 -17.43 2.970 2.970
DISS2 X tetagrot=18000. Hgardfou 0.5838 246.9 246.3 -41.48 -17.27 2.980 2.980
DISS3 X tetatemp=18000. 2400 mois && Completed 0.08606 246.5 246.4 -41.81 -17.91 2.963 2.964
DISS4 X tetagrot=10800. Hgardfou 0.1699 246.5 246.4 -41.76 -17.70 2.975 2.976
DISS5 X tetagrot=18000. Hgardfou 0.3514 246.4 246.0 -41.90 -17.58 2.970 2.971
DISS6 X dissip_factz=2. Hgardfou 0.3809 247.0 246.6 -41.41 -17.51 2.971 2.972
DISS7 X tau_top_bound=.5e-5 Hgardfou . . . . . . .
CLD01 X cld_tau_lsc=2700. cld_tau_con=2700. Completed -2.502 243.6 246.1 -44.69 -20.49 2.979 2.980
CLD02 X cld_lc_lsc=2.e-4 cld_lc_con=2.e-4 Completed 2.084 249.1 247.0 -39.36 -15.75 2.963 2.964
CLD03 X cld_lc_lsc=1.5e-4 cld_lc_con=1.5e-4 Completed 3.592 250.8 247.2 -37.61 -14.23 2.958 2.959
CLD04 X cld_lc_lsc=1.5e-4 cld_lc_con=1.5e-4 cld_tau_lsc=2700. cld_tau_con=2700. Completed 0.5260 247.3 246.8 -41.13 -17.37 2.973 2.973
CLD05 X ffallv_lsc=0.7 ffallv_con=0.7 Completed 1.130 243.9 242.8 -44.49 -16.81 2.914 2.914
CLD06 X coef_eva=1.e-5 Completed 1.612 248.7 247.1 -39.70 -16.23 2.979 2.979
CLD07 X ratqsbas=0.01 Completed 1.092 247.9 246.9 -40.45 -16.77 2.960 2.961
CLD08 X ratqshaut=0.15 Completed 1.073 239.2 238.2 -49.39 -18.42 2.906 2.907
CLD09 X ratqsbas=0.01 ratqshaut=0.15 Completed 2.447 240.8 238.4 -47.78 -16.96 2.892 2.892
CLD10 X rad_chaud1=14 rad_chaud2=14 Completed 0.6547 247.4 246.7 . . 2.974 2.975
CLD11 X epmax=0.9975 Completed 1.367 249.8 248.4 -38.57 -16.53 2.997 2.997
CLD12 X epmax=0.9975 ffallv_lsc=1.1 ffallv_con=1.1 Completed 1.104 251.7 250.6 -36.64 -16.75 3.029

Tests aerosols n°1

Tests aerosols n°2

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
  • 4 tests 10 ans, lev_histday=3, lev_histmth=4 pour avoir les diagnostics aerosols
  • Temps calcul /!\ Attention à l'impact de la charge SX9 sur les perfs /!\ :
    • REFERENCE (=CLD10 chez p86maf) : 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)
RUN .def Statusfnettopsadtopsaitopstoplevapr
REFERENCE Completed 0.2523 . . 247.0 246.7 2.973 2.974
LMDACT00 ok_ade=y, period=actuel Completed -0.55 -0.8216 . 246.2 246.7 2.974 2.975
LMDIACT00 ok_ade=y, ok_aie=y, period=actuel Completed -1.59 -0.8235 -1.557 245.1 246.7 2.964 2.965
LMDPRE00 ok_ade=y, period=preind Completed -0.01 -0.2568 . 246.7 246.7 2.973 2.974
LMDIPRE00 ok_ade=y, ok_aie=y, period=preind Completed -0.76 -0.2586 -1.288 245.9 246.7 2.970 2.971

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
  • 4 tests (10 ans), lev_histday=3, lev_histmth=4 pour avoir les diagnostics aerosols
  • Temps calcul :
    • LM6DACT : 1150 s avec diagnostics (elapsed time)
    • LM6DIACT : 1180 s avec diagnostics (elapsed time)
RUN .def Statusfnettopsadtopsaitopstoplevapr
REFERENCE Completed 0.2523 . . 247.0 246.7 2.973 2.974
LM6DACT00 ok_ade=y, period=actuel Completed -0.73 -1.003 . 246.0 246.8 2.970 2.971
LM6DIACT00 ok_ade=y, ok_aie=y, period=actuel Completed -1.84 -1.010 -1.688 244.9 246.8 2.960 2.962
LM6DPRE00 ok_ade=y, period=preind Completed -0.08 -0.4436 . 246.6 246.7 2.969 2.970
LM6DIPRE00 ok_ade=y, ok_aie=y, period=preind Completed -1.14 -0.4482 -1.449 245.7 246.8 2.967 2.968

Tests aerosols n°2 Bis en 19 niveaux

RUN .def Statusfnettopsadtopsaitopstoplevapr
AR4REFERENCE period=actuel Completed 1.319 . . 241.5 240.1 2.775 2.776
AR4DACTNEW ok_ade=y, new_aod=y Completed 0.4403 -0.8405 . 240.6 240.2 2.776 2.777
AR4IACTNEW ok_ade=y, ok_aie=y, new_aod=y Completed 0.3041 -0.8474 -1.565 240.5 240.2 2.775 2.776
AR4IACTNEW_OA ok_ade=y, ok_aie=y, new_aod=n Completed 0.3211 -1.030 -1.563 240.4 240.0 2.773 2.774
AR4DACTOLD ok_ade=y, new_aod=y Completed 0.6653 -0.6724 . 240.8 240.1 2.773 2.774
AR4IACTOLD ok_ade=y, ok_aie=y, new_aod=y Completed 0.3893 -0.6802 -1.831 240.5 240.1 2.772 2.773
AR4IACTOLD_OA ok_ade=y, ok_aie=y, new_aod=n Completed 0.1421 -0.8397 -1.826 240.2 240.1 2.772 2.773
RUN .def period=preind Statusfnettopsadtopsaitopstoplevapr
AR4DPRENEW ok_ade=y, new_aod=y Completed 1.103 -0.2647 . 241.2 240.1 2.774 2.775
AR4IPRENEW ok_ade=y, ok_aie=y, new_aod=y Completed 1.129 -0.2674 -1.315 241.3 240.2 2.774 2.775
AR4DPREOLD ok_ade=y, new_aod=y Completed 1.131 -0.1884 . 241.3 240.1 2.773 2.774
AR4IPREOLD ok_ade=y, ok_aie=y, new_aod=y Completed 1.109 -0.1915 -1.579 241.2 240.1 2.772 2.773

Tests aerosols n°3

  • Sources :
    • merge de LMDZ 1242 avec modifs récentes de Michael, Yves, Raffaella et moi.
      • aeropt_2bands.F90 : celle de Michael (25/08/09) avec tableaux sels marins et modifs pour faire SAVE sur les tableaux A1_ASSSM_b1
      • aeropt_5wv.F90 : celle de Yves (25/08/09) avec tableaux sels marins et modifs pour faire SAVE sur tableaux
      • hines_gwd.F : suppression de print
      • newmicro.F : déplacement de la colonne 73 du *1.E6 (!)
      • physiq.F : ajout de itap dans l'interface d'appel a readaerosol_optic + diag
      • readaerosol_interp.F90 : modification dans calcul de date
      • readaerosol_optic.F90 : modification dans calcul date pour interpolation
      • sw_aeroAR4.F90 : modification des calculs de diagnostics (2 fois)
      • radlwsw.F90 : calculs diags conditionnés
      • phys_output* diags
    • Orchidee 1 9 4
    • 96x95x39
    • C'est là : /work/cont003/p86denv/PARA_SX9_CM5_R1242_39L/modipsl/
RUN .def Statusfnettopsadtopsaitopstoplevapr
REFERENCE05 Completed 1.548 . . 247.6 246.0 2.950 2.951
LM6DIACT05 ALL aerosols, ok_ade=y, ok_aie=y, period=actuel Completed -1.265 -0.4860 -0.2503 245.2 246.5 2.924 2.925
LM6DIPRE05 ALL aerosols, ok_ade=y, ok_aie=y, period=preind Completed -0.4501 0.0000 2.459E-17 245.8 246.2 2.929 2.930
LMDIACT05 SO4 Only, ok_ade=y, ok_aie=y, period=actuel Completed -0.1158 -0.5631 -0.2782 245.8 245.9 2.939 2.941
LMDIPRE05 SO4 Only, ok_ade=y, ok_aie=y, period=preind Completed 0.7204 0.0000 1.082E-16 246.8 246.0 2.942 2.944

Tests aerosols n°4

  • Sources :
    • LMDZ 1246
    • Orchidee 1 9 4
    • 96x95x39
    • Réglage pour le milieu de la décennie 1995-2005 avec les ghg moyens suivants :

CO2 : 367.00, CH4 : 1760.00, N2O : 316.00, CFC11 : 741,20, CFC12 : 535.00.

  • C'est là : /work/cont003/p86denv/PARA_SX9_CM5_R1242_39L/modipsl/
RUN .def Status
GHGREF Completed
GHGREF1 epmax=0.997 Completed
GHGREF2 epmax=0.999 Completed
GHGREF3 epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed
GHGREF4 epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
GHGREF5 epmax=0.999 ffallv_lsc=1.3 ffallv_con=1.3 Completed
LM6DIACT06 ok_ade=y ok_aie=y epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LM6DIPRE06 ok_ade=y ok_aie=y aer_type=preind epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LMDIACT06 ok_ade=y ok_aie=y flag_aerosol=1 epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LMDIPRE06 ok_ade=y ok_aie=y flag_aerosol=1 aer_type=preind epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LM6DIACT07 ok_ade=y ok_aie=y epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LM6DIPRE07 ok_ade=y ok_aie=y aer_type=preind epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LMDIACT07 ok_ade=y ok_aie=y flag_aerosol=1 epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LMDIPRE07 ok_ade=y ok_aie=y flag_aerosol=1 aer_type=preind epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LM6DIACT08 ok_ade=y ok_aie=y epmax=0.999 ffallv_lsc=1.3 ffallv_con=1.3 Completed
LM6DIPRE08 ok_ade=y ok_aie=y aer_type=preind epmax=0.999 ffallv_lsc=1.3 ffallv_con=1.3 Completed
LMDIACT08 ok_ade=y ok_aie=y flag_aerosol=1 epmax=0.999 ffallv_lsc=1.3 ffallv_con=1.3 Completed
LMDIPRE08 ok_ade=y ok_aie=y flag_aerosol=1 aer_type=preind epmax=0.999 ffallv_lsc=1.3 ffallv_con=1.3 Completed

Tests aerosols n°5

  • Sources :
    • LMDZ 1246
    • Orchidee 1 9 4
    • 96x95x39
    • Réglage pour le milieu de la décennie 1995-2005 avec les ghg moyens suivants :

CO2 : 367.00, CH4 : 1760.00, N2O : 316.00, CFC11 : 741,20, CFC12 : 535.00.

  • C'est là : /work/cont003/p86denv/PARA_SX9_CM5_R1242_39L/modipsl/
  • Nouvelle carte d'aérosol actuel de Sophie Szopa
RUN .def Status
GHGREF3 epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed
GHGREF4 epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
GHGREF5 epmax=0.999 ffallv_lsc=1.3 ffallv_con=1.3 Completed
LM6DIACT06B GHGREF4 + ok_ade=y ok_aie=y epmax=0.999 ffallv_lsc=1.1 ffallv_con=1.1 Completed
LM6DIACT07B GHGREF3 + ok_ade=y ok_aie=y epmax=0.997 ffallv_lsc=1.1 ffallv_con=1.1 Completed

Tests sur les nuages en forcé : réglage finale

  • On est d'accord sur le diagnostic du réglage LM6DIACT08B:

la croissance de epmax (une décroissance de 1-epmax qui est en gros la fraction de l'eau montée par la convection qu'on relâche dans l'environnement en haut des tours convectives) et la croissance simultanée de fallv (facteur sur la vitesse de chute de la glace) on fait décroître déraisonnablement les nuages hauts et donc simultanément les CRFsw et CRFlw des nuages.

Y'a pas de doute qu'on a un malaise du cote des CRFs.

Dans les regimes convectifs, le LW CRF et le SW CRF sont tjrs trop faibles, d'autant plus que epmax est eleve (ce qui est attendu puisqu'on reduit les nuages hauts en augmentant epmax). La simul GHGREF (ou REFERENCE05) etait encore la meilleure.

Augmenter la vitesse de chute de la glace ne fait qu'agraver les choses puisqu'on reduit encore le contenu en eau des nuages hauts.

Si l'on croit aux obs HOAPS de flux latents en surface, alors l'evap sur oceans parait aussi systematiquement surestimee. Ca pourrait etre bien de regler les drags avec ça.

Les flux radiatifs ciel-clair TOA sont plutot une bonne surprise (surtout si l'on considere que le jeux CERES-EBAF est mieux que CERES-std).

La SST aussi (sauf dans la warm pool ou elle semble tjrs trop froide malgre le manque de SW CRF).

  • Sources :
    • LMDZ 1273
    • Orchidee 1 9 4_AR5
    • 96x95x39
    • Réglage pour le milieu de la décennie 1995-2005
    • C'est là : /work/cont003/p86denv/PARA_SX9_CM5_R1273_39L/modipsl/config/LMDZ4OR_v3
    • Nouvelle carte d'aérosol actuel de Sophie Szopa
    • Ozone jour/nuit avec les fichiers de forçages de David Cugnet
RUN .def Status
CLOUD01 epmax=0.999 ffallv_lsc=0.7 ffallv_con=0.7 Completed
CLOUD02 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 Completed
CLOUD03 epmax=0.999 ffallv_lsc=0.3 ffallv_con=0.3 Completed
CLOUD04 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 ratqshaut=0.15 Completed
CLOUD05 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 ratqshaut=0.66 Completed
CLOUD06 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 ratqshaut=0.15 cld_tau_lsc=2700. cld_tau_con=2700. Completed
CLOUD07 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 ratqshaut=0.66 cld_tau_lsc=2700. cld_tau_con=2700. Completed
CLOUD08 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 cld_tau_lsc=3600. cld_tau_con = 3600. Completed
CLOUD09 epmax=0.999 ffallv_lsc=0.5 ffallv_con=0.5 cld_lc_lsc=4.16e-4 cld_lc_con = 4.16e-4 Completed

Cf table CRFs : http://dodsp.idris.fr/rlmd093/RUNS_SEB/CLOUD/HTML/tab_simulations.html

Moyennes zonale : http://dodsp.idris.fr/rlmd093/RUNS_SEB/CLOUD/MoyZon/EXP3/MoyZonAn.pdf

monitoring : http://igcmg.lsce.ipsl.fr/monitoring/tmp/fegg_plot01_AN3Cn7_prod/

Atlas (avec les coupes méridiennes !) : http://dods.extra.cea.fr/data/p86denv/LMDZOR/CLOUD09/ATLAS/SE_1980_1989/ATM/ATM.html

En lien ci dessous, les diag en regime des dernieres simul CLOUD (pour gagner en lisibilite, j'ai retire des graphes les simuls qui etaient trop dans les choux). Pour completer ce que Ionela a fait sur les atlas, j'ai aussi ajoute les obs sur les fractions nuageuses en regimes.

Les simul CLOUD08 et CLOUD09 sont effectivement tres chouettes dans les tropiques. La plupart des variables sont ameliorees par rapport a CLOUD02.

CLOUD09 a l'air legerement meilleure que CLOUD08, mais les deux sont qd meme tres proches l'une de l'autre (par rapport au reste de la famille CLOUD).

Eventuellement, si on peut essayer de pousser encore un peu plus les potards des nuages bas (augmenter cld_lc_con/lsc a 2.6e-4 * 2?)... mais le mieux est souvent l'ennemi du bien et en plus on n'a peut-etre deja plus le temps...

En tout cas, ces dernieres simul ont l'air vraiment sympa.

Tests de stabilité HGARDFOU sur la dissipation en 144x142x39

ISSUE DU DEUXIÈME ROUND D'ÉQUILIBRAGE Résolutions Strato/Hines/NewDiffusion/TopBound? tetagdiv tetagrot tetatemp pmagic Flux net au sommet Durée
HX8L3925 144x142x39 Y/Y/Y/N 3600(it=1) 5400(it=2) 5400(it=2) 0.0 -6.2W/m² 2400 mois && COMPLETED
  • Sources :
    • LMDZ4-dev 1215
    • Orchidee 1_9_4
    • 144x142x39
    • C'est là : /work/cont003/p86denv/PARA_SX9_CM5_HI_R1215_39L/modipsl/config/LMDZOR
    • Synthèse en images
  • Une série de simulations pour tourner autour des paramètres de dissipation qui "tiennent", est ce que l'on a de la marge?
RUN Résolutions tetagdiv,tetagrot,tetatemp Status fnettopstoplevapr
HIREFERENCE 144x142x39 3600(it=1),5400(it=2),5400(it=2) 1182 mois && HGARDFOU -0.04514 247.2 247.3 2.986 2.987
HIDISS1 144x142x39 tetagdiv=7200. 384 mois && HGARDFOU 0.4852 . . . .
HIDISS2 144x142x39 tetagrot=18000. 475 mois && HGARDFOU 0.9409 . . . .
HIDISS3 144x142x39 tetatemp=18000. 1268 mois && HGARDFOU -0.2174 . . . .
HIDISS4 144x142x39 tetagrot=10800. tetatemp=10800. 92 mois && HGARDFOU 0.2337 . . . .
HIDISS5 144x142x39 tetagrot=18000. tetatemp=18000. 578 mois && HGARDFOU 0.09957 . . . .
HIDISS6 144x142x39 nitergdiv=2 tetagdiv=7200. 36 mois && HGARDFOU 2.004 . . . .
HIDISS7 144x142x39 nitergdiv=2 tetagrot=18000. 291 mois && PS négative 1.567 . . . .
HIDISS8 144x142x39 nitergdiv=2 tetatemp=18000. 535 mois && HGARDFOU 1.117 . . . .
HIDISS9 144x142x39 nitergdiv=2 tetagrot=10800. tetatemp=10800. 264 mois && DEXP -> ARG 1.441 . . . .
HIDISS10 144x142x39 nitergdiv=2 tetagrot=18000. tetatemp=18000. 464 mois && HGARDFOU 1.572 . . . .
HIDISS11 144x142x39 dissip_factz=2. 93 mois && HGARDFOU 0.09482 . . . .
HIDISS12 144x142x39 tau_top_bound=.5e-5 128 mois && HGARDFOU 0.06558 . . . .

Tests couplés sur le journal de bord lmdz

Tests en couplé

  • Sources :
    • IPSLCM5_v1 : branches/IPSLCM5_WORK 678
    • LMDZ 1196
    • Orchidee 1 9 4
    • NEMO 1340
  • Avec les mêmes réglages LMDZ que REFERENCE.

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à : 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/
  • En fait non superposable car LR04 est avec orchidee 1 9 2. Poursuivi au delà des 40 ans.
  • Plantage au mois 631, 52 ans, 1912-08

144x142x39

  • SX9 - C514143A
    • tetagdiv=5400. (it=1)
    • 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à : Monitoring C14143A
    • Plantage au mois 414 (34 ans) le 12 juillet
  • SX9 - C514143B
    • idem avec tetagdiv=3600. (it=1) comparable à HX8L3925
    • Voir : /work/cont003/p86maf/CMIP5/C514143/modipsl/config/IPSLCM5/C514143B
    • Lancé le 10 août 2008. Plantage 40 ans.
  • SX9 - C514143C
    • idem C514143B
    • tetagdiv=1800 et tetagrot=2700
    • plantage mois 283, 23 ans, 1883-07
  • SX9 - C514143D
    • idem C514143B
    • tau_top_bound=2.e-5
    • en machine, a fait 50 ans
  • SX9 - C514143E
    • idem C514143B
    • iflag_top_bound=0
    • plantage hgardfou au mois 60

Tests d'utilisation de COSP

Marche dans le couplé IPSLCM5_v2 depuis la revision 1332 de LMDZ.

Pour utiliser COSP, il faut :

  • le demander à la compilation. Ajouter -cosp true dans config/IPSLCM5/Makefile
  • avoir les 2 fichiers : PARAM/cosp_input_nl.txt et PARAM/cosp_output_nl.txt
  • avoir mis dans PARAM/physiq_L39 : ok_cosp=y
  • avoir ajouté dans COMP/lmdz.card les 3 fichiers (sorties) : histmthCOSP, histdayCOSP et histhfCOSP

Test fait sur brodie le 25 mars 2010 :

  • Taille des fichiers de sorties pour le mois de janvier :
    • 1786198596 : Output/DA/COSP3_18600101_18600131_1D_histdayCOSP.nc
    • 7144713516 : Output/HF/COSP3_18600101_18600131_HF_histhfCOSP.nc
    • 57645396 : Output/MO/COSP3_18600101_18600131_1M_histmthCOSP.nc
  • Quelques chiffrages :
    • 2000s de temps réel par mois sur brodie (3 procs atm) pour 1200 s sans COSP.
    • Mémoire 11GB (9 GB sans COSP).
Last modified 12 years ago Last modified on 03/26/10 15:39:50