filibert/GRADS..M-A Filiberti Note technique IPSL/ Pôle de modélisation Draft
Décembre 1999
Ce document décrit des programmes généraux pour traiter
les sorties de LMDZ.3.3, il peuvent àtre adaptés en fonction
des besoins de chacun.
Différentes étapes sont nécessaires pour créer
un atlas. Il faut tout d'abord créer un fichier au bon format
à partir des sorties du modèle.
A partir des sorties du modèle au format netcdf, générées
par IOIPSL dans physiq.F (fichiers de type histmth.nc), pour
obtenir des moyennes mensuelles sur plusieurs années de simulation,
et l'interpolation sur 15 niveaux verticaux de certains champs
pour comparer aux observations qui sont à cette résolution
verticale, il faut utiliser des programmes d'analyse issus de
ceux crées pour les sorties au format netcdf par L. Li (LMD)
et qui se trouvent sur Rhodes (IDRIS) sous /home/rech/hou/rhou909/GRADS
et sur Etoile au CCR (Jussieu) sous filibert/GRADS..
Il s'agit de compactage (histmth.nc) et de analyse (histmth.nc).
Pour ces deux programmes il faut juste modifier l'adresse de
retour du listing décrivant l'exécution du job en début
de programme.
Le premier programme compactage_histmth_nc récupère
des fichiers sur Gaya (IDRIS) dans le directory reper1. Il s'agit
de fichiers issus de LMDZT de la forme : histmth.y.m.nc oó y et m sont respectivement les années et les mois sur 2
chiffres.
Il faut modifier reper1 pour récuperer ses propres fichiers.
Puis on utilise le programme ncea de la librairie netcdf NCO qui fait la moyenne des
différents fichiers. Dans notre cas il s'agit de la moyenne
mensuelle sur plusieurs années d'expérience. En sortie
de ce programme on obtient des fichiers histmth.nc.mm (m étant le mois sur 1 ou 2 chiffres) qui sont stockés
directement sur Gaya (éventuellement dans un directory différent :
reper2 à modifier dans le programme).
Ce programme est à utiliser si on veut une moyenne mensuelle
sur plusieurs années de simulation, sinon on peut directement
utiliser le programme suivant, le programme Grads créant
l'atlas pouvant gérer plusieurs années.
Si on veut obtenir la moyenne de champs à l'intérieur
d'un màme fichier par exemple pour passer de sorties journalières
à des sorties mensuelles on peut de màme utiliser le
programme ncra de la librairie NCO. Ces deux programmes se trouvent
sur Rhodes sous rhou909/bin et sur Etoile sous filibert/GRADS.
Le second programme analyse_histmth_nc crée à partir
des fichiers mensuels des fichiers interpolés sur 15 niveaux
et contenant un nombre réduit de champs par rapport aux fichiers
initiaux.
Il faut modifier éventuellement les dimensions de la grille :
PARAMETER(iim=96, jjm=72, llm=19, jjp1=jjm+1) en début de
programme.
Ce programme récupère sur Gaya les fichiers mensuels
(format Netcdf) (eventuellement issus du programme précedent).
Les fichiers sont de la forme : histmth.nc.mm (m étant le mois sur 1 ou 2 chiffres) et sont stockés
sur Gaya dans le directory reper (à modifier en fonction
de la localisation de ses propres fichiers). En fin de programme
on copie le fichier annuel résultat histmth_champ3d.nc au format Netcdf dans le directory
reper.
Sorties dans la partie physique de LMDZ.3.3 (dans le programme
physiq.F)
Pour les sorties mensuelles : fichier histmth.nc.
Pour chaque champ stocké dans le fichier histmth_champ3d.nc
on indique son nom dans le fichier, son nom dans le modèle,
sa signification, son unité et s'il est oui ou non conservé
dans le fichier final.
On utilise le fichier annuel contenant des champs
mensuels crée à l'étape précédente. Deux
possibilités : soit sortir les champs horizontaux à la
surface ou à un niveau de pression donnée, soit des champs
2D verticaux en moyennes zonales.
Pour cela on utilise une version modifiée du programme Grads crée par S. Bony-Lena (LMD) et adaptée au format Netcdf par M. Forichon (LMD): n atlas_lmdz.gsqui permet de comparer les sorties aux observations.
Il faut àtre sur une machine sur laquelle les fonctions Grads associées à regrid2 sont installées par exemple sur piaf au LMD.
Il faut utiliser l'option portrait de Grads donc lancer : grads -pc ``run natlas_lmdz.gs``
ou
grads -p puis
run natlas_lmdz.gs
On peut éventuellement avoir plusieurs fichiers annuels le programme gérant les années.
Pour chaque champ stocké dans le fichier on indique les données correspondantes et la source de ses données.
Les données sont disponibles pour une certaine période (mois et année ou climatologie).
Le programme natlas_lmdz.gs pose différentes questions :
'Première année de simulation LMDZ? '
'Dernière année de simulation LMDZ? '
'Directory oó se trouvent les fichiers LMDZ (lisibles par Grads)?'
'Sorties sous forme de fichiers postscript? (oui/non)'
si rep='oui' => 'Directory oó seront écrits les fichiers postscripts?'
'Projection: lat/lon (latlon) ou north/south stéréo (nps/sps)?'
'Traitement de toutes les variables ou slt une seule? (tout/une)'
Si on ne veut traiter qu'une seule année il faut entrer la màme date pour les 2 premières questions. L'année est donnée sur 2 chiffres (à revoir pour des données ou des simulations postérieures à 1999).
Les comparaisons se font pour les mois de janvier et juillet sur la période de recouvrement maximum entre les données et le modèle. On trace le champ moyen LMDZ, le champ moyen dans les données et la différence des deux.
On peut le faire pour un champ donné ou pour tous.
On utilise le programme zonal_atlas_lmdz.gs écrit par M. Forichon (LMD)
qui trace des moyennes zonales pour janvier et juillet avec comparaison aux données NCEP.
Il faut àtre sur une machine sur laquelle les fonctions Grads associées à regrid2 sont installées par exemple sur piaf au LMD.
Il faut utiliser l'option portrait de Grads donc
lancer :
grads -pc ``run zonal_atlas_lmdz.gs``
ou
grads -p
puis
run zonal_atlas_lmdz.gs
Ce programme est plus simple que le précédent et ne sonne pas la possibilité de faire une partie seulement des sorties.
Pour visualiser les fichiers IOIPSL issus d'une simulation sans passer par l'interpolation verticale décrite précedemment : pour l'instant ces fichiers ne sont pas directement utilisables par les logiciels graphiques car ils ne contiennent pas les bonnes instructions sur le temps et les lat/lon. Cette démarche est inutile dans le cas des fichiers interpolés car ils sont réécris avec les données nécessaires au graphique.
Pour transformer les fichiers de sorties *.nc pour
qu'ils puissent àtre lus avec les logiciels graphiques, il
faut compiler le programme ncregu. F qui se situe dans mafipsl/UTILS
au LMD ou dans filibert/UTILS au CCR. Le README contient
les instructions de compilation de ce programme :
f77 -I/distrib/local/include ncregu.F -L/distrib/local/lib
-l netcdf
f90 -I/usr/local/netcdf-3.3.1/include ncregu.F -L/usr/local/netcdf-3.3.1/lib
-l netcdf
a.out fichier.nc
=> fichier.nc. Le fichier garde le màme nom.
Puis lancer grads (s'il existe sur le site oó on travaille sinon transférer les fichiers modifiés): (pas de *.ctl) et faire
sdfopen fichier.nc,
puis les commandes classiques de Grads.
Sorties graphiques avec LMDZT:
Comment créer un atlas des sorties avec LMDZ.3.3 avec
le logiciel Grads.
This document was generated using the LaTeX2HTML translator Version 96.1 (Feb 5, 1996) Copyright © 1993, 1994, 1995, 1996, Nikos Drakos, Computer Based Learning Unit, University of Leeds.
The command line arguments were:
latex2html -split 0 graflmdz.
The translation was initiated by Compte generique LMDZ on Thu Feb 17 17:11:19 MET 2000