17/02/2014
Compte Rendu du Point Hebdomadaire LMDZ du 17/02/2014
CR POIHL 2014/02/17
-------------------
Durée : 1h25
Joyeux anniversaire Lluis et merci pour les chouquettes.
Ehoaurn :
---------
Tests OK sauf compilo ada.
Ce serait bien de faire une testing.
Jean-Baptiste :
---------------
A continué les affaires Antartique.
Est tombé sur des problèmes avec grilles_gcm... Va faire du nettoyage.
A fait des simulations guidées sur un mois.
Veut passer sur un an.
Jean-Yves :
-----------
Travaille sur les corrections sur les traceurs. Changent les résultats météo.
* A ajouté un flag.
A passé du temps à faire des diagnostics de conservation.
* ferret est complètement 32 bits.
* a fait des scripts de conservation dans grads (64 bits).
* Seule fisrtilp conserve l'eau.
Propose une règle pour la signature des commentaires des corrections.
Jean-Yves fait une propostion pour la prochaine fois.
Alina :
-------
Post-doc. Au LMD jusuq'à fin juillet.
Marie-Alice :
-------------
Accelerateur sur IPSLCM6-RC0
Arnaud a commis sur LMDZ pour Oasis-MCT. Merci !
Marie-Pierre :
--------------
A travaillé avec Laurent sur RRTM pour commission.
A voulu faire un dernier test. Laurent a rajouté une clef CPP_RRTM
Jessica a donné ce qu'elle avait fait sur cas Cgils.
Ionela :
--------
A fait avancé les simulations à différentes résolution.
192x192x39 et 280x192x39 avancent mais se plantent systématiquement à cause du workdir de ada.
Rm d'Ehouarn : il existe une option qui permet de regrouper les messages des differents processeurs dans un seul listing, l'option "-labelio yes" de poe:
poe ./mon_beau_gcm.e -labelio yes
Rm Lionel: Que l'on peut combiner avec l'option stdoutmode et un script, cf message
au Poihl du 19 septembre 2013 :
https://manne.lmd.jussieu.fr/mailman/private/poihl/2013-September/003205.html
Tourne par an.
Lluis :
-------
WRF does not want to work bellow 10 Pa.
Frédéric
--------
- A réintégré une version des thermiques d'Arnaud Jam qui permet de faire des strato-cumulus sans trop casser les petits cumulus.
- A sorti le dsigmin contrôlant l'épaisseur de la première couche
- A ébauché une discrétisation à 59 niveaux avec première couche vers 8 m qui permet d'obtenir des bons résultats avec la version du code d'Arnaud sur le cas Fire.
Y compris avec un pas de temps de 15 minutes !!!
Premiers tests 3D en cours.
- A modifié le instal1d.sh pour pouvoir limité le nombre de cas à installer.
- A continué à nettoyer le 1D. Plus aucun fichier fortran ne doit être modifié.
Les routines et .h dans phy1d sont complémentaires de ceux de phylmd.
Un travail reste à faire si on veut éliminer les liens logiques vers dyn3d et dyn3d_common. Une partie des routines (comme pres2lev) pourraient
être dans bibio ou grid. Le infotrac devrait peut-êre être dans la physique, au moins en partie.
Base de réflexion/discussion pour la convergence avec les physiques planétaires.
- A ajouté des routines d'IOs.
LMDZ6
=====
On peut tourner par an.
Se mettre sur une configuration stochastique avec 59 niveaux et première couche à 10 m. Et sans doute avec le thermcell_plume d'Arnaud qui permet d'avoir les strato-cumulus.
Basculer au maximum sur dyn3dmem.
Faire un profiling du modele sur ada.
Faut-il rerégler les paramètres de l'orographie
Qui prend en charge XIOS ?
Rayonnement :
commettre les sources rapidement.
Réactiver les tests thermodynamique de la glace / déclenchement stochastique
Regarder le changement de densité des poches.
Jean-Yves fait remarquer qu'on pourrait mettre un klev_max.
Avoir 1+1=2 avec la nouvelle physique.
Suivi des actions en cours !
===========================!
Réduction et nettoyage des sorties
----------------------------------
2014/01/06 Réintroduction des sorties traceurs originalement dans histrac dans les
sorties standard.
2013/09/02 Nettoyage terminé à la fin du stage d'Ulysse Gérard.
2013/09/02 Abderrahmane a étendu à d'autres fichiers (hist*NMC* ...)
Interface physique/dynamique
----------------------------
2014/01/06 : Relancer une réunion.
2013/09/02 : compilation LMDZ/WRF reussie. Tests en cours
Réunion le mercredi 10 sur l'interface physique/dynamique.
Compilation de mac
------------------
2014/01/06 : . testé à l'occasion de la formation.
. problème avec awk, sed et autres.
. netcdf doit être installé sur la machine.
2012/../.. Installation réussie chez Aboudlaye Deme.
Développement du modèle.
-----------------------
2014/01/06 Lionel a passé physiq.F en physiq.F90
2012/../.. Réécrire physiq.F
Passer phylmd de *.F en *.F90
dyn3dem
-------
2014/01/06 : 1) l'opérateur non star n'est pas codé.
2) le géopotentiel n'est pas le même dans dyn3dpar et dyn3dmem.
Du coup pas de convergence avec phydev dont le rappel newtonien utilise
le geopotentiel.
3) ce0l ne fonctionne pas dans dyn3dmem.
2013/09/09 : tests encourageants.
2013/09/02 : importantes mises à jour par Yann. Tests en cours
2013/02/25 : Valider
Orchidee
--------
2013/07/* Nouvelle version disponible dans install.sh
Passer les fichiers de LMDZ sous licence Cecill
-----------------------------------------------
2014/01/06 : fait (Lionel)
ATTENTION : quand on ajoute un nouveau fichier, il faut faire une opération spécifique.
Conservation avec la nouvelle physique :
----------------------------------------
2013/04/13 : nouvelle routine pour la conservation de l'énergie
basée sur les termes source de l'équation pronostique de la TKE.
2012/03/05 : travail en cours
A regarder.
Lien avec les problèmes identifiés en 1D.
Problème identifié dans les sorties, en 1D.
1+1=2 :
-------
2014/01/06 : le problème est un peu différent. C'est day_step qui n'est pas un multiple
de idissip. Que fait-on ?
2013/09/02 : probleme identifié pour la dissipation : la période d'appel à la dissipation
effective idissip n'était pas automatiquement un diviseur de day_step
Résoudre la partie poches et convection de la nouvelle phys.
Ehouarn reprend la partie dissipation.
OK avec ancienne physique à partir de svn1626
ce0l
----
2014/02/14 : faire que ce0l marche en mode parallèle mais en séquentiel
Faire que grilles_gcm.nc sorte automatiquement.
2014/01/06 : problème de mémoire pour les grosses configs
2013/02/25 Tests en cours pour changer l'interpolation verticale
de splin vers linéaire en log-pression.
Nouveau rayonnement
-------------------
2013/10/07 : code prêt. A commettre de ce pas.
Inclusion des isotopes
----------------------
2012/../..
Identification d'une façon jolie d'inclure les isotopes dans la
dynamique.
Codé chez Camille.