!
! This example shows how to create datasets with region and object references to
! datasets located under the group /Data_Products/VIIRS-MOD-GTM-EDR-GEO
!   Main illustrative functions: H5LRcreate_ref_to_all_f
!
PROGRAM main
  
  USE ISO_C_BINDING
  USE HDF5 ! module of HDF5 library
  USE H5LT ! module of H5LT
  USE H5HL_REGION

  IMPLICIT NONE

  CHARACTER(LEN=78), PARAMETER :: filename = "GMGTO_npp_raw.h5"
  CHARACTER(LEN=65), PARAMETER :: Dataset_Gran_0 = "/Data_Products/VIIRS-MOD-GTM-EDR-GEO/VIIRS-MOD-GTM-EDR-GEO_Gran_0"
  CHARACTER(LEN=35), PARAMETER :: group = "/All_Data/VIIRS-MOD-GTM-EDR-GEO_All"
  CHARACTER(LEN=63), PARAMETER :: Dataset_Aggr = "/Data_Products/VIIRS-MOD-GTM-EDR-GEO/VIIRS-MOD-GTM-EDR-GEO_Aggr"

  INTEGER(hid_t) :: file_id  ! file identifier
  INTEGER :: status
  !
  !  Initialize FORTRAN interface.
  !
  CALL H5open_f(status)
  !
  ! Open the NPP file.
  !
  CALL H5Fopen_f(filename, H5F_ACC_RDWR_F, file_id, status)
  !
  !  Create a dataset with region references to all the datasets located under /All_Data/VIIRS-MOD-GTM-EDR-GEO_All 
  !  Store the region references in the /Data_Products/VIIRS-MOD-GTM-EDR-GEO/VIIRS-MOD-GTM-EDR-GEO_Gan_0 dataset.
  !
  CALL H5LRcreate_ref_to_all_f(file_id, group, Dataset_Gran_0, H5_INDEX_NAME_F, H5_ITER_INC_F, H5R_DATASET_REGION_F, status)
  !
  !  Create a dataset with object references to all the datasets located under /All_Data/VIIRS-MOD-GTM-EDR-GEO_All 
  !  Store the region references in the /Data_Products/VIIRS-MOD-GTM-EDR-GEO/VIIRS-MOD-GTM-EDR-GEO_Aggr dataset.
  !
  CALL H5LRcreate_ref_to_all_f(file_id, group, Dataset_Aggr, H5_INDEX_NAME_F, H5_ITER_INC_F, H5R_OBJECT_F, status);

  CALL H5Fclose_f(file_id, status)

END PROGRAM main