anchor.diff:          do i = 1, size(point_gauges%static_var, kind=ip)
               j = point_gauges%static_var(i)
   
  -            ! Initial stage 
  +            ! Initial stage
               if(j == STG) then
                   call check( nf90_def_var(iNcid, "stage0", NF90_REAL4, [iDimStation_ID], iVar_STG_static_ID), __LINE__)
                   call check(nf90_put_att(iNcid, iVar_STG_static_ID, "long_name", &
  
positive.code: #else
          ! Variables for writing out to netcdf
          character(len=charlen):: netcdf_gauge_output_file
          integer :: netcdf_gauge_output_file_ID 
          integer(ip) :: netcdf_time_var_ID
          integer(ip), allocatable:: time_series_ncdf_iVar_ID(:)
          ! Count how many steps have been output
  
negative.code: 
          ! This is the standard netcdf attribute for time-series
          call check(nf90_put_att(iNcid, nf90_global, "featureType", "timeSeries"), __LINE__)
  
          ! Here we add other attributes that might be useful (e.g. the name of the input stage and elevation)
          if((present(attribute_names)).and.(present(attribute_values))) then
              do i = 1, size(attribute_names, kind=ip)
                  call check(nf90_put_att(iNcid, nf90_global, attribute_names(i), attribute_values(i)), __LINE__)
  