105 | | NB: jusqu'en 2000, on garde les fichiers composites 3_RCP85. Au-delà, on utiliser les fichiers 3_${RUN}. |
106 | | |
107 | | STRATOSPHERE: |
108 | | Sorties utilisees: sur gaya, dans les repertoires /u/rech/hou/rhou079/LMDZ4/output/h2o-${RUN}-KE-8mpi-1omp/$Y/ |
109 | | * Partie historique (1960-2006): RUN=REF1, fichiers histday_tra_${Y}${M}.nc et histday_ins_${Y}${M}.nc |
110 | | * Partie a venir (2007-2098): RUN=REF2 et RUN=SCNB2c, memes noms de fichiers |
111 | | |
112 | | Pour chaque run (REF1, REF2 et SCNB2c): |
113 | | <A> ozone de jour+nuit (tro3A): meme demarche qu'en troposphere pour les fichiers histday_tra (moyennes journalieres). |
114 | | <B> ozone de jour: |
115 | | (1) Masquage des pixels de nuit (obtenu avec un indicateur de nuit, en l'occurrence l'espece photosensible NO) |
116 | | (2) Interpolation sur les niveaux standard. |
117 | | (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 |
118 | | les 10, 20 et 30 du mois) => moyennes zonales d'ozone de jour (tro3_daylightB) de 1961 a 2006. |
119 | | (4) Operations (2) et (3) (sans masquage) => moyennes zonales d'ozone de jour+nuit (tro3B) de 1961 a 2006. |
120 | | (5) creation d'une climatologie de tro3_daylightB et tro3B sur toute la duree de la simulation => determination du masque M |
121 | | indiquant la zone photosensible en filtrant le rapport (tro3B-tro3B_daylight)/tro3B |
122 | | (6) Reconstruction de tro3_daylight(y), ou y est l'annee, qui vaut: |
123 | | * tro3A(y) en-dehors du masque M |
124 | | * tro3A(y)*tro3_daylightB(y)/tro3B(y) dans le masque M |
| 109 | |
| 110 | |
| 111 | |
| 112 | == STRATOSPHERE: == |
| 113 | |
| 114 | Pour chaque scenario (REF1, REF2 et SCNB2c ; on donne en exemple le cas de REF1): |
| 115 | |
| 116 | (A) OZONE DE JOUR+NUIT (fichiers histday_tra, moyennes journalières): |
| 117 | 1) extraction des variables O3 et press |
| 118 | 2) interpolation sur les 50 niveaux standard (variable "presnivs") ; moyenne zonale. |
| 119 | 3) moyenne mensuelle, concaténation. |
| 120 | => 12 champs mensuels d'ozone de jour+nuit moyen tro3(y) pour y variant de 1960 à 2006 |
| 121 | |
| 122 | (B) OZONE DE JOUR (fichiers histins_tra, champs instantanés les 10, 20 et 30 de chaque mois): |
| 123 | 1) Masquage des pixels de nuit. |
| 124 | On utilise l'espèce photosensible NO sur le 45ième niveau (1.527329 hPa): une concentration supérieure à 1.0E-15 vmr indique la nuit. |
| 125 | 2) interpolation sur les 50 niveaux standard (variable "presnivs") ; moyenne zonale. |
| 126 | 3) moyenne mensuelle: O3m($M)=[O3(${M-1}30) + 2.O3(${M}10) + 2.O3(${M}20) + O3(${M}30)]/6 |
| 127 | => ozone de jour tro3_daylight_Ins(y) pour y variant de 1961 à 2006 |
| 128 | 4) étapes 2 et 3 sans masquage |
| 129 | => ozone de jour+nuit tro3_Ins(y) pour y variant de 1961 à 2006 |
| 130 | 5) création d'une climatologie de tro3_daylight(Ins) et tro3(Ins) sur toute la durée de la simulation => détermination du masque M |
| 131 | indiquant la zone photosensible en filtrant le rapport (tro3(Ins)-tro3_daylight(Ins))/tro3(Ins) |
| 132 | => 12 champs mensuels d'ozone de jour tro3_daylight(y) pour y variant de 1960 à 2006 |
| 133 | |
| 134 | 6) Reconstruction de tro3_daylight(y) (y est l'année), qui vaut: |
| 135 | * tro3(y) hors du masque M |
| 136 | * tro3(y)*tro3_daylight_Ins(y)/tro3_Ins(y) dans le masque M |
| 137 | |