540 likes | 641 Views
Cycle MASDEV4_7. Préparation de cas réels et idéaux. Spawning à 2 modèles et zoom des PGD Spawning des champs de surface Coordonnée SLEVE Filtrage de vortex: améliorations Chimie: schéma ReLACS avec OH pronostique Chimie et PREP_REAL_CASE. B. Simulation. DYNAMIQUE.
E N D
Cycle MASDEV4_7 • Préparation de cas réels et idéaux • Spawning à 2 modèles et zoom des PGD • Spawning des champs de surface • Coordonnée SLEVE • Filtrage de vortex: améliorations • Chimie: schéma ReLACS avec OH pronostique • Chimie et PREP_REAL_CASE
B. Simulation DYNAMIQUE • Schémas d’advection centré du 4ème ordre (CEN4TH) et schéma PPM pour les scalaires • Terme de Nudging (LNUDGING) • Diffusion numérique (LZDIFFU) • Nesting pour la surface (XWAY=3.) • Ajustement du nbre d’itérations du solveur de pression pendant le run • Microphysique chaude à 2 moments pour les Sc (KHKO) • Microphysique froide à 1 moment avec grêle (ICE4) • ICE3/ICE4 : introduction de sédimentation d’eau nuageuse C2R2/KHKO : introduction de l’activation par refroidissement radiatif • Nouvelle formulation de l’autoconversion Chaboureau et Pinty (2006) • Modifications du rayonnement (NOVLP) • Turbulence: moments du 3e ordre pour Theta (TOMs) • Amélioration de longueur de mélange de Deardorff • Chimie: aérosols (ORILAM) • Poussières désertiques (DUST) • Paramétrisation des éclairs dans le schéma de convection MICROPHYSIQUE RAYONNEMENT TURBULENCE CHIMIE
Cycle MASDEV4_7 C. Diagnostiques • Températures de brillance RTTOV • Calcul de rafales et de co-indices de réfraction • Nouveaux diagnostiques LES • Réinitilisation des traceurs lagrangiens D. Code • Nettoyage : • - suppression des $n (doc de Didier) • - suppression des branches mortes • 2. Bugs • 3. Suppression de l’option de compilation sur Fujitsu :-Wv -Ma
A.1 Spawning à 2 modèles Zoom des PGD 1: modèle père P 2: modèle fils F1 = interpolation champs P même résolution = interpolation champs P + recopie champs F1 3: modèle fils F2 • Intérêt: simulation à haute résolution de systèmes propagatifs
A.1.a. Spawning à 2 modèles &NAM_LUNIT2_SPA CINIFILE=‘Fichier_P’ YDOMAIN=‘PGD_zoomé_F2’ YSPAFILE=‘Fichier_F2’ / • Contribution: D. Barbary (CRC) • Recopie des champs (d’altitude) du frère dans le domaine commun • (les champs de surface restent interpolés à partir du père) • Mise en œuvre: dans SPAWN1.nam YSONFILE=‘Fichier_F1’ domaine_F2 = domaine_F1 • Cas particulier testé: spawning « classique » F1 P on a vérifié que champs_F2 = champs_F1 spawning à 2 modèles F2 P+F1
A.1.b. Zoom des PGD • Cycle MASDEV4_7 (surfex 1.3) • nouveau programme fichier_PGD programme ZOOM_PGD fichier_PGD_zoomé PRE_ZOOM1.nam • Cycles antérieurs et MASDEV4_5: • possibilité de zoomer lors de PREP_REAL_CASE • (&NAM_MESONH_DOM dans PRE_REAL1.nam) • Cycle MASDEV4_6 (surfex0.6): • zoom impossible • (domaine du PGD = domaine de simulation)
A.2. Spawning des champs de surface • Cycles antérieurs et MASDEV4_5: • SPAWNING: interpolation des champs d’altitude • et des champs de surface • Cycle MASDEV4_6 (surfex0.6): • SPAWNING: interpolation des champs d’altitude • PREP_REAL_CASE: interpolation des champs de surface • &NAM_PREP_SURF_ATM CFILE=‘input_spawn’ • Cycle MASDEV4_7(surfex1.3): • SPAWNING: interpolation des champs d’altitude et de surface • &NAM_SPAWN_SURF_ATM LSPAWN_SURF=.TRUE./(défaut)
A.3. Coordonnée SLEVE Rappel: coordonnée Gal-Chen, Sommerville Coordonnée SLEVE: Calcul d’une orographie lissée: ZSMT dansPREP_IDEAL_CASE, PREP_PGD, PREP_NEST_PGD) Ajout des variables: LSLEVE, XZSMT, XLEN1, XLEN2
A.3. Coordonnée SLEVE: cas test ReunionOrographie lissée SLEVE GC WM THM LCOLAREA=T LMARKER=T ZS ZSMT (moyenne sur 9 points)
A.3 Coordonnée SLEVE • Simulation: stretching (Sleve) > stretching(Gal-Chen) utiliser CRESI plutôt que RICHA pour le solveur (nouveau défaut) augmenter le nombre d’itérations NITR (LITRADJ: permet d’ajuster le nombre d’itérations au cours du run voir diapo « B.5 ») • Grid-nesting: PREP_NEST_PGD: • Remplacement du relief lissé du fils par celui du père interpolé • (plutôt que remplacement du champ du père par celui du fils moyenné)
A.4. Prep_real_case: filtrage de vortex • Problème dans les cas avec forts gradients de vent dans le domaine • (présence d’une tempête tropicale en plus du cyclone): • modification du vent en dehors de la zone de filtrage. • Correction dans removal_vortex.f90: • vent modifié (Schuman+filtrage+Schuman) uniquement • dans la zone de filtrage • Testée sur le cas ERICA • Passage dans PRE_REAL1.nam du nom du père du fichier bogussé: • &NAM_HURR_CONF CDADBOGFILE=‘fichier_père_bogussé’ • CDADATMFILE=‘fichier_père_filtré’ / • (évite de commenter un test dans ini_sizen lors du run multi-modèles…)
A.5. Chimie • Pas de contrôle sur les dates des 2 fichiers d’entrée • schéma ReLACS avec OH pronostique • Prep_real_case: possibilité d’initialiser champs atmosphériques • et champs chimiques à partir de 2 sources différentes. Arpege U,v,q,T Mocage ou MesoNH O3,NO,CO,… FM • Mise en œuvre: dans PRE_REAL1.nam &NAM_FILE_NAMES HATMFILE=‘Fichier_Arpege’, HATMFILETYPE=‘GRIBEX’, HCHEMFILE=‘Fichier_Mocage’,HCHEMFILETYPE=‘GRIBEX’/ OU &NAM_FILE_NAMES HATMFILE=‘Fichier_Arpege’, HATMFILETYPE=‘GRIBEX’, HCHEMFILE=‘Fichier_MesoNH’, HCHEMFILETYPE=‘MESONH’/
B.1 Nouveaux schémas d’advection • Schéma centré du 4ème ordre (CEN4TH) : Préconisé pour l’advection de U, V, W • PPM pour les scalaires : PPM_00 : sans contrainte • PPM_01 : monotone et défini positif • Mise en œuvre : dans EXSEGn.nam • &NAM_ADVn CUVW_ADV_SCHEME = ‘CEN4TH’ (défaut) pour u,v,w • ou ‘CEN2ND’ • CMET_ADV_SCHEME = ‘FCT2ND’ (défaut)pour q, r, TKE ou CEN2ND ou CEN4TH ou MPDATA ou PPM_00 ou PPM_01CSV_ADV_SCHEME = ‘FCT2ND’ (défaut)pour scalaires CEN2ND ou CEN4TH ou MPDATA ou PPM_00 ou PPM_01 • Suppression de CFV_ADV_SCHEME • Schéma temporel FIT pour l’advection des variables selon PPM et LF pour le reste du modèle Filtre d’asselin avec XASSELIN_SV (0.02 Par défaut) en plus de XASSELIN (0.2 Par défaut)
B.1 Réorganisation de la routine d’advection ADVECUVW HUVW_ADV_SCHEME ADVECUVW_4TH ADVECMET HMET_ADV_SCHEME ADVECMET_4TH FCT_MET ADVECTION MPDATA PPM_MET ADVECSCALAR HSV_ADV_SCHEME ADVECSCALAR_4TH Conseil sur le choix des schémas : - PPM_00 pour CMET_ADV_SCHEME - PPM_01 pour CSV_ADV_SCHEME FCT_SCALAR MPDATA_SCALAR PPM_SCALAR
B.2. Simulation: ajout d’un terme de NUDGING • Principe: rappel des variables pronostiques vers les valeurs LS • sur tout le domaine de simulation. • Ajout d’un terme source: • Mise en œuvre: dans EXSEGn.nam &NAM_NUDGINGn LNUDGING=.TRUE. XTNUDGING=21600. / • Bien qu’activable pour chaque modèle, l’intérêt est surtout • pour le modèle père (évite la dérive par rapport au modèle coupleur).
B.3 Simulation: diffusion numérique • Opérateur de diffusion adapté aux fortes pentes (Zängl, 2002) • Mise en œuvre: dans EXSEG1.nam &NAM_DYN LNUMDIFF=.TRUE., LZDIFFU=.TRUE. / XT4DIFF dans NAM_DYNn n’est pas utilisé THM_CV_ LZDIFFU=T LZDIFFU=F
B.3 Simulation: diffusion numérique Simulation MAP-POI8 (N.Asencio) : coupes verticales Sud-Nord Theta 12h LZDIFFU=T LZDIFFU=F
B.4 Simulation: Nesting pour la surface • Principe: Remontée d’information du fils vers le père pour les champs 2D en input de la surface : précipitations instantanées et cumulées explicites et convectives (INPRR, INPRS, INPRG, INPRH, PRCONV, PRSCONV, ACPRR, ACPRS, ACPRG, ACPRH, PACCONV), et les flux radiatifs SW (DIRFLASWD, SCAFLASWD, DIRSRFSWD) • But : Meilleur réalisme des champs de précipitation et d’humidité du sol du modèle père dans la zone de recouvrement • Mise en œuvre: dans EXSEG1.nam &NAM_NESTING XWAY(n° fils)=3. /
B.4 Simulation/ Nesting pour la surface Cas AMMA 30 Août 2004 Cumul 24h : 6utc-6utc 10km+5km XWAY=3 10km seul 10km+5km XWAY=2 Le 5km permet d’obtenir le système sur le Burkina Faso NOAA CPC
B.4 Nesting pour la surface Surestimation des pluies du père XWAY(3)=2 XWAY(3)=3 MAP POI2a 18/9/99 - 00TU Modèle 2 Cumul 12h Obs max=102mm max=53mm Modèle 3 max=93mm max=82mm
B.5 • Ajustement du nombre d’itérations du solveur de pression Possibilité d’ajuster le nombre d’itérations du solveur pendant le run (test sur la divergence résiduelle à chaque pas de temps: si div<1E-9 alors NITR=NITR-1 si div>1E-8 alors NITR=NITR+2 ) • Mise en œuvre: dans EXSEGn.nam • &NAM_DYNn LITRADJ= T
B.6 Microphysique chaude à 2 moments Khairoudinov-Khogan • Mise en œuvre: dans EXSEGn.nam &NAM_PARAMn CCLOUD= ‘KHKO’ / Mêmes réglages que pour C2R2 dansNAM_PARAM_C2R2 Dans le code, routines : ini_rain_c2r2.f90, rain_khko.f90, c2r2_adjust.f90 NCCN (Paramètre) Nact (m- 3) Activation : Cohard et al (1998) (C2R2) Autoconversion : KHKO qc (kg/kg) qr (kg/kg) qv (kg/kg) Cond / évap : Langlois (1973) (C2R2) Nc (m-3) Accretion : KHKO Nr (m-3) Sedimentation : KHKO Evaporation : KHKO Extrait de Geoffroy (2005) • Principe: Nouveau schéma uniquement pour les simulations LES de Sc • (rc max de 1 à 2g/kg) • Les processus relatifs à la pluie de C2R2 sont modifiés pour être adaptés à la pluie faible des Sc (les processus de formation nuageuse et l’ajustement sont ceux de C2R2)
B.7 Microphysique mixte à 1 moment avec grêle • Principe: 4 espèces de glace (ice, snow, graupel, hail) au lieu de 3 dans ICE3. • (champs supplémentaires RHM et RHT) • Mise en œuvre: dans EXSEGn.nam &NAM_PARAMn CCLOUD= ‘ICE4’ / • Dans le code, mêmes routines que ICE3 : ini_rain_ice.f90, rain_ice.f90, ice_adjust.f90 Voir Thèse F.Lascaux
B.8 ICE3/ICE4 et KESS : Modifications de la sédimentation • Sédimentation de l’eau nuageuse à partir d’une concentration en gouttelettes diagnostique (différence mer/continent/ville) : active uniquement pour nuages bas (stratus, brouillard, effet négligeable pour les Cu et Cb) • Mise en œuvre: dans EXSEGn.nam • &NAM_PARAM_ICE LSEDIC=T (LSEDIC=F par défaut) • Nouveau champ 2D INPRC fourni à la surface • Sédimentation dernier processus microphysique lent (et non plus premier) • .Algorithme de sédimentation :Distinction entre la pluie créée pendant le sous-pas de temps courant et celle restant des sous-pas de temps précédents • Avant : on appliquait n fois une chute pendant Dt/n à une quantité créée pendant Dt • Maintenant : on applique n fois une chute pendant Dt/n à une quantité créée pendant Dt/n
Cas 1D de brouillard – ICE3-Evolution temporelle sur 18h à partir de 18TU (schéma LW : RRTM, NOVLP=5) Sans sédimentation de rc Rc max=0.6g/kg q Rr max=0.007g/kg Avec sédimentation de rc q Rc max=0.5g/kg Rr max=0.004g/kg
B.8 C2R2/KHKO : Modifications • ACTIVATION : Introduction d’un terme de refroidissement radiatif dans l’équation d’évolution de la sursaturation pour les nuages bas (brouillards, stratus) • Mise en œuvre: dans EXSEGn.nam • &NAM_PARAM_C2R2 LACTIT=T (LACTIT=F par défaut) • Prise en compte d’une variation verticale de C fonction de la stabilité dans NCCN=C.Sk : C’=C.exp(-K.dq/dz). Par défaut, K=0. • SEDIMENTATION dernier processus microphysique lent (et non plus premier) • .SEDIMENTATION de l’eau nuageuse pour C2R2 : Introduction du facteur correctif de Cunningham pour la vitesse terminale de chute (comme dans ICE) • Nouveau champ 2D INPRC fourni à la surface
Cas 1D de brouillard – C2R2 LRAIN=F, LACTIT=F LRAIN=F, LACTIT=T LRAIN=T, LACTIT=T Rc max=0.5g/kg Rc max=0.5g/kg Rc max=0.9g/kg Nc max=453/cm3 Nc max=450/cm3 Rr max=0.02g/kg
Cas 1D de brouillard – KHKO C2R2 : LRAIN=T, LACTIT=T KHKO : LRAIN=T, LACTIT=F KHKO : LRAIN=T, LACTIT=T Rc max=0.9g/kg Rc max=0.4g/kg Rc max=0.4g/kg Rr max=0.02g/kg Rr max=0.005g/kg Rr max=0.003g/kg
B.9 RAYONNEMENT • Différentes modifications ont été apportées au code de transfert radiatif (Thouron, 2006) : • Corrections en SW sur le calcul des propriétés optiques (transformation de Delta-Eddington) et sur le calcul de l’efficacité de diffusion : routines swclr.f90, swr.f90, swde.f90, swni.f90 • Suppression des hypothèses de recouvrement pour le calcul de l’efficacité de diffusion : nouvelles options de recouvrement (paramètre NOVLP dans ini_radconf.f90): • NOVLP=5 : Recouvrement aléatoire pour la fraction de ciel clair : Imposé lorsque la condensation sous-maille n’est pas activée. • NOVLP=6 : Recouvrement maximum aléatoire pour la fraction de ciel clair : Adapté en présence de nuage multi-couches (valeur par défaut sauf si SUBGCOND, avec nécessité d’associer RRTM pour LW) . • NOVLP=7 : Recouvrement maximum pour la fraction de ciel clair : Adapté en l’absence de nuage multi-couches (faible résolution verticale). • NOVLP=8 : Recouvrement maximum aléatoire pour la fraction de ciel clair et pour l’angle zénithal solaire (option précédemment imposée d’office, non conseillée) • Modification de l’albédo de diffusion simple pour l’absorption SW par les inclusions d’aérosols dans les gouttelettes de nuage (Sandu et al., 2005) • Prise en compte des poussières (voir DUST)
B.10 Turbulence : Moments d’ordre 3 pour q • Principe: Meilleure représentation du contre-gradient dans les CL convectives • ( Tomas et Masson, 2006) • Les 2 moments d’ordre 3 relatifs à q sont ‘fittés’ à partir de simulations LES zi=Hauteur d’inversion, préalablement diagnostiquée w*=w convective q*= q convective • Mise en œuvre: dans EXSEGn.nam &NAM_TURBn CTOM = ‘TM06’ / CTOM=‘NONE’ par défaut • Modification de la longueur de mélange de Deardorff Critère de stabilité calculé en (ql,rt) et non plus qv
B.11 Chimie des aérosols: ORILAM SURFEX: Reconnaissance passage chimie gazeuse et aérosols via le nom des traceurs: - activation dépôt sec si WES89 - activation émission chimique / aérosols si spécifiée en namelist de surface • Représentation de type log-normale (Tulet et al., 2005) • Mise en oeuvre: dans EXSEG1.nam : nouvelle namelist NAM_ORILAM &NAM_ORILAM LORILAM=T avec&NAM_CH_MNHCnLUSECHEM=T - Nombre de modes: 2 (Aitken et mode d’accumulation) - Dispersion variable ou non (LVARSIGI pour le mode d’Aitken, LVARSIGJ pour l’accumulation) -Activation sédimention (LSEDIMAERO) -Activation lessivage par pluies convectives (LCH_CONV_SCAV dans NAM_CH_MNHCn) - Activation type de paramétrisation pour nucléation (CNUCLEATION) -Activation type d'équilibre thermodynamique (4 pour la fraction inorganique (CMINERAL), et 2 pour les organiques (CORGANIC)) - Fonctionne avec ReLACS / RACM (inorganique uniquement) et CACM / ReLACS2 (inorganique + organiques) (Pour une représentation correcte, il est recommandé d’avoir HNO3, H2SO4, NH3 et CO)
B.12 Poussières désertiques (DUST) Activation à partir de Méso-NH: - Nombre de modes: 1 à 3 (NMODE_DST) - Dispersion variable ou non (LVARSIG) -Activation sédimention (LSEDIMDUST) -Activation lessivage par pluies convectives (LCH_CONV_SCAV dans NAM_CH_MNHCn même si LUSECHEM =.F. ) Reconnaissance par SURFEX (via le nom des traceurs); de la demande d'émission de poussières et des caractéristiques (modes, dispersion ...etc) Activation du modèle DEAD dans ISBA: Emissions fonction: 1) Type de surface (sol nu ou Rock) 2) Contenu en eau du sol 3) Vitesse de friction (saltation) Retour Méso-NH Transport / sédimentation / lessivage Modification du schéma rayonnement ECMWF: modification de l'absorption/diffusion de la classe dust Interpolation d'un code de Mie en fonction de la représentation lognormale des modes N, SIGMA et Rg Si LDUST=T, l’effet radiatif direct est automatiquement pris en compte à partir d’une table d’interpolation de SHDOM Représentation de type log-normale (Grini et Tulet, 2005) Dans EXSEG1.nam : nouvelle namelist NAM_DUST (activable avec ou sans chimie)
B.13 Paramétrisation des éclairs dans le schéma de convection • Principe: calcul d’une source de NOx dans le schéma de convection profonde et des quantités d’éclairs associés. si la chimie est activée, ajout de cette source à la variable chimique NO, s’il n’y a pas de chimie, nouvelle variable scalaire LINOX. réf.: Mari, C., J.-P. Chaboureau, J.-P. Pinty, J. Duron, P. Mascart, J.-P. Cammas, F. Gheusi, T. Fehr, H. Schlager, A. Roiger, M. Lichtenstein, and P. Stock: Regional lightning NOx sources during the TROCCINOX experiment.Atmos. Chem. Phys. Discussion. • Mise en œuvre: dans EXSEGn.nam • &NAM_CH_MNHCn LCH_CONV_LINOX=T Eclairs intra-nuages LiNOx(300hPa)
Réinitialisation des traceurs Lagrangiens Possibilité de réinitialiser les traceurs Lagrangiens plusieurs fois pendant un segment (et non plus seulement en début de segment). • Mise en œuvre: dans EXSEG1.nam • &NAM_CONF CINIT_LG=‘FMOUT’ / • FMOUT: réinitialisation lors de l’écriture de chaque fichier de sortie
C.1. Diag: RTTOV • Calcul de température de brillance avec le code RTTOV (v8_7). • (pas de suppression du code de transfert radiatif ondes courtes CRAD_SAT) • Mise en œuvre: dans DIAG1.nam • &NAM_DIAG NRTTOVINFO(:,1)=3 6 20 0 /
C.2. Diag: autres • Calcul de rafales (comme en masdev4_5) • Mise en œuvre: dans DIAG1.nam • &NAM_DIAG_SURFn N2M=2 / • Calcul de co-indices de réfraction (M et N) (présentation Vivien ) • Mise en œuvre: dans DIAG1.nam • &NAM_DIAG LCOREF=T / • Nouveaux diagnostics LES (présentation Fleur)
D1. CODE : Nettoyage • Suppression des $n :Une documentation est disponible sur le site Web. • Suppression des “branches mortes” : • Suppression de CDRAG (&NAM_PARAMn) • Suppression du schéma de turbulence K-e : CTURBLEN=‘KEPS’ (&NAM_TURBn) suppression des variables pron. EPS et de LHORELAX_EPS dans NAM_DYNn • Suppression des CLL de type Davies: CLBCX,Y=‘DAVI’ • Suppression de l’interpolation de Clark et Farley • Suppression de XWAY( )=0 dans NAM_NESTING : pas d’interaction père-fils • Suppression de CCLOUD=‘KES2’ • Suppression de CFV_ADV_SCHEME dans NAM_ADVn • Fichiers de namelist: • EXSEG1.nam : • - LSTEADY_LS n’est plus variable de NAM_DYN • mis à .TRUE. si NAM_LUNITn contient des fichiers de couplage, • à .FALSE. sinon. • - Dans NAM_CONFn, les LUSERn ne sont plus à préciser, sauf LUSERV (notamment dans le cas où CCLOUD=‘NONE’) • - Suppression de CHEVRIMED_ICE dans NAM_PARAM_ICE • - Dans NAM_TURBn défauts modifiés : LSIGMAS=T et LSIG_CONV=F si LSUB_COND=T • - Suppression de LRAD_DIAG
D1. CODE : Nettoyage • Fichiers de namelist (suite): • PRE_IDEA1.nam • Suppression de L1D et L2D de NAM_CONF_PRE • Ajout d’une nouvelle namelist NAM_VER_GRID (la même que dans PRE_REAL1.nam) avec -NKMAX déplacé de NAM_DIMn_PRE à NAM_VER_GRID • -LTHINSHELL déplacé de NAM_CONF_PRE à NAM_VER_GRID • -CZGRID_TYPE, XDZGRD, XDZTOP, XZMAX_STRGRD, XSTRGRD, XSTRTOP déplacés de NAM_GRIDn_PRE à NAM_VER_GRID en changeant la première lettre (C->Y, X->Z) • -NAM_GRIDn_PRE renommé en NAM_GRIDH_PRE (avec comme variables restantes XLATCEN, XLONCEN, XDELTAX, XDELTAY, XHMAX, NEXPX, NEXPY, XAX, XAY, NIZS, NJZS) • Dans le cas CSTN, ajout de la date au format yyyy mm dd sec (ex: 2000 01 01 0.) • Fichiers FM: • Après PREP_REAL_CASE et PREP_IDEAL_CASE (CSTORAGE_TYPE=‘TT’), • stockage des variables en M uniquement dans le fichier FM • PREP_REAL_CASE : • read_grib_field.f90: Augmentation de la taille du buffer JPACK pour lire les fichiers CEP à partir du 01/02/2006 (cycle c30r1: T799L91)
4_5 FIRE 4_6 Max=0.591 g/kg rc rc Max=0.606 g/kg 4_7 24h 24h NOVLP pour AZE + KHKO + PPM_01 Max=0.649 g/kg 24h
IDF 04/08/1994 Pluies totales cumulées 18h Obs 4_4 4_5 4_6 Max=34mm Max=32mm Max=28mm
Cas LG sur IdF CUVW_ADV_SCHEME = CEN4TH CMET_ADV_SCHEME = PPM_01 PPM_00 CEN4TH FCT2ND
4_4 4_3 MAP POI2a 17/9/99 - 19TU Obs 4_5 4_6
MAP POI2a 17/9/99 - 19TU Obs 4_6/FCT/ICE3 4_7/FCT/ICE3 4_7/FCT/ICE4 4_7/PPM_00/ICE3
4_4 4_3 MAP POI2a 17/9/99 - 20TU Obs 4_5 4_6 max=53dB
Max=51.6dB Max=53.4dB MAP POI2a 17/9/99 - 20TU Max=52dB Obs 4_6/FCT/ICE3 4_7/FCT/ICE3 Max=62.7dB Max=52.0dB 4_7/FCT/ICE4 4_7/PPM_00/ICE3
4_4 4_3 MAP POI2a 17/9/99 - 21TU Obs 4_5 4_6
Max=54dB Max=54dB MAP POI2a 17/9/99 - 21TU Max=50dB Obs 4_6/FCT/ICE3 4_7/FCT/ICE3 Max=64dB Max=53.6dB 4_7/FCT/ICE4 4_7/PPM_00/ICE3
4_4 4_3 MAP POI2a 17/9/99 - 22TU Obs 4_5 4_6
Max=54dB Max=57dB MAP POI2a 17/9/99 - 22TU Max=51dB Obs 4_6/FCT/ICE3 4_7/FCT/ICE3 Max=67dB Max=55.6dB 4_7/FCT/ICE4 4_7/PPM_00/ICE3
4_4 4_3 MAP POI2a 17/9/99 - 23TU Obs 4_5 4_6