mfix_netcdf_mod.f

Go to the documentation of this file.

module mfix_netcdf

#ifdef NETCDF
  include 'netcdf.inc'
  USE netcdf ! ignore-depcomp
#else
  integer nf90_64bit_offset

  integer nf90_double
  integer nf90_int
  integer nf90_nowrite
  integer nf90_clobber

  parameter(nf90_64bit_offset = 0)
  parameter(nf90_double = 0)
  parameter(nf90_int = 0)
  parameter(nf90_nowrite = 0)
  parameter(nf90_clobber = 0)
#endif

  ! Overloaded variable functions
  interface mfix_nf90_def_var
     module procedure mfix_nf90_def_var_scalar, mfix_nf90_def_var_onedim, mfix_nf90_def_var_manydims
  end interface mfix_nf90_def_var ! MFIX_nf90_def_var

  interface mfix_nf90_put_var
     module procedure mfix_nf90_put_var_text,                                   &
          mfix_nf90_put_var_fourbyteint,  &
          mfix_nf90_put_var_fourbytereal, mfix_nf90_put_var_eightbytereal
     module procedure mfix_nf90_put_var_1d_text,                                      &
          mfix_nf90_put_var_1d_fourbyteint, &
          mfix_nf90_put_var_1d_fourbytereal, mfix_nf90_put_var_1d_eightbytereal
     module procedure mfix_nf90_put_var_2d_text,                                       &
          mfix_nf90_put_var_2d_fourbyteint,   &
          mfix_nf90_put_var_2d_fourbytereal, mfix_nf90_put_var_2d_eightbytereal
     module procedure mfix_nf90_put_var_3d_text,                                       &
          mfix_nf90_put_var_3d_fourbyteint, &
          mfix_nf90_put_var_3d_fourbytereal, mfix_nf90_put_var_3d_eightbytereal
     module procedure mfix_nf90_put_var_4d_fourbyteint , mfix_nf90_put_var_4d_eightbytereal
  end interface mfix_nf90_put_var ! MFIX_nf90_put_var

  interface mfix_nf90_get_var
     module procedure mfix_nf90_get_var_text,                                   &
          mfix_nf90_get_var_fourbyteint, & ! MFIX_nf90_get_var_EightByteInt, &
          mfix_nf90_get_var_fourbytereal, mfix_nf90_get_var_eightbytereal
     module procedure mfix_nf90_get_var_1d_text,                                      &
          mfix_nf90_get_var_1d_fourbyteint , &
          mfix_nf90_get_var_1d_fourbytereal, mfix_nf90_get_var_1d_eightbytereal
     module procedure mfix_nf90_get_var_2d_text,                                      &
          mfix_nf90_get_var_2d_fourbytereal, mfix_nf90_get_var_2d_eightbytereal
     module procedure mfix_nf90_get_var_3d_text,                                      &
          mfix_nf90_get_var_3d_fourbyteint, &
          mfix_nf90_get_var_3d_fourbytereal, mfix_nf90_get_var_3d_eightbytereal
     module procedure mfix_nf90_get_var_4d_text
     module procedure mfix_nf90_get_var_5d_text
     module procedure mfix_nf90_get_var_6d_text
     module procedure mfix_nf90_get_var_7d_text
  end interface mfix_nf90_get_var ! MFIX_nf90_get_var

contains
  function mfix_nf90_def_var_scalar(ncid, name, xtype, varid)
    integer,               intent( in) :: ncid
    character (len = *),   intent( in) :: name
    integer,               intent( in) :: xtype
    integer,               intent(inout) :: varid
    integer                            :: MFIX_nf90_def_var_Scalar

    ! Dummy - shouldn't get used
    integer, dimension(1) :: dimids

#ifdef NETCDF
    mfix_nf90_def_var_scalar = nf90_def_var(ncid, name, xtype, varid)
#else
    mfix_nf90_def_var_scalar = 1
