[Dart-dev] [4674] DART/trunk/diagnostics/threed_sphere/obs_diag.f90: Setting the value of 'ivert' so the write() statement

nancy at ucar.edu nancy at ucar.edu
Wed Jan 26 10:20:36 MST 2011


Revision: 4674
Author:   thoar
Date:     2011-01-26 10:20:36 -0700 (Wed, 26 Jan 2011)
Log Message:
-----------
Setting the value of 'ivert' so the write() statement
has something meaningful to write. The write() statement
logs all the observations that are more than 3 standard deviations 
to a file called LargeInnov.txt This code block is a remnant
from before the days of a useful DART error code and could easily
be removed. It has no bearing on any of the diagnostics being
calculated.

Modified Paths:
--------------
    DART/trunk/diagnostics/threed_sphere/obs_diag.f90

-------------- next part --------------
Modified: DART/trunk/diagnostics/threed_sphere/obs_diag.f90
===================================================================
--- DART/trunk/diagnostics/threed_sphere/obs_diag.f90	2011-01-26 15:37:44 UTC (rev 4673)
+++ DART/trunk/diagnostics/threed_sphere/obs_diag.f90	2011-01-26 17:20:36 UTC (rev 4674)
@@ -37,7 +37,7 @@
                              KIND_U_WIND_COMPONENT, KIND_V_WIND_COMPONENT
 use     location_mod, only : location_type, get_location, set_location_missing,   &
                              write_location, operator(/=), is_location_in_region, &
-                             set_location,                                        &
+                             set_location, query_location,                        &
                              vert_is_undef,    VERTISUNDEF,    &
                              vert_is_surface,  VERTISSURFACE,  &
                              vert_is_level,    VERTISLEVEL,    &
@@ -902,6 +902,7 @@
          obs_time = get_obs_def_time(obs_def)
          obs_loc  = get_obs_def_location(obs_def)
          obsloc3  = get_location(obs_loc)
+         ivert    = nint(query_location(obs_loc))
 
          obslon   = obsloc3(1) ! [  0, 360]
          obslat   = obsloc3(2) ! [-90,  90]


More information about the Dart-dev mailing list