various fixes in examples

examples/arkode/F2003_serial/ark_kpr_mri_f2003.f90

@@ -71,14 +71,13 @@
 ! at the end.
 !--------------------------------------------------------------------
 
-module ode_mod
+module ark_kpr_mri_mod
 
   use, intrinsic :: iso_c_binding
   use fsundials_core_mod
   use farkode_mod
   use farkode_arkstep_mod
   use farkode_mristep_mod
-  use fnvector_serial_mod
   use fsunmatrix_dense_mod
   use fsunlinsol_dense_mod
 
@@ -556,9 +555,17 @@ integer(c_int) function Ytrue(t, y) &
 
 program main
 
-  use ode_mod
-  implicit none
+  use, intrinsic :: iso_c_binding
+  use fsundials_core_mod
+  use farkode_mod
+  use farkode_arkstep_mod
+  use farkode_mristep_mod
+  use fnvector_serial_mod
+  use fsunmatrix_dense_mod
+  use fsunlinsol_dense_mod
+  use ark_kpr_mri_mod
 
+  implicit none
 
   ! general problem variables
   type(c_ptr) sunctx                       ! SUNDIALS simulation context
@@ -666,14 +673,14 @@ program main
       print *, "   solver: none/exp-3 (no slow, explicit fast)"
     case (2)
       reltol = max(hs*hs*hs, real(1e-10,8))
-      abstol = 1e-11
+      abstol = real(1e-11,8)
       print *, "   solver: none/dirk-3 (no slow, dirk fast)"
       print '(A,E12.4,A,E12.4)', "    reltol: ", reltol, " abstol: ", abstol
     case (3)
       print *, "   solver: exp-3/none (explicit slow, no fast)"
     case (4)
       reltol = max(hs*hs, real(1e-10,8))
-      abstol = 1e-11
+      abstol = real(1e-11,8)
       print *, "   solver: dirk-2/none (dirk slow, no fast)"
       print '(A,E12.4,A,E12.4)', "    reltol: ", reltol, " abstol: ", abstol
     case (5)
@@ -682,17 +689,17 @@ program main
       print *, "   solver: exp-4/exp-3 (MRI-GARK-ERK45a / ERK-3-3)"
     case (7)
       reltol = max(hs*hs*hs, real(1e-10,8))
-      abstol = 1e-11
+      abstol = real(1e-11,8)
       print *, "   solver: dirk-3/exp-3 (MRI-GARK-ESDIRK34a / ERK-3-3) -- solve decoupled"
       print '(A,E12.4,A,E12.4)', "    reltol: ", reltol, " abstol: ", abstol
     case (8)
       reltol = max(hs*hs*hs, real(1e-10,8))
-      abstol = 1e-11
+      abstol = real(1e-11,8)
       print *, "   solver: ars343/exp-3 (IMEX-MRI3b / ERK-3-3) -- solve decoupled"
       print '(A,E12.4,A,E12.4)', "    reltol: ", reltol, " abstol: ", abstol
     case (9)
       reltol = max(hs*hs*hs*hs, real(1e-14,8))
-      abstol = 1e-14
+      abstol = real(1e-14,8)
       print *, "   solver: imexark4/exp-4 (IMEX-MRI4 / ERK-4-4) -- solve decoupled"
       print '(A,E12.4,A,E12.4)', "    reltol: ", reltol, " abstol: ", abstol
   end select

examples/arkode/F2003_serial/test_ark_butcher_f2003.f90

@@ -28,7 +28,7 @@ integer function smoke_tests() result(ret)
     implicit none
     type(c_ptr) :: ERK, DIRK
     integer(c_int)  :: ierr, q(1), p(1)
-    integer(c_long_long) :: liw(1), lrw(1)
+    integer(c_long) :: liw(1), lrw(1)
     real(c_double)  :: b(2), c(2), d(2), A(4)
 
     !===== Setup ====

examples/idas/F2003_serial/idasAkzoNob_ASAi_dns_f2003.f90

@@ -408,7 +408,7 @@ program main
   retval = FIDASStolerancesB(mem, indexB(1), RTOLB, ATOLB)
   call check_retval(retval, "FIDASStolerancesB")
 
-  retval = FIDASetMaxNumStepsB(mem, indexB(1), 1000_8)
+  retval = FIDASetMaxNumStepsB(mem, indexB(1), 1000)
   call check_retval(retval, "FIDASetMaxNumStepsB")
 
   ! Create dense SUNMatrix for use in linear solves