#endif
  end function mfix_nf90_def_var_scalar
  ! -----
  function mfix_nf90_def_var_onedim(ncid, name, xtype, dimids, varid)
    integer,               intent( in) :: ncid
    character (len = *),   intent( in) :: name
    integer,               intent( in) :: xtype
    integer,               intent( in) :: dimids
    integer,               intent(inout) :: varid
    integer                            :: MFIX_nf90_def_var_oneDim

    integer, dimension(1) :: dimidsA
    dimidsa(1) = dimids
#ifdef NETCDF
    mfix_nf90_def_var_onedim = nf90_def_var(ncid, name, xtype, dimids, varid)
#else
    mfix_nf90_def_var_onedim = 1
#endif
  end function mfix_nf90_def_var_onedim
  ! -----
  function mfix_nf90_def_var_manydims(ncid, name, xtype, dimids, varid)
    integer,               intent( in) :: ncid
    character (len = *),   intent( in) :: name
    integer,               intent( in) :: xtype
    integer, dimension(:), intent( in) :: dimids
    integer,               intent(inout) :: varid
    integer                            :: MFIX_nf90_def_var_ManyDims

#ifdef NETCDF
    mfix_nf90_def_var_manydims = nf90_def_var(ncid, name, xtype, dimids, varid)
#else
    mfix_nf90_def_var_manydims = 1
#endif
  end function mfix_nf90_def_var_manydims


  function mfix_nf90_put_var_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *),             intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_text


#ifdef NETCDF
    mfix_nf90_put_var_text = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_text = 1
#endif
  end function mfix_nf90_put_var_text

  function mfix_nf90_get_var_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *),             intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_text


#ifdef NETCDF
    mfix_nf90_get_var_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_text = 1
#endif
  end function mfix_nf90_get_var_text


  function mfix_nf90_put_var_1d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_1D_text

    integer, parameter                    :: numDims = 1

#ifdef NETCDF
    mfix_nf90_put_var_1d_text = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_1d_text = 1
#endif

  end function mfix_nf90_put_var_1d_text


  function mfix_nf90_put_var_2d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_2D_text

    integer, parameter                    :: numDims = 2

#ifdef NETCDF
    mfix_nf90_put_var_2d_text = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_2d_text = 1
#endif

  end function mfix_nf90_put_var_2d_text


  function mfix_nf90_put_var_3d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_3D_text

    integer, parameter                    :: numDims = 3
#ifdef NETCDF
    mfix_nf90_put_var_3d_text = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_3d_text = 1
#endif

  end function mfix_nf90_put_var_3d_text



  function mfix_nf90_get_var_1d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_1D_text

    integer, parameter                  :: numDims = 1
#ifdef NETCDF
    mfix_nf90_get_var_1d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_1d_text = 1
#endif
  end function mfix_nf90_get_var_1d_text


  function mfix_nf90_get_var_2d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_2D_text

    integer, parameter                  :: numDims = 2
#ifdef NETCDF
    mfix_nf90_get_var_2d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_2d_text = 1
#endif
  end function mfix_nf90_get_var_2d_text


  function mfix_nf90_get_var_3d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_3D_text

    integer, parameter                  :: numDims = 3
#ifdef NETCDF
    mfix_nf90_get_var_3d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_3d_text = 1
#endif
  end function mfix_nf90_get_var_3d_text


  function mfix_nf90_get_var_4d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_4D_text

    integer, parameter                  :: numDims = 4
#ifdef NETCDF
    mfix_nf90_get_var_4d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_4d_text = 1
#endif
  end function mfix_nf90_get_var_4d_text


  function mfix_nf90_get_var_5d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :, :, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_5D_text

    integer, parameter                  :: numDims = 5

#ifdef NETCDF
    mfix_nf90_get_var_5d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_5d_text = 1
#endif
  end function mfix_nf90_get_var_5d_text


  function mfix_nf90_get_var_6d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :, :, :, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_6D_text

    integer, parameter                  :: numDims = 6
#ifdef NETCDF
    mfix_nf90_get_var_6d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_6d_text = 1
