[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 = 100
         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       
      EndSect  // METHOD
      
      [MODEL_ERROR_MODEL]         
        use_clock_based_seeds = False 
        number_of_model_errors = 1
        
        [MODEL_ERROR_1]
            include = 1
            type = 'adbc'
            fraction_id = 1
            n_bounds = 1 
            bound_codes = 3

            [Error_Formulation]	
                st_dev = 0.05	
                perturbation_type = 'additive'
                propagation_type = 'persistence' // 'AR(1)'
                propagation_parameter = 86400
                initialization_type = 1
            EndSect
        EndSect  // MODEL_ERROR_1
        
        [MODEL_ERROR_2]
          type = 'adsource'
          source_id = 1       // HD Source number
          fraction_id = 1      // = MT%SSC_FRACTION_1
          
          [ERROR_FORMULATION]
                st_dev = 0.1
                perturbation_type = 'multiplicative_non_negative'
                propagation_type = 'AR(1)'
                initialization_type = 1
          EndSect  // ERROR_FORMULATION
        EndSect  // MODEL_ERROR_2
        
        [MODEL_ERROR_3]
          type = 'MT_parameter'
          parameter_name = 'settling_velocity_coefficient'		  
          fraction_id = 2      // = MT%SSC_FRACTION_2
          
          [ERROR_FORMULATION]
                st_dev = 5
                perturbation_type = 'additive_non_negative'
                propagation_type = 'persistence'
                initialization_type = 1
          EndSect  // ERROR_FORMULATION
        EndSect  // MODEL_ERROR_3

      EndSect  // MODEL_ERROR_MODEL

      [MODEL_STATE_SPACE]      
         hd_state_space_type = 0 
         mt_state_space_type = 1      
      EndSect  // MODEL_STATE_SPACE

      [MEASUREMENTS]
         number_of_independent_measurements = 1

         [MEASUREMENT_1]
            name = 'SSC1'
            measured_variable = 'MT%SSC_FRACTION_1'
            type = 1                  
            file_name = |obs\obs_SSC1.dfs0|
            item_number = 1
            position = 3000, 3400
            st_dev = 0.01
            max_mean_innovation = 1
         EndSect  // MEASUREMENT_1

      EndSect  // MEASUREMENTS

      [ESTIMATOR]
         [REGULARIZATION]
            use_temporal_smoothing = false
            use_localization = false
         EndSect  // REGULARIZATION
      EndSect  // ESTIMATOR

      [ERRCOV_IO]     
        [INPUT]
            include = 0 
        EndSect  // INPUT
        
        [OUTPUTS]
          number_of_outputs = 0         
        EndSect  // OUTPUTS 
      EndSect  // ERRCOV_IO

      [ENSEMBLE_IO]     
         [INPUT]
            type = 0  // 0: not active, 1: Meanstate, 2: Ensemble   
         EndSect  // INPUT

         [OUTPUT]
            type = 2  // 0: not active, 1: Meanstate, 2: Ensemble			
            file_name_volume = 'State_Volume.dfsu'			
         EndSect  // OUTPUT
      EndSect  // ENSEMBLE_IO

      [DIAGNOSTICS]
        [OUTPUTS]
            number_of_outputs = 2
             
             [OUTPUT_1]
                type = 2
                position = 3500, 3500
                variable_name = 'MT%SSC_FRACTION_1'
                file_name = 'diagnostics_SSC1.dfs0'				
             EndSect  // OUTPUT_1
                          
             [OUTPUT_2]
                type = 2
                position = 3500, 3500
                variable_name = 'MT%SSC_FRACTION_2'
                file_name = 'diagnostics_SSC2.dfs0'				
             EndSect  // OUTPUT_2
             
        EndSect  // OUTPUTS
      EndSect  // DIAGNOSTICS	
      
   EndSect  // DATA_ASSIMILATION_MODULE   
   
EndSect  // FemEngine