anchor.diff:      !  wish to store once. e.g. [ELV], or [UH, VH, ELV]
       ! @param gauge_ids real array of size n_gauges giving an ID for each
       !  gauge. Make it REAL to avoid truncation issues for large IDs
  -    ! @param bounding_box If provided, only keep gauges inside a given bounding box
  -    ! @param priority_gauges If provided, a logical array with one entry per gauge. Gauges with .TRUE. are retained, others are
  -    ! removed.
  +    ! @param bounding_box If provided, only keep gauges inside a given bounding 
  +    ! box
  +    ! @param priority_gauges If provided, a logical array with one entry per 
  +    ! gauge. 
  +    ! Gauges with .TRUE. are retained, others are removed.
       !
       subroutine allocate_gauges(point_gauges, xy_coordinates, &
           time_series_var_indices, static_var_indices, gauge_ids, &
  
positive.code:     ! @param point_gauges variable of type point_gauge_type
      ! @param domain_dx cell size [dx,dy] in the domain
      ! @param domain_nx number of cells [nx,ny] in the domain
      ! @param domain_U the array with dimensions [nx, ny, :] in which we look up values at the gauges
      ! @param netcdf_gauge_output_file name of output file
      ! @param attribute_names character vector of names for additional attributes to add to the netcdf file
      ! @param attribute_values character vector of values for the additional attributes in the netcdf file
      !
      subroutine initialise_gauges(point_gauges, domain_lower_left, domain_dx, &
          domain_nx, domain_U,&
  
negative.code:         close(point_gauges%static_output_unit)
  
          ! Initialise time-varying output files
          open(newunit=point_gauges%time_series_output_unit, &
              file=point_gauges%time_series_output_file, &
              access='stream', form='unformatted')
  
          ! Initialise gauge metadata file
  