Post-traitements des simulations MesoNH: les outils annexes ou `` tools’’

1. Post-traitements des simulations MesoNH:les outils annexesou ``tools’’ Isabelle Mallet, Nicole Asencio, Gaëlle Tanguy CNRM/GMME

2. Plan • Les outils autour du format LFI dans la grille • du modèle: • Programme DIAG • Fichiers synchrones et diachroniques • Sorties graphiques : DIAPROG • Sorties ASCII après interpolation sur des iso-surfaces 2D: • fichier FICVAL dans DIAPROG • Les outils annexes (tools): • Traitement des observations • Vers des formats externes Netcdf Grib ASCII • Vers des grilles « observation » : variables colocalisées, • vent méridien et vent vent zonal • Traitement de sous-domaines/liste de variables • Post-traitements personnalisés phasés avec la Masdev

3. conv2dia Manipulation de fichiers LFI supercalculateur Synchrone(s) (x,y,z) .nnn Simulation MesoNH DIAG: + var dérivées CONVLFI Diachronique .000 (x,y,z,t,mask,p) Write ASCII Source perso. VIS5D GRIB ASCII-radar lfi2cdf lfi2v5d lfi2grb PC (x,y,z,1,1,1) diaprog Lecture+extraction 2D/1D Sortie ASCII-FICVAL Graphique gmeta Netcdf Vis5d Grib

4. avantages inconvénients Grille modèle et post-traitements Fichier FM (LFI) Vent : U,V grille Mesonh Uzonal,Vmeridien grille C • Les grilles décalées sont traitées • DIAG, CONVLFI, lfi2grb,lfi2v5d • Les grilles décalées et la projection • conforme sont traitées de façon • transparente par diaprog lors de la • constitution des graphiques. • Les sorties ASCII restent dans la grille modèle • (les grilles décalées et la projection conforme): • ASCII perso, ASCII-FICVAL, lfi2cdf • Les post-traitements de ces sorties ASCII doivent • tenir compte de ces spécificités • +les modifications de source personnelles sont • à phaser avec les versions de MesoNH

5. Avantages des outils annexes • Echanger des données (simulations, observations) pour les intercomparer. • == > interpolation des points de grille MesoNH • (projection conforme, grilles décalées, niveaux K) • vers une grille lat/lon (régulière ou non), • et une grille verticale en niveaux Z ou P • == > formats supplémentaires: NetCDF, ASCII lat-lon[-altitude]-valeur • Utiliser les ressources informatiques des postes de travail (PC) et réserver le • supercalculateur aux simulations et aux diagnostics « automatiques » • == > lecture /écriture d’un fichier FM diachronique en traitant • l’ensemble de la matrice 3D (x,y,z) + temps + processus • Limiter le nombre de logiciels graphiques à utiliser dans la communauté Mesonh • == > utilisation de l’outil standard diaprog pour traiter des observations • et les superposer aux données simulées • Un seul format d entrée aux outils: le format diachronique • == > utiliser CONV2DIA pour les fichiers synchrones

6. Les outils annexes (1): obs <->modèle MesoNH et son format FM diachronique Niveaux Z ou P ou K mesonh2obs 3D 2D 1 profil obs2mesonh 1 valeur OBSERVATIONS lat,lon, alt,valeur

7. Fichier FM AR11_mc2.06d.Z mesonh2obs: exemple obscoordlatlon 45 8 0 Fichier FM AR10_mc2.06d.Z Format du fichier de sortie ? llhv Nombre de niveaux verticaux pour l’interpolation ? 1 Nom du fichier de localisation géographique des obs.? obscoordlatlon Print de contrôle ? 1 Nom du fichier diachro en supprimant .lfi ? AR10_mc2.06d.Z Nom du champ MesoNH en majuscules? T2M END

8. obs2mesonh: ex. de comparaison à un champ 2D (T2M) G.Jaubert (CNRM/GMME) Directives diaprog: NIMNMX=1 XISOMIN=272 XISOMAX=284 XDIAINT=2 _file1_’fichiersim’ _file2_’fichierobs’ T2M(+273)_file1_ON_ LSPOT=T T2M11h_file2_ Directives obs2mesonh: fichiersim 3 ll fichierobs T2M11h K 2D END

9. obs2mesonh: ex. de comparaison à un champ 3D (réflectivité radar) Coupe verticale : modèle + obs radar Coupe horizontale à K=20 obs radar dBz Milan Ouest Est

