110 likes | 197 Views
MSDOL/PALM. D. Cugnet. Atelier ADOMOCA - Toulouse 21-22 novembre 2007. Géométrie et transport. Modèle stratosphérique à 36 niveaux: de 10 à 55,5 kms par pas de 1,3 kms (263,3 à 0,33 hPa) Coordonnées pressions constantes Grille horizontale: 2,5°x2,5°
E N D
MSDOL/PALM D. Cugnet Atelier ADOMOCA - Toulouse 21-22 novembre 2007
Géométrie et transport • Modèle stratosphérique à 36 niveaux: de 10 à 55,5 kms par pas de 1,3 kms (263,3 à 0,33 hPa) • Coordonnées pressions constantes • Grille horizontale: 2,5°x2,5° • Transport: semi-lagrangien, ∆t=30 min. ; trois méthodes d’interpolation, de précision et de coût croissants. Par défaut: ordre 1. Forçage: des fichiers météo ECMWF T106 (vent, T) interpolés sur 36 niveaux.
Chimie Module ROSE modifié, ∆t=30 min. ; séquence: • Angle zénithal, C.L. O3, O3 total -> photolyse • Taux de réactions fonctions de T (2 versions) • Taux de réactions pour les aérosols (optionnel) • Mise à jour des concentrations avec chimie polaire (optionnelle) • Conservation de la masse (optionnelle)
Implantation MAQUETTE ADOMOCA_V3 /MSDOL/src_sv (sources) /MSDOL_COMPILE_DIRECT/src_MSDOL (Makefile, Makefile.inc, inst_links) /MSDOL_COMPILE_ASSIMILATION/src_3DFGAT_MSDOL (Makefile, Makefile.inc, inst_links) /MSDOL_JOBS (jobMSDOL, msdolpalm.bash) /MSDOL_PALM/ (MSDOL.pil, etc…) /MSDOL_TOOLS (convsh, interp.tcl) LIBRAIRIES NECESSAIRES: * direct: hdf, NetCDF. * 3DFGAT: PALM, Hnov, lam, xrd, gribex, sparsekit
Utilisation • Compilation ; choix des options et répertoires dans jobMSDOL ; lancement (éventuellement en série avec msdolpalm.bash). • Version réécrite 100% conforme à l’original ; fonctionne correctement en 3DFGAT. • Documentation d’installation « à partir de rien » • PROBLEME: bug MSDOL (cf. tests décrits ci-après) => Routine de photolyse correcte?
Mode direct (moyennes zonales) MSDOL climatologie 1 heure 1 jour 9 jours MOCAGE
Mode direct (colonnes totales) Génération d’O3 près du terminateur -> routine de photolyse incréminée. Mais: Chaîne 3DFGAT montée et fonctionnelle. Résultats corrects dans l’hémisphère nord ; les incréments tendent à corriger l’excès d’ozone au sud.
MSDOL 3D-FGAT • Assimilation de données MIPAS: MOCAGE (4 jours) MSDOL (4 jours) MSDOL initialisé avec des C.I MOCAGE
Les enseignements • Inutile de perdre davantage de temps avec MSDOL: bug modèle • Mais: l’objectif était de savoir PALMer, acquérir les connaissances nécessaires en informatique (programmation, scripts, installateurs et compilateurs) et en modélisation de l’atmosphère (CTM). • Les (multiples) erreurs de tactiques seront évitées pour la suite.
Etapes utiles • Vérification du modèle (implantation éventuelle des sorties dûment vérifiées) • PALMérisation avec validation du modèle direct PALMé, vérification des redémarrages • 3DFGAT avec données vides • 3DFGAT avec données MIPAS ; comparaison avec d’autres modèles • Ajustement des paramètres d’assimilation (L corrélation, diffusion verticale, rapport B/R)
Plan de travail LMDZ (REPROBUS, puis INCA): • À tester • À PALMer • À Valider • Implémentation du filtre de Kalman suboptimal de MSDOL dans la chaîne de la MAQUETTE