add ierr to MPI_Finalize

test/H5testFpar.f90

@@ -12,9 +12,7 @@
       INTEGER*8 nstep
       INTEGER*8 ndata
       INTEGER*8 npoints,pointoffset
-      INTEGER*8 step
-      INTEGER*8 I
-      INTEGER*8 J
+      INTEGER*8 I,J,K,step
       REAL*8,ALLOCATABLE:: X(:),Y(:),Z(:),PX(:),PY(:),PZ(:)
       INTEGER*8,ALLOCATABLE:: ID(:)
       REAL*8 REALTIME(1)
@@ -23,6 +21,11 @@
       INTEGER*8 rc
       INTEGER:: comm
       
+      INTEGER*8:: offset
+      INTEGER*8, allocatable, dimension(:,:) :: Localnum1
+      INTEGER*8, allocatable, dimension(:,:,:) :: Localnum
+      INTEGER nprocrow, nproccol
+
       call MPI_Init(ierr)
       comm = MPI_COMM_WORLD
       call MPI_Comm_rank(comm,myproc,ierr)
@@ -51,11 +54,12 @@
 !     write an annotation to the file
       rc=h5pt_writefileattrib_string(file,"Annotation","Testing 1 2 3")
       print *,"and now write the steps"
-      do I=1,nstep
+ 
+      do step=1,nstep
 !        Set the step number
-         print *,"Write Step ",I
-         rc = h5pt_setstep(file,I)
-!        Now start writing the data arrays for this step
+         print *,"Write Step ",step
+         rc = h5pt_setstep(file,step)
+         !        Now start writing the data arrays for this step
          rc = h5pt_writedata_r8(file,"x",X)
          rc = h5pt_writedata_r8(file,"y",Y)
          rc = h5pt_writedata_r8(file,"z",Z)
@@ -66,14 +70,45 @@
          do J=1,npoints
             ID(J)=ID(J)+10
          enddo
-!        And write a simple floatingpoint attribute associated with this timestep
+         !        And write a simple floatingpoint attribute associated with this timestep
          REALTIME(1) = REAL(I)*0.1
          rc = h5pt_writestepattrib_r8(file,"RealTime",REALTIME,1_8)
+         print*,'Write REALTIME ',REALTIME
+         !
+         ! Test file arrtibutes
+         !
+         nprocrow=2
+         nproccol=2
+         print*,'Write localnum ',Localnum
+         allocate(Localnum1(1,nprocrow*nproccol))
+         allocate(Localnum(2,0:nprocrow-1,0:nproccol-1))
+         !
+         ! init testdata
+         !
+         do k = 0, 1
+            do i = 0,(nprocrow-1)
+               do j = 0,(nproccol-1)
+                  Localnum(k,i,j) = k+i+j
+               enddo
+            enddo
+         enddo
+         !
+         offset=0
+         do i = 0, (nprocrow-1)
+            do j = 0,(nproccol-1)
+               Localnum1(offset) = Localnum(1,i,j)
+               offset = offset + 1
+            enddo
+         enddo
+         rc = h5pt_writestepattrib_i8(file,"Localnum1",Localnum1(:),offset)   
+         print*,'rc (Localnum1)= ',rc
+         
       enddo
+      
       print *,"Done writing, now close the file"
       rc = h5pt_close(file)
-
-! **************** Clean out some variables ***************
+      
+      ! **************** Clean out some variables ***************
       nstep=0
 !      npoints=0
       do I=1,npoints
@@ -82,16 +117,25 @@
          Z(I)=-1.0
          ID(I)=0
       enddo
+      
+      do k = 0, 1
+         do i = 0,(nprocrow-1)
+            do j = 0,(nproccol-1)
+               Localnum(k,i,j) = 0
+            enddo
+         enddo
+      enddo
+      
+      REALTIME(1) = 0.0
 
-! *****************Now Reopen for Reading ******************
-
-
+      ! *****************Now Reopen for Reading ******************
+      
       print *,"Open file for reading"
       file = h5pt_openr_par("testfilef.h5",comm)
       print *,"  Opened testfilef.h5"
       nstep = h5pt_getnsteps(file)
       print *,"  Nsteps = ",nstep
-      rc = h5pt_setstep(file,1_8)
+       rc = h5pt_setstep(file,1_8)
       print *,"now get the number of datasets"
       ndata = h5pt_getndatasets(file)
       print *,"  Ndata=",ndata
@@ -109,15 +153,19 @@
 !        set the current step
          rc = h5pt_setstep(file,step)
          rc = h5pt_readdata_i8(file,"id",ID)
+         rc = h5pt_readstepattrib_r8(file,"RealTime",REALTIME)
+         print*,'Read REALTIME ',REALTIME
+         rc = h5pt_readstepattrib_i8(file,"Localnum1",Localnum1)   
+         print*,'Read Localnum ',Localnum1
 !        read the Z data from the current step
 !        rc = h5prt_readdata(file,step,X,Y,Z,PX,PY,PZ,ID)
-         do J=1,10
-            print *,"    ID(",J,")==",ID(J)
-         enddo
+!         do J=1,10
+!            print *,"    ID(",J,")==",ID(J)
+!         enddo
       enddo
 
       rc = h5pt_close(file)
       print *,"done"      
-      call MPI_Finalize()
+      call MPI_Finalize(ierr)
       
       end