Exercice de tuning automatique du modele couplé printemps 2023

Reunion du 14 Avril 2023

Point irene: difficultés openMP. Ca semble aller normalement depuis ce matin vendredi 14 avril. Feu vert à venir dans les prochaines heures.

Petit souci irene-skl avec des deadlock impromptus. Ils étaient la avant, et le sont toujours;

Besoin de prévenir le TGCC: Arnaud s'en occupe, cc Julie.

Discussion sur outils de communication: wiki, lien éventuels vers fichiers excel ou autre. On décide de ne pas faire de slack/mattermost/discord supplementaire.

1) Paramètres inclus dans l’exploration

Paramètres NEMO

RNALB RNCDN RNCE RNLC

nom raccourci	min	max	nominal exploration	nom complet
0.49	0.58	0.50 linear	rn_alb_imlt (albedo of melting bare sea ice)
0.10	0.50	0.31	linear	rn_cnd_s (thermal conductivity of the snow over sea ice, W/m/K)
0.06	0.08	0.06	linear	rn_ce
0.05	0.5	0.15	linear	rn_lc

Guillaume Gastineau insiste pour modifier le paramètre albédo de sorte qu'on affecte les 4 valeurs d'albédo et pas 1 seule - nécessite écriture d'un (bout de) script pour transformer le scaling générique en scaling pour chacun des 4 albédos et remplacer

RNALB

0.49

0.58

0.50

linear

rn_alb_imlt (albedo of melting bare sea ice)

par

RNALG

0

1

0.50

linear

facteur multiplicatif pour régler les 4 albédos en même temps

Paramètres LMDZ

FALLV OMEPMX DZ EVAP GKDRAG

CLC	1e-4	1e-3	6.5e-4	linear
0.3	2.	0.8	linear
0.0001	0.1	0.001	log
0.05	0.2	0.07	linear
5e-5	5e-4	1e-4	log
0.2	2	0.6	linear	sso_gkdrag

Paramètres ORCHIDEE

PCENT

0.3

1

0.8

linear

Pcent

Philippe Peylin : "Utiliser le paramètre Pcent avec les valeurs que nous t’avons données ; Ce paramètre aura une influence medium sur le flux de transpiration sur tout le globe. Catherine a aussi suggéré un paramètre du vieillissement de la neige qui aurait un impact fort mais uniquement sur les régions boréales : tcst_snowa qui est défini dans constantes.f90 (décroissance de l’albédo de la neige en fonction de son age) ; Défaut 10 jour, Range de variations plausibles: 5 - 15 jours (linear) ; Attention : il s'agit d'une constante définie avant compilation ; il faudrait donc modifier le code pour externaliser ce paramètre

Patricia vérifie le code et confirme qu'il n'y a pas grand chose à faire, mais préfère confirmer auprès de Josefine la manière de procéder -> commit à inclure dans IPSLCM6.5_work ? ou autrement ?

Attention implique aussi d'ajouter une nouvelle ligne au tableau :

ASNOW

5

15

10

linear

tcst_snowa

-> total 12 paramètres à inclure dans précondtionnement 1D et 3D

2) Préconditionnement 1D et 3D

Revoir total de membres depuis qu'on upgrade à 12 variables... 120 ?

Attention : préconditonnement 3D avec LMDZOR

Étienne : Aussi DYNAMICO-LMDZOR seul ? (souci lié à orographie sous-maille avec GKDRAG)

Julie : Non, LMDZ lonlat exclusif pour préconditionnement.

Étienne : On peut espérer que l'utilisation de métriques radiatives exclusivement, pour préconditionnement, vont rendre l'exercice quand même utile.

Julie : Importance scientifique pour sonder interactions OCE-ATM même avec DYNAMICO pas parfait.

3) Préparation des configurations et des états initiaux

IPSLCM6.5_work

LR

(VLR - Myriam only)

(LR-hybride - Myriam only)

IPSLCM7_work

LR

Préparer Experiments LR : pdcontrol_tuning avec états intiaux pd pour LMDZ et NEMO-ICE et pi pour ORCHIDEE et PISCES -> Christian et Patricia (COMP, PARAM, DRIVER) + vérifier paramètres par défaut :

