about SHUFFLE+GZIP

hi,
i have met a hdf5 file, which appears as follows:
   Compression:SHUFFLE:Nbytes=4,GIZP:level=6, Storage allocation time: Incremental

which function can i use to realize the above "SHUFFLE+GZIP" format in my file with fortran code ?

any help will be appreciated,
thanks

dyjbean