[Dart-dev] DART/branches Revision: 12729

[dart at ucar.edu]

jla at ucar.edu
2018-07-11 11:49:08 -0600 (Wed, 11 Jul 2018)
136
Modifications to model, input, and namelist to serve as base case for exploration of error distribution on 
concentration observations.



Modified: DART/branches/concentration_observations/models/simple_advection/model_mod.f90
===================================================================
--- DART/branches/concentration_observations/models/simple_advection/model_mod.f90	2018-07-11 15:51:51 UTC (rev 12728)
+++ DART/branches/concentration_observations/models/simple_advection/model_mod.f90	2018-07-11 17:49:08 UTC (rev 12729)
@@ -145,7 +145,7 @@
                      lagrangian_for_wind, destruction_rate, &
                      source_random_amp_frac, source_damping_rate, &
                      source_diurnal_rel_amp, source_phase_noise, &
-                     template_file, my_ens_size
+                     template_file
 
 
 ! Define the location of the state variables in module storage
@@ -196,7 +196,7 @@
    source_random_amp(num_grid_points))
 
 ! Set the base value for the mean source; This case has a single enhanced source
-mean_source    = 0.1_r8
+mean_source    = 0.0_r8
 mean_source(1) = 1.0_r8
 
 ! A rapidly varying source
@@ -653,10 +653,11 @@
    where(state_ens_handle%vars(i,:) < 0.0_r8) state_ens_handle%vars(i,:) = 0.0_r8 
 
    ! Perturb the source
+   ! Not perturbed for Elia's initial exploration
    do j=1,state_ens_handle%my_num_copies
-      state_ens_handle%vars(num_grid_points + i,j) = random_gaussian(random_seq, &
-                                                     state_ens_handle%vars(num_grid_points + i,j), &
-                                                     state_ens_handle%vars(num_grid_points + i,j))
+      !!!state_ens_handle%vars(num_grid_points + i,j) = random_gaussian(random_seq, &
+                                                     !!!state_ens_handle%vars(num_grid_points + i,j), &
+                                                     !!!state_ens_handle%vars(num_grid_points + i,j))
    enddo
    where(state_ens_handle%vars(num_grid_points + i,:) < 0.0_r8) state_ens_handle%vars(num_grid_points + i,:) = 0.0_r8 
 

Modified: DART/branches/concentration_observations/models/simple_advection/work/filter_input_list.txt
===================================================================
--- DART/branches/concentration_observations/models/simple_advection/work/filter_input_list.txt	2018-07-11 15:51:51 UTC (rev 12728)
+++ DART/branches/concentration_observations/models/simple_advection/work/filter_input_list.txt	2018-07-11 17:49:08 UTC (rev 12729)
@@ -1 +1 @@
-filter_input.nc
+perfect_input.nc

Modified: DART/branches/concentration_observations/models/simple_advection/work/input.nml
===================================================================
--- DART/branches/concentration_observations/models/simple_advection/work/input.nml	2018-07-11 15:51:51 UTC (rev 12728)
+++ DART/branches/concentration_observations/models/simple_advection/work/input.nml	2018-07-11 17:49:08 UTC (rev 12729)
@@ -1,5 +1,5 @@
 &perfect_model_obs_nml
-   read_input_state_from_file = .true.
+   read_input_state_from_file = .false.
    single_file_in             = .true.
    input_state_files          = "perfect_input.nc"
 
@@ -49,7 +49,7 @@
 
    ens_size                     = 80
    num_groups                   = 1
-   perturb_from_single_instance = .false.
+   perturb_from_single_instance = .true.
    perturbation_amplitude       = 0.2
    distributed_state            = .true.
 
@@ -69,7 +69,7 @@
    inf_initial_from_restart    = .false.,                 .false.
    inf_sd_initial_from_restart = .false.,                 .false.
    inf_deterministic           = .true.,                  .true.
-   inf_initial                 = 1.0,                     1.2
+   inf_initial                 = 1.00,                     1.2
    inf_lower_bound             = 1.0,                     1.0
    inf_upper_bound             = 1000000.0,               1000000.0
    inf_damping                 = 1.0,                     1.0
@@ -99,7 +99,7 @@
 # doing a successful assimilation, change cutoff to 0.02 and rerun.
 &assim_tools_nml
    filter_kind                     = 1
-   cutoff                          = 100000000.0
+   cutoff                          = 0.2
    sort_obs_inc                    = .false.
    spread_restoration              = .false.
    sampling_error_correction       = .false.
@@ -128,21 +128,21 @@
    /
 
 &obs_kind_nml
-   assimilate_these_obs_types = 'TRACER_CONCENTRATION',
-                                'VELOCITY'
-   evaluate_these_obs_types   = 'TRACER_SOURCE'
+   assimilate_these_obs_types =  'TRACER_CONCENTRATION', 'VELOCITY',
+                                
+   evaluate_these_obs_types   = ''
    /
 
 &model_nml
-   num_grid_points        = 10
+   num_grid_points        = 40
    grid_spacing_meters    = 100000
    time_step_days         = 0