anchor.diff:          class(point_gauge_type), intent(inout):: point_gauges
   
           if (point_gauges%n_gauges > 0) then
  -#ifndef NONETCDF      
  -            ! Close the netcdf file 
  +#ifndef NONETCDF
  +            ! Close the netcdf file
               call check(nf90_close(point_gauges%netcdf_gauge_output_file_ID), __LINE__)
   #else
               close(point_gauges%time_series_output_unit)
  
positive.code:         ! which can speed up gauge-merging for multidomains
          logical :: priority_gauges_only = .FALSE.
  
          ! Indices of variables 
          integer(ip), allocatable:: time_series_var(:), static_var(:)
  
          ! Coordinates of gauges
  
negative.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, &
  