1 files changed, 87 insertions, 0 deletions
diff --git a/Ganlib/src/KDIOP.f90 b/Ganlib/src/KDIOP.f90
new file mode 100644
index 0000000..0af8d6a
--- /dev/null
+++ b/Ganlib/src/KDIOP.f90
@@ -0,0 +1,87 @@
+!
+!-----------------------------------------------------------------------
+!
+!Purpose:
+! Fortran-2003 bindings for kdi.
+!
+!Copyright:
+! Copyright (C) 2009 Ecole Polytechnique de Montreal
+! This library is free software; you can redistribute it and/or
+! modify it under the terms of the GNU Lesser General Public
+! License as published by the Free Software Foundation; either
+! version 2.1 of the License, or (at your option) any later version.
+!
+!Author(s): A. Hebert
+!
+!-----------------------------------------------------------------------
+!
+function KDIOP(name, iactio)
+   ! open a KDI file
+   use, intrinsic :: iso_c_binding
+   use LCMAUX
+   type(c_ptr) KDIOP
+   character(len=*) :: name
+   integer :: iactio
+   character(kind=c_char), dimension(13) :: name13
+   interface
+      function kdiop_c (name_c, iactio) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) kdiop_c
+         character(kind=c_char), dimension(*) :: name_c
+         integer(c_int), value :: iactio
+      end function kdiop_c
+   end interface
+   call STRCUT(name13, name)
+   KDIOP=kdiop_c(name13, iactio)
+end function KDIOP
+!
+subroutine KDIPUT(my_file, idata, iofset, length)
+   ! store a data array in a KDI file at offset iofset
+   use, intrinsic :: iso_c_binding
+   type(c_ptr) :: my_file, pt_data
+   integer, target, dimension(*) :: idata
+   integer :: iofset, length
+   interface
+      subroutine kdiput_c (my_file, idata, iofset, length) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr), value :: my_file, idata
+         integer(c_int), value :: iofset, length
+      end subroutine kdiput_c
+   end interface
+   pt_data=c_loc(idata)
+   call kdiput_c(my_file, pt_data, iofset, length)
+end subroutine KDIPUT
+!
+subroutine KDIGET(my_file, idata, iofset, length)
+   ! read a data array from a KDI file at offset iofset
+   use, intrinsic :: iso_c_binding
+   type(c_ptr) :: my_file, pt_data
+   integer, target, dimension(*) :: idata
+   integer :: iofset, length
+   interface
+      subroutine kdiget_c (my_file, idata, iofset, length) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr), value :: my_file, idata
+         integer(c_int), value :: iofset, length
+      end subroutine kdiget_c
+   end interface
+   pt_data=c_loc(idata)
+   call kdiget_c(my_file, pt_data, iofset, length)
+end subroutine KDIGET
+!
+function KDICL(my_file, istatu)
+   ! close a KDI file
+   use, intrinsic :: iso_c_binding
+   integer(c_int) KDICL
+   type(c_ptr) :: my_file
+   integer :: istatu
+   interface
+      function kdicl_c (my_file, istatu) bind(c)
+         use, intrinsic :: iso_c_binding
+         integer(c_int) kdicl_c
+         type(c_ptr), value :: my_file
+         integer(c_int), value :: istatu
+      end function kdicl_c
+   end interface
+   KDICL=kdicl_c(my_file, istatu)
+end function KDICL