[Dart-dev] [6706] DART/trunk/models/mpas_atm/model_mod.f90: if the get_close_obs() routine has to do a vertical conversion

[nancy at ucar.edu]

Revision: 6706
Author:   nancy
Date:     2013-12-31 11:55:22 -0700 (Tue, 31 Dec 2013)
Log Message:
-----------
if the get_close_obs() routine has to do a vertical conversion
for vertical localization, have it store the converted location
in the list that gets passed back to filter.  with this single-line
change, each observation which is close to other obs will only
have to do a vertical conversion once instead of multiple times.
this change was tested by so-young and produced identical results
with a faster running speed.

Modified Paths:
--------------
    DART/trunk/models/mpas_atm/model_mod.f90

-------------- next part --------------
Modified: DART/trunk/models/mpas_atm/model_mod.f90
===================================================================
--- DART/trunk/models/mpas_atm/model_mod.f90	2013-12-30 21:22:40 UTC (rev 6705)
+++ DART/trunk/models/mpas_atm/model_mod.f90	2013-12-31 18:55:22 UTC (rev 6706)
@@ -2043,11 +2043,11 @@
   else if (base_which /= vert_localization_coord) then
       call vert_convert(ens_mean, base_obs_loc, base_obs_kind, ztypeout, istatus1)
       if(debug > 5) then
-      call write_location(0,base_obs_loc,charstring=string1)
-      call error_handler(E_MSG, 'get_close_obs: base_obs_loc',string1,source, revision, revdate)
-  endif
+         call write_location(0,base_obs_loc,charstring=string1)
+         call error_handler(E_MSG, 'get_close_obs: base_obs_loc',string1,source, revision, revdate)
+     endif
+   endif
 endif
-endif
 
 if (istatus1 == 0) then
 
@@ -2070,6 +2070,7 @@
       if (.not. horiz_dist_only) then
           if (local_obs_which /= vert_localization_coord) then
               call vert_convert(ens_mean, local_obs_loc, obs_kind(t_ind), ztypeout, istatus2)
+              obs_loc(t_ind) = local_obs_loc
           else
               istatus2 = 0
           endif