Post traitement de INCA
Atlas de comparaison
(Anne Cozic 05/11/2007) en cours de developpement
Pour plus de renseignements sur la visualisation des runs et les catalogues d'atlas disponibles à l'IPSL voir ici
Les scripts suivant sont utilisables avec des fichiers de sorties moyennés mensuellements.
Pour créer un tel fichier à partir d'une sortie journalière vous pouvez utiliser la commande nco suivante :
ncwa -a time_counter file.nc MM_file.nc
sur les machines du LSCE
Vous devez dans un premier temps mettre à jour votre fichier .cshrc en ajoutant les lignes suivantes :
source /home/users/brock/.atlas_env_asterix_csh setenv FER_ATLAS_INCA /home/users/brock/ATLAS_INCA setenv FER_GO "$FER_GO $FER_ATLAS_INCA" setenv FER_DATA "$FER_DATA $FER_ATLAS_INCA" setenv FER_PALETTE "$FER_PALETTE $FER_ATLAS_INCA" setenv FER_EXTERNAL_FUNCTIONS "$FER_EXTERNAL_FUNCTIONS $FER_ATLAS_INCA" set path=($FER_ATLAS_INCA $path)
Avant de créer un atlas, vous devez rappatrier sur asterix vos fichiers de moyennes mensuelles (créés avec la commande indiquée plus haut)
La commande pour créer un catalogue d'atlas est alors la suivante :
atlas_inca -o output_directory -t "title" file1.nc file2.nc (avec file1.nc et file2.nc les deux fichiers que vous voulez comparer) Pour connaître toutes les options : atlas_inca -h
ATTENTION : vous devez indiquer l'adresse complète des fichiers file1 et file2
Le catalogue créé répartie les images en plusieurs catégories :
emissions - dépots - wet - wet3d - dry - load - od550 - od550_3d - od443 - od443_3d - od 670 .... - odivisir - asyvisir -asyuvvis
sur les machines du CCRT
Post traitement idl et python
Mail de Sophie Szopa suite à une réunion sur le post-traitement des sorties de INCA (06/02/2008)
Suite à la réunion sur le post-traitement et la visualisation des sorties INCA: Pensez à vous inscrire sur la liste idl via https://sympa.lsce.ipsl.fr/wws/info/idl si vous etes interesses. Vous trouverez les outils présentés aux adresses suivantes: - une premiere version distribuable du package AEROCOM presenté par Jan pour les champs 2D sera archivée la semaine prochaine (je renverrai un mail). Cet ensemble de routines idl (cf presentation de Jan jointe) permet la lecture de sorties de modèle+ mesures et la préparation de différents graphiques (série temporelle, cartes avec bullet pour le obs, histogrammes, scatter plot). - les outils PYTHON d'extraction aux stations développés par Francois se trouvent à ~delage/lib/python ou sur mercure à /home/cont003/p24fdela/EXTRACT_LMDZ_DATA/ Vous pouvez consulter sa présentation à ~delage/pres_final.pdf - l'atlas ferret (Anne Cozic et Patrick Brockmann) pour une previsualisation simple de comparaisons entre deux runs INCA se trouve pour le moment sur asterix, vous trouverez toutes les infos afferentes à cet atlas mais plus généralement à INCA3 sur la page wiki http://forge.ipsl.jussieu.fr/inca Suite à la discussion sur LMDz, j'ai placé (temporairement) le document préparé par l'équipe invsat dans le cadre de GEMS qui évalue les performances de LMDz pour le transport du CO2 à sur asterix à ~laval/LMDZ_GEMS.pdf. Merci à ceux qui ont participé à cette petite réunion. Sophie
