Merge pull request #636 from rjleveque/gauges_at_t0

modify so gauge values always written at start time regardless of min_time_increment

Author: mandli

src/2d/shallow/gauges_module.f90

@@ -141,7 +141,6 @@ subroutine set_gauges(restart, num_eqn, num_aux, fname)
                 ! note that for lagrangian gauges, the x,y values read here 
                 ! might be overwritten if this is a restart
                 gauges(i)%buffer_index = 1
-                gauges(i)%last_time = gauges(i)%t_start
                 ! keep track of last position for lagrangian gauges,
                 ! initialize here in case checkpoint happens before 
                 ! ever writing this gauge:
@@ -164,10 +163,17 @@ subroutine set_gauges(restart, num_eqn, num_aux, fname)
             read(UNIT, *)
             read(UNIT, *) (gauges(i)%gtype, i=1, num_gauges)
 
+
             ! Read in q fields
             read(UNIT, *)
             read(UNIT, *)
             do i = 1, num_gauges
+
+                ! initialize last_time so that first gauge output will be
+                ! at time gauges(i)%t_start regardless of min_time_increment:
+                gauges(i)%last_time = gauges(i)%t_start - 1.d0 &
+                                      - gauges(i)%min_time_increment
+
                 allocate(gauges(i)%q_out_vars(num_eqn))
                 read(UNIT, *) gauges(i)%q_out_vars