Changes between Version 52 and Version 53 of InfosOzone


Ignore:
Timestamp:
01/19/11 19:15:44 (14 years ago)
Author:
cugnet
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • InfosOzone

    v52 v53  
    5252 * Passage par défaut à la version v2 des fichiers climoz. Voir : [1150] [1151] 
    5353 
    54 == Génération des fichiers hybrides bruts == 
    55  
    56 Ce sont des fichiers annuels de 12 champs (un par mois). Ils contiennent quatre champs: 
    57   * tro3 (ozone de jour + nuit) 
    58   * tro3_daylight (ozone de jour) 
    59   * Ptrop (pression à la tropopause, utilisée pour l'hybridation) 
    60   * night (0: nuit ; 1: jour ; s'applique à l'ozone de jour) 
     54 
     55 
     56 
     57== GENERATION DES FICHIERS HYBRIDES == 
    6158 
    6259 
     
    7875  * Pour RUN=RCP26, RCP45 ou RCP85: LOI_IPCC_${RUN}_${Y}${M}01_${Y}${M}30_1M_inca_avgr.nc[[BR]] 
    7976 
     77 
     78 
    8079=== Principe résumé === 
    8180 
     
    8988 
    9089 
     90 
    9191=== Principe détaillé === 
    9292 
    93 ETAPE 1: Obtention de moyennes zonales lissées temporellement. 
    94  
    95 TROPOSPHERE: 
    96  
    97 Pour chaque scenario: 
     93 
     94== ETAPE 1: Obtention de moyennes zonales lissées temporellement. == 
     95 
     96 
     97     == TROPOSPHERE: == 
     98 
     99Pour chaque scenario (3, RCP26, RCP45, RCP60, RCP85): 
    98100  1) extraction des variables O3 et pmid 
    99101  2) interpolation sur les 19 niveaux standard (variable "presnivs") ; moyenne zonale. 
    100102  3) pour RUN=RCP26, RCP45 ou RCP60: moyenne mensuelle, concaténation. 
    101103=> fichiers annuels à 12 champs (un par mois) 
     104 
     105Pour chaque paire de scenarii 3 (de 1850 à 2000), ${RUN} (de 2001 à 2100): 
    102106  4) moyenne temporelle glissante (+/- 5 ans) 
    103107 
    104108'''PRODUITS''': moyennes zonales lissees tro3_${Y}.nc de 1855 a 2095. 
    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 
     114Pour 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. 
     124On 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 
    125138  <C> moyenne temporelle glissante (+/- 5 ans) pour les deux associations REF1/REF2 et REF1/SCNB2c 
    126139      => moyennes zonales lissees tro3_${Y}.nc de 1961 a 2095 pour les variables tro3 et tro3_daylight.