ln_mle = true (par défaut dans IPSLCM6.5_work et IPSLCM7_work)

do_airsol = False (a vérifier)

nn_etau = 0 ou 1 en fonction de la famille d'ensemble

+ modif pour Carbone (pd et pi en meme temps) (Patricia : rajouter + de details ici please)

4) Scripts de lancement des ensembles de simulation

ins_job -e : produit simulation

sles.sh (en préparation par David Niezgoda) : script qui prend donnée tableau préconditonnement et modifie fichiers dans PARAM de chaque simulation - insérer ici modif pour albédos NEMO-ICE

Action Julie : transmettre à David liste de noms des variables correspondant aux params sortis du precondtionnement - verifier au passage que pas de doublons dans PARAM !

Qsub_ensemble_ppe1.sh : boucle sur toutes les simulations d'une vague, va dans répertoire de chaque simu et lance qsub

+ modifs dans libIGCM pour mettre sorties dans SCRATCH (1 répertoire par ensemble) et non STORE

Comment partager ces scripts et modifs ?

Arnaud : Attention à mettre à jour avec la dernière version de libIGCM (trunk)

+ création d'une branche libIGCM pour tuning_ensemble ? OK dans un premier temps

Documentation pour utiliser ces scripts : en cours de rédaction pour sles.sh ; pour les autres scripts c'est déjà inclus dans doc libIGCM.

5) Gestion des sorties, packing et archivage dans STORE

NEMO

on peut ne pas sortir 3D dans l'océan en modifiant field_def (et EXP/config.xml) pour rajouter des diagnostiques

Christian : ca va poser des problèmes pour le monitoring tel qu'il est ; il faudrait le modifier pour qu'il prenne les nouveaux diagnostiques (t200...) ; est-ce qu'on veut avoir le monitoring pour suivre en live les simulations ? a priori non car simulations très courtes

Priorité : réduire au maximum les sorties 3D de NEMO -> modifications en cours de test par Guillaume Gastineau, à intégrer dans finalisation des configurations par Christian

ATM

histday precip nécessaires pour diagnostiques MJO

OK mais bien supprimer les autres histday pour réduire poids des sorties

discussion sur différentes solutions pour modifier... xml ?

Guillaume Gastineau fait les modifications nécessaires et envoie info à Christian pour finalisation configs

en même temps validation par Juliette et Fred - inclure conclusion du point (6) ci-dessous [u,v,t,q en 3D en mensuel + 2D classiques (inclus flux, precip, nuages...)]

ORCHIDEE

Patricia souligne que si diagsHF (3h) sortis alors utile pour finaliser développements ORCHIDEE (simus offline)

utile surtout pour LR a priori

Attention : préparer un pack à la volée pour réduire nombre d'inodes

PISCES ?

Patricia et Christian réfléchissent aux variables à conserver... sachant qu'on ne fait tourner que 20 ans en tout -> moyennes annuelles ?

6) Script de post-processing pour calculer les métriques

scripts dépendent des sorties des simulations (mais déjà modulable) ; pour l'instant surtout metriques pour ocean (inclus variables atmosphériques), en grand nombre

au contraire, coté atmosphère, il y a deja eu une forte réduction des métriques pour se concentrer sur l'essentiel (ie surtout métriques radiatives), ce qui exclut des variables pertinentes pour l'océan (vent...)

-> être prudent coté sorties atmosphériques pour conserver variables nécessaires pour bien qualifier / discriminer les simus -> u,v,t,q en 3D en mensuel + 2D classiques (inclus flux, précip, nuages...)

7) Qui fait quoi ?

Nom	mail	production des vagues ?
Myriam Khodri	myriam.khodri@…	VLR et LR hybride, sur Irene ROME
Juliette Mignot	juliette.mignot@…	LR sur Irene SKL
Guillaume Gastineau	guillaume.gastineau@…	LR sur Irene SKL
Arnaud Caubel	arnaud.caubel@…	ICO sur SKL
Sebastien Fromang	sebastien.fromang@…	ICO sur SKL