anchor.diff:          logical :: all_gauges
           logical, allocatable :: points_inside(:)
           integer(ip), allocatable :: indices_inside(:)
  -            
  -            ! At the moment the space dimension must be 2 
  +
  +            ! At the moment the space dimension must be 2
               if (point_gauges%space_dim /= size(xy_coordinates, 1)) then
                   write(log_output_unit,*) 'gauge xy dimension is not equal ', &
                       point_gauges%space_dim
                   flush(log_output_unit)
               end if
   
  -            ! Make space for all gauges. Depending on optional arguments 
  +            ! Make space for all gauges. Depending on optional arguments
               ! passed, this may be updated again below
               all_gauges = .TRUE.
               n_gauges = size(xy_coordinates, 2, kind=ip)
   
               if(present(bounding_box) .or. present(priority_gauges)) then
  -                ! It is possible that not all the gauges will be in the 
  +                ! It is possible that not all the gauges will be in the
                   ! bounding box. Identify those which are inside
                   all_gauges = .FALSE.
   
  
positive.code:     subroutine setup_gauge_netcdf_output(point_gauges, netcdf_gauge_output_file, &
          attribute_names, attribute_values)
  
          type(point_gauge_type), intent(inout):: point_gauges 
          character(charlen), intent(in):: netcdf_gauge_output_file
          character(charlen), optional, intent(in):: attribute_names(:)
          character(charlen), optional, intent(in):: attribute_values(:)
  
negative.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
  