SETTINGS_co2_config 3.29 KiB
# ==================================================
#
# FLEXINVERT CONFIGURATION SETTINGS
#
# comment lines start with '#'
# each parameter line starts with 'parameter name:'
#
# ==================================================
# Run mode:
# 0 = run forward model
# 1 = run optimization
# 2 = randomly perturb for MC
run_mode: 1
# Random number seed
# (only used if run_mode = 2)
seed: 100
# Run dates: format yyyymmdd
# datei = start date
# datef = end date
datei: 20120101
datef: 20120131
# Inversion method ('analytic', 'congrad' or 'm1qn3')
method: m1qn3
# Average/select flexpart releases (true or false)
average_fp: .false.
# Number of iterations
# only used if method is 'congrad' or 'm1qn3'
maxiter: 2
# Include ocean boxes in inversion (true or false)
# currently only for ghg species
inc_ocean: .false.
# Optimize initial mixing ratios (true or false)
opt_cini: .true.
# Use spatial correlation in error covariance matrix (true or false)
# if use regions based ecosystems then should be false
spa_corr: .true.
# Use best guess estimate from previous inversion (congrad only)
# 0 = no best guess
# 1 = use best guess (file must be specified in SETTINGS_files)
prior_bg: 0
# Restart a crashed run
# for congrad/m1qn3 will pick-up from last iteration
# for analytic will only use pre-calculated covariance matrix and boundary conditions
# 0 = new run
# 1 = restart crashed run
restart: 1
# Verbose output
# only use for debugging small runs
verbose: .true.
# Species ("co2" or "ghg")
spec: co2
# Molar mass (in flux files, e.g. C=12, CH4=16)
molarmass: 12.
# Coefficient to convert from grid units of ppt to observation unit (e.g. ppm)
coeff: 1.e-6
# Use nested flexpart output (true or false)
nested: .false.
# Inversion domain:
# if using nested output it must match the inversion domain
# w_edge_lon = lon of western edge of inversion grid
# s_edge_lat = lat of southern edge inversion grid
# e_edge_lon = lon of eastern edge of inversion grid
# n_edge_lat = lat of northern inversion grid
# xres = longitudinal grid resolution
# yres = latitudinal grid resolution
w_edge_lon: -15
s_edge_lat: 30
e_edge_lon: 35
n_edge_lat: 75
xres: 1.0
yres: 1.0
# Spatial aggregation of grid (true or false)
regions: .false.
# State vector time resolution:
# time resolution at which NEE fluxes are optimized
# statres determines averaging interval over 1 or more days (given in days)
# statres_hr determines time resolution within 1 day (given in hours)
# example: statres = 5, statres_hr = 6
# gives 4 time intervals per day that are averaged over 5 days
statres: 5
statres_hr: 12
# non-NEE prior fluxes:
# number of flux fields per year in input
nt_flx: 12
# NEE prior fluxes:
# num_nee_day = number of NEE values per 24h
num_nee_day: 8
# Measurement error: unit same as obs
# used if error in obs input <= zero
measerr: 0.5
# Scalar of initial conc error: fraction
cinierr: 0.01
# Prior flux error: fraction
flxerr: 0.25
# Fossil fuel error: fraction
ffferr: 0.05
# Spatial correlation length for land: unit (km)
sigma_land: 20.
# Spatial correlation length for ocean: unit (km)
sigma_ocean: 1000.
# Temporal correlation length: unit (days)
sigmatime: 7.
# Total error for domain (Tg/y)
globerr: 300.