examples/sunlinsol/test_sunlinsol.f90

@@ -83,7 +83,7 @@ integer(C_INT) function Test_FSUNLinSolSpace(S, myid) result(failure)
 
     type(SUNLinearSolver), pointer :: S
     integer(C_INT)                 :: myid
-    integer(c_int64_t)                :: lenrw(1), leniw(1)
+    integer(C_LONG)                :: lenrw(1), leniw(1)
 
     failure = 0
 

examples/sunmatrix/band/test_fsunmatrix_band_mod.f90

@@ -42,7 +42,7 @@ integer(C_INT) function smoke_tests() result(fails)
     type(SUNMatrix), pointer :: A, B               ! SUNMatrix
     type(N_Vector),  pointer :: x, y               ! NVectors
     real(C_DOUBLE),  pointer :: matdat(:)          ! matrix data pointer
-    integer(c_int64_t)          :: lenrw(1), leniw(1) ! matrix real and int work space size
+    integer(C_LONG)          :: lenrw(1), leniw(1) ! matrix real and int work space size
     integer(c_int64_t)          :: val
     type(C_PTR),     pointer :: cptr
 

examples/sunmatrix/dense/test_fsunmatrix_dense_mod.f90

@@ -40,8 +40,8 @@ integer(C_INT) function smoke_tests() result(fails)
     type(SUNMatrix), pointer :: A, B               ! SUNMatrix
     type(N_Vector),  pointer :: x, y               ! NVectors
     real(C_DOUBLE),  pointer :: matdat(:)          ! matrix data pointer
-    integer(c_int64_t)          :: lenrw(1), leniw(1) ! matrix real and int work space size
-    integer(c_int64_t)          :: val
+    integer(C_LONG)          :: lenrw(1), leniw(1) ! matrix real and int work space size
+    integer(c_int64_t)       :: val
 
     fails = 0
 

examples/sunmatrix/sparse/test_fsunmatrix_sparse_mod.f90

@@ -41,7 +41,7 @@ integer(C_INT) function smoke_tests() result(fails)
     type(N_Vector),  pointer      :: x, y
     real(C_DOUBLE),  pointer      :: matdat(:)          ! matrix data pointer
     integer(c_int64_t), pointer :: inddat(:)          ! indices data pointer
-    integer(c_int64_t)               :: lenrw(1), leniw(1) ! matrix real and int work space size
+    integer(C_LONG)               :: lenrw(1), leniw(1) ! matrix real and int work space size
 
     integer(c_int64_t) :: tmp1
     integer(C_INT)  :: tmp2

examples/sunmatrix/test_sunmatrix.f90

@@ -473,7 +473,7 @@ integer(C_INT) function Test_FSUNMatSpace(A, myid) result(failure)
 
     integer(C_INT)  :: myid
     type(SUNMatrix) :: A
-    integer(c_int64_t) :: lenrw(1), leniw(1)
+    integer(C_LONG) :: lenrw(1), leniw(1)
 
     failure = 0
 

examples/sunnonlinsol/fixedpoint/test_fsunnonlinsol_fixedpoint_mod.f90

@@ -45,7 +45,7 @@ integer(C_INT) function unit_tests() result(retval)
     type(SUNNonlinearSolver), pointer :: NLS         ! test nonlinear solver
     type(N_Vector),           pointer :: x, y0, y, w ! test vectors
     real(C_DOUBLE),           pointer :: ydata(:)
-    integer(c_int64_t)                   :: niters(1)
+    integer(C_LONG)                   :: niters(1)
     integer(C_INT)                    :: tmp
 
     x  => FN_VNew_Serial(NEQ, sunctx)

examples/sunnonlinsol/newton/test_fsunnonlinsol_newton_mod.f90

@@ -57,7 +57,7 @@ integer(C_INT) function unit_tests() result(retval)
 
     type(SUNNonlinearSolver), pointer :: NLS        ! test nonlinear solver
     real(C_DOUBLE),           pointer :: ydata(:)
-    integer(c_int64_t)                   :: niters(1)
+    integer(C_LONG)                   :: niters(1)
     integer(C_INT)                    :: tmp
     type(IntegratorMem),      pointer :: Imem