To set up additional information required by constituent obs and
not provided in obsSpaceData. Also to assign observation layer
top and bottom levels (and averaging kernel matrices).
See ‘oopc_CHobsoperators’.
Arguments
kmode [integer ,in] :: Mode of observation operator
assim_* arrays could instead be made available to all families
by moving them to a different input namelist (and changing its
dimensions settings).
genOperConstraintType
Reference profile type for weighted integration
or layer averaging (generalized observation)
operator.
Relevant for operatorSubType(2,i)=’genOper’.
================================================
‘Trial’ use trial field xb for mass weighted
increment distribution
‘Diff’ use a combination of the difference of
an external reference xc and the trial
field xb, i.e. mass weighted increment
distribution as a(xc-xb) + b*xc where a
and b depend on the size of
sum[(xc-xb)/sig(xb)]^2 over the profile
genOperHCorrlenExpnt
Used with operatorSubType(2,i) =’genOper’
Exponent for partially mitigating the effect of
the influence of neighbouring column amonunt obs
from background error correlations.
Emperically obtained exponent value.
Not optimal for all possible local horizontal data densities.
genOperOmAStatsFactor
OmA RMS (or std dev) conservation factor for
operatorSubType(2,i) =’genOper’.
assim_fam
List of families to which filt_diagnOnly is to
apply.
assim_exclude_flag
Array specifying bits for identifying
diagnostic-only observations for observations
that would otherwise be assimilated according to
the other assim_* arrays
assim_exclude_nflag
Number of bit flags to specify in
assim_exclude_flag array
assim_all
Logical indicating if all assimilatable obs of
the specified family will be assimilated
(default is .true.)
When assim_all is .false., account for the setttings
of assim_num, assim_varno, assim_stnid, assim_nlev.
assim_num
Relevant when assim_all = ,false.
Number combinations (stnid, bufr element,
multi/uni-level) identified for assimilation.
All others will not be assimilated. OmP and OmA
diagnostics and output will still be produced
for non-assimilated datasets.
0
none are to be assimilated when assim_all
is .false. (default)
>0
sets of (stnid, bufr varno,
multi/uni-levels) to be assimilated
assim_varno
Bufr elements of obs sets for assimilation. A
value of 0 implies that all are to be used.
assim_stnid
Stnids of obs sets for assimilation. ‘*’ denote
wild cards
assim_nlev
0
multi-level and uni-level
1
uni_level
>1
multi-level
tropo_mode
Integer indicating if special treatment is to be
given to the troposphere when assimilating total
column measurements. Values indicate
0
No special treatment given (default)
1
Values of the adjoint model above
obsoper%columnBound set to zero. If
specified, generalized innovation
operator only applied below
obsoper%columnBound in the tangent
linear model.
2
Values of tangent linear model and
adjoint model above obsoper%columnBound
set to zero.
Array index refers to BUFR code element of Table
08046 (iconstituentId) identifying the
constituent. Relevant for total column
measurements only.
tropo_bound
Integer indicating which column top value to use
if tropo_mode is non-zero.
0
Use fixed value of tropo_column_top
1
Use model determination of tropopause
2
Use model determination of PBL
Options 1 and 2 will default to the value set
in tropo_column_top if the model derived column
top could not be determined. Relevant for total
column measurements only.
tropo_column_top
Default value to use for the column boundary
(in Pa). Array index refers to BUFR code element
of Table 08046 (iconstituentId) identifying
the constituent. Relevant for total column
measurements only.
obsdata_maxsize
Max allowed size of work arrays (in terms of
number of obs) associated to ordered observation
indices
modelName
Identifier of forecast model
Default: ‘GEM-MACH’
Set to ‘GEM’ for varNames of ‘O3L’, ‘CH4L’, and ‘N2OL’
:operatorSubType(2,i):Operator sub-type name.
Index (2,i) for sub-type to apply for stnid in element (1,i)
See related “obsoper@operatorCategory” automatically assigned
based on obs BUFR element and obsinfo_chm content.
Operator Sub-type name Description
Category
============= =====================================================================
‘Interp’ ‘default’ Piecewise linear interpolation (default)
‘genOper’ Weighted vertical layer averaging - see ‘genOper*’ parameters
Notes:
‘genOper’ requires NAMBCHM namelist parameter settings
getPhysSpaceStats=.true. and
getPhysSpaceHCorrel=.true.
Application of averaging kernels is directed only
by the content of the obsinfo_chm file ‘SECTION III’
storeOperators
Logical indicating if linear operators are stored for re-use in TL and AD calc.
If so, the linear operators will not be re-calculated at different iterations.
Not used when tropo_mode>=1
To return reference model levels or layer boundaries for an observation.
Combination of STNID, element variable number and number of vertical
levels to determine association to the observations. Default values
for top and bottom layers for total column measurements are to be
provided.
To find the averaging kernel for an observation if one is
specified. Returns 0 if either not found or not specified.
Combination of STNID, BUFR element and number of vertical levels
to determine association to the observations.
Arguments
cstnid [character ,in] :: station id
varno [integer ,in] :: BUFR descriptor element
nlev [integer ,in] :: number of levels in the observation
Return
istnid [integer ] :: Index of averaging kernel in oopc_avgkern. Zero if not found.
To identify whether or not the obs set identified by the
combination of (cstnid,varno,nobslev) will be assimilated or
else used for independent verifications after
assimilation/minimization
