diff --git a/examples/H5Part/Bench.c b/examples/H5Part/Bench.c
deleted file mode 100644
index c5e4499..0000000
--- a/examples/H5Part/Bench.c
+++ /dev/null
@@ -1,263 +0,0 @@
-#include <stdio.h>
-#include <stdlib.h>
-#include <string.h>
-#include <mpi.h>
-/* #include <mpio.h> */
-#include <unistd.h>
-#include <sys/types.h>
-#ifndef PARALLEL_IO
-#define PARALLEL_IO
-#endif
-
-#ifndef DISABLE_H5PART
-#include "H5hut.h"
-#endif
-
-#define FILENAME "testio"
-/* normally 64 steps for real benchmark */
-/* #define NSTEPS 5 */
- 
-/* normally 51e6 for real benchmark */
-#define NPARTICLES 51e4
-#define NTRIALS 3
-
-/*
-
-bench <nParticles> <nSteps>
-
-*/
-
-
-int main(int argc,char *argv[]){
-
-  if (argc < 3) {
-    printf("Usage: bench <nParticles> <nSteps> \n");
-    exit(-1);
-  }
-  else {
-    printf("nparticles: %d, nsteps: %d \n", atoi(argv[1]), atoi(argv[2]));
-  }
-
-  MPI_Info info;
-  int nprocs,rank;
-  int trial;
-  int i,j,n; /* iteration variables */
-  double starttime,curtime, endtime;
-
-  int nparticles = atoi(argv[1]);
-  int nsteps = atoi(argv[2]);
-
-  double *x,*y,*z,*px,*py,*pz;
-  typedef double *ddouble;
-  ddouble data[6];
-  MPI_Datatype chunktype;
-  int offset;
-  int localnp;
-  char filename[128]; /*= FILENAME; */
-#ifndef DISABLE_H5PART
-  h5_file_t f;
-#endif
-  char newfilename[128];
-  FILE *fd;
-  MPI_File file;
-  MPI_Offset foffset;
-
-  MPI_Comm dcomm = MPI_COMM_WORLD;
-  
-  MPI_Init(&argc,&argv);
-  MPI_Comm_rank(dcomm,&rank);
-  MPI_Comm_size(dcomm,&nprocs);
-
-  localnp=nparticles/(int64_t)nprocs;
-  for(offset=0,i=0;i<rank;i++){
-    offset+=localnp;
-  }
-  
-  data[0]=x=(double*)malloc(sizeof(double)*(size_t)localnp);
-  data[1]=y=(double*)malloc(sizeof(double)*(size_t)localnp);
-  data[2]=z=(double*)malloc(sizeof(double)*(size_t)localnp);
-  data[3]=px=(double*)malloc(sizeof(double)*(size_t)localnp);
-  data[4]=py=(double*)malloc(sizeof(double)*(size_t)localnp);
-  data[5]=pz=(double*)malloc(sizeof(double)*(size_t)localnp);
-
-  
-  /* printf("about to call create subarray with nparticles=%u localnp=%u offset=%u\n",
-     nparticles,localnp,offset); */
-  MPI_Type_create_subarray(1, /* rank */
-			   &nparticles, /* size of the global array */
-			   &localnp, /* size of my local chunk */
-			   &offset, /* offset of this chunk in global */
-			   MPI_ORDER_FORTRAN, /* fortran storage order */
-			   MPI_DOUBLE,
-			   &chunktype);
-  MPI_Type_commit(&chunktype);
-  MPI_Info_create(&info);
-  MPI_Info_set(info, "IBM_largeblock_io", "true" );
-
-  if(rank==0) printf("Nprocs=%u Particles=%u*6attribs*sizeof(double) Particles/proc=%u Nsteps=%u Ntrials=%u\n",
-		     nprocs,nparticles,localnp,nsteps,NTRIALS);
-	
-
-  for(trial=0;trial<NTRIALS;trial++){
-    if(rank==0) printf("---------------------- Trial %u of %u ---------------------\n",trial+1,NTRIALS);
-
-
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-    sprintf(filename,"%s.%u.mpio.dat",FILENAME,nprocs);
-
-    if(rank==0) unlink(filename);
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-
-    MPI_File_open(MPI_COMM_WORLD,filename,
-		  MPI_MODE_CREATE | MPI_MODE_RDWR,
-		  info,&file);
-	
-    MPI_File_set_view(file,0,MPI_DOUBLE,chunktype,"native",info);
-    /* now a barrier to get the start timers roughly synced*/
-    MPI_Barrier(MPI_COMM_WORLD);
-    curtime = starttime = MPI_Wtime();
-    endtime = starttime+5.0*60.0; /* end in 5 minutes */
-    MPI_Bcast(&endtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-    /* must touch the entire array after each write */
-    /* ensures cache-invalidation */
-    foffset=0;
-    i=0;
-    curtime=starttime;
-    for(i=0;i<nsteps;i++){
-      int n;
-      MPI_Status status;
-      for(j=0;j<6;j++){
-	/* touch data */
-	for(n=0;n<localnp;n++)
-	  (data[j])[n]=(double)rank;
-	/* write to that file */
-	/*  MPI_File_set_view(file,foffset,MPI_DOUBLE,chunktype,"native",info);*/
-
-#ifdef COLLECTIVE_IO
-	MPI_File_write_at_all(file,
-			      foffset,
-			      data[j],
-			      localnp,
-			      MPI_DOUBLE,&status);
-#else
-	MPI_File_write_at(file,
-			  foffset,
-			  data[j],
-			  localnp,
-			  MPI_DOUBLE,&status);
-#endif
-
-	foffset+=nparticles/nprocs;
-      }
-      curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
-      MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-    }
-
-
-    MPI_File_close(&file);
-    MPI_Barrier(MPI_COMM_WORLD);
-    endtime=MPI_Wtime();
-
-
-    /* foffset*=nprocs; if we want total megabytes written */
-    if(rank==0){
-      puts("*");
-      unlink(filename); 
-      puts("======================================================");
-      printf("Raw MPI-IO Total Duration %lf seconds, iterations=%u %lf Megabytes written per processor Nprocs= %u \n",
-	     (endtime-starttime),i,((double)foffset)/(1024.0*1024.0),nprocs);
-      printf("Raw MPI-IO Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
-	     (double)(nprocs*localnp*sizeof(double))*((double)nsteps)*6.0/((endtime-starttime)*1024.0*1024.0),
-	     (double)(localnp*sizeof(double))*((double)nsteps)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
-      puts("======================================================");
-    }
-
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-
-    /* OK, now we do this using POSIX IO */
-    sprintf(newfilename,"testio%u.%u.dat",rank,nprocs);
-    unlink(newfilename);
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-    fd = fopen(newfilename,"w");
-    /* start the timer */
-    starttime=endtime=MPI_Wtime();
-    for(i=0;i<nsteps;i++){
-      for(j=0;j<6;j++){
-	/* touch data */
-	for(n=0;n<localnp;n++)
-	  (data[j])[n]=(double)rank;
-	fwrite(data[j],sizeof(double),localnp,fd);
-      }
-      curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
-      MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-    }
-    fclose(fd);
-    MPI_Barrier(MPI_COMM_WORLD);
-    endtime=MPI_Wtime();
-    if(rank==0) puts("*");
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-    unlink(newfilename); 
-    MPI_Barrier(MPI_COMM_WORLD);
-    if(rank==0){
-      puts("======================================================");
-      printf("Raw 1-file-per-proc Total Duration %lf seconds, iterations=%u %lf Megabytes written Nprocs= %u \n",
-	     (endtime-starttime),nsteps,((double)foffset)/(1024.0*1024.0),nprocs);
-      printf("Raw 1-file-per-proc Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
-	     (double)(nprocs*localnp*sizeof(double))*((double)nsteps)*6.0/((endtime-starttime)*1024.0*1024.0),
-	     (double)(localnp*sizeof(double))*((double)nsteps)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
-      puts("======================================================");
-    }
-
-#ifndef DISABLE_H5PART
-
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-    /* OK, now we do this using H5Part */
-
-    sprintf(filename,"%s.%u.h5.dat",FILENAME,nprocs); 
-    if(rank==0) unlink(filename);
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-
-    f = H5OpenFile(filename,H5_O_WRONLY,MPI_COMM_WORLD);
-    MPI_Barrier(MPI_COMM_WORLD); /* to prevent unlink from interfering with file open */
-    /* start the timer */
-    starttime=endtime=MPI_Wtime();
-    H5PartSetNumParticles(f,localnp);
-    for(i=0;i<nsteps;i++){
-      for(j=0;j<6;j++){
-	/* touch data */
-	for(n=0;n<localnp;n++)
-	  (data[j])[n]=(double)rank;
-      }
-      H5SetStep(f,i);
-      H5PartWriteDataFloat64(f,"x",x);
-      H5PartWriteDataFloat64(f,"y",y);
-      H5PartWriteDataFloat64(f,"z",z);
-      H5PartWriteDataFloat64(f,"px",px);
-      H5PartWriteDataFloat64(f,"py",py);
-      H5PartWriteDataFloat64(f,"pz",pz);
-
-      curtime=MPI_Wtime(); /* ensure no race condition by broadcasting time */
-      MPI_Bcast(&curtime,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
-    }
-    H5CloseFile(f);
-    MPI_Barrier(MPI_COMM_WORLD);
-    endtime=MPI_Wtime();
-    if(rank==0){
-      puts("*");
-      unlink(filename); 
-      puts("======================================================");
-      printf("H5Part Total Duration %lf seconds, iterations=%u %lf Megabytes written Nprocs= %u \n",
-	     (endtime-starttime),nsteps,((double)foffset)/(1024.0*1024.0),nprocs);
-      printf("H5Part Effective Data Rate = %lf Megabytes/sec global and %lf Megabytes/sec per task Nprocs= %u \n",
-	     (double)(nprocs*localnp*sizeof(double))*((double)nsteps)*6.0/((endtime-starttime)*1024.0*1024.0),
-	     (double)(localnp*sizeof(double))*((double)nsteps)*6.0/((endtime-starttime)*1024.0*1024.0),nprocs);
-      puts("======================================================");
-    }
-    MPI_Barrier(MPI_COMM_WORLD);
-#endif
-  } /* trials */
-
-  MPI_Finalize();
-
-  return 0;
-}
diff --git a/examples/H5Part/H5test.c b/examples/H5Part/H5test.c
deleted file mode 100644
index e3de2e4..0000000
--- a/examples/H5Part/H5test.c
+++ /dev/null
@@ -1,263 +0,0 @@
-
-#include <stdlib.h>
-#include "H5hut.h"
-
-/*
-  A simple regression test that shows how you use this API
-  to write and read multi-timestep files of particle data.
-*/
-
-
-#ifdef PARALLEL_IO
-
-int main(int argc,char *argv[]){
-        int sz=5;
-        double *x,*y,*z;
-        h5_int64_t *id;
-        h5_file_t file;
-        int i,t,nt,nds;
-
-        int nprocs = 1;
-        int myproc = 0;
-        MPI_Comm comm = MPI_COMM_WORLD;
-
-#if PARALLEL_IO
-        MPI_Init (&argc,&argv);
-        MPI_Comm_size (comm,&nprocs);
-        MPI_Comm_rank (comm,&myproc);
-#endif
-        x = (double*)malloc (sz*nprocs*sizeof(double));
-        y = (double*)malloc (sz*nprocs*sizeof(double));
-        z = (double*)malloc (sz*nprocs*sizeof(double));
-        id=(h5_int64_t*)malloc (sz*nprocs*sizeof(h5_int64_t));
-
-        file=H5OpenFile ("parttest.h5",H5_O_WRONLY,comm);
-        if(!file) {
-                perror("File open failed:  exiting!");
-                exit(0);
-        }
-
-        H5PartWriteFileAttribString (
-                file,
-                "File Description",
-                "This file is created by H5PartTest.cc. "
-                "Simple H5Part file for testing purpose...");
-        char* FileAttrib = "Created by H5PartTest.cc";
-        H5PartWriteFileAttrib (
-                file,
-                "Origin",
-                H5T_NATIVE_CHAR,
-                FileAttrib ,strlen(FileAttrib));
-
-        for(t=0;t<5;t++){
-                MPI_Barrier(comm);
-                for(i=0;i<sz;i++) {
-                        x[i]=(double)(i+t)+10.0*(double)myproc;
-                        y[i]=0.1 + (double)(i+t);
-                        z[i]=0.2 + (double)(i+t*10);
-                        id[i]=i+sz*myproc;
-                }
-                printf("Proc[%u] Writing timestep %u \n",myproc,t);
-                H5SetStep(file,t); /* must set the current timestep in file */
-                H5PartSetNumParticles(file,sz); /* then set number of particles to store */
-                /* now write different tuples of data into this timestep of the file */
-                H5PartWriteDataFloat64(file,"x",x); 
-                H5PartWriteDataFloat64(file,"y",y);
-                H5PartWriteDataFloat64(file,"z",z);
-                
-                H5PartWriteDataFloat64(file,"px",x); 
-                H5PartWriteDataFloat64(file,"py",y);
-                H5PartWriteDataFloat64(file,"pz",z);
-
-                H5PartWriteDataInt64(file,"id",id);
-                H5PartWriteStepAttribString (
-                        file,
-                        "Step Description",
-                        "STEP STEP STEP");
-                char* StepAttrib = "STEP";
-                H5PartWriteStepAttrib (
-                        file,
-                        "Step",
-                        H5T_NATIVE_CHAR,
-                        StepAttrib ,strlen(StepAttrib));
-
-        }
-
-        printf("Done writing p[%u]\n",myproc);
-        H5CloseFile(file);
-        fprintf(stderr,"Closed files p[%u]\n",myproc);
-        MPI_Barrier(comm);
-  
-        fprintf(stderr,
-                "p[%u:%u] : OK, close file and reopen for reading idStart %u  idEnd %u \n",
-                myproc,nprocs,idStart,idEnd);
-
-        unsigned int idStart = 0+sz*myproc;
-        unsigned int idEnd = (sz-1)+sz*myproc;
-
-
-        file=H5OpenFile("parttest.h5",H5_O_RDONLY,comm);
-        H5SetStep(file,0);
-        // unsigned int np = 0;
-        unsigned int np = (int)H5PartGetNumParticles(file);
-        nt=H5GetNumSteps(file); /* get number of steps in file */
-        nds=H5PartGetNumDatasets(file); /* get number of datasets in timestep 0 */
-
-        MPI_Barrier(comm);
-  
-        H5PartSetView(file,idStart,idEnd);
-
-        np = (int)H5PartGetNumParticles(file);
-        printf("After SetView(%d,%d): steps= %u  datasets= %u particles= %u\n",
-               (int)idStart,(int)idEnd,
-               nt,nds,np);
-        
-        free(x); 
-        free(y);
-        free(z);
-        free(id);
-
-        H5CloseFile(file);
-        MPI_Barrier(comm);
-        fprintf(stderr,"proc[%u]:  done\n",myproc);
-        return MPI_Finalize();
-}
-
-#else
- 
-  /*+++++++++++++ Reopen File for Reading +++H5PartSetStep(h5partFile,0)++++++++*/
-  file=H5PartOpenFile("parttest.h5",H5_O_RDONLY);
-
-  
-  /********************************************/
-  H5PartSetStep(file,0);
-  nt=H5GetNumSteps(file); /* get number of steps in file */
-  nds=H5PartGetNumDatasets(file); /* get number of datasets in timestep 0 */
-  np=H5PartGetNumParticles(file);
-  
-
-  fprintf(stdout,"OK, close file and reopen for reading\n");
-  fprintf(stdout,"steps= %u\tdatasets=%u\tparticles= %u\n",
-	 nt,nds,np);
-  
-  // clear the particles
-  for(i=0;i<np;i++){
-    x[i]=y[i]=z[i]=0.0;
-    id[i]=0;
-  }
-
-  H5PartReadDataFloat64(file,"x",x);
-  H5PartReadDataFloat64(file,"y",y);
-  H5PartReadDataFloat64(file,"z",z);
-  H5PartReadDataInt64(file,"id",id);
-
-  for(i=0;i<np;i++){
-    fprintf(stdout,
-	    "\tp[%3u] x=%lf y=%lf z=%lf id=%lld\n",
-	    i,x[i],y[i],z[i],(long long)(id[i]));
-  }
-  /************************  std::cout << "nParticles: " << nParticles << std::endl;
-********************/
-  printf("Set to last step and reload data\n");
-  H5PartSetStep(file,nt-1);
-  H5PartReadDataFloat64(file,"x",x);
-  H5PartReadDataFloat64(file,"y",y);
-  H5PartReadDataFloat64(file,"z",z);
-  H5PartReadDataInt64(file,"id",id);
-  for(i=0;i<np;i++){
-    fprintf(stdout,"\tp[%3u] x=%lf y=%lf z=%lf id=%lld\n",
-	    i,x[i],y[i],z[i],(long long) (id[i]));
-  }
-  
-  /********************************************/
-  idEnd=np;
-  printf("Old View is %d:%d\n",(int)idStart,(int)idEnd);
-  H5PartSetView(file,idStart,idEnd>>1);
-  printf("Set new view = %d:%d\n",(int)idStart,(int)(idEnd>>1));
-  H5PartGetView(file,&idStart,&idEnd);
-  np=H5PartGetNumParticles(file);
-  printf("steps= %u  datasets= %u particles= %d with view %d:%d\n",
-	 nt,nds,(int)np,(int)idStart,(int)idEnd);
-  H5PartSetStep(file,nt-1); // set to last step
-  printf("Setting to last step = %u\n",nt-1);
-  for(i=0;i<10;i++){ x[i]=y[i]=z[i]=0.0; id[i]=0; } /* clear the arrays */
-  H5PartReadDataFloat64(file,"x",x);
-  H5PartReadDataFloat64(file,"y",y);
-  H5PartReadDataFloat64(file,"z",z);
-  H5PartReadDataInt64(file,"id",id);
-
-  for(i=0;i<np;i++){
-    fprintf(stdout,
-	    "\tp[%3u] x=%lf y=%lf z=%lf id=%lld\n",
-	    i,x[i],y[i],z[i],(long long)id[i]);
-  }
-
-  /********************************************/
-  printf("Now set the view to the latter half of the data in step #%u\n",nt-1);
-  H5PartResetView(file);
-  H5PartGetView(file,&idStart,&idEnd);
-  printf("Reset view = %d:%d\nSetting to %u:%u\n",
-	 (int)idStart,(int)idEnd,
-	 (int)idEnd>>1,(int)idEnd);
-  H5PartSetView(file,(idEnd>>1),idEnd);
-  np=H5PartGetNumParticles(file);
-  printf("Now particles in selection are %d\n",np);
-  printf("doubleCheck=%lld\n", (long long)H5PartGetView(file,0,0));
-  
-  for(i=0;i<10;i++){ x[i]=y[i]=z[i]=0.0; id[i]=0; } /* clear the arrays */
- 
-  H5PartReadDataFloat64(file,"x",x);
-  H5PartReadDataFloat64(file,"y",y);
-  H5PartReadDataFloat64(file,"z",z);
-  H5PartReadDataInt64(file,"id",id);
-  for(i=0;i<np;i++){
-    fprintf(stdout,
-	    "\tp[%3u] x=%lf y=%lf z=%lf id=%lld\n",
-	    i,x[i],y[i],z[i],(long long)id[i]);
-  }
-
-  // read dataset names
-  h5_int64_t status = H5_SUCCESS;
-  const h5_int64_t lenName = 64;
-  char datasetName[lenName];
-  h5_int64_t datasetType;
-  h5_int64_t datasetNElems;
-
-  H5PartSetStep(file,0);
-  for (h5_int64_t i=0; i < nds; i++) {
-    status = H5PartGetDatasetInfo(file, i, datasetName, lenName,
-				  &datasetType, &datasetNElems);
-
-    if (status != H5_SUCCESS) {
-      perror("Could not retrieve dataset names!");
-    }
-    else {
-      printf("datasetName: %s, type: %lld, nElements: %lld   ", 
-	     datasetName, datasetType, datasetNElems);
-      if (datasetType ==  H5_INT64_T) {
-	printf("H5PPART_INT64 \n");
-      }
-      else {
-	printf("H5PPART_FLOAT64 \n");
-      }
-    }
-  }
-
-  if(x) 
-    free(x); 
-  if(y) 
-    free(y);
-  if(z) 
-    free(z);
-  if(id) 
-    free(id);
-
-  H5PartCloseFile(file);
-
-  fprintf(stderr,"done\n");
-}
-
-#endif
-
-
-
diff --git a/examples/H5Part/H5testF.f90 b/examples/H5Part/H5testF.f90
deleted file mode 100644
index 3ba08e4..0000000
--- a/examples/H5Part/H5testF.f90
+++ /dev/null
@@ -1,924 +0,0 @@
-PROGRAM H5PartExampleF
-   IMPLICIT   NONE
-#ifdef PARALLEL_IO
-   INCLUDE 'mpif.h'
-#endif
-   INCLUDE 'H5hutF.h'
-
-   INTEGER*8 :: h5_err = 0
-   INTEGER :: mpi_err
-   INTEGER :: num_args = 0
-   INTEGER :: i
-   CHARACTER(LEN=32) :: arg_str
-   INTEGER :: opt_read = 0
-   INTEGER :: opt_write = 0
-   INTEGER*8 :: opt_debug = 2
-   INTEGER :: comm = 0
-   INTEGER :: myproc = 0
-   INTEGER :: num_procs = 0
-   CHARACTER(LEN=128) :: fname = "testfile.h5"
-   INTEGER*8 :: file
-
-   REAL*8    :: r8_attrib (2)
-   REAL*4    :: r4_attrib (2)
-   INTEGER*8 :: i8_attrib (2)
-   INTEGER*4 :: i4_attrib (2)
-   CHARACTER (LEN=H5_MAX_NAME_LEN) :: string_attrib
-   LOGICAL skip
-
-   DATA r8_attrib / 42.0, 43.0 /
-   DATA r4_attrib / 42.0, 43.0 /
-   DATA i8_attrib / 42, 43 /
-   DATA i4_attrib / 42, 43 /
-   DATA string_attrib / "The answer is 42." /
-
-   num_args = IARGC ()
-   IF (num_args == 0) THEN
-      PRINT *, "Usage: H5PartExampleF -w | -r [-g]"      
-      CALL EXIT (1)
-   END IF
-   DO i = 1, num_args
-      PRINT *, "Loop..."
-      if (skip .EQV. .TRUE.) THEN
-         CYCLE
-      END IF
-      CALL GETARG (i, arg_str)
-      IF (arg_str == "-r") THEN
-         opt_read = 1
-      ELSE IF (arg_str == "-w") THEN
-         opt_write = 1
-      ELSE IF (arg_str == "-d") THEN
-         CALL GETARG (i+1, arg_str)
-         arg_str = TRIM (arg_str)
-         READ (arg_str, *, ERR=999) opt_debug
-         skip = .TRUE.
-         CYCLE
-999      PRINT *, "Debug level must be an integer: ", arg_str
-         CALL EXIT (1)
-      ELSE
-         PRINT *, "Illegal option ", arg_str, "\n"
-         PRINT *, "Usage: H5BlockExampleF -w | -r [-d LEVEL]"
-         CALL EXIT (1)
-      END IF
-   END DO
-
-#ifdef PARALLEL_IO
-   CALL MPI_Init (mpi_err)
-   comm = MPI_COMM_WORLD
-   CALL MPI_Comm_rank (comm, myproc, mpi_err)
-   CALL MPI_Comm_size (comm, num_procs, mpi_err)
-
-   h5_err = h5_set_verbosity_level (opt_debug)
-
-   IF (opt_write == 1) THEN
-     PRINT "('[proc ', I3, ']: Open file for writing ...')", myproc
-     file = h5_openw_par (fname, comm)
-     IF (file == 0) THEN
-        PRINT "('[proc ', I3, ']: Error opening file ...')", myproc
-        h5_err = -2
-        GOTO 911
-     ENDIF
-     h5_err = write_file (file, myproc, num_procs)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Failed to write file ', A, '!')", myproc, fname
-        GOTO 911
-     END IF
-  ELSE IF (opt_read == 1) THEN
-     PRINT "('[proc ', I3, ']: Open file for reading ...')", myproc
-     file = h5_openr_par (fname, comm)
-     IF (file == 0) THEN
-        PRINT "('[proc ', I3, ']: Error opening file ...')", myproc
-        GOTO 911
-     ENDIF
-     h5_err = read_file (file, myproc, num_procs)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Failed to read file ', A, '!')", myproc, fname
-        GOTO 911
-     END IF
-  ENDIF
-
-  PRINT "('[proc ', I3, ']: Done.')", myproc
-911 CONTINUE  
-
-  h5_err = h5_close (file)
-  IF (h5_err < 0) THEN
-     PRINT "('[proc ', I3, ']: Error closing file ...')", myproc
-  ENDIF
-  PRINT "('[proc ', I3, ']: File closed.')", myproc
-
-  h5_err = h5_finalize ()
-  IF (h5_err < 0) THEN
-     PRINT "('[proc ', I3, ']: Error closing H5hut library ...')", myproc
-  ENDIF
-  PRINT "('[proc ', I3, ']: H5hut finalized.')", myproc
-
-  CALL MPI_Finalize (mpi_error)
-  PRINT "('[proc ', I3, ']: MPI finalized.')", myproc
-
-#else
-  !! SERIAL CODE
-   h5_err = h5_set_verbosity_level (opt_debug)
-  IF (opt_write == 1) THEN
-     PRINT "('[proc ', I3, ']: Open file for writing ...')", myproc
-     file = h5_openw (fname)
-     IF (file == 0) THEN
-        PRINT "('[proc ', I3, ']: Error opening file ...')", myproc
-        h5_err = -2
-        GOTO 911
-     ENDIF
-     h5_err = write_file (file, myproc, num_procs)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Failed to write file ', A, '!')", myproc, fname
-        GOTO 911
-     END IF
-
-  ELSE IF (opt_read == 1) THEN
-     file = h5_openr (fname)
-     IF (file == 0) THEN
-        PRINT "('[proc ', I3, ']: Error opening file ...')", myproc
-        GOTO 911
-     ENDIF
-     h5_err = read_file (file, myproc, num_procs)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Failed to read file ', A, '!')", myproc, fname
-        GOTO 911
-     END IF
-  ENDIF
-  PRINT "('[proc ', I3, ']: Done.')", myproc
-911 CONTINUE
-
-  h5_err = h5_close (file)
-  IF (h5_err < 0) THEN
-     PRINT "('[proc ', I3, ']: Error closing file ...')", myproc
-  ENDIF
-
-#endif
-
-  CALL EXIT (h5_err)
-
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!!
-   CONTAINS
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION write_file (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) ::        file
-     INTEGER, INTENT(IN) ::          myproc
-     INTEGER, INTENT(IN) ::          num_procs
-
-     INCLUDE 'H5hutF.h'
-
-     INTEGER   :: mpi_err
-     INTEGER*8 :: step = 1
-     INTEGER*8 :: num_steps = 2
-
-     h5_err = write_file_attribs (file, myproc)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing file attributes ...')", myproc
-        write_file = h5_err
-        RETURN
-     ENDIF
-
-     DO step = 1, num_steps
-        h5_err = h5_setstep (file, step)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Error setting step to ', I2)", myproc, step
-           write_file = h5_err
-           RETURN
-        ENDIF
-
-        h5_err = write_step_attribs (file, myproc, num_procs)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Error writing step attributes ...')", myproc
-           write_file = h5_err
-           RETURN
-        ENDIF
-
-        h5_err = write_data (file, myproc, num_procs)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Error writing step data ...')", myproc
-           write_file = h5_err
-           RETURN
-        ENDIF
-
-     ENDDO
-
-     write_file = 0
-   END FUNCTION write_file
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION write_file_attribs (file, myproc)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-
-     INCLUDE 'H5hutF.h'
-
-     PRINT "('[proc ', I3, ']: Writing file attribute of type H5_FLOAT64_T ...')", myproc
-     h5_err = h5_writefileattrib_r8 (file, "r8_attrib", r8_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing file attribute of type H5_FLOAT64_T')", myproc
-        write_file_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing file attribute of type H5_FLOAT32_T ...')", myproc
-     h5_err = h5_writefileattrib_r4 (file, "r4_attrib", r4_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing file attribute of type H5_FLOAT32_T')", myproc
-        write_file_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing file attribute of type H5_INT64_T ...')", myproc
-     h5_err = h5_writefileattrib_i8 (file, "i8_attrib", i8_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing file attribute of type H5_INT64_T')", myproc
-        write_file_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing file attribute of type H5_INT32_T ...')", myproc
-     h5_err = h5_writefileattrib_i4 (file, "i4_attrib", i4_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing file attribute of type H5_INT32_T')", myproc
-        write_file_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing file attribute of type H5_STRING_T ...')", myproc
-     h5_err = h5_writefileattrib_string (file, "string_attrib", string_attrib)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing file attribute of type H5_STRING_T')", myproc
-        write_file_attribs = h5_err
-        RETURN
-     ENDIF
-
-     write_file_attribs = 0
-   END FUNCTION write_file_attribs
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION write_step_attribs (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-     INTEGER, INTENT(IN) :: num_procs
-
-     INCLUDE 'H5hutF.h'
-
-     PRINT "('[proc ', I3, ']: Writing step attribute of type H5_FLOAT64_T ...')", myproc
-     h5_err = h5_writestepattrib_r8 (file, "r8_attrib", r8_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing step attribute of type H5_FLOAT64_T')", myproc
-        write_step_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing step attribute of type H5_FLOAT32_T ...')", myproc
-     h5_err = h5_writestepattrib_r4 (file, "r4_attrib", r4_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing step attribute of type H5_FLOAT32_T')", myproc
-        write_step_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing step attribute of type H5_INT64_T ...')", myproc
-     h5_err = h5_writestepattrib_i8 (file, "i8_attrib", i8_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing step attribute of type H5_INT64_T')", myproc
-        write_step_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing step attribute of type H5_INT32_T ...')", myproc
-     h5_err = h5_writestepattrib_i4 (file, "i4_attrib", i4_attrib, 2_8)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing step attribute of type H5_INT32_T')", myproc
-        write_step_attribs = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing step attribute of type H5_STRING_T ...')", myproc
-     h5_err = h5_writestepattrib_string (file, "string_attrib", string_attrib)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing step attribute of type H5_STRING_T')", myproc
-        write_step_attribs = h5_err
-        RETURN
-     ENDIF
-
-     write_step_attribs = 0
-   END FUNCTION write_step_attribs
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION write_data (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-     INTEGER, INTENT(IN) :: num_procs
-
-     INCLUDE 'H5hutF.h'
-
-     INTEGER*8 :: h5_err = 0
-
-     REAL*8,    ALLOCATABLE :: r8_data (:)
-     REAL*4,    ALLOCATABLE :: r4_data (:)
-     INTEGER*8, ALLOCATABLE :: i8_data (:)
-     INTEGER*4, ALLOCATABLE :: i4_data (:)
-
-     INTEGER*8 :: num_points = 1024
-     INTEGER*8 :: offset = 0
-     INTEGER*8 :: step = 1
-
-     ALLOCATE (r8_data (num_points))
-     ALLOCATE (r4_data (num_points))
-     ALLOCATE (i8_data (num_points))
-     ALLOCATE (i4_data (num_points))
-
-     step = h5_getstep (file)
-
-     offset = num_points * myproc
-     DO i = 1, num_points
-        r8_data (i) = REAL (i + offset + step)
-        r4_data (i) = REAL (i + offset + step)
-        i8_data (i) = i + offset + step
-        i4_data (i) = i + offset + step
-     ENDDO
-
-     h5_err = h5pt_setnpoints (file, num_points)
-
-     PRINT "('[proc ', I3, ']: Writing file data of type H5_FLOAT64_T ...')", myproc
-     h5_err = h5pt_writedata_r8 (file,"r8_data", r8_data)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing r8 data')", myproc
-        write_data = h5_err
-        RETURN
-     ENDIF
-
-     PRINT "('[proc ', I3, ']: Writing file data of type H5_FLOAT32_T ...')", myproc
-     h5_err = h5pt_writedata_r4 (file,"r4_data", r4_data)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing r4 data')", myproc
-        write_data = h5_err
-        RETURN
-     ENDIF
-     
-     PRINT "('[proc ', I3, ']: Writing file data of type H5_INT64_T ...')", myproc
-     h5_err = h5pt_writedata_i8 (file,"i8_data", i8_data)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing i8 data')", myproc
-        write_data = h5_err
-        RETURN
-     ENDIF
-     
-     PRINT "('[proc ', I3, ']: Writing file data of type H5_INT32_T ...')", myproc
-     h5_err = h5pt_writedata_i4 (file,"i4_data", i4_data)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error writing i4 data')", myproc
-        write_data = h5_err
-        RETURN
-     ENDIF
-
-     write_data = 0
-   END FUNCTION write_data
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION read_file (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-     INTEGER, INTENT(IN) :: num_procs
-
-     INCLUDE 'H5hutF.h'
-
-     INTEGER   :: mpi_err
-     INTEGER*8 :: step = 1
-     INTEGER*8 :: num_steps = 0
-
-     h5_err = read_file_attribs (file, myproc, num_procs)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Error reading file attributes ...')", myproc
-        read_file = h5_err
-        RETURN
-     ENDIF
-
-     num_steps = h5_getnsteps (file)
-     PRINT "('[proc ', I3, ']: Number of steps = ', I3)", myproc, num_steps
-     DO step = 1, num_steps
-        PRINT "('[proc ', I3, ']: Set step = ', I3)", myproc, step
-        h5_err = h5_setstep (file, step)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Error setting step to ', I2)", myproc, step
-           read_file = h5_err
-           RETURN
-        ENDIF
-
-        h5_err = read_step_attribs (file, myproc, num_procs)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Error reading step data ...')", myproc
-           read_file = h5_err
-           RETURN
-        ENDIF
-
-        h5_err = read_data (file, myproc, num_procs)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Error reading step attributes ...')", myproc
-           read_file = h5_err
-           RETURN
-        ENDIF
-
-     ENDDO
-
-     read_file = 0
-911  CONTINUE
-   END FUNCTION read_file
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION cmp_i4 (x, y, size)
-     IMPLICIT   NONE
-     INTEGER*4, DIMENSION (1:size), INTENT (IN) :: x, y
-     INTEGER*8, INTENT (IN) :: size
-     INTEGER*8  i
-     DO i = 1, size
-        IF (x (i) /= y (i)) THEN
-           PRINT "('[proc ', I3, ']: Value error ', I6, ' /= ', I6)", &
-                myproc, x (i), y (i)
-           cmp_i4 = -2
-           RETURN
-        END IF
-     END DO
-     cmp_i4 = 0
-   END FUNCTION cmp_i4
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION cmp_i8 (x, y, size)
-     IMPLICIT   NONE
-     INTEGER*8, DIMENSION (1:size), INTENT (IN) :: x, y
-     INTEGER*8, INTENT (IN) :: size
-     INTEGER*8  i
-     DO i = 1, size
-        IF (x (i) /= y (i)) THEN
-           PRINT "('[proc ', I3, ']: Value error ', I6, ' /= ', I6)", &
-                myproc, x (i), y (i)
-           cmp_i8 = -2
-           RETURN
-        END IF
-     END DO
-     cmp_i8 = 0
-   END FUNCTION cmp_i8
-
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION cmp_r4 (x, y, size)
-     IMPLICIT   NONE
-     REAL*4, DIMENSION (1:size), INTENT (IN) :: x, y
-     INTEGER*8, INTENT (IN) :: size
-     INTEGER*8  i
-     DO i = 1, size
-        IF (x (i) /= y (i)) THEN
-           PRINT "('[proc ', I3, ']: Value error ', F6.2, ' /= ', F6.2)", &
-                myproc, x (i), y (i)
-           cmp_r4 = -2
-           RETURN
-        END IF
-     END DO
-     cmp_r4 = 0
-   END FUNCTION cmp_r4
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION cmp_r8 (x, y, size)
-     IMPLICIT   NONE
-     REAL*8, DIMENSION (1:size), INTENT (IN) :: x, y
-     INTEGER*8, INTENT (IN) :: size
-     INTEGER*8  i
-     DO i = 1, size
-        IF (x (i) /= y (i)) THEN
-           PRINT "('[proc ', I3, ']: Value error ', F6.2, ' /= ', F6.2)", &
-                myproc, x (i), y (i)
-           cmp_r8 = -2
-           RETURN
-        END IF
-     END DO
-     cmp_r8 = 0
-   END FUNCTION cmp_r8
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION cmp_string (x, y, size)
-     IMPLICIT   NONE
-     CHARACTER (LEN=size), INTENT (IN) :: x, y
-     INTEGER*8, INTENT (IN) :: size
-     INTEGER*8  i
-     DO i = 1, size
-        IF (x /= y) THEN
-           PRINT "('[proc ', I3, ']: Value error ', A, ' /= ', A)", &
-                myproc, x, y
-           cmp_string = -2
-           RETURN
-        END IF
-     END DO
-     cmp_string = 0
-   END FUNCTION cmp_string
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION read_file_attribs (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-     INTEGER, INTENT(IN) :: num_procs
-
-     INCLUDE 'H5hutF.h'
-
-     INTEGER*8 :: num_file_attribs
-     INTEGER*8 :: idx
-     CHARACTER(LEN=128) :: name
-     INTEGER*8 :: type
-     INTEGER*8 :: num_elems
-     INTEGER*8 :: h5_err
-
-     INTEGER*4, ALLOCATABLE :: i4 (:)
-     INTEGER*8, ALLOCATABLE :: i8 (:)
-     REAL*4,    ALLOCATABLE :: r4 (:)
-     REAL*8,    ALLOCATABLE :: r8 (:)
-     CHARACTER (LEN=H5_MAX_NAME_LEN) :: str
-
-     num_file_attribs = h5_getnfileattribs (file)
-     IF (num_file_attribs == H5_FAILURE) THEN
-        PRINT "('[proc ', I3, ']: Cannot read number of file attributes ...')", myproc
-        read_file_attribs = h5_err
-        RETURN
-     END IF
-     PRINT "('[proc ', I3, ']: Number of file attributes: ', I3)", myproc, num_file_attribs
-
-     DO idx = 1, num_file_attribs
-        h5_err = h5_getfileattribinfo (file, idx, name, type, num_elems)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Cannot get info about step attribute ', I2)", myproc, idx
-           RETURN
-        END IF
-        SELECT CASE (type)
-        CASE (H5_INT32_T)
-           PRINT "('[proc ', I3, ']: Reading file attribute of type H5_INT32_T ...')", myproc
-           ALLOCATE (i4 (num_elems))
-           h5_err = h5_readfileattrib_i4 (file, name, i4)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-           h5_err = cmp_i4 (i4, i4_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-           
-        CASE (H5_INT64_T)
-           PRINT "('[proc ', I3, ']: Reading file attribute of type H5_INT64_T ...')", myproc
-           ALLOCATE (i8 (num_elems))
-           h5_err = h5_readfileattrib_i8 (file, name, i8)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-           h5_err = cmp_i8 (i8, i8_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-
-        CASE (H5_FLOAT32_T)
-           PRINT "('[proc ', I3, ']: Reading file attribute of type H5_FLOAT32_T ...')", myproc
-           ALLOCATE (r4 (num_elems))
-           h5_err = h5_readfileattrib_r4 (file, name, r4)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-           h5_err = cmp_r4 (r4, r4_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-
-        CASE (H5_FLOAT64_T)
-           PRINT "('[proc ', I3, ']: Reading file attribute of type H5_FLOAT64_T ...')", myproc
-           ALLOCATE (r8 (num_elems))
-           h5_err = h5_readfileattrib_r8 (file, name, r8)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_file_attribs = h5_err
-              RETURN
-            END IF
-           h5_err = cmp_r8 (r8, r8_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_file_attribs = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_STRING_T)
-           PRINT "('[proc ', I3, ']: Reading file attribute of type H5_STRING_T ...')", myproc
-           h5_err = h5_readfileattrib_string (file, name, str)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_file_attribs = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_string (str, string_attrib, H5_MAX_NAME_LEN)
-           IF (h5_err < 0) THEN
-              read_file_attribs = h5_err
-              RETURN
-           END IF
-        END SELECT
-     END DO
-     read_file_attribs = 0
-   END FUNCTION read_file_attribs
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION read_step_attribs (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-     INTEGER, INTENT(IN) :: num_procs
-
-     INCLUDE 'H5hutF.h'
-
-
-     INTEGER*8 :: num_step_attribs
-     INTEGER*8 :: idx
-     CHARACTER(LEN=128) :: name
-     INTEGER*8 :: type
-     INTEGER*8 :: num_elems
-     INTEGER*8 :: h5_err
-
-     INTEGER*4, ALLOCATABLE :: i4 (:)
-     INTEGER*8, ALLOCATABLE :: i8 (:)
-     REAL*4,    ALLOCATABLE :: r4 (:)
-     REAL*8,    ALLOCATABLE :: r8 (:)
-     CHARACTER (LEN=H5_MAX_NAME_LEN) :: str
-
-     num_step_attribs = h5_getnstepattribs (file)
-     IF (num_step_attribs == H5_FAILURE) THEN
-        PRINT "('[proc ', I3, ']: Cannot read number of step attributes ...')", myproc
-        read_step_attribs = h5_err
-     END IF
-     PRINT "('[proc ', I3, ']: Number of step attributes: ', I3)", myproc, num_step_attribs
-
-     DO idx = 1, num_step_attribs
-        h5_err = h5_getstepattribinfo (file, idx, name, type, num_elems)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Cannot get info about step attribute ', I2)", myproc, idx
-           RETURN
-        END IF
-        SELECT CASE (type)
-        CASE (H5_INT32_T)
-           PRINT "('[proc ', I3, ']: Reading step attribute of type H5_INT32_T ...')", myproc
-           ALLOCATE (i4 (num_elems))
-           h5_err = h5_readstepattrib_i4 (file, name, i4)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_i4 (i4, i4_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-           
-        CASE (H5_INT64_T)
-           PRINT "('[proc ', I3, ']: Reading step attribute of type H5_INT64_T ...')", myproc
-           ALLOCATE (i8 (num_elems))
-           h5_err = h5_readstepattrib_i8 (file, name, i8)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_i8 (i8, i8_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_FLOAT32_T)
-           PRINT "('[proc ', I3, ']: Reading step attribute of type H5_FLOAT32_T ...')", myproc
-           ALLOCATE (r4 (num_elems))
-           h5_err = h5_readstepattrib_r4 (file, name, r4)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_r4 (r4, r4_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_FLOAT64_T)
-           PRINT "('[proc ', I3, ']: Reading step attribute of type H5_FLOAT64_T ...')", myproc
-           ALLOCATE (r8 (num_elems))
-           h5_err = h5_readstepattrib_r8 (file, name, r8)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_r8 (r8, r8_attrib, num_elems)
-           IF (h5_err < 0) THEN
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_STRING_T)
-           PRINT "('[proc ', I3, ']: Reading step attribute of type H5_STRING_T ...')", myproc
-           h5_err = h5_readstepattrib_string (file, name, str)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading attribute')", myproc
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_string (str, string_attrib, H5_MAX_NAME_LEN)
-           IF (h5_err < 0) THEN
-              read_step_attribs = h5_err
-              RETURN
-           END IF
-
-        END SELECT
-     END DO
-
-     read_step_attribs = 0
-   END FUNCTION read_step_attribs
-
-   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-   !!
-   INTEGER*8 FUNCTION read_data (file, myproc, num_procs)
-     IMPLICIT   NONE
-
-     INTEGER*8, INTENT(IN) :: file
-     INTEGER, INTENT(IN) :: myproc
-     INTEGER, INTENT(IN) :: num_procs
-
-     INCLUDE 'H5hutF.h'
-
-     INTEGER*8 :: h5_err = 0
-     INTEGER*8 :: num_points = 0
-     INTEGER*8 :: num = 0
-     INTEGER*8 :: num_points_total = 0
-     INTEGER*8 :: num_datasets = 0
-     INTEGER*8 :: offset = 0
-     INTEGER*8 :: idx = 0
-     INTEGER*8 :: step = 0
-
-     CHARACTER(LEN=128) :: name
-     INTEGER*8 :: type
-
-     REAL*8,    ALLOCATABLE :: r8_data (:)
-     REAL*4,    ALLOCATABLE :: r4_data (:)
-     INTEGER*8, ALLOCATABLE :: i8_data (:)
-     INTEGER*4, ALLOCATABLE :: i4_data (:)
-     REAL*8,    ALLOCATABLE :: r8 (:)
-     REAL*4,    ALLOCATABLE :: r4 (:)
-     INTEGER*8, ALLOCATABLE :: i8 (:)
-     INTEGER*4, ALLOCATABLE :: i4 (:)
-
-     num_datasets = h5pt_getndatasets (file)
-     IF (num_datasets < 0) THEN
-        PRINT "('[proc ', I3, ']: Cannot read number of datasets ...')", myproc
-        read_data = h5_err
-        RETURN
-     END IF
-     PRINT "('[proc ', I3, ']: Number of datasets: ', I3)", myproc, num_datasets
-
-     h5_err = h5pt_resetview (file)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Cannot reset view ...')", myproc
-        read_data = h5_err
-        RETURN
-     END IF
-
-     num_points_total = h5pt_getnpoints (file)
-     IF (num_points_total < 0) THEN
-        PRINT "('[proc ', I3, ']: Cannot read total number of points ...')", myproc
-        read_data = h5_err
-        RETURN
-     END IF
-     PRINT "('[proc ', I3, ']: Total number of points: ', I6)", myproc, num_points_total
-
-     ! select subset
-     num_points = num_points_total / num_procs
-     offset = num_points * myproc + 1;
-     h5_err = h5pt_setview (file, offset, offset + num_points - 1)
-     IF (h5_err < 0) THEN
-        PRINT "('[proc ', I3, ']: Cannot set view ...')", myproc
-        read_data = h5_err
-     END IF
-     PRINT "('[proc ', I3, ']: View: ', I4, ':', I4)", myproc, offset, offset+num_points-1
-
-     ALLOCATE (r8_data (num_points))
-     ALLOCATE (r4_data (num_points))
-     ALLOCATE (i8_data (num_points))
-     ALLOCATE (i4_data (num_points))
-     ALLOCATE (r8 (num_points))
-     ALLOCATE (r4 (num_points))
-     ALLOCATE (i8 (num_points))
-     ALLOCATE (i4 (num_points))
-
-     step = h5_getstep (file)
-
-     offset = num_points * myproc
-     DO i = 1, num_points
-        r8_data (i) = REAL (i + offset + step)
-        r4_data (i) = REAL (i + offset + step)
-        i8_data (i) = i + offset + step
-        i4_data (i) = i + offset + step
-     ENDDO
-
-     DO idx = 1, num_datasets
-        h5_err = h5pt_getdatasetinfo (file, idx, name, type, num)
-        IF (h5_err < 0) THEN
-           PRINT "('[proc ', I3, ']: Cannot get info about datset ', I2)", myproc, idx
-        END IF
-        SELECT CASE (type)
-        CASE (H5_INT32_T)
-           PRINT "('[proc ', I3, ']: Reading dataset of type H5_INT32_T ...')", myproc
-           h5_err = h5pt_readdata_i4 (file, name, i4)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading datset')", myproc
-              read_data = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_i4 (i4, i4_data, num_points)
-           IF (h5_err < 0) THEN
-              read_data = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_INT64_T)
-           PRINT "('[proc ', I3, ']: Reading dataset of type H5_INT64_T ...')", myproc
-           h5_err = h5pt_readdata_i8 (file, name, i8)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading datset')", myproc
-              read_data = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_i8 (i8, i8_data, num_points)
-           IF (h5_err < 0) THEN
-              read_data = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_FLOAT32_T)
-           PRINT "('[proc ', I3, ']: Reading dataset of type H5_FLOAT32_T ...')", myproc
-           h5_err = h5pt_readdata_r4 (file, name, r4)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading datset')", myproc
-              read_data = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_r4 (r4, r4_data, num_points)
-           IF (h5_err < 0) THEN
-              read_data = h5_err
-              RETURN
-           END IF
-
-        CASE (H5_FLOAT64_T)
-           PRINT "('[proc ', I3, ']: Reading dataset of type H5_FLOAT64_T ...')", myproc
-           h5_err = h5pt_readdata_r8 (file, name, r8)
-           IF (h5_err < 0) THEN
-              PRINT "('[proc ', I3, ']: Error reading datset')", myproc
-              read_data = h5_err
-              RETURN
-           END IF
-           h5_err = cmp_r8 (r8, r8_data, num_points)
-           IF (h5_err < 0) THEN
-              read_data = h5_err
-              RETURN
-           END IF
-           
-        END SELECT
-     END DO
-
-     read_data = 0
-   END FUNCTION read_data
-END PROGRAM