[FemEngine] [MODULE_SELECTION] // other modules... mode_of_data_assimilation_module = 2 EndSect // MODULE_SELECTION [DATA_ASSIMILATION_MODULE] // https://dhi-developer-documentation.azurewebsites.net/mike_fm_da/overview/ [TIME] start_time_step = 0 start_time_step_assimilation = 1200 time_step_factor_assimilation = 5 EndSect // TIME [METHOD] type = 1 // 0=Free, 1=EnKF (ensemble), 2=Simple, 3=Steady, 4=EnOI ensemble_size = 10 algorithm = 1 // (for type=1) 1=serialESRF, 2=DEnKF, 3=ETKF Rfactor_anomalies = 2.0 // inflation only where model was updated (e.g. 2.0) Rfactor_all = 1.0 // factor on st_dev for *all* measurements EndSect // METHOD [MODEL_ERROR_MODEL] use_clock_based_seeds = false // false=same seq. of random numbers every time number_of_model_errors = 2 [MODEL_ERROR_1] type = 'wlbc' n_bounds = 2 bound_codes = 2, 3 [Error_Formulation] st_dev = 0.05 // 5cm perturbation_type= 'additive' propagation_type = 'AR(1)' propagation_parameter = 3600 initialization_type = 0 horizontal_discretization_type = 1 // 0: constant, 1: (piecewise) linear, 2: equidistant grid horizontal_corr_function = 1 // 1: gaussian, 2: exponential horizontal_corr = 70000 // correlation length in meters EndSect EndSect // MODEL_ERROR_1 [MODEL_ERROR_2] type = 'wind' [Error_Formulation] st_dev = 0.1 // 10% since multiplicative perturbation_type = 'multiplicative' propagation_type = 'AR(1)' propagation_parameter = 3600 initialization_type = 0 horizontal_discretization_type = 2 // 0: constant, 1: (piecewise) linear, 2: equidistant grid horizontal_grid_spacing = 40000 horizontal_corr_function = 1 // 1: gaussian, 2: exponential horizontal_corr = 120000 // correlation length in meters EndSect EndSect // MODEL_ERROR_2 EndSect // MODEL_ERROR_MODEL [MEASUREMENTS] number_of_independent_measurements = 4 [MEASUREMENT_1] include = 1 // 0: inactive, 1: active used for DA, 2: active but only for validation name = 'ts_WL_NW' measured_variable = 'water level' // 'water level' (default), 'u', 'v', 'ua', 'va' (3d), 'temperature' type = 1 // 1=fixed station (dfs0); 2=track (dfs0) file_name = |.\..\Measurements\ts_WL_NW.dfs0| item_number = 1 data_offset = 0.0 // add this amount to all data in file (default=0.0) type_time_interpolation = 1 // 0: discrete/no interp, 1: piecewise linear, 2: cubic spline // time_window_in_seconds = 300 // in case type_time_interpolation=0 or type=2 (use AD time step size) position = 251, 752 coordinate_type = 'NON-UTM' // default: same as model group = 1 // measurements may belong to group (used by localization) st_dev = 0.01 type_of_space_interpolation = 1 // 0=constant e.g. no interpolation, 1=linear reconstruction by gradient (default), 2=inverse distance weighing to neighbors neighbors // error_temporal_corr = 1 // if type>1: correlated measurement errors along track [in seconds] // error_horizontal_corr = 1 // if type>1: correlated measurement errors along track [in meters] EndSect // MEASUREMENT_1 [MEASUREMENT_2] name = 'ts_WL_SW' file_name = |.\..\Measurements\ts_WL_SW.dfs0| position = 251, 252 group = 1 st_dev = 0.01 EndSect // MEASUREMENT_2 [MEASUREMENT_3] include = 2 name = 'ts_WL_SE' file_name = |.\..\Measurements\ts_WL_SE.dfs0| position = 751, 252 group = 1 st_dev = 0.01 EndSect // MEASUREMENT_3 [MEASUREMENT_4] name = 'ts_WL_NE' file_name = |.\..\Measurements\ts_WL_NE.dfs0| position = 751, 752 group = 1 st_dev = 0.01 EndSect // MEASUREMENT_4 EndSect // MEASUREMENTS [ESTIMATOR] type = 1 [REGULARIZATION] use_temporal_smoothing = false smoothing_halftime = 7200 use_localization = false [LOCALIZATION] number_of_groups = 1 [GROUP_1] horizontal_corr_function = 3 // 1:gauss, 2:exp, 3:Gasperi&Cohn(default) horizontal_corr = 500 // correlation length in meters EndSect // GROUP_1 EndSect // LOCALIZATION EndSect // REGULARIZATION EndSect // ESTIMATOR [ERRCOV_IO] // only for DA%algorithm=1 [INPUT] include = 0 EndSect // INPUT [OUTPUTS] number_of_outputs = 0 [OUTPUT_1] include = 1 first_time_step = 0 last_time_step = 360 time_step_frequency = 360 // like time_step_factor time_average = true // true=only a single step output smoothed_out = false // true=output time-smoothed error_covariance (instead of normal error_cov) file_name_area = 'ErrCovIO_avg_Area.dfsu' file_name_volume = 'ErrCovIO_avg_Volume.dfsu' file_name_err01_wlbc = 'ErrCovIO_avg_wlbc_err.dfs1' file_name_err02_wind_u = 'ErrCovIO_avg_Wind_err_u.dfs2' file_name_err02_wind_v = 'ErrCovIO_avg_Wind_err_v.dfs2' EndSect // OUTPUT_1 EndSect // OUTPUTS EndSect // ERRCOV_IO [ENSEMBLE_IO] [INPUT] // can be used for hot-in for ensemble run type = 0 // 0: not active, 1: Meanstate, 2: Ensemble file_name_area = |.\DA\State_Area.dfsu| file_name_volume = |.\DA\State_Volume.dfsu| file_name_err01_wlbc = |.\DA\State_wlbc_err.dfs1| file_name_err02_wind_u = |.\DA\State_Wind_err_u.dfs2| file_name_err02_wind_v = |.\DA\State_Wind_err_v.dfs2| EndSect // INPUT [OUTPUT] // can be used for hot-out for ensemble run type = 2 // 0: not active, 1: Meanstate, 2: Ensemble modelerror_only = false file_name_area = 'State_Area_out.dfsu' file_name_volume = 'State_Volume_out.dfsu' file_name_err01_wlbc = 'State_wlbc_err.dfs1' file_name_err02_wind_u = 'State_Wind_err_u.dfs2' file_name_err02_wind_v = 'State_Wind_err_v.dfs2' first_time_step = 0 last_time_step = 3600 time_step_frequency = 100 EndSect // OUTPUT EndSect // ENSEMBLE_IO [DIAGNOSTICS] [OUTPUTS] number_of_outputs = 4 [OUTPUT_1] type = 2 // 1: measurement point, 2: validation point position = 500, 500 coordinate_type = 'NON-UTM' variable_name = 'water level' file_name = 'Diagnostics_TS2.dfs0' EndSect // OUTPUT_1 [OUTPUT_2] type = 2 position = 1, 500 coordinate_type = 'NON-UTM' variable_name = 'water level' file_name = 'Diagnostics_TS_x1_y500.dfs0' EndSect // OUTPUT_2 [OUTPUT_3] include = 1 type = 1 measurement_id = 2 file_name = 'Diagnostics_measTS2_SW.dfs0' EndSect // OUTPUT_3 [OUTPUT_4] type = 1 measurement_id = 1 file_name = 'Diagnostics_measTS2_NW.dfs0' EndSect // OUTPUT_4 EndSect // OUTPUTS EndSect // DIAGNOSTICS EndSect // DATA_ASSIMILATION_MODULE EndSect // FemEngine