From 4f35ba1356c2b9f0a0421f24e5dbbdf803b8aff9 Mon Sep 17 00:00:00 2001 From: Achim Gsell Date: Tue, 12 Sep 2006 22:18:58 +0000 Subject: [PATCH] test/H5testFpar.f - renamed to "H5testFpar.f90" --- .gitattributes | 1 - test/H5testFpar.f | 123 ---------------------------------------------- 2 files changed, 124 deletions(-) delete mode 100644 test/H5testFpar.f diff --git a/.gitattributes b/.gitattributes index 7660043..f4d6a35 100644 --- a/.gitattributes +++ b/.gitattributes @@ -106,7 +106,6 @@ test/H5PartTestParallel.cc -text test/H5ecloudExample.cc -text test/H5test.cc -text test/H5testF.f -text -test/H5testFpar.f -text test/H5testFpar.f90 -text test/Makefile.am -text tools/Makefile.am -text diff --git a/test/H5testFpar.f b/test/H5testFpar.f deleted file mode 100644 index 700a2e6..0000000 --- a/test/H5testFpar.f +++ /dev/null @@ -1,123 +0,0 @@ -! ============== -! -! Sample Fortran program that uses HDF5 bindings -! -! ============== - - program H5testFpar - implicit none - include 'mpif.h' - include 'H5PartF90.inc' - INTEGER*8 file - INTEGER*8 nstep - INTEGER*8 ndata - INTEGER*8 npoints,pointoffset - INTEGER*8 step - INTEGER*8 I - INTEGER*8 J - REAL*8,ALLOCATABLE:: X(:),Y(:),Z(:),PX(:),PY(:),PZ(:) - INTEGER*8,ALLOCATABLE:: ID(:) - REAL*8 REALTIME(1) - INTEGER myproc,nprocs - INTEGER ierr - INTEGER*8 rc - INTEGER:: comm - - call MPI_Init(ierr) - comm = MPI_COMM_WORLD - call MPI_Comm_rank(comm,myproc,ierr) - call MPI_Comm_size(comm,nprocs,ierr) - - file = h5pt_openw_par("testfilef.h5",comm) - print *,"Opened file testfilef.h5 for writing" - npoints = 1024 - nstep = 10 - ALLOCATE(X(npoints),Y(npoints),Z(npoints)) - ALLOCATE(PX(npoints),PY(npoints),PZ(npoints)) - ALLOCATE(ID(npoints)) - print *," Npoints=",npoints," nsteps=",nstep - print *," writing X,Y,Z,PX,PY,PZ,ID" - print *," ... initialize the data arrays" - do I=1,npoints - X(INT(I))=0.0 - Y(I)=1.0+REAL(I) - Z(I)=100.0+REAL(I)*2.0 - ID(I)=I - enddo - print *,"Tell h5pt how many particles are stored in the file" -! set the number of points - rc = h5pt_setnpoints(file,npoints) - print *,"write an attribute string" -! 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 -! Set the step number - print *,"Write Step ",I - rc = h5pt_setstep(file,I) -! 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) - rc = h5pt_writedata_r8(file,"px",PX) - rc = h5pt_writedata_r8(file,"py",PY) - rc = h5pt_writedata_r8(file,"pz",PZ) - rc = h5pt_writedata_i8(file,"id",ID) - do J=1,npoints - ID(J)=ID(J)+10 - enddo -! 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) - enddo - print *,"Done writing, now close the file" - rc = h5pt_close(file) - -! **************** Clean out some variables *************** - nstep=0 -! npoints=0 - do I=1,npoints - X(I)=-1.0 - Y(I)=-1.0 - Z(I)=-1.0 - ID(I)=0 - enddo - -! *****************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) - print *,"now get the number of datasets" - ndata = h5pt_getndatasets(file) - print *," Ndata=",ndata - npoints = h5pt_getnpoints(file) - print *," Total NP=",npoints - npoints = npoints/nprocs - pointoffset = npoints*myproc; - rc = h5pt_setview(file,pointoffset,pointoffset+npoints) - print *," Pointoffset=",pointoffset," Local NP=",npoints -! Now we need to set our view to read correct points - - - do step=1,nstep - print *,"Read step ",step -! set the current step - rc = h5pt_setstep(file,step) - rc = h5pt_readdata_i8(file,"id",ID) -! 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 - enddo - - rc = h5pt_close(file) - print *,"done" - call MPI_Finalize() - - end