[Dart-dev] [5698] DART/branches/development/obs_sequence/obs_seq_coverage.f90: Fixed code segment writing the observation definitions.

nancy at ucar.edu nancy at ucar.edu
Fri Apr 13 16:46:02 MDT 2012


Revision: 5698
Author:   thoar
Date:     2012-04-13 16:46:02 -0600 (Fri, 13 Apr 2012)
Log Message:
-----------
Fixed code segment writing the observation definitions.
Incorrectly looping over the first N observation times for
each 'station' instead of all N times and ignoring the 
time slots that were/are empty.

Modified Paths:
--------------
    DART/branches/development/obs_sequence/obs_seq_coverage.f90

-------------- next part --------------
Modified: DART/branches/development/obs_sequence/obs_seq_coverage.f90
===================================================================
--- DART/branches/development/obs_sequence/obs_seq_coverage.f90	2012-04-13 22:29:52 UTC (rev 5697)
+++ DART/branches/development/obs_sequence/obs_seq_coverage.f90	2012-04-13 22:46:02 UTC (rev 5698)
@@ -1417,9 +1417,12 @@
 
    call set_obs_def_kind(    obs_def, stations(i)%obs_type)
    call set_obs_def_location(obs_def, stations(i)%location)
-   TimeLoop : do j = 1,stations(i)%ntimes
-      call set_obs_def_time( obs_def, stations(i)%times(j))
-      call write_obs_def(iunit, obs_def, i, 'formatted')
+
+   TimeLoop : do j = 1,num_verification_times
+      if (stations(i)%times(j) /= no_time) then
+         call set_obs_def_time( obs_def, stations(i)%times(j))
+         call write_obs_def(iunit, obs_def, i, 'formatted')
+      endif
    enddo TimeLoop
 
    if (verbose) then


More information about the Dart-dev mailing list