******************************************************************
*******************************************************************
VERSION 4.5 DE LA PROGRAMMATHEQUE DES PARAMETRAGES
LE 8 JUIN 2007
*******************************************************************
*******************************************************************
-------------------------------------------------------------------
DESCRIPTION GENERALE DE LA VERSION 4.5
-------------------------------------------------------------------
Avertissements importants :
-------------------------
1. La refonte de l'interface dynamique-physique a entraîné des
changements importants dans le namelist de la physique. Il est
fortement recommandé de consulter le fichier suivant :
http://iweb.cmc.ec.gc.ca/~armnphy/v_4.5/phy_namelist
2. Les usagers du modèle GEM qui désirent DÉSACTIVER la physique
doivent s'assurer que la clé "phy_pck_version", qui définit
la version de la physique à utiliser, prenne la valeur suivante
dans le namelist physics_cfgs :
&physics_cfgs
phy_pck_version = 'NIL'
/
On peut aussi obtenir le meme resultat en eliminant le namelist physics_cfgs.
3. Les options correspondant aux clés P_pbl_monta_L et P_pbl_corts_L ne sont plus
disponibles. Pour ceux qui avaient l'habitude d'utiliser
P_pbl_monta_L = .false.
avec la physique 4.4 (et les versions antérieures) devront avoir
cortm = .false.
dans le namelist physics_cfgs pour obtenir l'effet équivalent.
4. Un nouvel outil de vérification de la syntaxe des namelists est disponible.
Pour faciliter votre travail avec la nouvelle version, vous avez la
possibilite de vous servir du script "check_nml.sh".
a) Placez vous dans le repertoire ou vous avez ouvert l'experience
b) Assurez-vous que le fichier gem_settings.nml se trouve dans le repertoire
c) Lancez le script
~armnphy/bin/check_nml.sh
d) Suivez les instructions qui apparaissent a l'écran.
A) Nouvelles options :
-----------------
1. Le rayonnement cccmarad, auteurs Li et Baker. Le schéma s´initialise
dans le namelist avec la clef
radia = 'CCCMARAD'
Le sous-programme principal "cccmarad" est appelé dans phy_exe.
2. Le paramétrage des ondes de gravité non orographiques, auteur Hines;
ce code a été importé de l´Institut Max-Planck, en Allemagne. Le schéma
s´active avec la clé (par défaut .false.)
NON_ORO = .true.
Pour spécifier si on calcule la tendance de température générée par les
ondes de gravité non orographiques, on utilise l´option (par défaut 0)
IHEATCAL = 1
et on peut initialiser l´écart-type de vitesse du vent avec
(par défaut 1.0)
RMSCON = la_valeur
Le sous-programme principal "gwspectrum" est appelé dans gwd5.
3. Le schéma de condensation à moments d'ordre multiple de Milbrandt-Yau.
Il y a deux nouvelles valeurs possibles pour la clé STCOND :
MULTIMOM_SM moment de premier ordre
MULTIMOM_DM moment d'ordre deux, alpha fixé
Dans une prochaine version de la programmathèque, deux nouvelles options
seront rendues disponibles :
MULTIMOM_DM2 moment d'ordre deux, alpha diagnostiqué
MULTIMOM_TM moment d'ordre trois
Les sous-programmes "multimoment_**" sont appelés dans vkuokon6.
4. Une nouvelle source d´humiditité due à l´oxidation du méthane
dans la stratosphère a été conçue pour les experiences climatiques
et stratosphériques. Le schéma (par défaut .false.) s´active avec
LMETOX=.true.
Le sous-programme principal "metox" est appelé dans phy_exe.
B) Autres modifications importantes :
--------------------------------
1. Nouveau traitement des tendances
L'ajout de nouveaux paramétrages nous a incité à simplifier le
traitement des tendances physiques. Ces dernières sont accumulées
dans la physique et seule la tendance totale est passée à la dynamique
pour chacune des variables prédictives du modèle. Les nouveaux noms
suivants ont été créés pour ces tendances :
uphytd vitesse du vent, composante U
vphytd vitesse du vent, composante V
tphytd température virtuelle
huphytd humidité spécifique
qcphytd rapport de mélange de l'eau nuageuse
qiphytd rapport de mélange de la glace ou de la neige
qgphytd rapport de mélange de la neige roulée ("graupel")
qrphytd rapport de mélange de la pluie
Ce traitement des tendances ne permet pas le lissage des tendances
convectives et celui-ci a été eliminé dans la version 3.3.0 de la
dynamique.
2. Révision de l'interface dynamique-physique :
L'interface entre la dynamique et la physique a été refondue pour
faciliter le branchement d'autres ensembles de paramétrages physiques
et chimiques. Les sous-programmes d'interface ont été renommmés ainsi :
phyexe1 --> phy_exe
phydebu4 --> phy_debu
phyopt --> phy_opt
getbus1 --> phy_getbus
zeracc --> phy_zeracc
De plus, les sous-programmes phy_init (pour la lecture des namelists)
et check_options (pour vérifier la cohérence des options choisies)
ont été créés.
L'impact le plus important sur les usagers est le changement au namelist
"physics". Les noms de variables du namelist correspondent désormais
aux noms qui étaient déja utilisés dans la physique (voir le common
options.cdk, le sous-programme phy_opt.ftn et le fichier de documentation
http://iweb.cmc.ec.gc.ca/~armnphy/v_4.5/phy_namelist).
De plus, le namelist physics_2 a été créé pour définir les options
d'un deuxième jeu de paramétrages utilisé sur un sous-domaine. Notez
que le parametre npeOMP du fichier configexp.dot.cfg, qui definit
le nombre de fils d'execution ("threads") OpenMP, doit alors etre
egal a 1.
Enfin, le namelist "physics_simplified" permet de définir les options
de la physique linéarisée.
3. Progiciel des séries temporelles, version 2 :
Voici les nouveautés :
* noms de variables de 4 caracteres
* étiquette de 12 caractères
* heure définie en real*8 : utile pour les très longues intégrations
* prise en compte des accumulateurs
* prise en compte de la coordonnée verticale hybride de GEM
La compatibilité arrière avec la version 1 est assurée.
4. Mosaïque pour CLASS :
Le schéma de surface CLASS peut dorénavant etre utilisé en mode
"mosaïque", c'est-à-dire que chaque pavé du modèle peut être
subdivisé en sous-pavés ayant des caractéristiques spécifiques.
Comment définit-on une variable mosaïque dans phy_ini? Par exemple,
la température du sol, qui a 2 niveaux (température superficielle
et température profonde), pourrait etre définie sur 3 niveaux de
mosaïque en modifiant l'appel à gesdict de la façon suivante :
VS=row*2
devra être remplacé par
VS=row*2@3
Ainsi, chacun des 3 sous-pavés sera caractérisé par une température
superficielle et une température profonde distincte des autres sous-
pavés.
Pour utiliser CLASS en mode mosaïque, des modifications doivent être
apportées au code de CLASS ainsi qu'à l'entrée et à la sortie du modèle
GEM (voir Vincent Fortin pour informations supplémentaires).