!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
#include <petsc/finclude/petscsys.h>
subroutine F90Array1dCreateScalar(array, start, len1, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start, len1
  PetscScalar, target :: array(start:start + len1 - 1)
  PetscScalar, pointer :: ptr(:)

  ptr => array
end subroutine

subroutine F90Array1dCreateReal(array, start, len1, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start, len1
  PetscReal, target :: array(start:start + len1 - 1)
  PetscReal, pointer :: ptr(:)

  ptr => array
end subroutine

subroutine F90Array1dCreateInt(array, start, len1, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start, len1
  PetscInt, target :: array(start:start + len1 - 1)
  PetscInt, pointer :: ptr(:)

  ptr => array
end subroutine

subroutine F90Array1dCreateMPIInt(array, start, len1, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start, len1
  PetscMPIInt, target :: array(start:start + len1 - 1)
  PetscMPIInt, pointer :: ptr(:)

  ptr => array
end subroutine

subroutine F90Array1dCreateFortranAddr(array, start, len1, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start, len1
  PetscFortranAddr, target :: array(start:start + len1 - 1)
  PetscFortranAddr, pointer :: ptr(:)

  ptr => array
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array1dAccessScalar(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:)
  PetscFortranAddr address
  PetscInt start

  if (.not. associated(ptr)) then
    address = 0
  else
    start = lbound(ptr, 1)
    call F90Array1dGetAddrScalar(ptr(start), address)
  end if
end subroutine

subroutine F90Array1dAccessReal(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:)
  PetscFortranAddr address
  PetscInt start

  if (.not. associated(ptr)) then
    address = 0
  else
    start = lbound(ptr, 1)
    call F90Array1dGetAddrReal(ptr(start), address)
  end if
end subroutine

subroutine F90Array1dAccessInt(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:)
  PetscFortranAddr address
  PetscInt start

  if (.not. associated(ptr)) then
    address = 0
  else
    start = lbound(ptr, 1)
    call F90Array1dGetAddrInt(ptr(start), address)
  end if
end subroutine

subroutine F90Array1dAccessMPIInt(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscMPIInt, pointer :: ptr(:)
  PetscFortranAddr address
  PetscInt start

  if (.not. associated(ptr)) then
    address = 0
  else
    start = lbound(ptr, 1)
    call F90Array1dGetAddrMPIInt(ptr(start), address)
  end if
end subroutine

subroutine F90Array1dAccessFortranAddr(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:)
  PetscFortranAddr address
  PetscInt start

  if (.not. associated(ptr)) then
    address = 0
  else
    start = lbound(ptr, 1)
    call F90Array1dGetAddrFortranAddr(ptr(start), address)
  end if
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array1dDestroyScalar(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:)

  nullify (ptr)
end subroutine

subroutine F90Array1dDestroyReal(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:)

  nullify (ptr)
end subroutine

subroutine F90Array1dDestroyInt(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:)

  nullify (ptr)
end subroutine

subroutine F90Array1dDestroyMPIInt(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscMPIInt, pointer :: ptr(:)

  nullify (ptr)
end subroutine

subroutine F90Array1dDestroyFortranAddr(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:)

  nullify (ptr)
end subroutine
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array2dCreateScalar(array, start1, len1, start2, len2, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscScalar, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1)
  PetscScalar, pointer :: ptr(:, :)

  ptr => array
end subroutine

subroutine F90Array2dCreateReal(array, start1, len1, start2, len2, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscReal, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1)
  PetscReal, pointer :: ptr(:, :)

  ptr => array
end subroutine

subroutine F90Array2dCreateInt(array, start1, len1, start2, len2, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1)
  PetscInt, pointer :: ptr(:, :)

  ptr => array
end subroutine

subroutine F90Array2dCreateFortranAddr(array, start1, len1, start2, len2, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscFortranAddr, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1)
  PetscFortranAddr, pointer :: ptr(:, :)

  ptr => array
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array2dAccessScalar(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :)
  PetscFortranAddr address
  PetscInt start1, start2

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  call F90Array2dGetAddrScalar(ptr(start1, start2), address)
end subroutine

subroutine F90Array2dAccessReal(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:, :)
  PetscFortranAddr address
  PetscInt start1, start2

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  call F90Array2dGetAddrReal(ptr(start1, start2), address)
end subroutine

subroutine F90Array2dAccessInt(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:, :)
  PetscFortranAddr address
  PetscInt start1, start2

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  call F90Array2dGetAddrInt(ptr(start1, start2), address)
end subroutine

subroutine F90Array2dAccessFortranAddr(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:, :)
  PetscFortranAddr address
  PetscInt start1, start2

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  call F90Array2dGetAddrFortranAddr(ptr(start1, start2), address)
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array2dDestroyScalar(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :)

  nullify (ptr)
end subroutine

subroutine F90Array2dDestroyReal(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:, :)

  nullify (ptr)
end subroutine

subroutine F90Array2dDestroyInt(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:, :)

  nullify (ptr)
end subroutine

subroutine F90Array2dDestroyFortranAddr(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:, :)

  nullify (ptr)
end subroutine
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array3dCreateScalar(array, start1, len1, start2, len2, start3, len3, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscScalar, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1)
  PetscScalar, pointer :: ptr(:, :, :)

  ptr => array
end subroutine

subroutine F90Array3dCreateReal(array, start1, len1, start2, len2, start3, len3, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscReal, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1)
  PetscReal, pointer :: ptr(:, :, :)

  ptr => array
end subroutine

subroutine F90Array3dCreateInt(array, start1, len1, start2, len2, start3, len3, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscInt, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1)
  PetscInt, pointer :: ptr(:, :, :)

  ptr => array
end subroutine

subroutine F90Array3dCreateFortranAddr(array, start1, len1, start2, len2, start3, len3, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscFortranAddr, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1)
  PetscFortranAddr, pointer :: ptr(:, :, :)

  ptr => array
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array3dAccessScalar(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  call F90Array3dGetAddrScalar(ptr(start1, start2, start3), address)
end subroutine

subroutine F90Array3dAccessReal(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  call F90Array3dGetAddrReal(ptr(start1, start2, start3), address)
end subroutine

subroutine F90Array3dAccessInt(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  call F90Array3dGetAddrInt(ptr(start1, start2, start3), address)
end subroutine

subroutine F90Array3dAccessFortranAddr(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  call F90Array3dGetAddrFortranAddr(ptr(start1, start2, start3), address)
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array3dDestroyScalar(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :, :)

  nullify (ptr)
end subroutine

subroutine F90Array3dDestroyReal(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:, :, :)

  nullify (ptr)
end subroutine

subroutine F90Array3dDestroyInt(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:, :, :)

  nullify (ptr)
end subroutine

subroutine F90Array3dDestroyFortranAddr(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:, :, :)

  nullify (ptr)
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
subroutine F90Array4dCreateScalar(array, start1, len1, start2, len2, start3, len3, start4, len4, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscInt start4, len4
  PetscScalar, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1, start4:start4 + len4 - 1)
  PetscScalar, pointer :: ptr(:, :, :, :)

  ptr => array
end subroutine

subroutine F90Array4dCreateReal(array, start1, len1, start2, len2, start3, len3, start4, len4, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscInt start4, len4
  PetscReal, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1, start4:start4 + len4 - 1)
  PetscReal, pointer :: ptr(:, :, :, :)

  ptr => array
end subroutine

subroutine F90Array4dCreateInt(array, start1, len1, start2, len2, start3, len3, start4, len4, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscInt start4, len4
  PetscInt, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1, start4:start4 + len4 - 1)
  PetscInt, pointer :: ptr(:, :, :, :)

  ptr => array
end subroutine

subroutine F90Array4dCreateFortranAddr(array, start1, len1, start2, len2, start3, len3, start4, len4, ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt start1, len1
  PetscInt start2, len2
  PetscInt start3, len3
  PetscInt start4, len4
  PetscFortranAddr, target :: array(start1:start1 + len1 - 1, start2:start2 + len2 - 1, start3:start3 + len3 - 1, start4:start4 + len4 - 1)
  PetscFortranAddr, pointer :: ptr(:, :, :, :)

  ptr => array
end subroutine

subroutine F90Array4dAccessScalar(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3, start4

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  start4 = lbound(ptr, 4)
  call F90Array4dGetAddrScalar(ptr(start1, start2, start3, start4), address)
end subroutine

subroutine F90Array4dAccessReal(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:, :, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3, start4

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  start4 = lbound(ptr, 4)
  call F90Array4dGetAddrReal(ptr(start1, start2, start3, start4), address)
end subroutine

subroutine F90Array4dAccessInt(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:, :, :, :)
  PetscFortranAddr address
  PetscInt start1, start2, start3, start4

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  start4 = lbound(ptr, 4)
  call F90Array4dGetAddrInt(ptr(start1, start2, start3, start4), address)
end subroutine

subroutine F90Array4dAccessFortranAddr(ptr, address)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :, :, :)
  PetscFortranAddr address
  PetscFortranAddr start1, start2, start3, start4

  start1 = lbound(ptr, 1)
  start2 = lbound(ptr, 2)
  start3 = lbound(ptr, 3)
  start4 = lbound(ptr, 4)
  call F90Array4dGetAddrFortranAddr(ptr(start1, start2, start3, start4), address)
end subroutine

subroutine F90Array4dDestroyScalar(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscScalar, pointer :: ptr(:, :, :, :)

  nullify (ptr)
end subroutine

subroutine F90Array4dDestroyReal(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscReal, pointer :: ptr(:, :, :, :)

  nullify (ptr)
end subroutine

subroutine F90Array4dDestroyInt(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscInt, pointer :: ptr(:, :, :, :)

  nullify (ptr)
end subroutine

subroutine F90Array4dDestroyFortranAddr(ptr)
  use, intrinsic :: ISO_C_binding
  implicit none
  PetscFortranAddr, pointer :: ptr(:, :, :, :)

  nullify (ptr)
end subroutine

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!