Changes between Version 1 and Version 2 of InfosOzone

Timestamp:

06/25/10 18:40:51 (14 years ago)

Author:

cugnet

Comment:

--

InfosOzone

@@ -1 +1 @@
 = Information sur les fichiers de forçage Ozone =
+
+Ce sont des fichiers annuels de 12 champs (un par mois). Ils contiennent quatre champs:
+  tro3:          ozone de jour + nuit
+  tro3_daylight: ozone de jour
+  Ptrop:         pression à la tropopause (utilisée pour l'hybridation)
+  night:         0: nuit ; 1: jour (s'applique à l'ozone de jour)
+
+STRATOSPHERE:
+  Sont utilisées les sorties du run REF-B1 du projet CCMVal, produites par LMDZ-REPROBUS, résolution 96x72x50.
+  Emplacement à l'IDRIS: '''/u/rech/hou/rhou355/LMDZ4/output/h2o-REF1-KE-8mpi-1omp/'''
+
+TROPOSHPERE:
+  Sont utilisées les sorties du run RCP prooduites par LMDZ-INCA, résolution 96x95x19.
+  Emplacement au CCRT:  '''/dmnfs11/cont003/p24lava/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_IPCC_3/CHM/Output/MO'''
+
+
+
+== PRINCIPE RÉSUMÉ: ==
+
+  * détection des tropopauses chimique (premier point renconté en partant du sol dont la concentration en ozone vaut 0.15 ppmV) dans les champs stratosphériques et troposhpériques.
+  * interpolation horizontale des deux champs à la résolution maximale (ici: 96x95)
+  * Le champ hybride est composé de trois domaines:
+    1) au-dessus  de la tropopause + 2 kms: champ troposphérique
+    2) au-dessous de la tropopause - 1 km : champ stratosphérique
+    3) au voisinage de la tropopause (-1 -> +2 kms): mélange des deux champs précédents, 
+on conserve le champ strato, et en troposphère le champ tropo dilaté (on lui imposer la tropopause stratosphérique), mélange des deux champs autour de la tropopause (entre 2kms au-dessus et 1 km en-dessous) avec un profil de pondération en sinus.
+  * Climatologie par moyenne glissante sur 11 ans.
+
+
+== PRINCIPE DÉTAILLÉ: ==
+
+1) Moyennes zonales d'ozone stratosphérique de jour+nuit à partir de sorties journalières moyennées (fréquence: 1 jour) -> précis
+   * interpolation des champs 3D sur les niveaux standard (presnivs)
+   * moyenne zonale
+2) Moyennes zonales d'ozone stratosphérique de jour+nuit et jour à partir de sorties instantanées (fréquence: 10 jours) -> très imprécis
+   * utilisation d'un masque jour/nuit (ici la concentration en NO sur le 45ième niveau standard du modèle: 1.527329 hPa, indique le jour lorsqu'elle est inférieure à 1.0E-15 vmr).
+   * même méthode qu'en étape 1 pour obtenir les deux champs
+
+Interlude: nécessité d'un lissage.
+
+Les raisons:
+  * L'échantillonnage temporel des champs instantanés est très faible: pour une climatologie mensuelle, trois fichiers sont utilisés (en réalité, 4, dont deux avec une semi-pondération), alors qu'une moyenne basée sur des fichiers mensuels fait intervenir 1440 champs (un parpas de temps).
+  * L'échantillonage spatial des champs instantanés est plus faible: seule la moitié des champs est utilisés (le jour ou la nuit)
+  * L'échantillonage spatial des champs instantanés est inhomogène: les champs sont tous pris à 0h UTC. Si une structure zonalement dissymétrique présente une stabilité géographique (particulièrement un vortex polaire non zonal), elle biaise la climatologie de jour.
+
+Les remèdes:
+  * Aucun cycle diurne n'est attendu, hormis près de la stratopause.
+=> Il est légitime de moduler la climatologie d'ozone obtenue en point 1) avec un coefficient correctif basé sur les champs produits en 2), restreint près de la stratopause.
+  * Ce cycle diurne n'a pas de raison physique valable de avrier significativement dans la période disponible (1960-2006).
+=> Il est légitime d'utiliser une climatologie sur toute la période disponible des champs produits en 2) pour construire le coefficient correctif.
+
+3) Ajustement des champs précédents.
+  * obtention d'une climatologie sur toute la période disponible (1961-2006)
+  * construction du rapport rO3nd=(O3-O3d)/O3  (d signifie 'day' ; l'absence d'indice sous-entend 'day+night')
+  * lissage de rO3nd en partant du sommet du modèle (niveau 50):
+    - près des pôles  (~[-90°,-72°] et [+72°,+90°]): rapport mis à 0 dès qu'il s'annule ou que deux valeurs consécutives sont inférieures à 0.2
+    - entre les pôles (~[-72° 72°]): rapport mis à zéro dès qu'il se remet à augmenter.
+  Le coefficient correctif Corr utilisé pour construire l'ozone de jour à partir de l'ozone jour+nuit vaut 1. si rO3nd=0., O3d/O3 sinon.
+  * génération de fichiers: tro3 vaut O3(calculé en 1) ; tro3_daylight vaut Corr*O3(calculé en 1)
+
+Interlude: l'interpolation temporelle.
+  Normalement, les valeurs de O3d situées dans la nuit polaire ne sont pas utilisées, justement car il n'y fait pas jour...
+  Toutefois, certains cas limite obligent à s'en préoccuper. Considérons un point P du globe près du terminateur ; nous voulons interpoler l'ozone de jour en P à la date ti. Soient t1 et t2 les dates adjacentes du fichier d'ozone. En t1 et ti, P est dans la nuit polaire ; en t2, il est dans le jour. Il faut une valeur d'ozone de jour en P afin de pouvoir interpoler la concentration en ti. Il est alors préférable de prendre pour cela l'ozone de jour+nuit, puisqu'adopter l'ozone de jour conduirait à surévaluer le résultat. C'est pour cela que les fichiers d'ozone de jour ont été complétés par l'ozone de 
+, tel qu'au 15 du par exemple, l'interpolation temporelle en un point qui est le mois précédent dans la nuit polaire, le mois
+
+Les raisons poussant à réaliser un lissage: 
+  * l'échantillonnage temporel est assez faible ; les champs sont tous pris à 0h UTC. Si une structure présentant une dissymétrie zonale (notamment les vortex polaires) subsiste près de la transition jour-nuit, elle nuit grandement à la qualité des moyennes zonale jour et jour+nuit. 
+
+  L'hybridation des fichiers passe apr une détection de la tropopause chimique (concentration en ozone
+
+Les fichiers s'obtiennent en 
+**** ANCIENNE METHODE ****
+
+
 
 Les fichiers d'ozone sont produits.