28/07/2010

Point sur la nouvelle physique du 28/07/2010

Le point sur la nouvelle physique fin Juillet 2010, et discussion autour des priorités pour la suite.

Climatologie

Points positifs : cycle diurne de la convection, nuages bas et moyens, moussons (plutôt un peu meilleur).
Moyennant réglage : bonne répartition en latitude et en régime des forçages radiatifs des nuages.
Rm : même avec alp_offset=0, les résultats sont très raisonnables avec les dernières versions (cf NPE1 et CMP28). Ca devrait permettre des
tests passionnants.

Problèmes en forcé : pluies trop fortes, humidité relative à 600 hPa trop forte

Problèmes en couplé : circulation thermohaline atlantique quasi inexistante
El Nino faible, mais peut-être pour de bonne raisons.

La haute résolution améliore beaucoup la position des jets, avec une forte diminution des biais en T et U sur les moyennes zonales.

Elle pourrait sans doute aussi améliorer El Nino.
Sans doute pas la THC ...

Informatique

1+1=2 non testé.
Non reproductibilité des simulations sur vargas, peut-être due simplement au fait que adjust=y (mais je n'ai pas eu le temps de tester avec n)
QUI PEUT FAIRE CA ?
Modèle sans doute insuffisamment optimisé.
Performance à contrôler.
En particulier, le modèle semble plus rapide sur vargas sur 32 procs que sur SX9 sur 8.

Travail et simulations en cours

Tests sur les paramètres contrôlant la convection profonde en 1D (Toga long) et
en parallèle en 3D (à faire, et ne pas se limiter aux quelques tests qui
semblent optimum pour le 1D).
Test en cours avec les .def récents sur les cas Eurocs_cu, bomex et eurocs_FG
pour voir si on retrouve bien les résultats précédemment publiés.
Sorties Cosp et comparaison aux obs Calipso sur les dernières simulations
(Abderrahmane).
Comprendre pourquoi la THC s'effondre à ce point (Jean-Louis et al.).

Quelques paramètres aident à rééquilibrer le modèle

Quand on change des choses, il reste toujours à la fin une phase de rééquilibrage. En général pour quelques W/m2.
Ca sera surement le cas, une fois choisi un novueau réglage pour les paramètres de la convection.
Il va falloir à nouveau refroidir le système quand on passe à la haute résolution par exemple.

Si on veut faire décroître le flux net au sommet (refroidir l'atmsophère), on peut
- diminuer epmax (l'éloigner de 1). Ca rajoute des nuages, donc du CRFSW
et LW mais ca fait globalement diminuer le net. Sandrine voudrait
qu'on soit dans les environ de 0.999. On est plutôt à 0.995 dans les
derniers réglages.
- augmenter fallv_... pour faire tomber plus vite les cristaux. Pas
évident comme conclusion mais efficace.
- augmenter ratqshaut.
- augmenter coef_eva (mais il est sans doute déjà très haut)
- augmenter les cld_tau... et les cld_lc... mais ils sont déjà très haut
aussi.
- diminuer wdens (mettre une moins grande densité de poches froides).
En général, les tests passés sur une autre config, permettent de "viser" assez juste en termes de rééquilibrage. Si passer de ratqshaut=0.1 à 0.2 changeait le net de -2W/m2, ca restera proche dans d'autres configs.

Rm : Tous ces tests ont  finalement assez peu d'effet sur le biais humide à 600hPa dans les tropiques.

Il faut quand on fait ca absolument contrôler qu'on ne détruit pas la distribution des CRFSW et LW en régime et en latitude (les X 1 et 2 dans le tableau).
Le reste bouge moins en général.

Dans les cartons

Couplage Nicolas Rochetin
Relief/convection de Jingmei

Discussion des priorités pour l'été

Par rapport à MH07 qui est toujours en machine, et qui correspondait en forcé à NNP40, on a
- des configurations réglées avec les nuages d'Arnaud (NNP40 était juste
avant)  avec des forçages plutôt meilleurs (cf. NNPE10 et CMP32 pour
comparer les différentes versions réglées de la nouvelle physiuqe,
successivement ANP00, NNP40 et NPE10, et la physique ESM/AR4, AP00).
- des configurations réglées de même en 144x143.

En en discutant avec Jean-Louis, il ne semble pas utile à ce stade de lancer autre chose que des contrôles actuels et des 1 % CO2.

On propose donc de laisser tourner MH07 jusqu'au retour de Laurent, et qu'il lance quand il sera là un 1% CO2 en partant 70 ans en arrière par rapport au point où en sera cette simulation.

On peut donc arrêter la simulation préindustrielle sur SX9.

Que met-on à la place ?
La haute résolution ancienne physique ?
L'équivalent de HNE10 (en basse résolution) ?
La même re-réglée en 142x143 ?

Avant de relancer sur SX9 la nouvelle physique, ce serait bien de regarder sérieusement les perfs.
Sinon on peut aussi penser tourner sur vargas ou on passe un an / 10 heures elapsed en 142x143x39 actuellement.
Est-ce que le couplé peut tourner dessus ?
Est-ce qu'on peut avoir de l'aide sur la mise en place d'un couplé nouvelle physique sur vargas. Ca peut être une chose que Laurent fait quand il passe au mois d'août. Mais ca sera peut-être un peu chaud.

Je suis après re-réflexion convaincu qu'il faut commencer à faire tourner du 142x143 nouvelle physique en couplé d'une façon ou d'une autre cet été.
Ne serait ce que pour voir si ça améliore bien El Nino comme espéré.
Et les biais en latitudes aussi. Ca devrait pouvoir commencer à être vraiment pas mal de ce point de vue.
Ca n'améliorera sans doute pas la convection en atlantique nord, et c'est très important de bosser là dessus aussi.
Et ce sera un élément important pour savoir si on privilégie du 96x95 ou du 142x143 pour la dernière ligne droite.

Peut être faut il lancer en fait, avec la même physique, un couplé 96x95x39 et un 144x143x39. Sur vargas si c'est possible.
Même si n'a qu'une 30 aine d'années en septembre, ca devrait être extrêmement utile.