anchor.diff:  !!------------------------------------------------------------
   module io_routines
       use netcdf
  +    
       implicit none
       ! maximum number of dimensions for a netCDF file
       integer,parameter::io_maxDims=10
  -    
  +    ! maximum length of a dimension name not this could be in data_structures, but that simply adds
  +    ! one more dependency on data_structures and prevents this being compiled in parallel
  +    integer,parameter::MAXDIMLENGTH=255
       !>------------------------------------------------------------
       !! Generic interface to the netcdf read routines
       !!------------------------------------------------------------
  
positive.code:         ! Close the file, freeing all resources.
          call check( nf90_close(ncid) )
      end subroutine io_write2d
      
      !>------------------------------------------------------------
      !! Simple error handling for common netcdf file errors
  
negative.code:         integer, parameter :: ndims = 6
          ! This will be the netCDF ID for the file and data variable.
          integer :: ncid, varid,temp_dimid,dimids(ndims)
          character(len=MAXDIMLENGTH), dimension(6) :: dims
          
          if (present(dimnames)) then
              dims = dimnames
  