[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
      EndSect  // TIME

      [METHOD]
          type = 1
          ensemble_size = 5
          algorithm = 1
      EndSect  // METHOD
      
      [MODEL_ERROR_MODEL]                 
        use_clock_based_seeds = False 
        number_of_model_errors = 1
        
        [MODEL_ERROR_1]
            type = 'TR_parameter'
            component_id = 2
            [Error_Formulation]	
                st_dev = 2e-6	
                propagation_type = 'AR(1)'  // 'AR(1)' or 'persistence'				
                propagation_parameter = 3600
                initialization_type = 1
            EndSect // Error_Formulation
            
        EndSect  // TRPAR_ERROR

        [MODEL_ERROR_2]
          type = 'adsource'
          module_name = 'TR'
          source_id = 1       // HD Source number
          component_id = 1      // = COMPONENT_1
          
          [ERROR_FORMULATION]
                st_dev = 0.1
                perturbation_type = 'multiplicative_non_negative'
                propagation_type = 'persistence'  // 'AR(1)' or 'persistence'								
                initialization_type = 1
          EndSect  // ERROR_FORMULATION
        EndSect  // MODEL_ERROR_2
        
      EndSect  // MODEL_ERROR_MODEL

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

      [MEASUREMENTS]
         number_of_independent_measurements = 1

        [MEASUREMENT_1]
            include = 1
            name = 'ts_WL'
            measured_variable = 'TR%COMPONENT_2'
            type = 1                  
            file_name = |.\obs\tracer.dfs0|
            item_number = 1            
            position = 606300, 6905800
            st_dev = 0.005      
            group = 1
         EndSect  // MEASUREMENT_1
      EndSect  // MEASUREMENTS

      [ESTIMATOR]
         [REGULARIZATION]
            use_temporal_smoothing = false
            smoothing_halftime = 600
            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 = 2  // 0: not active, 1: Meanstate, 2: Ensemble
            file_name_area = |.\Initial\State_Area_input.dfsu|     
            file_name_err01_trpar	= |.\Inital\TRpar_input.dfs0|		
         EndSect  // INPUT

         [OUTPUT]
            type = 2  // 0: not active, 1: Meanstate, 2: Ensemble
            include_model_variables = true
            include_model_errors = true
            
            file_name_area = 'DAState_Area_out.dfsu'         
            file_name_err01_trpar = 'DATRpar_out.dfs0'
            
            first_time_step = 480
            last_time_step = 480
            time_step_frequency = 1
         EndSect  // OUTPUT
      EndSect  // ENSEMBLE_IO

      [DIAGNOSTICS]
        [OUTPUTS]
            number_of_outputs = 2
             
             [OUTPUT_1]
                type = 2
                position = 606300, 6905800
                variable_name = 'TR%COMPONENT_2'
                file_name = 'Diagnostics_TS.dfs0'
             EndSect  // OUTPUT_1
             
             [OUTPUT_2]
                type = 1
                measurement_id = 1
                variable_name = 'TR%COMPONENT_2'
                file_name = 'Diagnostics_TS_Meas.dfs0'			
             EndSect  // OUTPUT_1
             
        EndSect  // OUTPUTS
      EndSect  // DIAGNOSTICS		
      
   EndSect  // DATA_ASSIMILATION_MODULE	
   
EndSect  // FemEngine