[Dart-dev] DART/branches Revision: 12181

[dart at ucar.edu]

nancy at ucar.edu
2017-12-05 14:37:25 -0700 (Tue, 05 Dec 2017)
157
we already had a namelist item for throwing out obs
above a given model level.  consolidate them.
also add the new namelist items to the default input.nml.




Modified: DART/branches/recam/models/cam-fv/model_mod.f90
===================================================================
--- DART/branches/recam/models/cam-fv/model_mod.f90	2017-12-05 21:30:14 UTC (rev 12180)
+++ DART/branches/recam/models/cam-fv/model_mod.f90	2017-12-05 21:37:25 UTC (rev 12181)
@@ -106,11 +106,10 @@
 integer            :: assimilation_period_days        = 0
 integer            :: assimilation_period_seconds     = 21600
 logical            :: use_log_vertical_scale          = .false.
-integer            :: no_assim_above_this_model_level = 5
+integer            :: no_assim_above_this_model_level = 4
 logical            :: use_damping_ramp_at_model_top   = .false.  
 integer            :: debug_level                     = 0
 logical            :: suppress_grid_info_in_output    = .false.
-integer            :: highest_obs_limit_in_state      = -1
 logical            :: custom_routine_to_generate_ensemble = .false.
 character(len=32)  :: fields_to_perturb(MAX_PERT)     = "QTY_TEMPERATURE"
 real(r8)           :: perturbation_amplitude(MAX_PERT)= 0.00001_r8
@@ -140,7 +139,6 @@
    no_assim_above_this_model_level, &
    use_damping_ramp_at_model_top,   &
    suppress_grid_info_in_output,    &
-   highest_obs_limit_in_state,      &
    custom_routine_to_generate_ensemble, &
    fields_to_perturb,               &
    perturbation_amplitude,          &
@@ -735,8 +733,8 @@
       ! level 1 is top, so test that the level numbers are *smaller* than the limit.
       ! (meaning the obs is above the given limit in at least one ensemble member)
 
-      if (highest_obs_limit_in_state > 0) then
-         if (any(four_bot_levs(icorner,:) <= highest_obs_limit_in_state)) then
+      if (no_assim_above_this_model_level > 0) then
+         if (any(four_bot_levs(icorner,:) <= no_assim_above_this_model_level)) then
             istatus(:) = 14
             return
          endif

Modified: DART/branches/recam/models/cam-fv/work/input.nml
===================================================================
--- DART/branches/recam/models/cam-fv/work/input.nml	2017-12-05 21:30:14 UTC (rev 12180)
+++ DART/branches/recam/models/cam-fv/work/input.nml	2017-12-05 21:37:25 UTC (rev 12181)
@@ -143,9 +143,12 @@
    assimilation_period_seconds = 21600
    assimilation_period_days    = 0
    no_assim_above_this_model_level = 5
-   use_damping_ramp_at_model_top = .false.
-   suppress_grid_info_in_output  = .false.
-   debug_level                 = 13
+   use_damping_ramp_at_model_top   = .false.
+   suppress_grid_info_in_output    = .false.
+   custom_routine_to_generate_ensemble = .true.
+   fields_to_perturb               = 'QTY_TEMPERATURE'
+   perturbation_amplitude,         = 0.0001
+   debug_level                     = 13
    /
 
 !#!    vert_coord          = 'log_invP'