Version 52 (modified by cugnet, 14 years ago) (diff) |
---|
Information sur les fichiers de forçage Ozone
Interpolation spatiale et temporelle des fichiers bruts AVEC "cycle saisonnier" ; dits v2
Les fichiers d'ozone hybrides bruts "avec cycle saisonnier" sont stockés sur les comptes communs là : .../IGCM/INIT/ATM/LMDZ/Ozone/HYBRIDE/v2.clim/tro3_????.new.nc
Ces fichiers sont ensuite interpolés spatialement (pour les différentes résolutions) et temporellement via la configuration CREATE :
- En input de la config CREATE : tro3_${year}.new.nc
- En ouptut de la config CREATE : climoz_LMDZ_${year}_v2.nc
Ces fichiers obtenus sont ensuite stockés sur les comptes communs :
- pour la resolution 9695 : .../IGCM/BC/ATM/IPSLCM5A/LMD9695/AR5/HISTORIQUE/climoz_LMDZ_????_v2.nc
- pour la résolution 144142 : .../IGCM/BC/ATM/IPSLCM5A/LMD144142/AR5/HISTORIQUE/climoz_LMDZ_????_v2.nc
Ces fichiers climoz_LMDZ_${year}_v2.nc sont les fichiers de forçage Ozone utilisés par défaut maintenant par le modèle.
Un nouveau jeu de fichier a été fourni par David Cugnet, ils diffèrent du jeu fourni le 24 septembre 2010 à partir de l'année 2003. Les fichiers différents ont remplacé les anciens dans:
- /dmnfs/cont003/p86ipsl/IGCM/STORAGE/INIT/ATM/Ozone/HYBRIDE/v2.clim
Les fichiers fournis le 24 septembre 2010 sont dans:
- /dmnfs/cont003/p86ipsl/IGCM/STORAGE/INIT/ATM/Ozone/HYBRIDE/v2.clim.bad
Les fichiers interpollés ont remplacé les précédents fichiers suffixés -v2
- IGCM/BC/ATM/IPSLCM5A/LMD9695/AR5/HISTORIQUE/climoz_LMDZ_????_v2.nc
- IGCM/BC/ATM/IPSLCM5A/LMD144142/AR5/HISTORIQUE/climoz_LMDZ_????_v2.nc
Interpolation spatiale et temporelle des fichiers bruts SANS "cycle saisonnier"
Les fichiers d'ozone hybrides bruts "sans cycle saisonnier" sont stockés sur les comptes communs là : .../IGCM/INIT/ATM/LMDZ/Ozone/HYBRIDE/clim/climO3_LMDZORINCAREPRO_????.nc
Ces fichiers sont ensuite interpolés spatialement (pour les différentes résolutions) et temporellement via la configuration CREATE :
- En input de la config CREATE : climO3_LMDZORINCAREPRO_${year}.nc
- En ouptut de la config CREATE : climoz_LMDZ_${year}.nc
Ces fichiers obtenus sont ensuite stockés sur les comptes communs :
- pour la resolution 9695 : .../IGCM/BC/ATM/IPSLCM5A/LMD9695/AR5/HISTORIQUE/climoz_LMDZ_????.nc
- pour la résolution 144142 : .../IGCM/BC/ATM/IPSLCM5A/LMD144142/AR5/HISTORIQUE/climoz_LMDZ_????.nc
Ces fichiers climoz_LMDZ_${year}.nc sont les fichiers de forçage Ozone utilisés pour faire les simulations ayant démarré avec un ozone sans "cycle saisonnier".
- historical[234]
- amip[12345]
Historique
22 octobre 2010
- Génération des fichiers ozone "sans cycle saisonnier" pour la période 1996-2008. Les précédents fichiers pour cette période sont renommés *_v0bis.nc
30 septembre 2010
Génération des fichiers hybrides bruts
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 des runs REF-B1, REF-B2 et SCN-B2c du projet CCMVal, produites par LMDZ-REPROBUS, résolution 96x72x50.
Emplacement à l'IDRIS (avec RUN=REF1, REF2 ou SCNB2c):
/u/rech/hou/rhou079/LMDZ4/output/h2o-${RUN}-KE-8mpi-1omp/${Y}/
Fichiers utiles:
- Pour l'ozone de jour et de jour+nuit: histday_tra_${Y}${M}.nc
- Pour l'ozone de jour uniquement: histday_ins_${Y}${M}.nc
TROPOSPHERE:
Sont utilisées les sorties des runs LOI_IPCC_3, LOI_IPCC_RCP26, LOI_IPCC_RCP45, LOI_IPCC_RCP60 et LOI_IPCC_RCP85 produites par LMDZ-INCA, résolution 96x95x19.
Emplacement au CCRT (avec RUN=3, RCP26, RCP45, RCP60 ou RCP85sstA2):
/dmnfs11/cont003/p24lava/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_IPCC_${RUN}/CHM/Output/MO
Fichiers utiles:
- Pour RUN=3 ou RCP85sstA2: clim_${Y}.nc
- Pour RUN=RCP26, RCP45 ou RCP85: LOI_IPCC_${RUN}_${Y}${M}01_${Y}${M}30_1M_inca_avgr.nc
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 troposphé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
- Climatologie par moyenne glissante sur 11 ans.
Principe détaillé
ETAPE 1: Obtention de moyennes zonales lissées temporellement.
TROPOSPHERE:
Pour chaque scenario:
1) extraction des variables O3 et pmid 2) interpolation sur les 19 niveaux standard (variable "presnivs") ; moyenne zonale. 3) pour RUN=RCP26, RCP45 ou RCP60: moyenne mensuelle, concaténation.
=> fichiers annuels à 12 champs (un par mois)
4) moyenne temporelle glissante (+/- 5 ans)
PRODUITS: moyennes zonales lissees tro3_${Y}.nc de 1855 a 2095. NB: jusqu'en 2000, on garde les fichiers composites 3_RCP85. Au-delà, on utiliser les fichiers 3_${RUN}.
STRATOSPHERE: Sorties utilisees: sur gaya, dans les repertoires /u/rech/hou/rhou079/LMDZ4/output/h2o-${RUN}-KE-8mpi-1omp/$Y/
- Partie historique (1960-2006): RUN=REF1, fichiers histday_tra_${Y}${M}.nc et histday_ins_${Y}${M}.nc
- Partie a venir (2007-2098): RUN=REF2 et RUN=SCNB2c, memes noms de fichiers
Pour chaque run (REF1, REF2 et SCNB2c):
<A> ozone de jour+nuit (tro3A): meme demarche qu'en troposphere pour les fichiers histday_tra (moyennes journalieres). <B> ozone de jour:
(1) Masquage des pixels de nuit (obtenu avec un indicateur de nuit, en l'occurrence l'espece photosensible NO) (2) Interpolation sur les niveaux standard. (3) Moyennes mensuelles: [O3(${M-1}30) + 2.O3(${M}10) + 2.O3(${M}20) + O3(${M}30)]/6 (Les fichiers histins_tra donnent un champ
les 10, 20 et 30 du mois) => moyennes zonales d'ozone de jour (tro3_daylightB) de 1961 a 2006.
(4) Operations (2) et (3) (sans masquage) => moyennes zonales d'ozone de jour+nuit (tro3B) de 1961 a 2006. (5) creation d'une climatologie de tro3_daylightB et tro3B sur toute la duree de la simulation => determination du masque M
indiquant la zone photosensible en filtrant le rapport (tro3B-tro3B_daylight)/tro3B
(6) Reconstruction de tro3_daylight(y), ou y est l'annee, qui vaut:
- tro3A(y) en-dehors du masque M
- tro3A(y)*tro3_daylightB(y)/tro3B(y) dans le masque M
<C> moyenne temporelle glissante (+/- 5 ans) pour les deux associations REF1/REF2 et REF1/SCNB2c
=> moyennes zonales lissees tro3_${Y}.nc de 1961 a 2095 pour les variables tro3 et tro3_daylight.
<D> construction de fichiers synthetiques pour tous les scenarii (meme methode pour tro3 que pour tro3_daylight).
Dans ce qui suit, '*' designe le scenario considere ; 'S' et 'R' designent respectivement les scenarii REF1_SCN-B2c et REF1_REF2.
tro3* = tro3S + (tro3R-tro3S)*(CO2*-CO2S)/(CO2R-CO2S)
Les grandeurs CO2 representent les concentrations en CO2 dans les scenarii consideres. N.B.: REF1_REF2 représente le scenario RCP6.0.
=> moyennes zonales lissees tro3_${Y}.nc de 1966 a 2093 pour tous les scenarii.
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 champs sont utilisés (en réalité, 4, dont deux avec une semi-pondération), alors qu'une moyenne basée sur des fichiers mensuels en fait intervenir 1440 (un par pas de temps).
- L'échantillonnage spatial des champs instantanés est plus faible: seule la moitié des champs est utilisée pour l'ozone de jour.
- L'échantillonnage 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 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 varier 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 nécessitent de 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 (milieux de deux mois consécutifs). En t1, P est dans la nuit polaire, alors qu'en ti et t2, il est dans une zone insolée. Il faut une valeur d'ozone de jour en t1 afin de pouvoir interpoler la concentration en ti. La routine de lecture des fichiers d'ozone complète la nuit polaire du champ d'ozone de jour par duplication vers les pôles des valeurs adjacentes les plus proches (près du terminateur). En l'occurrence, cela conduirait à surévaluer la concentration interpolée. Il est préférable de choisir dans la nuit polaire l'ozone de jour+nuit. C'est ainsi qu'ont été complétés les champs d'ozone de jour dans la nuit polaire, en association avec la variable "night".
4) Hybridation des fichiers:
- Interpolation horizontale (96x95), puis détection de la tropopause des champs stratosphériques (cette tropopause est prise comme référence dans la suite, car la résolution verticale de ces champs est supérieure).
- Interpolation verticale (50 niv), puis détection de la tropopause des champs troposphériques.
- Étirement (par interpolation) du champ troposphérique de manière à "tirer" la tropopause troposphérique vers la tropopause stratosphérique et conserver le même champ de pression au sol. Cela revient à considérer que les points troposhpériques sont non pas en Pt, mais en Pt*, avec:
( Pt ) LOG(Ptrops/PtS) Pt* = PtS ( ---- ) ^ ----------------- ("S": "Sol", "trop": "tropopause") ( PtS ) LOG(Ptropt/PtS)
- Construction des champs globaux (12 par an, fichiers annuels): a) 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 note P1 (resp. P2) la pression côté strato: tropopause + 2 kms (resp. côté tropo: tropopause - 1 km). En clair:
O3t(P)-O3s(P) PI lnP -lnP1 O3*(P) = O3s(P)+ ------------- [ 1 + SIN( -- (1 + --------- )) ] 2 2 lnP2-lnP1
DIFFICULTÉ: les périodes des deux runs diffèrent. tropo: 1850-2000 ; strato: 1961-2006.
ANCIENNE MÉTHODE:
- 1850<=y<=1960: tropo(y), strato(1961-1971)
- 1961<=y<=2000: tropo(y), strato(y)
- 2001<=y<=2006: tropo(1990-2000), strato(y)
NOUVELLE MÉTHODE:
- 1850<=y<=1960: tropo(y), strato(y-11*E[(y-1961-11+1)/11])
- 1961<=y<=2000: tropo(y), strato(y)
- 2001<=y<=2006: tropo(y-11*E[(y-2000+11-1)/11]), strato(y)
Les 11 ans correspondent à la largeur de moyennage temporel (étape suivante). Les fichiers annuels de la nouvelle méthode n'ont aucun sens physique: seules leurs climatologies sont utilisables.
5) Moyenne glissante sur 11 ans.
COMMENTAIRE SUR LES DEUX MÉTHODES: La première méthode produit des fichiers annuels utilisables, mais la climatologie tropo passe continûement de la moyenne 1961-1971 (notée SC66) à celle même moyenne, par l'intermédiaire de champs différents, de 1955 à 1966 (dans la suite, une grandeur dSYY désigne une anomalie par rapport à la moyenne SC66, ie: dSYY=SYY-SC66 ; on remarque que la somme des dSYY pour YY variant de 61 à 71 est nulle):
Sclim55 =SC66 Sclim56=(10*SC66+S61 )/11=SC66 +dS61/11 Sclim57=(9 *SC66+S61+S62 )/11=Sclim56+dS62/11 Sclim58=(8 *SC66+S61+S62+S63 )/11=Sclim57+dS63/11 ... Sclim60=(6 *SC66+S61+S62+S63+S64+S65 )/11=Sclim59+dS65/11=Sclim61-dS66/11 Sclim61=(5 *SC66+S61+S62+S63+S64+S65+S66 )/11=Sclim60+dS66/11=Sclim62-dS67/11 ... Sclim63=(3 *SC66+S61+S62+S63+S64+S65+S66+S67+S68 )/11=Sclim64-dS69/11 Sclim64=(2 *SC66+S61+S62+S63+S64+S65+S66+S67+S68+S69 )/11=Sclim65-dS70/11 Sclim65=( SC66+S61+S62+S63+S64+S65+S66+S67+S68+S69+S70 )/11=SC66 -dS71/11 Sclim66=( S61+S62+S63+S64+S65+S66+S67+S68+S69+S70+S71)/11=SC66
Cela n'est pas très satisfaisant...
En revanche, avec la seconde méthode, la climatologie strato reste toujours identique à SC66: grâce au petit stratagème "cyclique", les strato de 1961 à 1971 se trouvent toujours comptabilisées dans la moyenne.
Notes:
Contrairement à ce que l'on pourrait penser, les fichiers avant 1956, par exemple, produits par les méthodes 1 et 2 ne sont pas identiques. La stratosphère est égale à SC66 (tout au moins jusqu'à l'enveloppe supérieure des tropopauses des fichiers impliqués), par construction, mais la troposphère diffère. En effet, dans le premier cas, la tropopause est toujours celle du champ SC66, et celle de la moyenne temporelle également. Dans le second cas, les tropopauses successives sont celles des champs S61 à S71, dans un ordre particulier. Les troposphères sont donc, avant moyennage, dilatées différemment, et la moyenne globale n'est pas la même qu'avec la méthode 1.
Le fait d'imposer, d'une manière ou une autre, une stratosphère "pré-industrielle" n'a en fait pas grande réalité physique, mais n'est pas moins artificiel que de combiner des champs d'ozone issus de deux réalisations et modèles différents (et donc deux histoires météorologiques différentes). Partant du principe qu'il n'y a pas d'autre alternative (faute d'un modèle chimique stratosphère-troposphère unifié), on peut considérer comme raisonnables les hypothèses qui sous-tendent ce qui est décrit sur la présente page.
DC
Attachments (5)
- FonctionChapeau.png (17.2 KB) - added by cugnet 7 years ago.
- TropopauseLMDZ.png (136.7 KB) - added by cugnet 7 years ago.
- TropopauseLMDZ.2.png (483.6 KB) - added by cugnet 7 years ago.
- TropopauseLMDZ2.png (88.4 KB) - added by cugnet 7 years ago.
- TropopauseLMDZ.3.png (259.8 KB) - added by cugnet 7 years ago.
Download all attachments as: .zip