PROGRAM MYCHUNK USE HDF5 ! This module contains all necessary modules IMPLICIT NONE INTEGER(HID_T) :: file_id ! File identifier INTEGER(HID_T) :: dset_id ! Dataset identifier INTEGER(HID_T) :: filespace ! Dataspace identifier in file INTEGER(HID_T) :: dataspace ! Dataspace identifier in file INTEGER(HID_T) :: crp_list ! dataset creatation property identifier INTEGER :: error ! Error flag !Maximum dimensions ! INTEGER(HSIZE_T), DIMENSION(2) :: maxdims INTEGER(HSIZE_T), DIMENSION(2) :: dims ! !Variables used in reading data back ! INTEGER(HSIZE_T), DIMENSION(2) :: chunk_dims = (/12,12/) INTEGER :: RANK = 2 !dataset name is "ExtendibleArray" ! CHARACTER(LEN=25), PARAMETER :: dsetname = "ExtendibleVector" INTEGER, allocatable, dimension(:,:) :: vampnew integer :: i1reach = 100 allocate(vampnew(100,100)) vampnew = 1 CALL h5open_f(error) CALL h5fcreate_f("smwf.h5", H5F_ACC_TRUNC_F, file_id, error) ! Create the data space for the dataset. maxdims = (/H5S_UNLIMITED_F, H5S_UNLIMITED_F/) dims(1)=i1reach dims(2)=i1reach CALL h5screate_simple_f(rank,dims,dataspace,error,maxdims) CALL h5pcreate_f(H5P_DATASET_CREATE_F,crp_list,error) CALL h5pset_chunk_f(crp_list,2,chunk_dims,error) CALL h5dcreate_f(file_id,dsetname,H5T_NATIVE_INTEGER,& dataspace,dset_id,error,crp_list) !CALL h5dextend_f(dset_id, dims, error) CALL h5dwrite_f(dset_id, H5T_NATIVE_INTEGER, vampnew,& dims,error,H5S_ALL_F,H5S_ALL_F,H5P_DEFAULT_F) CALL h5sclose_f(dataspace, error) CALL h5dclose_f(dset_id, error) CALL h5pclose_f(crp_list, error) CALL h5fclose_f(file_id, error) END PROGRAM MYCHUNK