Commit 4ac42c7f authored by Espen Sollum's avatar Espen Sollum
Browse files

Updated config file

parent 9ac663d4
......@@ -27,15 +27,15 @@ verbose : 2
logfile: pyCIF.logtest
# Execution directory
workdir : /homevip/espen/FLEX_INV/CIF/flextest/
workdir : &workdir /homevip/espen/FLEX_INV/CIF/flextest/
# Initial date
# Use the following compatible format for the date:
# YYYY-mm-dd or YYYY-mm-dd HH:mm
datei : 2012-01-01
datei : 2012-02-01
# End date (same as for the initial date)
datef : 2012-01-31
datef : 2012-02-11 00:00
#####################################################################
#####################################################################
# Transport model
......@@ -67,12 +67,47 @@ model :
# How to build your observation vector and observation uncertainties if needed
# Also projects information from the observation to the model space
# - fic_obsvect: observation vector from previous simulations
# - dump_debug: write out extra information (for debugging)
# For FLEXPART the background plugin takes the following parameters
# - dir_initconc: directory containing concentration files
# - file_initconc: name of concentration files
# - varname_initconc: name of fata variable in the netCDF file
# - optimize_cini: optimize initial concentration or not
# - cini_lat: northern edges of latitude bands for cini optimization
# - cini_err: Initial concentration error (unit: fraction). Default 28.
# - cini_res : Time resolution for initial mixing ratio scalars (days)
obsvect:
plugin:
name: standard
version: std
# fic_obsvect : /homevip/espen/FLEX_INV/CIF/flextest/monitor_obsvect.nc
# file_obsvect : /homevip/espen/FLEX_INV/CIF/flextest/monitor_obsvect.nc
dump: True
dump_debug: True
dump_type: nc
sort_order: ['station', 'index']
background:
plugin:
name: FLEXPART
version: nc
# parameters:
# CH4 :
dir_netcdf : *workdir
# Initial concentrations
dir_initconc : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/INIT_CONC/GHG/
file_initconc: ch4_noaa_%Y%m.nc
varname_conc: CH4
varname_pressure: pressure
optimize_cini: True
# Time resolution for initial mixing ratio scalars (days)
cini_res : 28.
cini_err: 0.03
cini_lat: [-30., 0., 30., 90.]
#####################################################################
#####################################################################
# Measurements to account for
......@@ -85,6 +120,8 @@ measurements :
version: std
# File where to save data, if does not exist. Reads from there if exists
# fic_monitor : /homevip/espen/FLEX_INV/CIF/flextest/monitor_allobs.nc
# file_monit : /homevip/espen/FLEX_INV/CIF//monitor_allobs.nc
# file_monitor : /homevip/espen/FLEX_INV/CIF/monitor_allobs.nc
dump_type : nc
species :
......@@ -93,6 +130,7 @@ measurements :
provider : WDCGG
format : std
dir_obs : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/OBS/GHG/*.cn.*
# rescale : false
rescale : true
na_values : -999999.99
default_unit : ppb
......@@ -137,6 +175,8 @@ domain :
nlat_glob : 180
dx_glob : 1.
dy_glob : 1.
heights : [100., 500., 1000., 2000., 3000., 5000., 7000., 9000.,
12000., 15000., 20000., 50000.]
......@@ -155,7 +195,10 @@ mode:
version: flexpart
reload_from_previous: False
maxiter: 10
df1: 0.01
df1: 1.0
# df1: 0.01
nsim: 100
m: 1
save_out_netcdf: True
####################################################################
......@@ -165,6 +208,8 @@ mode:
# Please refer to the documentation to know what to include here
# For the standard FLEXPART, include the following:
# - ficlsm: land-sea mask (must be consistent with LMDZ grid)
# - regions_lsm : if true, use land/sea mask encoded in fileregions (sea: < 0, land: > 0)
# This overrides usage of lsm file, if this is also provided
# - correl: Use correlation or not in B
# - dircorrel: path to pre-computed correlations
# - sigma_land: spatial correlation length for prior errors over land (km)
......@@ -177,10 +222,8 @@ mode:
# if bands, define a list of latitudes as band limits (n+1 for n bands)
# - inc_ocean: for 'regions', optimize ocean boxes (default true)
# - periodflux: period of variation for increments within a month (days)
# - glob_err (optional) = uncertainty on global budget
# - flxerr_ll = Lower limit flux error (kg/m2/h)
# - flxerr_ll: lower limit flux error: unit (kg/m2/h)
# - glob_err (optional): uncertainty on global budget. Units (Tg/y).
controlvect:
plugin:
name: standard
......@@ -192,13 +235,16 @@ controlvect:
plugin:
name: 'FLEXPART'
version: 'nc'
# resol : hpixels
type : flexpart
to_netcdf : True
dir_netcdf : /homevip/espen/FLEX_INV/CIF/flextest/
dir_netcdf : *workdir
hresol : regions
inc_ocean : true
fileregions : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/regions_ghg.nc
errtype : avg
# errtype : max
# errtype : avg
regions_lsm : True
err : 0.5
# Lower limit flux error: unit (kg/m2/h)
flxerr_ll: 1.e-8
......@@ -206,6 +252,7 @@ controlvect:
glob_err: 10.
numscale : 1.E12
xb_scale : 1.
periodflux : 10D
period : 10D
dir : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/FLUXES/GHG/
file : CH4_TOTAL_%Y_05x05.nc
......@@ -227,14 +274,19 @@ controlvect:
landsea: True
# landsea: False
filelsm : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/lsm_0.5x0.5_VERIFY.nc
dircorrel : /homevip/espen/FLEX_INV/CIF/flextest/
dircorrel : *workdir
dump_hcorr : True
sigma_land: 250
sigma_sea: 1000
sigma_land: 250.
sigma_sea: 1000.
# sigma_land: 1.
# sigma_sea: 1.
evalmin : 1.e-6
# evalmin : 1.e-15
tcorrelations :
sigma_t : 30
# sigma_t : 30
sigma_t : 10
dump_tcorr : True
dircorrel : /homevip/espen/FLEX_INV/CIF/flextest/
dircorrel : *workdir
evalmin : 1.e-5
####################################################################
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment