Skip to content

GitLab

Commit 0668a26c authored by Espen Sollum's avatar Espen Sollum
Browse files

Updated config file for flexpart/nilu test

parent b87fe6a2
Hide whitespace changes
Inline Side-by-side
examples/config_inv_flexpart_ch4.yml
View file @ 0668a26c
......@@ -27,7 +27,7 @@ verbose : 2
logfile: pyCIF.logtest
# Execution directory
workdir : /home/eso/repos/CIF/flextest/
workdir : &workdir /homevip/espen/FLEX_INV/CIF/flextest/
# Initial date
# Use the following compatible format for the date:
......@@ -53,46 +53,14 @@ model :
name : FLEXPART
version : std
periods : 1MS
# periods : 3MS
nested : yes
run_dir_nest : /home/eso/FLEX_INV/CH4/TEST_OUTPUT/FLEXOUT/GHG/NEST/
run_dir_glob : /home/eso/FLEX_INV/CH4/TEST_OUTPUT/FLEXOUT/GHG/NEST/
run_dir_nest : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/FLEXOUT/GHG/NEST/
run_dir_glob : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/FLEXOUT/GHG/NEST/
coeff : 1.e-3
mmair : 28.97
molarmass : 16.
numscale : 1.E12
####################################################################
####################################################################
# Measurements to account for
# Main keys are:
# - infos: infos on tracers to include in the simulation
# - fic_monitor: standard txt file including all observations
measurements :
plugin:
name: standard
version: std
# File where to save data, if does not exist. Reads from there if exists
# file_monit : /home/eso/repos/CIF/flextest/monitor_allobs.nc
# file_monitor : /home/eso/repos/CIF/flextest/monitor_allobs.nc
dump_type : nc
species :
CH4 :
dailyPM :
provider : WDCGG
format : std
dir_obs : /home/eso/FLEX_INV/CH4/TEST_INPUT/OBS/GHG/*.cn.*
# rescale : false
rescale : true
na_values : -999999.99
default_unit : ppb
default_duration: 1
dump : True
# Minimum measeurement error
measerr : 5.0
#####################################################################
#####################################################################
......@@ -100,19 +68,20 @@ measurements :
# 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)
# - sort_order: (optional) datastore sort order, ['index', 'station'] is default
# 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: same as observations, or fraction???)
# - 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
# file_obsvect : /home/eso/repos/CIF/flextest/monitor_obsvect.nc
# 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
......@@ -124,20 +93,54 @@ obsvect:
version: nc
# parameters:
# CH4 :
dir_netcdf : /home/eso/repos/CIF/flextest/
dir_netcdf : *workdir
# Initial concentrations
dir_initconc : /home/eso/FLEX_INV/CH4/TEST_INPUT/INIT_CONC/GHG/
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.]
cini_err: .01
#####################################################################
#####################################################################
# Measurements to account for
# Main keys are:
# - infos: infos on tracers to include in the simulation
# - fic_monitor: standard txt file including all observations
measurements :
plugin:
name: standard
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 :
CH4 :
dailyPM :
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
default_duration: 1
dump : True
# Minimum measeurement error
measerr : 5.0
####################################################################
####################################################################
# spaces to control/observation spaces
# - read_only: do not run the transport model, read output from existing run
......@@ -166,7 +169,7 @@ domain :
nlat : 90
nlev : 1
type : deg
xmin_glob : -180.
xmin_glob : -179.
ymin_glob : -90.
nlon_glob : 360
nlat_glob : 180
......@@ -232,21 +235,26 @@ controlvect:
plugin:
name: 'FLEXPART'
version: 'nc'
# resol : hpixels
type : flexpart
to_netcdf : True
dir_netcdf : /home/eso/repos/CIF/flextest/
dir_netcdf : *workdir
hresol : regions
inc_ocean : true
fileregions : /home/eso/FLEX_INV/CH4/TEST_OUTPUT/regions_ghg.nc
fileregions : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/regions_ghg.nc
# errtype : max
# errtype : avg
regions_lsm : True
err : 0.5
# Lower limit flux error: unit (kg/m2/h)
flxerr_ll: 1.e-8
# glob_err: 50.
# Total error inversion domain: unit (Tg/y) -- todo: not implemented
glob_err: 10.
numscale : 1.E12
xb_scale : 1.
period : 1D
# period : 10D
dir : /home/eso/FLEX_INV/CH4/TEST_OUTPUT/FLUXES/GHG/
periodflux : 10D
period : 10D
dir : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/FLUXES/GHG/
file : CH4_TOTAL_%Y_05x05.nc
file_glob : CH4_TOTAL_%Y_10x10.nc
......@@ -257,12 +265,16 @@ controlvect:
offsets : True
# offsets : False
# Initial concentrations
dir_initconc : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/INIT_CONC/GHG/
file_initconc: ch4_noaa_%Y%M.nc
varname_init: CH4
hcorrelations :
landsea: True
# landsea: False
filelsm : /home/eso/FLEX_INV/CH4/TEST_INPUT/lsm_0.5x0.5_VERIFY.nc
dircorrel : /home/eso/repos/CIF/flextest/
filelsm : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/lsm_0.5x0.5_VERIFY.nc
dircorrel : *workdir
dump_hcorr : True
sigma_land: 250.
sigma_sea: 1000.
......@@ -274,8 +286,7 @@ controlvect:
# sigma_t : 30
sigma_t : 10
dump_tcorr : True
dircorrel : /home/eso/repos/CIF/flextest/
# evalmin : 1.e-6
evalmin : 1.e-99
dircorrel : *workdir
evalmin : 1.e-5
####################################################################
examples/config_inv_flexpart_ch4_nilu.yml 0 → 100755
View file @ 0668a26c
####################
# pyCIF config file #
#####################
# Define here all parameters for pyCIF following Yaml syntax
# For details on Yaml syntax, please see:
# http://docs.ansible.com/ansible/latest/YAMLSyntax.html
# For non-specified parameters, one can either comment the corresponding line
# or give the value "null"
# Some parameters are switches that need to be activated or not.
# pyCIF expects a boolean True/False
# Yaml accepts the following syntax to be converted to booleans:
# y|Y|yes|Yes|YES|n|N|no|No|NO|true|True|TRUE|false|False|FALSE|on|On|ON|off|Off|OFF
# This offers flexibility, but one should stick to a style for consistency
#####################################################################
# pyCIF parameters
# Verbose level
# 1 = Basic information
# 2 = Debugging
verbose : 2
# Log file (to be saved in $wordkir)
logfile: pyCIF.logtest
# Execution directory
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-02-01
# End date (same as for the initial date)
datef : 2012-02-11 00:00
#####################################################################
#####################################################################
# Transport model
# Accepts any model registered in pyvar.models
# fp_run_dir - directory from which to read FLEXPART footprints
# nested - yes if using nested output domain
# coeff - coefficient to convert from grid units of ppt to ppbv/ppmv (1e-3/1e-6)
# (flexinvert uses ppb, pycif use 1e-6)
# mmair - standard molecular mass of air
# molarmass - Molar mass (in flux files, e.g. CO2-C=12, CH4=16)
# numscale - numeric scaling
model :
plugin:
name : FLEXPART
version : std
periods : 1MS
nested : yes
run_dir_nest : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/FLEXOUT/GHG/NEST/
run_dir_glob : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/FLEXOUT/GHG/NEST/
coeff : 1.e-3
mmair : 28.97
molarmass : 16.
numscale : 1.E12
#####################################################################
#####################################################################
# 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
# Main keys are:
# - infos: infos on tracers to include in the simulation
# - fic_monitor: standard txt file including all observations
measurements :
plugin:
name: standard
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 :
CH4 :
dailyPM :
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
default_duration: 1
dump : True
# Minimum measeurement error
measerr : 5.0
####################################################################
####################################################################
# spaces to control/observation spaces
# - read_only: do not run the transport model, read output from existing run
obsoperator:
plugin:
name: fp
version: std
autorestart: false
read_only: true
####################################################################
#####################################################################
# Domain definition
# coordinates refers to edge of inversion domain
# For FLEXPART, include definition of global domain
domain :
plugin :
name : FLEXPART
version : std
xmin : -15.
xmax : 35.
ymin : 30.
ymax : 75.
nlon : 100
nlat : 90
nlev : 1
type : deg
xmin_glob : -179.
ymin_glob : -90.
nlon_glob : 360
nlat_glob : 180
dx_glob : 1.
dy_glob : 1.
heights : [100., 500., 1000., 2000., 3000., 5000., 7000., 9000.,
12000., 15000., 20000., 50000.]
mode:
plugin:
name: 4dvar
version: std
minimizer:
plugin:
# name: congrad
name: M1QN3
version: std
simulator:
plugin:
name: gausscost
version: flexpart
reload_from_previous: False
maxiter: 10
df1: 1.0
# df1: 0.01
nsim: 100
m: 1
save_out_netcdf: True
####################################################################
#####################################################################
# Arguments to define the state vector
# These are specific to the state vector and inversion method you chose
# 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)
# - sigma_sea: spatial correlation length for prior errors over ocean (km)
# - tracers: list of tracers to put in the state vector (with definition arguments):
# - calcstd: calculate global standard deviation
# - resol: resolution at which fields are scaled
# (choice = bands,regions,pixels;
# if regions, provide a netcdf file ficregion
# 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)
# - flxerr_ll: lower limit flux error: unit (kg/m2/h)
# - glob_err (optional): uncertainty on global budget. Units (Tg/y).
controlvect:
plugin:
name: standard
version: std
components:
fluxes:
parameters:
CH4 :
plugin:
name: 'FLEXPART'
version: 'nc'
# resol : hpixels
type : flexpart
to_netcdf : True
dir_netcdf : *workdir
hresol : regions
inc_ocean : true
fileregions : /homevip/espen/FLEX_INV/CH4/TEST_OUTPUT/regions_ghg.nc
# errtype : max
# errtype : avg
regions_lsm : True
err : 0.5
# Lower limit flux error: unit (kg/m2/h)
flxerr_ll: 1.e-8
# Total error inversion domain: unit (Tg/y) -- todo: not implemented
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
file_glob : CH4_TOTAL_%Y_10x10.nc
varname_flx: emisch4
lonname_flx: longitude
latname_flx: latitude
timename_flx: time
offsets : True
# offsets : False
# Initial concentrations
dir_initconc : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/INIT_CONC/GHG/
file_initconc: ch4_noaa_%Y%M.nc
varname_init: CH4
hcorrelations :
landsea: True
# landsea: False
filelsm : /homevip/espen/FLEX_INV/CH4/TEST_INPUT/lsm_0.5x0.5_VERIFY.nc
dircorrel : *workdir
dump_hcorr : True
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 : 10
dump_tcorr : True
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