- lupdstat
update the statistics
défaut : lupdstat = .false.
contrôles : ERR si merge=.true. ou lplot=.true.
Olive : lupdstat=.true. si OBSTAT_PRECALC (Use pre-calculated obs-model) ="no" ou
="yes" mais le fichier de statistiques n'existe pas
- carchformat
format du fichier
valeurs : netcdf (défaut) ou netepy
Olive : Statistics definition/type netcdf
- cbestwind
for scatt wind selection
détection bestwind par rapport à andep ou fgdep (for scatterometer datas)
défaut : '2an'
remarque : ajout de l'information dans sta()%comment
Olive : gnam.calc
- cexpver
experiment version
→ indication dans le fichier (kind ≠ 8)
kind = 8 : initialisation du nom de fichier, attribut expver dans le fichier, ligne 4 du titre du tracé
défaut : xxxx (calcul)
remarque : cexpver='xxxx' (globvar), cexpver(maxfiles,3)='' (globvar1 - plot kind=8)
Olive : cexpver=${OBSTAT_EXPERIMENT_XPID} (Experiment ID)
${OBSTAT_REFERENCE_XPID} (Reference ID)
- cobsformat
selects the input obs format
défaut : cobsformat = 'undef'
valeurs : 'ODB' 'USER'
pour odb2 (cobsformat='ODB') initaliser également
odbfiletype (MFP ou OFB)
Olive : cobsformat = 'ODB' (calcul des statistiques à partir de bases de données odb)
- cstream
${modele}${cutoff}.${vconf}
→ indication dans le fichier (kind ≠ 8)
kind = 8 : initialisation du nom de fichier
défaut : xxxx
Olive : cstream='${OBSTAT_EXPERIMENT_MODEL}${OBSTAT_EXPERIMENT_CUTOFF}'
OBSTAT_EXPERIMENT_MODEL = 'arp','aro','ald' (Experiment model)
OBSTAT_EXPERIMENT_CUTOFF = 'A', 'P', 'C' (Experiment cut-off)
OBSTAT_EXPERIMENT_VCONF = '4dvarfr', '3dvarfr', 'pifrance', 'aearp', 'courtfr', ...
- datecycle
date réseau (yyyymmddhhmn)
défaut : datecycle = '999999999999'
Olive : Date and time
- hsris_calc
type de calcul (airs/iasi/cris)
défaut : hsris_calc = 3 (channel number)
remarque : valeurs entre 1 et 5 (voir
airs_channels,
iasi_channels,
cris_channels,
description)
1 → presure (Pa)
2 → channel index
3 → channel number
4 → wave number (cm-1)
5 → frequency (microns)
Olive : gnam.calc
- imaxread
stop obs decoding after that many reads
défaut : imaxread = 999999999
remarque : limitation du nombre d'observations lues par pool
Olive : gnam.calc
- laddglobal calcul global des statistiques
défaut : laddglobal = .false.
remarque : action si
sta()%areatag ='NH.Tr.SH'
→ calcul pour "N.Hemis", "Tropics", "S.Hemis" et "GLOBE"
- lcalcdiag
calculs des diagnostics Desroziers
défaut : lcalcdiag = .true.
remarque : initialisation de %ldiag(item)
calcul des items 11, 28, 12 et 40 (cf
diagnostics)
Olive : Statistics definition/Calc. Desroziers diag.
- lgeogrid_norm
normalisation des calculs en haute résolution
défaut : lgeogrid_norm = .true.
remarque : information archivée dans le fichier netcdf GEO* ("normalised :" True/False)
Olive : gnam.calc
- lminute
use time steps as minutes instead of hours
défaut : lminute = .false.
remarque : information archivée dans l'entête du fichier TXT*
Olive : lminute = .true. si Date and time/step = 15mn
- lnanqcheck
with this option unrealistic values will be filtered
défaut : lnanqcheck = .false. (cf odbread.F90, odb2read.F90)
- lobstat_ecmwf
pas top - différences avec cep
défaut : lobstat_ecmwf = .true.
remarque : calcul différent dans dataqc.F90
Olive : lobstat_ecmwf = .false.
- lprthard
print hard statistics
défaut : lprthard = .true.
Olive : lprthard = .false. pour
lmerge = .true. et
lplot =.true.
remarque : écriture sur unité 6 ; chunk==pool
Référence de toutes les données dans la partie
"== BEGIN TPI STATISTICS =="
- lscaling (à partir du cycle 37)
activate scaling of values
défaut : lscaling = .true. (appel odbscaling.F90)
remarques :
pour codetype 1,29,57,58,181,182,183,184,185,186,188,206,193 et 250
(cf odbcode.cfg
PARAMS)
types 1, 57 : géopotential → mètre géopotentiel
types 29, 58 : humidité relative → %
types 285, 206 : ozone kg/m2 → Dobson
types 181 : NO2 → NO2 * 1308900 (1015 molecules/cm2)
types 182 : SO2 kg/m2 → Dobson
types 183 : CO → CO * 2150 (1018 molecules/cm2)
types 184 : HCHO → HCHO * 2004008 (1015 molecules/cm2)
types 186 : CO2 → CO2 * 10000000
types 188 : CH4 → CH4 * 10**9
types 250 : gps → gps * 10**6
- lwrtempty
write out missing data
défaut : lwrtempty = .false.
limitations : kind = 1,
sta()%nbin /= 0 et
sta()%nbin2 = 0
→ Ecriture des niveaux de valeurs manquantes (-9999.) sur
cstatoutFname (statout.txt)
- odbfiletype
ECMA,CCMA pour odb ; format MFB ou OFB pour odb2
- obstat_nbslot
nombre de slots pour obstat
nombre de divisions dans window_length_4d pour le calcul obstat → définir
timeslot (calcul sur une seule tranche horaire)
défaut : obstat_nbslot = 1
- specialsql
get specific values to cater for studid Graeme fix
défaut : specialsql = 'none'
Olive : specialsql = '$OBSTAT_SQL' (Build parameters/SQL request)
fonction : conditionne l'appel à la requête obstat_xxx.sql
xxx définis : 'conv', 'fcdep', 'fdcdep_gpsro', 'geos', 'gpsro', 'mwimg', 'radar',
'resat', 'satob', 'scatt', 'smos', 'smos_land', 'tovs'
- timeslot
numéro du timeslot utilisé
défaut : timeslot = 999
remarques : timeslot=0, 999 ou nbslots=1 → calculs sur tous les slots → même résultats
différence dans la première
ligne
exemple : window_length_4d=6, obstat_nblot=6, timelost=2 → calcul sur deuxième heure de la période
Olive : passage argument initialisation par gnam.calc
- window_length_4d
4d-var total window length
défaut : window_length_4d = 999
remarque : → sta()%nbr_tslot (window_length/tslot)
Olive : définition automatique à partir du choix du modèle
- window_offset_4d
4d-var offset time
défaut : window_offset_4d = 3
remarque : → calcul des pas de temps
- lplot
to generate plots
contrôles : if (nproc > 1) lplot=.false.
Peut être complété par les options
loverlay,
lplotdiff
ou
lploterrorbars
- plotmode
plotmode=0 1 graphe par fichier sauf airs,iasi,cris
plotmode=1 tous les graphes dans le même fichier
plotmode=2 1 graphe par fichier
- cgroupparam
liste des paramètres à grouper 'x1,x2'
défaut : cgroupparam = 'none'
x1[,x2] : numéro indiqué dans
cstat2Fname params
- cgroupstat
liste des statistiques ` grouper 'x1,x2,...'
défaut : cgroupstat = 'none'
x1,x2,... : numéro d'ordre dans
cstat2Fname ; peut être différent de celui de "
BEGIN STATDEF xx"
- clayout
plot layout option
défaut : clayout = 'A4'
valeurs : clayout = A4, A3, EPS
Olive : $OBSTAT_LAYOUT
- coeff_plt
pour gérer l'axe des X avec lploterrorbars
Bornes du tracé = coeff_plt * sigma(répartition des valeurs) si coeff_plt > 0.
défaut : 0.
- ctitle
titre du tracé
Première ligne du titre pour kind=8
kind=1 à 4 : utilisation des informations contenues dans le fichier de statistiques (
statout.txt)
dans les procédures olive ou d'exécutions directes
défaut : ctitle = 'untitled' → "." sur le plot (Pb magics)
- cmax,cmin,cint
bornes du tracé
défaut : -99.
- hsris_band_split (à partir du cycle 37)
split hsris into 4 spectral bands (LW,O3,WV and SW)
défaut : hsris_band_split = .true.
remarque : valeurs entre 1 et 5
- hsris_plot (à partir du cycle 37)
hsris (airs/iasi/cris) plot type
défaut : hsris_plot = 3
remarque : valeurs entre 1 et 5 (voir
airs_channels,
iasi_channels,
cris_channels,
description)
1 → presure (Pa)
2 → channel index
3 → channel number
4 → wave number (cm-1)
5 → frequency (microns)
Olive : gnam.plot
- humminlev
limit TEMP humidity from Surface to humminlev hPa
défaut : humminlev = 70
remarque : uniquement pour kind = 1 (
cstatinFname -
stat.def.odb)
utilisation pour subtype (odb) ≥ 35 et ≤ 40 , = 135
(
subtypes)
paramètre 7 (odb) = humidité spécifique (
params)
paramètres 2 et 7 (odb) = température et humidité spécifique (
params)
limitation du tracé aux pressions > humminlev
- lautoscale
échelle des axes à partir des valeurs
défaut : lautoscale = .true.
lautoscale = .false. → échelle à partir des valeurs définies
dans cmin/cmax.cint ou dans
general.cfg
- lfixYaxis
keep axes scales constant en Y airs/iasi/cris
défaut : lfixYaxis = .false.
- loghistogram
use log scaling for histogram bar plots
défaut : loghistogram = .false.
remarque : uniquement pour kind = 2 (
cstatinFname -
stat.def.odb)
- loverlay
overlay plot with 2ndary stats
défaut : loverlay = .false.
Olive : OBSTAT_Oplotype
remarques :
loverlay =.true. →
cstat2Fname existe
lectures et calculs en "sommes" sur
sta2()
tracé en rouge , noté "ref"
- lsymxaxis
keep x-axes for bias from min to max for most cases
défaut : lsymxaxis = .true.
remarque : uniquement pour kind = 1,2 (
cstatinFname -
stat.def.odb)
pas pour lploterrorbars/rms
lsymxaxis = .false. → échelle asymétrique
- minpop
min. population required to plot rms/bias
défaut : minpop = 5
Olive : $
OBSTAT_MINPOP
remarque : uniquement pour kind = 1 (
cstatinFname -
stat.def.odb)
comparaison avec la valeur de "number" dans
statout.txt"
- lplotdiff
plot différence entre deux expériences
défaut : lplotdiff = .false.
- lplotdiag
plot diagnostics Desroziers (pour lploterrorbars = .false.)
défaut : lplotdiag = .false.
- lplotempty
plot dataless plots avec plotmode=1
défaut : lplotempty = .false.
- lploterrorbars
plot avec barre d'erreur
défaut : lploterrorbars = .false.
remarques :
pour kind=1 uniquement
→
lplotdiag=.true.
→ utilisation de fichier TXTD_* pour la référence et TXT_* pour l'expérience
TXT_* expérience et référence → TXTD_* avec un merge
lcalcdiff=.true.
- lplotisol
plot isolated levels for airs, iasi, cris
défaut : lplotisol = .false.
remarque : uniquement pour kind = 1 (
cstatinFname -
stat.def.odb)
- lplotnorm
avec deux expériences
défaut : lplotnorm = .false.
- lplotobs
plot curves with observations
plot item idval/obsvalue@body (8) s'il est présent
défaut : lplotobs = .false.
remarque : pour kind = 1,2,3 (
cstatinFname -
stat.def.odb)
- lstdev
plot standard deviation instead of rms
défaut : lstdev = .false.
Olive : $
OBSTAT_STDEV
remarque : uniquement pour kind = 1,5 (
cstatinFname -
stat.def.odb)
- lwarning
change axis colour to red if data outside bounds
défaut : lwarning = .true.
- lwhiskers
plot box-and-whiskers over histograms
défaut : lwhiskers = .true.
