[Dart-dev] DART/branches Revision: 12348

[dart at ucar.edu]

nancy at ucar.edu
2018-01-17 14:54:45 -0700 (Wed, 17 Jan 2018)
72
move the initialization for high obs limit into
a separate subroutine.




Modified: DART/branches/recam/models/cam-fv/model_mod.f90
===================================================================
--- DART/branches/recam/models/cam-fv/model_mod.f90	2018-01-17 21:45:56 UTC (rev 12347)
+++ DART/branches/recam/models/cam-fv/model_mod.f90	2018-01-17 21:54:45 UTC (rev 12348)
@@ -290,16 +290,9 @@
    are_damping = .true.
 endif
 
-! set top limit where obs are discarded
+! set top limit where obs are discarded.  -1 to disable.
 if (no_assim_above_pressure > 0.0_r8) then
-   write(string1, *) 'discarding observations above a pressure level of ', &
-                      no_assim_above_pressure, ' Pascals' 
-   call error_handler(E_MSG, 'static_init_model', string1, source, revision, revdate)
-
-   ! compute both height and pressure columns once, based on a surface
-   ! pressure of 1010 mb.  use for quick conversions when absolute accuracy 
-   ! isn't a primary concern.
-   call store_generic_columns()
+   call init_discard_high_obs()
 endif
 
 
@@ -3381,6 +3374,24 @@
 end subroutine free_generic_columns
 
 !--------------------------------------------------------------------
+
+subroutine init_discard_high_obs()
+
+! compute a conversion table between height and pressure based on
+! a surface pressure of 1010 mb.  this is a fixed table and does not
+! vary with temperature, humidity or surface elevation. 
+! use only for quick conversions when absolute accuracy 
+! isn't a primary concern.
+
+call store_generic_columns()
+
+write(string1, *) 'Discarding observations above a pressure level of ', &
+                   no_assim_above_pressure, ' Pascals' 
+call error_handler(E_MSG, 'init_discard_high_obs', string1, source, revision, revdate)
+
+end subroutine init_discard_high_obs
+
+!--------------------------------------------------------------------
 ! initialize the info needed to damp the assimilation increments at
 ! the model top.  compute where the damping ramp starts and save that
 ! info in 2 globals: ramp_start and ramp_start_loc