diff --git a/examples/arkode/F2003_serial/ark_bruss1D_FEM_klu_f2003.f90 b/examples/arkode/F2003_serial/ark_bruss1D_FEM_klu_f2003.f90
index b6799607ec..31bdd7cbdd 100644
--- a/examples/arkode/F2003_serial/ark_bruss1D_FEM_klu_f2003.f90
+++ b/examples/arkode/F2003_serial/ark_bruss1D_FEM_klu_f2003.f90
@@ -1091,28 +1091,28 @@ program main
   nout   = ceiling(tend/dtout)
 
   ! create and assign SUNDIALS N_Vectors
-  sunvec_y => FN_VNew_Serial(neq, ctx)
+  sunvec_y => FN_VNew_Serial(int(neq,c_long), ctx)
   if (.not. associated(sunvec_y)) then
      print *, 'ERROR: sunvec = NULL'
      stop 1
   end if
   yvec(1:neqreal,1:N)  => FN_VGetArrayPointer(sunvec_y)
 
-  sunvec_u => FN_VNew_Serial(neq, ctx)
+  sunvec_u => FN_VNew_Serial(int(neq,c_long), ctx)
   if (.not. associated(sunvec_u)) then
      print *, 'ERROR: sunvec = NULL'
      stop 1
   end if
   umask(1:neqreal,1:N) => FN_VGetArrayPointer(sunvec_u)
 
-  sunvec_v => FN_VNew_Serial(neq, ctx)
+  sunvec_v => FN_VNew_Serial(int(neq,c_long), ctx)
   if (.not. associated(sunvec_v)) then
      print *, 'ERROR: sunvec = NULL'
      stop 1
   end if
   vmask(1:neqreal,1:N) => FN_VGetArrayPointer(sunvec_v)
 
-  sunvec_w => FN_VNew_Serial(neq, ctx)
+  sunvec_w => FN_VNew_Serial(int(neq,c_long), ctx)
   if (.not. associated(sunvec_w)) then
      print *, 'ERROR: sunvec = NULL'
      stop 1
@@ -1158,13 +1158,13 @@ program main
 
   ! Tell ARKODE to use a sparse linear solver for both Newton and mass matrix systems.
   sparsetype = 1
-  sunmat_A => FSUNSparseMatrix(neq, neq, nnz, sparsetype, ctx)
+  sunmat_A => FSUNSparseMatrix(int(neq,c_long), int(neq,c_long), int(nnz,c_long), sparsetype, ctx)
   if (.not. associated(sunmat_A)) then
      print *, 'ERROR: sunmat_A = NULL'
      stop 1
   end if
 
-  sunmat_M => FSUNSparseMatrix(neq, neq, nnz, sparsetype, ctx)
+  sunmat_M => FSUNSparseMatrix(int(neq,c_long), int(neq,c_long), int(nnz,c_long), sparsetype, ctx)
   if (.not. associated(sunmat_M)) then
      print *, 'ERROR: sunmat_M = NULL'
      stop 1
diff --git a/examples/arkode/F2003_serial/ark_diurnal_kry_bp_f2003.f90 b/examples/arkode/F2003_serial/ark_diurnal_kry_bp_f2003.f90
index ba295a786a..3bc9da629f 100644
--- a/examples/arkode/F2003_serial/ark_diurnal_kry_bp_f2003.f90
+++ b/examples/arkode/F2003_serial/ark_diurnal_kry_bp_f2003.f90
@@ -86,11 +86,11 @@ module DiurnalKryBP_mod
   integer(c_long), parameter :: mxsteps = 10000
 
   ! ODE non-constant parameters
-  real(c_double)  :: q3
-  real(c_double)  :: q4
-  real(c_double)  :: c1
-  real(c_double)  :: c2
-  integer(c_long) :: jx, jy
+  real(c_double) :: q3
+  real(c_double) :: q4
+  real(c_double) :: c1
+  real(c_double) :: c2
+  integer(c_int) :: jx, jy
 
 contains
 
diff --git a/examples/arkode/F2003_serial/ark_roberts_dnsL_f2003.f90 b/examples/arkode/F2003_serial/ark_roberts_dnsL_f2003.f90
index 2aedc669ec..fc10c55300 100644
--- a/examples/arkode/F2003_serial/ark_roberts_dnsL_f2003.f90
+++ b/examples/arkode/F2003_serial/ark_roberts_dnsL_f2003.f90
@@ -508,7 +508,7 @@ subroutine PrintOutput(arkode_mem, t, y)
   real(c_double) :: t, y(neq)
 
   ! internal variables
-  integer(c_int)  :: retval, kused(1)
+  integer(c_int)  :: retval
   integer(c_long) :: nst(1)
   real(c_double)  :: hused(1)