[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)				
      EndSect  // METHOD
      
      [MODEL_ERROR_MODEL]         
        use_clock_based_seeds = false // false=same seq. of random numbers every time
        number_of_model_errors = 1
                
        [MODEL_ERROR_1]
            type = 'windspeed'
            [Error_Formulation]	
                st_dev = 1.5
                perturbation_type= 'additive'
                propagation_type = 'AR(1)'
                propagation_parameter = 7200
                initialization_type = 1
                
                horizontal_discretization_type = 2
                horizontal_grid_spacing = 300000.0
                horizontal_corr_function = 1
                horizontal_corr = 2000000.0
            EndSect
        EndSect  // MODEL_ERROR_1
        
        [MODEL_ERROR_2]
            type = 'SW_parameter'
            parameter_name = 'background_Charnock_parameter'  // global constant
            [Error_Formulation]	
                st_dev = 0.05
                perturbation_type= 'additive'
                propagation_type = 'AR(1)'
                propagation_parameter = 360000
                initialization_type = 1
            EndSect
        EndSect  // MODEL_ERROR_2
        
        [MODEL_ERROR_3]
            type = 'SW_parameter'
            parameter_name = 'Nikuradse_roughness'
            [Error_Formulation]	
                st_dev = 0.05
                perturbation_type= 'additive'
                propagation_type = 'AR(1)'
                propagation_parameter = 360000
                initialization_type = 1
                
                horizontal_discretization_type = 2
                horizontal_grid_spacing = 1000000.0
                horizontal_corr_function = 1
                horizontal_corr = 10000000.0
            EndSect
        EndSect  // MODEL_ERROR_3
        
      EndSect  // MODEL_ERROR_MODEL

      [MODEL_STATE_SPACE]  
        sw_state_space_type = 1	
        type_Hm0 = 1
        type_Tp = 0
        type_windspeed = 0
      EndSect  // MODEL_STATE_SPACE

      [MEASUREMENTS]
         number_of_independent_measurements = 2

         [MEASUREMENT_1]
            include = 2              // 0: inactive, 1: active used for DA, 2: active but only for validation
            name = 'Hm0_1'
            measured_variable = 'SW%Hm0' 
            
            type = 1                      // 1=fixed station (dfs0); 2=track (dfs0)
            file_name = |.\..\Measurements\Hm0_1.dfs0|			
            item_number = 1 
            
            position = 3.276, 51.999
            
            st_dev = 0.1
         EndSect  // MEASUREMENT_1

         [MEASUREMENT_2]
            include = 2            
            name = 'Hm0_2'
            measured_variable = 'SW%Hm0' 
            file_name = |.\..\Measurements\Hm0_2.dfs0|
            item_number = 1
            position = 1.476, 51.218
            
            st_dev = 0.1
         EndSect  // MEASUREMENT_2	 
         
      EndSect  // MEASUREMENTS

      [ESTIMATOR]
         type = 1
         [REGULARIZATION]
            use_temporal_smoothing = false
            smoothing_halftime = 3600
            use_localization = false
            [LOCALIZATION]                 
               number_of_groups = 1
               [GROUP_1]
                  horizontal_corr_function = 3  // 1:gauss, 2:exp, 3:Gasperi&Cohn(default) 
                  horizontal_corr = 50000		// 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 = 1200
                    last_time_step = 3600
                    time_step_frequency = 2400 // make one time averaged output
                    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_spectrum = 'ErrCovIO_avg_Spectrum.dfsu'
                    file_name_err01_windspeed = 'ErrCovIO_avg_Windspeed_err.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_spectrum = |.\DA\State_Spectrum.dfsu|
            file_name_err01_windspeed = |.\DA\State_Windspeed_err.dfs2|
            
         EndSect  // INPUT

         [OUTPUT]  // can be used for hot-out for ensemble run
            type = 2  // 0: not active, 1: Meanstate, 2: Ensemble
            include_model_variables = true
            include_model_errors = true
            
            file_name_area = 'State_Area_out.dfsu'            
            file_name_spectrum = ''
            file_name_err01_windspeed = 'State_windspeed_err.dfs2'         
            
            first_time_step = 0
            // last_time_step = 3600    // till the end
            time_step_frequency = 100
         EndSect  // OUTPUT
      EndSect  // ENSEMBLE_IO

      [DIAGNOSTICS]
        [OUTPUTS]
            number_of_outputs = 2
             
             [OUTPUT_1]
                type = 2    // 1: measurement point, 2: validation point
                position = 1.476, 51.218		
                variable_name = 'SW%Hm0'
                file_name = 'Diagnostics_x_1p476_y_51p218.dfs0'
             EndSect  // OUTPUT_1
             
             [OUTPUT_2]
                type = 1
                measurement_id = 1
                file_name = 'Diagnostics_meas1.dfs0'
             EndSect  // OUTPUT_2
             
        EndSect  // OUTPUTS
      EndSect  // DIAGNOSTICS	
      
   EndSect  // DATA_ASSIMILATION_MODULE   
   
EndSect  // FemEngine