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: !!------------------------------------------------------------
  module io_routines
      use netcdf
      implicit none
      ! maximum number of dimensions for a netCDF file
      integer,parameter::io_maxDims=10
      
      !>------------------------------------------------------------
      !! Generic interface to the netcdf read routines
      !!------------------------------------------------------------
  
negative.code:         ! End define mode. This tells netCDF we are done defining metadata.
          call check( nf90_enddef(ncid) )
          
          ! write the actual data to the file
          call check( nf90_put_var(ncid, varid, data_out), trim(filename)//":"//trim(varname))
      
          ! Close the file, freeing all resources.
  