10. Les outils annexes (2): formats externes &sous-ensemble (variables,domaines) MesoNH et son format FM diachronique extraction d’un sous-domaine x,y,z,t,mask,p extraction d’une liste de variables vent: + calcul de dd(0-360), ff(module en m/s) avec/sans changement de grille horizontale (lat/lon), verticale (Z,P) + changement de format extractdia + Rhodref + Orographie, Uzonal,Vmerid nouveau fichier FM diachronique OBSERVATIONS lat,lon, alt,valeur Format libre à préciser NetCDF DIAPROG

11. Les outils annexes (3): 2 solutions pour modifier la référence temporelle: variable TIME MesoNH et son format FM synchrone MesoNH et son format FM synchrone MesoNH et son format FM synchrone MesoNH et son format FM synchrone MesoNH et son format FM diachronique Directives pour conv2dia: 4 C102nophop1_20060725.00.Z C102nophop1_20060725.06.Z C102nophop1_20060725.12.Z C102nophop1_20060726.00.Z C102_20060725.4ech n n 2 e n OU Utile pour traiter en evolution temporelle: • plusieurs fichiers de couplage • des séries diachroniques issues de segments de run différents conv2dia concat_time_diafile Options identiques à celles d’extractdia + TIME_orig= en secondes 1 seul fichier FM diachronique: Plusieurs écheances,toutes les variables, tout le domaine Voir les sorties possibles d’extractdia DIAPROG

12. Les outils annexes (4): post-traitements personnalisés en Fortran MesoNH et son format FM diachronique voir $MESONH/MAKE/tools/diachro/src/EXTRACTDIA/exrwdia.f90 Utilisateur averti sur les grilles MesoNH! programme FORTRAN utilisateur (exemple exrwdia.f90 à adapter) + makefile gmaketools: lecture de fichier(s) FM donc grille modèle traitements personnels écriture au format … + Rhodref + Orographie, Uzonal,Vmerid nouveau fichier FM diachronique OBSERVATIONS lat,lon, alt,valeur Format libre à préciser NetCDF DIAPROG

13. Les outils annexes (5): compute_r00_pc sur PC ou DIAG sur supercalculateur MesoNH et son format FM diachronique MesoNH et son format FM diachronique MesoNH et son format FM diachronique MesoNH et son format FM diachronique MesoNH et son format FM diachronique Fichiers consécutifs reconstitution des trajectoires lagrangiennes sur l’ensemble de la période et stockage de paramètres le long de ces trajectoires: cf documentation sur les traceurs lagrangiens compute_r00_pc nouveau fichier FM diachronique DIAPROG notamment les lachers de particules

14. Les outils annexes (6): radarascii2llv ASCII-radar issu de DIAG radarascii2llv nouveau fichier ASCII : latitude longitude valeur ( nécessite obs2mesonh pour le tracer avec DIAPROG)

15. Aide à l’utilisation: documentation & liste de diffusion • Documentation • sur le site web MesoNH: • Tools related to Meso-NH (version Mar 21, 2005) • Sur les disques CNRM: • $MESONH/ doc/tools.ps • sous $MESONH/MAKE/readme • mesonh2obs.LISEZMOI • obs2mesonh .LISEZMOI • extractdia.LISEZMOI • exrwdia.LISEZMOI • compute_r00.LISEZMOI • Liste de diffusion pour partager ses expériences, ses problèmes d’utilisation: • mesonhtools@meteo.fr

16. Aide à l’utilisation: FAQ Au CNRM (version March 14 2007): file:///miroir/mesonh/doc/faq_MesonhTools/index.html sur le site Web Mesonh: FAQ on Meso-NH tools (version Nov 20 2006)

17. 1. Des outils interactifsqui répondent aux besoins exprimés par de nombreux utilisateurs: mesonh2obs , obs2mesonh, extractdia, cocat_time_diafile,conv2dia, compute_r00_pc, radarascii2llv Conclusion: partage d’outils • Ces outils ont été testés par de nombreuses utilisations • Les modules et routines utilisés suivent les phasages des versions de MesoNH Outils en constante évolution pour répondre aux besoins des utilisateurs. • 2. Des développements personnels basés sur exrwdia.f90 utilisant des fonctions mises à la disposition de tous les utilisateurs: • zmoy …………….. Moyenne sur une couche définie en K, Z, P • dd …………….. Calcul de la direction du vent 0-360 • ff …………….. Calcul du module du vent m/s • readvar ……………….Lecture d une variable dans un fichier LFI • writevar ………………Écriture d une variable dans un fichier diachronique LFI • change_a_grid ………changement de grille Mesonh • zinter ……………. ..Interpolation au niveau Altitude=val • pinter …………….. .Interpolation au niveau Pression=val • writellhv ……………..Ecriture lon lat alt val • writecdl ………………Ecriture en Netcdf • routines pour basculer dans un repère cylindrique (cyclones, …) …. intégration prévue • par DIRRE/CRC