10/02/2014
Compte Rendu du Point Hebdomadaire LMDZ du 10/02/2014
Durée : 1h15
Jean-Yves
---------
S'est apperçu que la précision interne de calcul de ferret est de 32 bits. Et que ça ne suffit sans doute pas pour regarder la conservation.
Se bat aussi avec Grads
Marie-Pierre
------------
A terminé de refaire tourner RRTM sur la version du 24 janvier.
A trouvé une bug dans le cas Sandu.
A fait un instal1d de la version du 5 février.
Catherine
---------
C'est finalemement plutôt la convection profonde qui ne conserve pas l'eau.
Les TP de M2 se sont bien passées.
Abderrahmane
------------
Retour sur la réunion XIOS.
Pour l'instant XIOS marche pour les 5 fichiers histoire mais pas pour les fichiers stations ni pour les fichiers de sorties.
Essaie de faire des tests sur ada.
Laurent va mettre à disposition une version XIOS.
Marie-Alice
----------- - le CR de la réunion XIOS est visible là: http://wiki.ipsl.jussieu.fr/Pole/ESCI/GT/IPSLCM6#A06.2BAC8-02.2BAC8-2014_:_r.2BAOk-union_XIOS_dans_LMDZ_et_ORCHIDEE
dans la page qui contient les CR concernant la mise au point de IPSLCM6.
- Autres points sur XIOS
avancées sur IPSLCM6-RC0 : Nemo3.6 OK (a trouvé un bug. Retrouve les
mêmes résultats qu'avant) ; OASIS MCT (qui n'a pas besoin d'un
executable. Arnaud va demander à faire des commissions sur LMDZ ; voir
avec Laurent) ; LMDZ (pour le moment dyn3dpar. Veut passer sur
dyn3dmem).
- curie : c'est le souk. Passage à netcdf4. Changement de résultats
numériques. embêtés par le freinage sur les inodes. Jean-Louis va
faire remonter des demandes groupées pour mettre la pression.
- Réunion Genci/CEA/CNRS pour Cmip6.
- Sur le site web de l'Insu, prospective AA : parle de besoin calcul et le
multi-disciplinaire.
Lluis
-----
Première simulation WRF/LMDZ qui tourne avec les conditions réalistes.
Essaie de tourner une aquaplanète avec la dernière version.
Va essayer de travailler sur Mesu, la nouvelle machine de l'UPMC.
Demain : réunion remember.
Frédéric
--------
- Retour sur la réunion physique/dynamique
WRF/LMDZ tourne. LMDZ/MAR compile et commence à tourner.
Motivation forte pour continuer le nettoyage du code physique.
Faire en sorte que la physique tourne de façon vraiment autonome.
La physique est sensée pouvoir tourner sans dimensions.h
De facon connexe, faire en sorte que le 1D réintègre le répertoire
phylmd.
- Motivé à la fois par cette réunion et par la reprise en main d'un TD
basé sur une physique 1D vidée de son contenu dans laquelle les
étudiants codent, travail sur le 1D.
- Modification du 1D. Changement de l'ordre real/entier dans les common
pour éviter des flopper de warning. Elimination de declarations inutiles
pour les memes raisons.
Le code 1D mérite pas mal de nettoyage.
Les common en dur dans plusieurs routines sont dangereux. Il faudrait
les sortir dans des .h. Peut-etre regrouper tous les commons dans un
seul .h.
Les routines qui ne sont pas redondantes avec celles de dyn3d peuvent
être dans un .F plutôt qu'un .h.
- Elimination de ocean_force_mod du 1D en mettant la routine de lecture
1D dans ocean_force_mod, appelé si klon_glo=1 ce qui était déjà le cas
dans phy1d.
- modification de create_make_gcm pour prendre en main le couplé.
- On a discuté à la réunion physique/dynamique de l'opportunité de passer
toutes les variables de la physique dans des modules avec des allocates.
C'est coûteux.
LMDZ6
=====
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/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.