#endif
  end function mfix_nf90_get_var_6d_text


  function mfix_nf90_get_var_7d_text(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    character (len = *), dimension(:, :, :, :, :, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_7D_text

    integer, parameter                  :: numDims = 7
#ifdef NETCDF
    mfix_nf90_get_var_7d_text = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_7d_text = 1
#endif
  end function mfix_nf90_get_var_7d_text

  function mfix_nf90_put_var_fourbyteint(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    integer , intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_put_var_FourByteInt

#ifdef NETCDF
    mfix_nf90_put_var_fourbyteint = nf90_put_var(ncid, varid, values, start)
#else
    mfix_nf90_put_var_fourbyteint = 1
#endif
  end function mfix_nf90_put_var_fourbyteint

  function mfix_nf90_put_var_fourbytereal(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    real , intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_put_var_FourByteReal


#ifdef NETCDF
    mfix_nf90_put_var_fourbytereal = nf90_put_var(ncid, varid, values, start)
#else
    mfix_nf90_put_var_fourbytereal = 1
#endif
  end function mfix_nf90_put_var_fourbytereal


  function mfix_nf90_put_var_eightbytereal(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    double precision, intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_put_var_EightByteReal


#ifdef NETCDF
    mfix_nf90_put_var_eightbytereal = nf90_put_var(ncid, varid, values, start)
#else
    mfix_nf90_put_var_eightbytereal = 1
#endif
  end function mfix_nf90_put_var_eightbytereal



  function mfix_nf90_get_var_twobyteint(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    integer , intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_get_var_TwoByteInt

    integer                               :: counter


#ifdef NETCDF
    mfix_nf90_get_var_twobyteint = nf90_get_var(ncid, varid, values, start)
#else
    mfix_nf90_get_var_twobyteint = 1
#endif
  end function mfix_nf90_get_var_twobyteint


  function mfix_nf90_get_var_fourbyteint(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    integer , intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_get_var_FourByteInt

    integer                               :: counter
    integer                               :: defaultInteger

#ifdef NETCDF
    mfix_nf90_get_var_fourbyteint = nf90_get_var(ncid, varid, values, start)
#else
    mfix_nf90_get_var_fourbyteint = 1
#endif
  end function mfix_nf90_get_var_fourbyteint


  !   function MFIX_nf90_get_var_EightByteInt(ncid, varid, values, start)
  !     integer,                         intent( in) :: ncid, varid
  !     integer , intent(inout) :: values
  !     integer, dimension(:), optional, intent( in) :: start
  !     integer                                      :: MFIX_nf90_get_var_EightByteInt
  !
  !
#ifdef NETCDF
  !     MFIX_nf90_get_var_EightByteInt = nf90_get_var(ncid, varid, values, start)
#else
  !     MFIX_nf90_get_var_EightByteInt = 1
#endif
  !   end function MFIX_nf90_get_var_EightByteInt


  function mfix_nf90_get_var_fourbytereal(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    real , intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_get_var_FourByteReal


#ifdef NETCDF
    mfix_nf90_get_var_fourbytereal = nf90_get_var(ncid, varid, values, start)
#else
    mfix_nf90_get_var_fourbytereal = 1
#endif
  end function mfix_nf90_get_var_fourbytereal


  function mfix_nf90_get_var_eightbytereal(ncid, varid, values, start)
    integer,                         intent( in) :: ncid, varid
    double precision, intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start
    integer                                      :: MFIX_nf90_get_var_EightByteReal

#ifdef NETCDF
    mfix_nf90_get_var_eightbytereal = nf90_get_var(ncid, varid, values, start)
#else
    mfix_nf90_get_var_eightbytereal = 1
#endif
  end function mfix_nf90_get_var_eightbytereal


  function mfix_nf90_get_var_1d_twobyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_1D_TwoByteInt

#ifdef NETCDF
    mfix_nf90_get_var_1d_twobyteint = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_1d_twobyteint = 1
#endif

  end function mfix_nf90_get_var_1d_twobyteint


  function mfix_nf90_get_var_2d_twobyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_2D_TwoByteInt

#ifdef NETCDF
    mfix_nf90_get_var_2d_twobyteint = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_2d_twobyteint = 1
#endif
  end function mfix_nf90_get_var_2d_twobyteint


  function mfix_nf90_get_var_3d_twobyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_3D_TwoByteInt

#ifdef NETCDF
    mfix_nf90_get_var_3d_twobyteint = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_3d_twobyteint = 1
#endif
  end function mfix_nf90_get_var_3d_twobyteint

  function mfix_nf90_get_var_1d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_1D_FourByteInt
#ifdef NETCDF
    mfix_nf90_get_var_1d_fourbyteint = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_1d_fourbyteint = 1
#endif
  end function mfix_nf90_get_var_1d_fourbyteint


  function mfix_nf90_get_var_2d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_2D_FourByteInt

#ifdef NETCDF
    mfix_nf90_get_var_2d_fourbyteint = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_2d_fourbyteint = 1
#endif
  end function mfix_nf90_get_var_2d_fourbyteint


  function mfix_nf90_get_var_3d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_3D_FourByteInt
#ifdef NETCDF
    mfix_nf90_get_var_3d_fourbyteint = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_3d_fourbyteint = 1
#endif
  end function mfix_nf90_get_var_3d_fourbyteint



  function mfix_nf90_get_var_1d_fourbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    real , dimension(:), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_1D_FourByteReal
#ifdef NETCDF
    mfix_nf90_get_var_1d_fourbytereal = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_1d_fourbytereal = 1
#endif
  end function mfix_nf90_get_var_1d_fourbytereal


  function mfix_nf90_get_var_2d_fourbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    real , dimension(:, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_2D_FourByteReal

#ifdef NETCDF
    mfix_nf90_get_var_2d_fourbytereal = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_2d_fourbytereal = 1
#endif
  end function mfix_nf90_get_var_2d_fourbytereal


  function mfix_nf90_get_var_3d_fourbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    real , dimension(:, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_3D_FourByteReal
#ifdef NETCDF
    mfix_nf90_get_var_3d_fourbytereal = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_3d_fourbytereal = 1
#endif
  end function mfix_nf90_get_var_3d_fourbytereal

  function mfix_nf90_get_var_1d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_1D_EightByteReal

#ifdef NETCDF
    mfix_nf90_get_var_1d_eightbytereal = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_1d_eightbytereal = 1
#endif
  end function mfix_nf90_get_var_1d_eightbytereal


  function mfix_nf90_get_var_2d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_2D_EightByteReal

#ifdef NETCDF
    mfix_nf90_get_var_2d_eightbytereal = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_2d_eightbytereal = 1
#endif
  end function mfix_nf90_get_var_2d_eightbytereal


  function mfix_nf90_get_var_3d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:, :, :), &
         intent(inout) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_get_var_3D_EightByteReal

#ifdef NETCDF
    mfix_nf90_get_var_3d_eightbytereal = nf90_get_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_get_var_3d_eightbytereal = 1
#endif
  end function mfix_nf90_get_var_3d_eightbytereal

  function mfix_nf90_inquire_dimension(ncid, dimid, name, len)
    integer,                       intent( in) :: ncid, dimid
    character (len = *), optional, intent(inout) :: name
    integer,             optional, intent(inout) :: len
    integer                                    :: MFIX_nf90_inquire_dimension

#ifdef NETCDF
    mfix_nf90_inquire_dimension = nf90_inquire_dimension(ncid, dimid, name, len)
#else
    mfix_nf90_inquire_dimension = 1
#endif
  end function mfix_nf90_inquire_dimension

  function mfix_nf90_put_var_1d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_1D_FourByteInt

#ifdef NETCDF
    mfix_nf90_put_var_1d_fourbyteint = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_1d_fourbyteint = 1
#endif
  end function mfix_nf90_put_var_1d_fourbyteint


  function mfix_nf90_put_var_2d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_2D_FourByteInt

#ifdef NETCDF
    mfix_nf90_put_var_2d_fourbyteint = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_2d_fourbyteint = 1
#endif
  end function mfix_nf90_put_var_2d_fourbyteint


  function mfix_nf90_put_var_3d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_3D_FourByteInt
#ifdef NETCDF
    mfix_nf90_put_var_3d_fourbyteint = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_3d_fourbyteint = 1
#endif
  end function mfix_nf90_put_var_3d_fourbyteint

  function mfix_nf90_put_var_4d_fourbyteint(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    integer , dimension(:, :, : , :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_4D_FourByteInt
#ifdef NETCDF
    mfix_nf90_put_var_4d_fourbyteint = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_4d_fourbyteint = 1
#endif
  end function mfix_nf90_put_var_4d_fourbyteint

  function mfix_nf90_put_var_1d_fourbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    real , dimension(:), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_1D_FourByteReal

#ifdef NETCDF
    mfix_nf90_put_var_1d_fourbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_1d_fourbytereal = 1
#endif
  end function mfix_nf90_put_var_1d_fourbytereal


  function mfix_nf90_put_var_2d_fourbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    real , dimension(:, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_2D_FourByteReal

#ifdef NETCDF
    mfix_nf90_put_var_2d_fourbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_2d_fourbytereal = 1
#endif
  end function mfix_nf90_put_var_2d_fourbytereal


  function mfix_nf90_put_var_3d_fourbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    real , dimension(:, :, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_3D_FourByteReal
#ifdef NETCDF
    mfix_nf90_put_var_3d_fourbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_3d_fourbytereal = 1
#endif
  end function mfix_nf90_put_var_3d_fourbytereal


  function mfix_nf90_put_var_1d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_1D_EightByteReal

#ifdef NETCDF
    mfix_nf90_put_var_1d_eightbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_1d_eightbytereal = 1
#endif
  end function mfix_nf90_put_var_1d_eightbytereal


  function mfix_nf90_put_var_2d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_2D_EightByteReal
#ifdef NETCDF
    mfix_nf90_put_var_2d_eightbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_2d_eightbytereal = 1
#endif
  end function mfix_nf90_put_var_2d_eightbytereal


  function mfix_nf90_put_var_3d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:, :, :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_3D_EightByteReal
#ifdef NETCDF
    mfix_nf90_put_var_3d_eightbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_3d_eightbytereal = 1
#endif
  end function mfix_nf90_put_var_3d_eightbytereal


  function mfix_nf90_put_var_4d_eightbytereal(ncid, varid, values, start, count, stride, map)
    integer,                         intent( in) :: ncid, varid
    double precision, dimension(:, :, : , :), &
         intent( in) :: values
    integer, dimension(:), optional, intent( in) :: start, count, stride, map
    integer                                      :: MFIX_nf90_put_var_4D_EightByteReal
#ifdef NETCDF
    mfix_nf90_put_var_4d_eightbytereal = nf90_put_var(ncid, varid, values, start, count, stride, map)
#else
    mfix_nf90_put_var_4d_eightbytereal = 1
#endif
  end function mfix_nf90_put_var_4d_eightbytereal

  function mfix_usingnetcdf()
    logical :: MFIX_usingNETCDF
#ifdef NETCDF
    mfix_usingnetcdf = .true.
#else
    mfix_usingnetcdf = .false.
#endif
  end function mfix_usingnetcdf

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  !                          check_netcdf                             !
  !
  subroutine mfix_check_netcdf( status )
    implicit none
    integer, intent ( in) :: status

#ifdef NETCDF
    if (status /= nf90_noerr) then
       write (*,*) ' ******************************************'
       write (*,*) trim(nf90_strerror(status))
       write (*,*) ' ******************************************'
    end if
#endif
  end subroutine mfix_check_netcdf

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  function mfix_nf90_create(path, cmode, ncid, initialsize, chunksize)
    character (len = *), intent(in   ) :: path
    integer,             intent(in   ) :: cmode
    integer,             intent(inout) :: ncid
    integer, optional,   intent(in   ) :: initialsize
    integer, optional,   intent(inout) :: chunksize
    integer                            :: MFIX_nf90_create

    integer :: fileSize

#ifdef NETCDF
    mfix_nf90_create = nf90_create(path, cmode, ncid, initialsize, chunksize)
#else
    mfix_nf90_create = 0
#endif
    return
  end function mfix_nf90_create


  function mfix_nf90_def_dim(ncid, name, len, dimid)
    integer,             intent( in) :: ncid
    character (len = *), intent( in) :: name
    integer,             intent( in) :: len
    integer,             intent(inout) :: dimid
    integer                          :: MFIX_nf90_def_dim

#ifdef NETCDF
    mfix_nf90_def_dim = nf90_def_dim(ncid, name, len, dimid)
#else
    mfix_nf90_def_dim = 1
#endif
  end function mfix_nf90_def_dim


  function mfix_nf90_enddef(ncid, h_minfree, v_align, v_minfree, r_align)
    integer,           intent( in) :: ncid
    integer, optional, intent( in) :: h_minfree, v_align, v_minfree, r_align
    integer                        :: MFIX_nf90_enddef

#ifdef NETCDF
    mfix_nf90_enddef = nf90_enddef(ncid, h_minfree, v_align, v_minfree, r_align)
#else
    mfix_nf90_enddef = 1
#endif
  end function mfix_nf90_enddef


  function mfix_nf90_open(path, mode, ncid, chunksize)
    character (len = *), intent(in   ) :: path
    integer,             intent(in   ) :: mode
    integer,             intent(inout) :: ncid
    integer, optional,   intent(inout) :: chunksize
    integer                            :: MFIX_nf90_open
#ifdef NETCDF
    mfix_nf90_open = nf90_open(path, mode, ncid, chunksize)
#else
    mfix_nf90_open = 1
#endif
  end function mfix_nf90_open


  function mfix_nf90_inquire(ncid, nDimensions, nVariables, nAttributes, unlimitedDimId, formatNum)
    integer,           intent( in) :: ncid
    integer, optional, intent(inout) :: nDimensions, nVariables, nAttributes, unlimitedDimId, formatNum
    integer                        :: MFIX_nf90_inquire

    integer :: nDims, nVars, nGAtts, unlimDimId, frmt

#ifdef NETCDF
    mfix_nf90_inquire = nf90_inquire(ncid, ndimensions, nvariables, nattributes, unlimiteddimid, formatnum)
#else
    mfix_nf90_inquire = 1
#endif
  end function mfix_nf90_inquire


  function mfix_nf90_inq_dimid(ncid, name, dimid)
    integer,             intent( in) :: ncid
    character (len = *), intent( in) :: name
    integer,             intent(inout) :: dimid
    integer                          :: MFIX_nf90_inq_dimid

#ifdef NETCDF
    mfix_nf90_inq_dimid = nf90_inq_dimid(ncid, name, dimid)
#else
    mfix_nf90_inq_dimid = 1
#endif
  end function mfix_nf90_inq_dimid

  function mfix_nf90_inq_varid(ncid, name, varid)
    integer,             intent( in) :: ncid
    character (len = *), intent( in) :: name
    integer,             intent(inout) :: varid
    integer                          :: MFIX_nf90_inq_varid

#ifdef NETCDF
    mfix_nf90_inq_varid = nf90_inq_varid(ncid, name, varid)
#else
    mfix_nf90_inq_varid = 1
#endif
  end function mfix_nf90_inq_varid

  function mfix_nf90_close(ncid)
    integer,             intent( in) :: ncid
    integer                          :: MFIX_nf90_close

#ifdef NETCDF
    mfix_nf90_close = nf90_close(ncid)
#else
    mfix_nf90_close = 1
#endif
  end function mfix_nf90_close

  subroutine mfix_ncvinq(ncid,varid,varnam,vartyp,nvdims,vdims,nvatts,rcode)
    implicit none

    integer       :: ncid , varid , vartyp , nvdims , nvatts , rcode
    integer       :: vdims(*)
    character(len=*)    varnam

#ifdef NETCDF
    call ncvinq(ncid,varid,varnam,vartyp,nvdims,vdims,nvatts,rcode)
#endif
    return
  end subroutine mfix_ncvinq


end module mfix_netcdf