(network Common Data Form) is a set of software libraries and machine-independent data formats that support the creation, access, and sharing of array-oriented scientific data (Fortran library).
| f90tst_nc4_par.f90 (netcdf-fortran-4.4.4) | : | f90tst_nc4_par.f90 (netcdf-fortran-4.4.5) |
|---|
| ! This parallel test was contributed by Jim Edwards at UCAR. Thanks Jim! | | ! This parallel test was contributed by Jim Edwards at UCAR. Thanks Jim! |
| program f90tst_nc4_par | | program f90tst_nc4_par |
| use netcdf | | use netcdf |
| use mpi | | use mpi |
| implicit none | | implicit none |
| | |
|
| character (len = *), parameter :: FILE_NAME = "f90tst_nc4_par.nc" | | integer :: ierr, my_task, nprocs, nmode |
| integer :: nmode, ierr, fh, my_task, nprocs, i, varid | | |
| integer :: dimid(3), start(3), count(3) | | |
| real :: f(3) | | |
| | |
| call MPI_INIT(ierr) | | call MPI_INIT(ierr) |
| call MPI_COMM_RANK(MPI_COMM_WORLD, my_task, ierr) | | call MPI_COMM_RANK(MPI_COMM_WORLD, my_task, ierr) |
| call MPI_COMM_SIZE(MPI_COMM_WORLD, nprocs, ierr) | | call MPI_COMM_SIZE(MPI_COMM_WORLD, nprocs, ierr) |
| | |
| if(nprocs/=8)then | | if(nprocs/=8)then |
| stop 'requires 8 tasks' | | stop 'requires 8 tasks' |
| end if | | end if |
| | |
|
| nmode = ior(NF90_CLOBBER,NF90_NETCDF4) | | #ifdef NF_HAS_PNETCDF |
| | nmode = IOR(nf90_clobber, nf90_mpiio) |
| | call parallel_io(nmode) |
| | #endif |
| | |
| | #ifdef NF_HAS_PARALLEL4 |
| | nmode = ior(NF90_CLOBBER, NF90_NETCDF4) |
| nmode = IOR(nmode, nf90_mpiio) | | nmode = IOR(nmode, nf90_mpiio) |
|
| | call parallel_io(nmode) |
| | #endif |
| | |
| | call MPI_Finalize(ierr) |
| | |
| | contains |
| | ! This subroutine handles errors by printing an error message and |
| | ! exiting with a non-zero status. |
| | subroutine handle_err(errcode) |
| | use netcdf |
| | implicit none |
| | integer, intent(in) :: errcode |
| | |
| | if(errcode /= nf90_noerr) then |
| | print *, 'Error: ', trim(nf90_strerror(errcode)) |
| | stop 2 |
| | endif |
| | end subroutine handle_err |
| | |
| | subroutine parallel_io(nmode) |
| | use netcdf |
| | use mpi |
| | implicit none |
| | |
| | character (len = *), parameter :: FILE_NAME = "f90tst_nc4_par.nc" |
| | integer :: nmode, fill_mode, ierr, fh, my_task, nprocs, i, varid |
| | integer :: dimid(3), start(3), count(3) |
| | real :: f(3) |
| | |
| | call MPI_COMM_RANK(MPI_COMM_WORLD, my_task, ierr) |
| | |
| call handle_err(nf90_create(FILE_NAME, nmode, fh, & | | call handle_err(nf90_create(FILE_NAME, nmode, fh, & |
| comm = MPI_COMM_WORLD, info = MPI_INFO_NULL)) | | comm = MPI_COMM_WORLD, info = MPI_INFO_NULL)) |
| | |
|
| call handle_err(nf90_set_fill(fh, NF90_NOFILL, nmode)) | | call handle_err(nf90_set_fill(fh, NF90_NOFILL, fill_mode)) |
| | |
| call handle_err(nf90_def_dim(fh, 'dim1', 6, dimid(1))) | | call handle_err(nf90_def_dim(fh, 'dim1', 6, dimid(1))) |
| call handle_err(nf90_def_dim(fh, 'dim2', 4, dimid(2))) | | call handle_err(nf90_def_dim(fh, 'dim2', 4, dimid(2))) |
| call handle_err(nf90_def_dim(fh, 'dim3', 1, dimid(3))) | | call handle_err(nf90_def_dim(fh, 'dim3', 1, dimid(3))) |
| | |
| call handle_err(nf90_def_var(fh, 'var1', NF90_DOUBLE, dimid, varid)) | | call handle_err(nf90_def_var(fh, 'var1', NF90_DOUBLE, dimid, varid)) |
| call handle_err(nf90_enddef(fh)) | | call handle_err(nf90_enddef(fh)) |
| | |
| do i=1,3 | | do i=1,3 |
| f(i) = my_task*3+i | | f(i) = my_task*3+i |
| | |
| skipping to change at line 60 | | skipping to change at line 93 |
|---|
| ! Reopen the file and check it. | | ! Reopen the file and check it. |
| call handle_err(nf90_open(FILE_NAME, NF90_MPIIO, fh, & | | call handle_err(nf90_open(FILE_NAME, NF90_MPIIO, fh, & |
| comm = MPI_COMM_WORLD, info = MPI_INFO_NULL)) | | comm = MPI_COMM_WORLD, info = MPI_INFO_NULL)) |
| | |
| call handle_err(nf90_get_var(fh, varid, f, start=start, count=count)) | | call handle_err(nf90_get_var(fh, varid, f, start=start, count=count)) |
| do i=1,3 | | do i=1,3 |
| if (f(i) .ne. my_task*3+i) stop 3 | | if (f(i) .ne. my_task*3+i) stop 3 |
| end do | | end do |
| | |
| call handle_err(nf90_close(fh)) | | call handle_err(nf90_close(fh)) |
|
| call MPI_Finalize(ierr) | | end subroutine parallel_io |
| | |
| contains | | |
| ! This subroutine handles errors by printing an error message and | | |
| ! exiting with a non-zero status. | | |
| subroutine handle_err(errcode) | | |
| use netcdf | | |
| implicit none | | |
| integer, intent(in) :: errcode | | |
| | |
| if(errcode /= nf90_noerr) then | | |
| print *, 'Error: ', trim(nf90_strerror(errcode)) | | |
| stop 2 | | |
| endif | | |
| end subroutine handle_err | | |
| | | |
| end program f90tst_nc4_par | | end program f90tst_nc4_par |
| | | |
| End of changes. 5 change blocks. |
|---|
| 22 lines changed or deleted | | 40 lines changed or added |
|---|
"Fossies" - the Fresh Open Source Software Archive 