anchor.diff:      interface io_write
           module procedure io_write6d, io_write3d, io_write2d, io_write3di
       end interface
  +    
  +    !>------------------------------------------------------------
  +    !! Generic interface to the netcdf read_attribute_TYPE routines
  +    !!------------------------------------------------------------
  +    interface io_read_attribute
  +        module procedure io_read_attribute_r
  +    end interface
  +    !, io_read_attribute_d, &
  +    !  io_read_attribute_i, io_read_attribute_c
   
   !   All routines are public
   contains
  
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
  