anchor.diff:              end if
   
           end do
  -        
  +
           ! Get static values -- could potentially write them here, and the
           ! deallocate the array, if conserving memory is very important
           call update_gauge_var(point_gauges, point_gauges%static_values, &
  
positive.code:     ! Convenience subroutine to lookup variables in domain_U at point_gauges,
      ! and populate time_series_values with the values of quantities var_inds
      !
      ! This can be used to populate either point_gauges%static_values or 
      ! point_gauges%time_series_values
      !
      ! @param point_gauges point_gauge_type
      ! @param time_series_values array with size [n_gauges, size(var_inds)] 
      ! that holds the time series data
      ! @param var_inds Array giving indices in the 3rd dimension of domain_U 
      ! that we write out. e.g. [STG, UH, VH], or [STG]
      !
      subroutine update_gauge_var(point_gauges, time_series_values, domain_U, &
  
negative.code: 
          ! Time
          call check( nf90_def_var(iNcid, 'time', NF90_REAL4, (/ iDimTime_ID /), iVarTIME_ID) , __LINE__)
          point_gauges%netcdf_time_var_ID = iVarTIME_ID
          call check(nf90_put_att(iNcid, iVarTIME_ID, "standard_name", &
              "time"), __LINE__)
          call check(nf90_put_att(iNcid, iVarTIME_ID, "long_name", &
  