[svn-r2576] Purpose:

	Adding the Fortran interface to the HDF5 library
Description:
	Fortran is now a subdirectory of the HDF5 library tree.
Platforms tested:
	Solaris and IRIX (O2K)

1 files changed, 159 insertions, 0 deletions

diff --git a/fortran/examples/refregexample.f90 b/fortran/examples/refregexample.f90
new file mode 100644
index 0000000..47ebaac
--- /dev/null
+++ b/fortran/examples/refregexample.f90
@@ -0,0 +1,159 @@
+!
+!    This program shows how to create, store and dereference references 
+!    to the dataset regions.
+!    Program creates a file and writes two dimensional integer dataset
+!    to it. Then program creates and stores references to the hyperslab
+!    and 3 points selected in the integer dataset, in the second dataset. 
+!    Program reopens the second dataset, reads and dereferences region
+!    references, and then reads and displays selected data from the 
+!    integer dataset.
+!    
+     PROGRAM REG_REFERENCE
+
+        USE HDF5 ! This module contains all necessary modules 
+        
+     IMPLICIT NONE
+     CHARACTER(LEN=10), PARAMETER :: filename = "FORTRAN.h5"
+     CHARACTER(LEN=6), PARAMETER :: dsetnamev = "MATRIX"
+     CHARACTER(LEN=17), PARAMETER :: dsetnamer = "REGION_REFERENCES"
+      
+     INTEGER(HID_T) :: file_id       ! File identifier 
+     INTEGER(HID_T) :: space_id      ! Dataspace identifier 
+     INTEGER(HID_T) :: spacer_id     ! Dataspace identifier 
+     INTEGER(HID_T) :: dsetv_id      ! Dataset identifier 
+     INTEGER(HID_T) :: dsetr_id      ! Dataset identifier 
+     INTEGER     ::   error
+     TYPE(hdset_reg_ref_t_f) , DIMENSION(2) :: ref     ! Buffers to store references
+     TYPE(hdset_reg_ref_t_f) , DIMENSION(2) :: ref_out !
+     INTEGER(HSIZE_T), DIMENSION(2) :: dims = (/2,9/)  ! Datasets dimensions
+     INTEGER(HSIZE_T), DIMENSION(1) :: dimsr = (/2/)   ! 
+     INTEGER(HSSIZE_T), DIMENSION(2) :: start
+     INTEGER(HSIZE_T), DIMENSION(2) :: count
+     INTEGER :: rankr = 1 
+     INTEGER :: rank = 2
+     INTEGER , DIMENSION(2,9) ::  data 
+     INTEGER , DIMENSION(2,9) ::  data_out = 0 
+     INTEGER(HSSIZE_T) , DIMENSION(2,3) :: coord
+     INTEGER(SIZE_T) ::num_points = 3  ! Number of selected points
+     INTEGER :: i, j
+     coord = reshape((/1,1,2,7,1,9/), (/2,3/))   ! Coordinates of selected points
+     data = reshape ((/1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6/), (/2,9/))
+     !
+     !  Initialize FORTRAN predefined datatypes.
+     !
+     CALL h5init_types_f(error)
+     !
+     !  Create a new file.
+     !
+     CALL h5fcreate_f(filename, H5F_ACC_TRUNC_F, file_id, error)
+                                         ! Default file access and file creation
+                                         ! properties are used. 
+     ! 
+     ! Create  dataspaces:
+     ! 
+     ! for dataset with references to dataset regions 
+     !
+     CALL h5screate_simple_f(rankr, dimsr, spacer_id, error)
+     !
+     ! for integer dataset 
+     !
+     CALL h5screate_simple_f(rank, dims, space_id, error)
+     !
+     ! Create  and write datasets:
+     !
+     ! Integer dataset 
+     !
+     CALL h5dcreate_f(file_id, dsetnamev, H5T_NATIVE_INTEGER, space_id, &
+                      dsetv_id, error)
+     CALL h5dwrite_f(dsetv_id, H5T_NATIVE_INTEGER, data, error)
+     CALL h5dclose_f(dsetv_id, error)
+     !
+     ! Dataset with references
+     !
+     CALL h5dcreate_f(file_id, dsetnamer, H5T_STD_REF_DSETREG, spacer_id, &
+                      dsetr_id, error)
+     !
+     ! Create a reference to the hyperslab selection.
+     !
+     start(1) = 0 
+     start(2) = 3 
+     count(1) = 2
+     count(2) = 3
+     CALL h5sselect_hyperslab_f(space_id, H5S_SELECT_SET_F, &
+                                start, count, error) 
+     CALL h5rcreate_f(file_id, dsetnamev, space_id, ref(1), error) 
+     !
+     ! Create a reference to elements selection.
+     !
+     CALL h5sselect_none_f(space_id, error)
+     CALL h5sselect_elements_f(space_id, H5S_SELECT_SET_F, rank, num_points,&
+                               coord, error) 
+     CALL h5rcreate_f(file_id, dsetnamev, space_id, ref(2), error) 
+     !
+     ! Write dataset with the references. 
+     !
+     CALL h5dwrite_f(dsetr_id, H5T_STD_REF_DSETREG, ref, error) 
+     !
+     ! Close all objects.
+     !
+     CALL h5sclose_f(space_id, error)
+     CALL h5sclose_f(spacer_id, error)
+     CALL h5dclose_f(dsetr_id, error)
+     CALL h5fclose_f(file_id, error)
+     !
+     ! Reopen the file to test selections.
+     !
+     CALL h5fopen_f (filename, H5F_ACC_RDWR_F, file_id, error)
+     CALL h5dopen_f(file_id, dsetnamer, dsetr_id, error)
+     !
+     ! Read references to the dataset regions.
+     !
+     CALL h5dread_f(dsetr_id, H5T_STD_REF_DSETREG, ref_out, error) 
+     ! 
+     ! Dereference the first reference.
+     ! 
+     CALL H5rdereference_f(dsetr_id, ref_out(1), dsetv_id, error)
+     CALL H5rget_region_f(dsetr_id, ref_out(1), space_id, error) 
+     !
+     ! Read selected data from the dataset.
+     !
+     CALL h5dread_f(dsetv_id, H5T_NATIVE_INTEGER, data_out, error, &
+                     mem_space_id = space_id, file_space_id = space_id)
+          write(*,*) "Hypeslab selection"
+          write(*,*)
+          do i = 1,2 
+          write(*,*) (data_out (i,j), j = 1,9)
+          enddo
+          write(*,*)
+     CALL h5sclose_f(space_id, error)
+     CALL h5dclose_f(dsetv_id, error)
+     data_out = 0
+     !  
+     ! Dereference the second reference.
+     ! 
+     CALL H5rdereference_f(dsetr_id, ref_out(2), dsetv_id, error)
+     CALL H5rget_region_f(dsetr_id, ref_out(2), space_id, error) 
+     !
+     ! Read selected data from the dataset.
+     !
+     CALL h5dread_f(dsetv_id, H5T_NATIVE_INTEGER, data_out, error, &
+                     mem_space_id = space_id, file_space_id = space_id)
+          write(*,*) "Point selection"
+          write(*,*)
+          do i = 1,2 
+          write(*,*) (data_out (i,j), j = 1,9)
+          enddo
+     !
+     ! Close all objects
+     !
+     CALL h5sclose_f(space_id, error)
+     CALL h5dclose_f(dsetv_id, error)
+     CALL h5dclose_f(dsetr_id, error)
+     !
+     ! Close FORTRAN predefined datatypes.
+     !
+     CALL h5close_types_f(error) 
+
+     END PROGRAM REG_REFERENCE 
+
+