Création des fichiers de guidage et adaptation de l'exécution : Différence entre versions

Version actuelle en date du 19 mai 2022 à 14:41

Créer les fichiers de guidage

1 => Créer le fichier grilles_gcm.nc avec ce0l. Pour cela, avant de lancer ce0l, ajouter dans run.def :

grilles_gcm_netcdf=y

(Note : avec libIGCM, ce fichier est créé automatiquement dans $STORE/IGCM_OUT/LMDZ/EL*-$resol/ATM/Output/Grid , lorsqu'on fait tourner l'expérience CREATE pour créer les fichiers initiaux et limites pour le grid choisi - régulier ou zoomé, cf gcm.def)

2 => Extraire les scripts pour interpolation de fichiers de guidage :

wget http://forge.ipsl.jussieu.fr/igcmg/browser/TOOLS/INTERP_NUDGE
ou (avec ID et mdp) :
svn co http://forge.ipsl.jussieu.fr/igcmg/svn/TOOLS/INTERP_NUDGE

On utilisera interp_from_era.ksh pour fichiers ERA*, ou interp_from_TS.ksh pour séries de temps sorties du modèle. Les choix de l'utilisateur se font dans ces scripts, qui lancent automatiquement era2gcm.ksh pour l'interpolation.

(un exemple d'utilisation de CDO pour interpoler les champs se trouve [ici])

3 => Vérifier l'emplacement, la disponibilité et le droit d'accès aux fichiers non-interpolés (ERA* ou TS). Télécharger les fichiers si besoin.

4 => Interpoler avec interp_from_*.ksh :

 4.a => Adapter les variables en tête du script, dans la section "USER DEFINITION":
gridfile=chemin pour grilles_gcm.nc
indir= ou se trouver les fichier era-i
varlist="u v" etc : les variables a interpoler
outdir= répertoire de resultats
first_year=2011
last_year=2011
rundir= répertoire temporaire
(for ERA*) : OutFreq='monthly' ou 'yearly' : fichiers de sortie par mois ou par année

 4.b => Lancer le script ; par ex :
  ./interp_from_era.ksh

Adapter le run

5 => Créer/adapter le fichier guide.def, par exemple :

## debut guide.def
# ok_guide=y : guidage active
ok_guide=y
(Note : avec libIGCM, mettre ok_guide=y in lmdz.card ; c'est lmdz.driver qui va changer automatiquement guide.def )

guide_u= y
guide_v= y
guide_T= n
guide_P= n
guide_Q= n

tau_min_u=0.0208333
tau_max_u=0.125
tau_min_v=0.0208333
tau_max_v=0.125
## fin guide.def

Remarque :
Ici, dans ce guide.def, uniquement le guidage des variables u et v sont activés.
Temps de relaxation, pour une grille régulière, la valeur tau_max est utilisé, ici tau_max=0.125 => 3H.
Temps de relaxation, pour une grille zoomée, la valeur tau_max est utilisé à l'intérieur de la zone du zoom et tau_min à l'extérieur de la zone du zoom, ici tau_min=0.0208333 => 30min.
Plus le temps de relaxation(tau) est petit, plus le guidage est fort.

6 => Ajouter en tête de run.def la prise en compte du fichier guide.def :

INCLUDEDEF=guide.def

7 => Copier (ou créer des liens vers) les fichier de guidage interpolés dans le répertoire où le modèle sera lancé, avec des noms génériques : u.nc, v.nc, T.nc, hur.nc ; Exemple :

ln -s ERAI_interp/2011/u_201101_erai.nc u.nc
ln -s ERAI_interp/2011/v_201101_erai.nc v.nc

(Note : avec libIGCM, les fichiers doivent être indiqués dans lmdz.card dans la section BoundaryFiles )

8 => Lancer l’exécution ; exemple :

./gcm.e > out_guide 2>&1

2022 mai 16