[Dart-dev] DART/branches Revision: 11087

[dart at ucar.edu]

hendric at ucar.edu
2017-02-09 16:50:03 -0700 (Thu, 09 Feb 2017)
61
setting observation time for regression factor diagnostics.




Modified: DART/branches/rma_single_file/assim_tools/assim_tools_mod.f90
===================================================================
--- DART/branches/rma_single_file/assim_tools/assim_tools_mod.f90	2017-02-09 23:18:42 UTC (rev 11086)
+++ DART/branches/rma_single_file/assim_tools/assim_tools_mod.f90	2017-02-09 23:50:03 UTC (rev 11087)
@@ -455,7 +455,7 @@
 ! HK I would like to move this to before the calculation of the forward operator so you could
 ! overwrite the vertical location with the required localization vertical coordinate when you 
 ! do the forward operator calculation
-call get_my_obs_loc(ens_handle, obs_ens_handle, obs_seq, keys, my_obs_loc, my_obs_kind, my_obs_type)
+call get_my_obs_loc(ens_handle, obs_ens_handle, obs_seq, keys, my_obs_loc, my_obs_kind, my_obs_type, obs_time)
 
 if (.not. lanai_bitwise) then
    ! convert the verical of all my observations to the localization coordinate
@@ -2682,7 +2682,7 @@
 
 !----------------------------------------------------------------------
 !> gets the location of of all my observations
-subroutine get_my_obs_loc(state_ens_handle, obs_ens_handle, obs_seq, keys, my_obs_loc, my_obs_kind, my_obs_type)
+subroutine get_my_obs_loc(state_ens_handle, obs_ens_handle, obs_seq, keys, my_obs_loc, my_obs_kind, my_obs_type, my_obs_time)
 
 type(ensemble_type),      intent(in)  :: state_ens_handle
 type(ensemble_type),      intent(in)  :: obs_ens_handle
@@ -2690,6 +2690,7 @@
 integer,                  intent(in)  :: keys(:)
 type(location_type),      intent(out) :: my_obs_loc(:)
 integer,                  intent(out) :: my_obs_type(:), my_obs_kind(:)
+type(time_type),          intent(out) :: my_obs_time
 
 type(obs_type) :: observation
 type(obs_def_type)   :: obs_def
@@ -2716,6 +2717,9 @@
    endif
 end do Get_Obs_Locations
 
+! Need the time for regression diagnostics potentially; get from first observation
+my_obs_time = get_obs_def_time(obs_def)
+
 end subroutine get_my_obs_loc
 
 !--------------------------------------------------------------------