[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