"implement option and code for running a sensitivity analysis of N2O GC simulations (MR !12)
issue
- we need simulations with quasi-randomised selection of sampling times, i.e., plusminus 30 sec around pre-selected sample times
- this may help providing an estimate on uncertainties associated with the GC simulation
solution and changes
- add new parameter and options to
preprocessing/additional_info/simulate_gc_n2o_data.txt
-
simGCsampleTimeSensiAnDefSec
- default sample times for sensitivity analysis -
modifySimGCsampleTimeSec
- the range in seconds the default sample timessimGCsampleTimeSensiAnDefSec
can be modified
-
- add to
preprocessing/run_all_main.R
the new optionsimulateGCdataSensiAn
- to run sensitivity analysis for "using the PGA data to simulate GC data"
- in
preprocessing/run_all_main.R
the new varrandMod_simGCsampleTimeSensiAnDefSec
is created holding the changes of the simGC sample times- this var is passed to the
f_main_data_preprocessing()
and used there in section "7. Sample N2O gas concentration data for GC simulation" ifsimulateGCdataSensiAn == "T"
- this var is passed to the
- in this section 7 if
simulateGCdataSensiAn == "T"
:- the new simGC sample times are calculated and stored in the var
simGCsampleTimeSensiAnModSec
- this var is passed as argument to
f_sample_gas_concentrations_for_gc_simulation_per_mp()
where they are used for simGC calculations
- this var is passed as argument to
- a logfile is created and data is written to it, i.e., modifications of original sample times and the new sample times
- the new simGC sample times are calculated and stored in the var
note
- for the sensitivity analysis the script needs to be run multiple times (manually) and the output data form the script need to be moved/stored & analysed manually
- i.e., there is no script that does this automatically
testing
- preprocessing is running fine with meas campaign 202307
- output and logfile show that the modified sample times are used
Closes #23 (closed)
Edited by Jan Engel