| 90 | |
| 91 | ETAPE 1: Obtention de moyennes zonales lissées temporellement. |
| 92 | |
| 93 | TROPOSPHERE: |
| 94 | |
| 95 | Pour chaque scenario: |
| 96 | 1) extraction des variables O3 et pmid |
| 97 | 2) interpolation sur les 19 niveaux standard (variable "presnivs") ; moyenne zonale. |
| 98 | 3) pour RUN=RCP26, RCP45 ou RCP60: moyenne mensuelle, concaténation. |
| 99 | => fichiers annuels à 12 champs (un par mois) |
| 100 | 4) moyenne temporelle glissante (+/- 5 ans) |
| 101 | |
| 102 | '''PRODUITS''': moyennes zonales lissees tro3_${Y}.nc de 1855 a 2095. |
| 103 | NB: jusqu'en 2000, on garde les fichiers composites 3_RCP85. Au-delà, on utiliser les fichiers 3_${RUN}. |
| 104 | |
| 105 | STRATOSPHERE: |
| 106 | Sorties utilisees: sur gaya, dans les repertoires /u/rech/hou/rhou079/LMDZ4/output/h2o-${RUN}-KE-8mpi-1omp/$Y/ |
| 107 | * Partie historique (1960-2006): RUN=REF1, fichiers histday_tra_${Y}${M}.nc et histday_ins_${Y}${M}.nc |
| 108 | * Partie a venir (2007-2098): RUN=REF2 et RUN=SCNB2c, memes noms de fichiers |
| 109 | |
| 110 | Pour chaque run (REF1, REF2 et SCNB2c): |
| 111 | <A> ozone de jour+nuit (tro3A): meme demarche qu'en troposphere pour les fichiers histday_tra (moyennes journalieres). |
| 112 | <B> ozone de jour: |
| 113 | (1) Masquage des pixels de nuit (obtenu avec un indicateur de nuit, en l'occurrence l'espece photosensible NO) |
| 114 | (2) Interpolation sur les niveaux standard. |
| 115 | (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 |
| 116 | les 10, 20 et 30 du mois) => moyennes zonales d'ozone de jour (tro3_daylightB) de 1961 a 2006. |
| 117 | (4) Operations (2) et (3) (sans masquage) => moyennes zonales d'ozone de jour+nuit (tro3B) de 1961 a 2006. |
| 118 | (5) creation d'une climatologie de tro3_daylightB et tro3B sur toute la duree de la simulation => determination du masque M |
| 119 | indiquant la zone photosensible en filtrant le rapport (tro3B-tro3B_daylight)/tro3B |
| 120 | (6) Reconstruction de tro3_daylight(y), ou y est l'annee, qui vaut: |
| 121 | * tro3A(y) en-dehors du masque M |
| 122 | * tro3A(y)*tro3_daylightB(y)/tro3B(y) dans le masque M |
| 123 | <C> moyenne temporelle glissante (+/- 5 ans) pour les deux associations REF1/REF2 et REF1/SCNB2c |
| 124 | => moyennes zonales lissees tro3_${Y}.nc de 1961 a 2095 pour les variables tro3 et tro3_daylight. |
| 125 | <D> construction de fichiers synthetiques pour tous les scenarii (meme methode pour tro3 que pour tro3_daylight). |
| 126 | Dans ce qui suit, '*' designe le scenario considere ; 'S' et 'R' designent respectivement les scenarii REF1_SCN-B2c et REF1_REF2. |
| 127 | tro3* = tro3S + (tro3R-tro3S)*(CO2*-CO2S)/(CO2R-CO2S) |
| 128 | Les grandeurs CO2 representent les concentrations en CO2 dans les scenarii consideres. |
| 129 | N.B.: REF1_REF2 représente le scenario RCP6.0. |
| 130 | |
| 131 | => moyennes zonales lissees tro3_${Y}.nc de 1966 a 2093 pour tous les scenarii. |
| 132 | |
| 133 | |