#include "ESMFPIO.h" !=================================================== ! DO NOT EDIT THIS FILE, it was generated using genf90.pl ! Any changes you make to this file may be lost !=================================================== #define __PIO_FILE__ 'piodarray' module piodarray use pio_types, only : file_desc_t, io_desc_t, var_desc_t, pio_noerr, iosystem_desc_t, & pio_iotype_pbinary, pio_iotype_binary, pio_iotype_direct_pbinary, & pio_iotype_netcdf, pio_iotype_pnetcdf, pio_iotype_netcdf4p, pio_iotype_netcdf4c, & PIO_MAX_VAR_DIMS use pio_kinds use pio_support use pionfwrite_mod, only : write_nf use pionfread_mod, only : read_nf use iompi_mod use alloc_mod use rearrange #ifdef TIMING use perf_mod, only : t_startf, t_stopf !_EXTERNAL #endif implicit none include 'mpif.h' !_EXTERNAL private public :: pio_read_darray, pio_write_darray !> !! @defgroup PIO_write_darray PIO_write_darray !! @brief The overloaded PIO_write_darray writes a distributed array to disk. !< # 27 "piodarray.F90.in" interface PIO_write_darray ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_1d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_2d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_3d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_4d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_5d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_1d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_2d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_3d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_4d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_5d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_1d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_2d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_3d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_4d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure write_darray_5d_double end interface !> !! @defgroup PIO_read_darray PIO_read_darray # 36 "piodarray.F90.in" !! @brief The overloaded PIO_read_darray function reads a distributed array from disk. !< # 38 "piodarray.F90.in" interface PIO_read_darray ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_1d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_2d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_3d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_4d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_5d_real ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_1d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_2d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_3d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_4d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_5d_int ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_1d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_2d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_3d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_4d_double ! TYPE real,int,double ! DIMS 1,2,3,4,5 module procedure read_darray_5d_double end interface !> !! @private !< character(len=*), parameter, private :: modName='piodarray' # 49 "piodarray.F90.in" contains ! TYPE real,int,double !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 1D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 64 "piodarray.F90.in" subroutine write_darray_1d_real (File,varDesc,ioDesc, array, iostat, fillval) use pio_msg_mod, only : pio_msg_writedarray ! !DESCRIPTION: ! Writes a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r4), dimension(:), target, intent(in) :: & array ! array to be written real(r4), optional, intent(in) :: fillval ! rearrange receiver fill value type(iosystem_desc_t), pointer :: ios integer(i4), intent(out) :: iostat integer :: msg, ierr character(len=*), parameter :: subName=modName//'::write_darray_real' ios => file%iosystem if(ios%async_interface .and. .not. ios%ioproc) then msg = PIO_MSG_WRITEDARRAY if(debugasync) print *,__PIO_FILE__,__LINE__ if(ios%comp_rank==0) call mpi_send(msg, 1, mpi_integer, ios%ioroot, 1, ios%union_comm, ierr) call mpi_bcast(file%fh, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%varid, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%rec, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%ndims, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(iodesc%async_id, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(MPI_REAL4, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) if(debugasync) print *,__PIO_FILE__,__LINE__, MPI_REAL4 if(present(fillval)) then call mpi_bcast(1, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(fillval, 1, MPI_REAL4, ios%compmaster, ios%intercomm, ierr) else call mpi_bcast(0, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) end if if(debugasync) print *,__PIO_FILE__,__LINE__ endif select case(File%iotype) case(pio_iotype_pbinary, pio_iotype_direct_pbinary) if (present(fillval)) then call write_darray_bin_real(File,varDesc,iodesc, array, iostat, fillval) else call write_darray_bin_real(File,varDesc,iodesc, array, iostat) endif case(pio_iotype_pnetcdf, pio_iotype_netcdf, pio_iotype_netcdf4c, pio_iotype_netcdf4p) if (present(fillval)) then call write_darray_nf_real(File,varDesc,iodesc, array, iostat, fillval) else call write_darray_nf_real(File,varDesc,iodesc, array, iostat) endif case(pio_iotype_binary) print *, subName,': IO type not supported' end select end subroutine write_darray_1d_real ! TYPE real,int,double !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 1D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 64 "piodarray.F90.in" subroutine write_darray_1d_int (File,varDesc,ioDesc, array, iostat, fillval) use pio_msg_mod, only : pio_msg_writedarray ! !DESCRIPTION: ! Writes a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor integer(i4), dimension(:), target, intent(in) :: & array ! array to be written integer(i4), optional, intent(in) :: fillval ! rearrange receiver fill value type(iosystem_desc_t), pointer :: ios integer(i4), intent(out) :: iostat integer :: msg, ierr character(len=*), parameter :: subName=modName//'::write_darray_int' ios => file%iosystem if(ios%async_interface .and. .not. ios%ioproc) then msg = PIO_MSG_WRITEDARRAY if(debugasync) print *,__PIO_FILE__,__LINE__ if(ios%comp_rank==0) call mpi_send(msg, 1, mpi_integer, ios%ioroot, 1, ios%union_comm, ierr) call mpi_bcast(file%fh, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%varid, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%rec, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%ndims, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(iodesc%async_id, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(MPI_INTEGER, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) if(debugasync) print *,__PIO_FILE__,__LINE__, MPI_INTEGER if(present(fillval)) then call mpi_bcast(1, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(fillval, 1, MPI_INTEGER, ios%compmaster, ios%intercomm, ierr) else call mpi_bcast(0, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) end if if(debugasync) print *,__PIO_FILE__,__LINE__ endif select case(File%iotype) case(pio_iotype_pbinary, pio_iotype_direct_pbinary) if (present(fillval)) then call write_darray_bin_int(File,varDesc,iodesc, array, iostat, fillval) else call write_darray_bin_int(File,varDesc,iodesc, array, iostat) endif case(pio_iotype_pnetcdf, pio_iotype_netcdf, pio_iotype_netcdf4c, pio_iotype_netcdf4p) if (present(fillval)) then call write_darray_nf_int(File,varDesc,iodesc, array, iostat, fillval) else call write_darray_nf_int(File,varDesc,iodesc, array, iostat) endif case(pio_iotype_binary) print *, subName,': IO type not supported' end select end subroutine write_darray_1d_int ! TYPE real,int,double !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 1D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 64 "piodarray.F90.in" subroutine write_darray_1d_double (File,varDesc,ioDesc, array, iostat, fillval) use pio_msg_mod, only : pio_msg_writedarray ! !DESCRIPTION: ! Writes a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r8), dimension(:), target, intent(in) :: & array ! array to be written real(r8), optional, intent(in) :: fillval ! rearrange receiver fill value type(iosystem_desc_t), pointer :: ios integer(i4), intent(out) :: iostat integer :: msg, ierr character(len=*), parameter :: subName=modName//'::write_darray_double' ios => file%iosystem if(ios%async_interface .and. .not. ios%ioproc) then msg = PIO_MSG_WRITEDARRAY if(debugasync) print *,__PIO_FILE__,__LINE__ if(ios%comp_rank==0) call mpi_send(msg, 1, mpi_integer, ios%ioroot, 1, ios%union_comm, ierr) call mpi_bcast(file%fh, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%varid, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%rec, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%ndims, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(iodesc%async_id, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(MPI_REAL8, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) if(debugasync) print *,__PIO_FILE__,__LINE__, MPI_REAL8 if(present(fillval)) then call mpi_bcast(1, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(fillval, 1, MPI_REAL8, ios%compmaster, ios%intercomm, ierr) else call mpi_bcast(0, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) end if if(debugasync) print *,__PIO_FILE__,__LINE__ endif select case(File%iotype) case(pio_iotype_pbinary, pio_iotype_direct_pbinary) if (present(fillval)) then call write_darray_bin_double(File,varDesc,iodesc, array, iostat, fillval) else call write_darray_bin_double(File,varDesc,iodesc, array, iostat) endif case(pio_iotype_pnetcdf, pio_iotype_netcdf, pio_iotype_netcdf4c, pio_iotype_netcdf4p) if (present(fillval)) then call write_darray_nf_double(File,varDesc,iodesc, array, iostat, fillval) else call write_darray_nf_double(File,varDesc,iodesc, array, iostat) endif case(pio_iotype_binary) print *, subName,': IO type not supported' end select end subroutine write_darray_1d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 2D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_2d_real (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r4), intent(in) :: & array(:,:) ! array to be written real(r4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_real (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_2d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 3D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_3d_real (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r4), intent(in) :: & array(:,:,:) ! array to be written real(r4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_real (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_3d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 4D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_4d_real (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r4), intent(in) :: & array(:,:,:,:) ! array to be written real(r4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_real (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_4d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 5D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_5d_real (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r4), intent(in) :: & array(:,:,:,:,:) ! array to be written real(r4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_real (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_real (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_5d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 2D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_2d_int (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor integer(i4), intent(in) :: & array(:,:) ! array to be written integer(i4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat integer(i4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_int (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_2d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 3D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_3d_int (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor integer(i4), intent(in) :: & array(:,:,:) ! array to be written integer(i4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat integer(i4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_int (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_3d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 4D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_4d_int (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor integer(i4), intent(in) :: & array(:,:,:,:) ! array to be written integer(i4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat integer(i4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_int (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_4d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 5D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_5d_int (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor integer(i4), intent(in) :: & array(:,:,:,:,:) ! array to be written integer(i4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat integer(i4) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_int (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_int (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_5d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 2D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_2d_double (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r8), intent(in) :: & array(:,:) ! array to be written real(r8), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r8) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_double (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_2d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 3D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_3d_double (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r8), intent(in) :: & array(:,:,:) ! array to be written real(r8), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r8) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_double (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_3d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 4D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_4d_double (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r8), intent(in) :: & array(:,:,:,:) ! array to be written real(r8), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r8) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_double (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_4d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_write_darray !! @brief Writes a 5D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 150 "piodarray.F90.in" subroutine write_darray_5d_double (File,varDesc,ioDesc, array, iostat, fillval) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! variable descriptor real(r8), intent(in) :: & array(:,:,:,:,:) ! array to be written real(r8), optional, intent(in) :: fillval ! rearrange receiver fill value integer(i4), intent(out) :: iostat real(r8) :: transvar(1), dumbvar(0) # 170 "piodarray.F90.in" ! cannot call transfer function with a 0 sized array if(size(array)==0) then call write_darray_1d_double (File, varDesc, iodesc, dumbvar, iostat) else if(present(fillval)) then call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat, fillval) else call write_darray_1d_double (File, varDesc, iodesc, transfer(array,transvar), iostat) end if end subroutine write_darray_5d_double ! TYPE real,int,double !> !! @public !! @ingroup PIO_read_darray !! @brief !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 193 "piodarray.F90.in" subroutine read_darray_1d_real (File,varDesc, ioDesc, array, iostat) use pio_msg_mod, only : pio_msg_readdarray ! !DESCRIPTION: ! Reads a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r4), dimension(:), intent(out) :: & array ! array to be read integer(i4), intent(out) :: iostat character(len=*), parameter :: subName=modName//'::read_darray_real' type(iosystem_desc_t), pointer :: ios integer :: ierr, msg array = 0 ios => File%iosystem if(ios%async_interface .and. .not. ios%ioproc) then msg = PIO_MSG_READDARRAY if(DebugAsync) print *,__PIO_FILE__,__LINE__ if(ios%comp_rank==0) call mpi_send(msg, 1, mpi_integer, ios%ioroot, 1, ios%union_comm, ierr) call mpi_bcast(file%fh, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%varid, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%rec, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(iodesc%async_id, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(MPI_REAL4, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) if(DebugAsync) print *,__PIO_FILE__,__LINE__, MPI_REAL4 endif select case(File%iotype) case(pio_iotype_pbinary, pio_iotype_direct_pbinary) call read_darray_bin_real (File,varDesc,iodesc,array, iostat) case(pio_iotype_pnetcdf, pio_iotype_netcdf, pio_iotype_netcdf4c, pio_iotype_netcdf4p) call read_darray_nf_real (File,varDesc,iodesc,array, iostat) case(pio_iotype_binary) print *, subName,': IO type not supported' end select end subroutine read_darray_1d_real ! TYPE real,int,double !> !! @public !! @ingroup PIO_read_darray !! @brief !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 193 "piodarray.F90.in" subroutine read_darray_1d_int (File,varDesc, ioDesc, array, iostat) use pio_msg_mod, only : pio_msg_readdarray ! !DESCRIPTION: ! Reads a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor integer(i4), dimension(:), intent(out) :: & array ! array to be read integer(i4), intent(out) :: iostat character(len=*), parameter :: subName=modName//'::read_darray_int' type(iosystem_desc_t), pointer :: ios integer :: ierr, msg array = 0 ios => File%iosystem if(ios%async_interface .and. .not. ios%ioproc) then msg = PIO_MSG_READDARRAY if(DebugAsync) print *,__PIO_FILE__,__LINE__ if(ios%comp_rank==0) call mpi_send(msg, 1, mpi_integer, ios%ioroot, 1, ios%union_comm, ierr) call mpi_bcast(file%fh, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%varid, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%rec, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(iodesc%async_id, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(MPI_INTEGER, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) if(DebugAsync) print *,__PIO_FILE__,__LINE__, MPI_INTEGER endif select case(File%iotype) case(pio_iotype_pbinary, pio_iotype_direct_pbinary) call read_darray_bin_int (File,varDesc,iodesc,array, iostat) case(pio_iotype_pnetcdf, pio_iotype_netcdf, pio_iotype_netcdf4c, pio_iotype_netcdf4p) call read_darray_nf_int (File,varDesc,iodesc,array, iostat) case(pio_iotype_binary) print *, subName,': IO type not supported' end select end subroutine read_darray_1d_int ! TYPE real,int,double !> !! @public !! @ingroup PIO_read_darray !! @brief !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 193 "piodarray.F90.in" subroutine read_darray_1d_double (File,varDesc, ioDesc, array, iostat) use pio_msg_mod, only : pio_msg_readdarray ! !DESCRIPTION: ! Reads a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r8), dimension(:), intent(out) :: & array ! array to be read integer(i4), intent(out) :: iostat character(len=*), parameter :: subName=modName//'::read_darray_double' type(iosystem_desc_t), pointer :: ios integer :: ierr, msg array = 0 ios => File%iosystem if(ios%async_interface .and. .not. ios%ioproc) then msg = PIO_MSG_READDARRAY if(DebugAsync) print *,__PIO_FILE__,__LINE__ if(ios%comp_rank==0) call mpi_send(msg, 1, mpi_integer, ios%ioroot, 1, ios%union_comm, ierr) call mpi_bcast(file%fh, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%varid, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(vardesc%rec, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(iodesc%async_id, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) call mpi_bcast(MPI_REAL8, 1, mpi_integer, ios%compmaster, ios%intercomm, ierr) if(DebugAsync) print *,__PIO_FILE__,__LINE__, MPI_REAL8 endif select case(File%iotype) case(pio_iotype_pbinary, pio_iotype_direct_pbinary) call read_darray_bin_double (File,varDesc,iodesc,array, iostat) case(pio_iotype_pnetcdf, pio_iotype_netcdf, pio_iotype_netcdf4c, pio_iotype_netcdf4p) call read_darray_nf_double (File,varDesc,iodesc,array, iostat) case(pio_iotype_binary) print *, subName,': IO type not supported' end select end subroutine read_darray_1d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 2D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_2d_real (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r4), intent(out) :: array(:,:) ! array to be read integer(i4), intent(out) :: iostat real(r4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_real (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_2d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 3D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_3d_real (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r4), intent(out) :: array(:,:,:) ! array to be read integer(i4), intent(out) :: iostat real(r4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_real (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_3d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 4D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_4d_real (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r4), intent(out) :: array(:,:,:,:) ! array to be read integer(i4), intent(out) :: iostat real(r4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_real (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_4d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 5D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_5d_real (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r4), intent(out) :: array(:,:,:,:,:) ! array to be read integer(i4), intent(out) :: iostat real(r4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_real (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_5d_real ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 2D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_2d_int (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor integer(i4), intent(out) :: array(:,:) ! array to be read integer(i4), intent(out) :: iostat integer(i4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_int (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_2d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 3D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_3d_int (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor integer(i4), intent(out) :: array(:,:,:) ! array to be read integer(i4), intent(out) :: iostat integer(i4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_int (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_3d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 4D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_4d_int (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor integer(i4), intent(out) :: array(:,:,:,:) ! array to be read integer(i4), intent(out) :: iostat integer(i4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_int (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_4d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 5D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_5d_int (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor integer(i4), intent(out) :: array(:,:,:,:,:) ! array to be read integer(i4), intent(out) :: iostat integer(i4), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_int (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_5d_int ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 2D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_2d_double (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r8), intent(out) :: array(:,:) ! array to be read integer(i4), intent(out) :: iostat real(r8), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_double (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_2d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 3D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_3d_double (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r8), intent(out) :: array(:,:,:) ! array to be read integer(i4), intent(out) :: iostat real(r8), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_double (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_3d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 4D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_4d_double (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r8), intent(out) :: array(:,:,:,:) ! array to be read integer(i4), intent(out) :: iostat real(r8), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_double (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_4d_double ! TYPE real,int,double ! DIMS 2,3,4,5 !> !! @public !! @ingroup PIO_read_darray !! @brief Reads a 5D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 263 "piodarray.F90.in" subroutine read_darray_5d_double (File,varDesc,ioDesc, array, iostat) ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! iodecomp descriptor real(r8), intent(out) :: array(:,:,:,:,:) ! array to be read integer(i4), intent(out) :: iostat real(r8), pointer :: tmpvar(:) call alloc_check(tmpvar,size(array)) call read_darray_1d_double (File, varDesc, iodesc, tmpvar, iostat) array = reshape(tmpvar,shape(array)) call dealloc_check(tmpvar) end subroutine read_darray_5d_double ! TYPE real,int,double !> !! @private !! @brief Write a 1D array of type real defined by varDesc using the decomposition described in iodesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 300 "piodarray.F90.in" subroutine write_darray_nf_real (File,varDesc,ioDesc,array, iostat, fillval) ! !DESCRIPTION: ! Writes a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor real(r4), target, intent(in) :: & array(:) ! array to be written real(r4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(pio_offset), pointer :: start(:), count(:) integer(i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::write_darray_nf_real' real(r4), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if IO processor integer (i4) :: len, &! length of IO decomp segmap glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform ndims ! number of variable dimensions integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind) :: UseRearranger #if DEBUG_REARR real(r4), dimension(:), pointer :: array2 integer i #endif real(r4) :: rsum integer(i4) :: ierr #ifdef TIMING call t_startf("pio_write_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! --------------------------------------------------------- ! pull information out of the decomposition data structure ! --------------------------------------------------------- glen = iodesc%glen offset = iodesc%IOmap%start len = iodesc%IOmap%length if(Debug) & print *,__PIO_FILE__,__LINE__,' NAME : IAM: ',File%iosystem%comp_rank, & ' UseRearranger: ',UseRearranger, glen, offset, len #ifdef TIMING call t_startf("pio_rearrange_write") #endif if(UseRearranger) then if (IOproc) then if(Debug) print *, subName,': IAM: ',File%iosystem%comp_rank,'Before call to allocate(IOBUF): ',len call alloc_check(IOBUF,len,' TYPE :IOBUF') if (present(fillval)) then IOBUF=fillval else IOBUF= -1.0_r8 endif !------------------------------------------------ ! set the IO buffer to a particular test pattern !------------------------------------------------ !JMD IOBUF(:) = real(File%iosystem%io_rank,kind=r8) if(Debug) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len !,' IOBUF: ',IOBUF else call alloc_check(IOBUF,0) IOBUF= -1.0_r8 endif !------------------------------------ ! Rearrange data from comp->io decomp !------------------------------------ ! "array" is comp data call rearrange_comp2io(File%iosystem,iodesc, array, iobuf) #if DEBUG_REARR call alloc_check(array2,size(array),'array2') call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array2) do i=1,size(array) if (array(i) /= array2(i)) then print *, subName,': error: write ping-pong test failed on index',i end if end do print *, subName,': passed write ping-pong test' call dealloc_check(array2) !!!!!!! end debug #endif !-------------------------------------------- ! End data rearrange !-------------------------------------------- else iobuf=>array endif ! if(UseRearranger) #ifdef TIMING call t_stopf("pio_rearrange_write") #endif if (IOproc) then !---------------------------------------------- ! write the global 2-d slice from IO processors !---------------------------------------------- if (DebugIO.and.userearranger.and.len>1) then print *,__PIO_FILE__,__LINE__, & File%iosystem%comp_rank,': write IOBUF r8', & IOBUF(1:2),' ...',IOBUF(len-1:len), & glen,len,offset ! write an ascii version write(10+File%iosystem%comp_rank,*) IOBUF(1:len) close(10+File%iosystem%comp_rank) endif ! this is a time dependent multidimensional array if(vardesc%rec>=0 .and. iodesc%start(1)>=0) then ndims = size(iodesc%start)+1 call alloc_check(start,ndims) call alloc_check(count,ndims) start(1:ndims-1)=iodesc%start count(1:ndims-1)=iodesc%count start(ndims:ndims)=vardesc%rec count(ndims:ndims)=1 if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count,' ndims:',ndims, minval(iobuf),maxval(iobuf) ! this is a timedependent single value else if(vardesc%rec>=0) then call alloc_check(start,1) call alloc_check(count,1) start(1) = int(vardesc%rec,kind=PIO_Offset) count(1) = 1_PIO_Offset if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count ! this is a non-timedependent array else ndims = size(iodesc%start) call alloc_check(start,ndims) call alloc_check(count,ndims) start=iodesc%start count=iodesc%count if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count,' ndims:',ndims end if else ! some compilers have problems passing ! unassociated pointers when they are intent in call alloc_check(start, 0) call alloc_check(count, 0) endif #ifdef TIMING call t_startf("pio_write_nf") #endif ierr = write_nf(File,IOBUF,varDesc,iodesc,start,count) #ifdef TIMING call t_stopf("pio_write_nf") #endif call dealloc_check(start) call dealloc_check(count) if(UseRearranger) call dealloc_check(IOBUF) ! call MPI_Barrier(File%iosystem%comp_comm,ierr) !-------------------------- ! set the error return code !-------------------------- iostat=ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_write_darray") #endif end subroutine write_darray_nf_real ! TYPE real,int,double !> !! @private !! @brief Write a 1D array of type int defined by varDesc using the decomposition described in iodesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 300 "piodarray.F90.in" subroutine write_darray_nf_int (File,varDesc,ioDesc,array, iostat, fillval) ! !DESCRIPTION: ! Writes a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor integer(i4), target, intent(in) :: & array(:) ! array to be written integer(i4), optional, intent(in) :: fillval ! rearrange receiver fill value integer(pio_offset), pointer :: start(:), count(:) integer(i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::write_darray_nf_int' integer(i4), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if IO processor integer (i4) :: len, &! length of IO decomp segmap glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform ndims ! number of variable dimensions integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind) :: UseRearranger #if DEBUG_REARR integer(i4), dimension(:), pointer :: array2 integer i #endif integer(i4) :: rsum integer(i4) :: ierr #ifdef TIMING call t_startf("pio_write_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! --------------------------------------------------------- ! pull information out of the decomposition data structure ! --------------------------------------------------------- glen = iodesc%glen offset = iodesc%IOmap%start len = iodesc%IOmap%length if(Debug) & print *,__PIO_FILE__,__LINE__,' NAME : IAM: ',File%iosystem%comp_rank, & ' UseRearranger: ',UseRearranger, glen, offset, len #ifdef TIMING call t_startf("pio_rearrange_write") #endif if(UseRearranger) then if (IOproc) then if(Debug) print *, subName,': IAM: ',File%iosystem%comp_rank,'Before call to allocate(IOBUF): ',len call alloc_check(IOBUF,len,' TYPE :IOBUF') if (present(fillval)) then IOBUF=fillval else IOBUF= -1.0_r8 endif !------------------------------------------------ ! set the IO buffer to a particular test pattern !------------------------------------------------ !JMD IOBUF(:) = real(File%iosystem%io_rank,kind=r8) if(Debug) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len !,' IOBUF: ',IOBUF else call alloc_check(IOBUF,0) IOBUF= -1.0_r8 endif !------------------------------------ ! Rearrange data from comp->io decomp !------------------------------------ ! "array" is comp data call rearrange_comp2io(File%iosystem,iodesc, array, iobuf) #if DEBUG_REARR call alloc_check(array2,size(array),'array2') call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array2) do i=1,size(array) if (array(i) /= array2(i)) then print *, subName,': error: write ping-pong test failed on index',i end if end do print *, subName,': passed write ping-pong test' call dealloc_check(array2) !!!!!!! end debug #endif !-------------------------------------------- ! End data rearrange !-------------------------------------------- else iobuf=>array endif ! if(UseRearranger) #ifdef TIMING call t_stopf("pio_rearrange_write") #endif if (IOproc) then !---------------------------------------------- ! write the global 2-d slice from IO processors !---------------------------------------------- if (DebugIO.and.userearranger.and.len>1) then print *,__PIO_FILE__,__LINE__, & File%iosystem%comp_rank,': write IOBUF r8', & IOBUF(1:2),' ...',IOBUF(len-1:len), & glen,len,offset ! write an ascii version write(10+File%iosystem%comp_rank,*) IOBUF(1:len) close(10+File%iosystem%comp_rank) endif ! this is a time dependent multidimensional array if(vardesc%rec>=0 .and. iodesc%start(1)>=0) then ndims = size(iodesc%start)+1 call alloc_check(start,ndims) call alloc_check(count,ndims) start(1:ndims-1)=iodesc%start count(1:ndims-1)=iodesc%count start(ndims:ndims)=vardesc%rec count(ndims:ndims)=1 if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count,' ndims:',ndims, minval(iobuf),maxval(iobuf) ! this is a timedependent single value else if(vardesc%rec>=0) then call alloc_check(start,1) call alloc_check(count,1) start(1) = int(vardesc%rec,kind=PIO_Offset) count(1) = 1_PIO_Offset if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count ! this is a non-timedependent array else ndims = size(iodesc%start) call alloc_check(start,ndims) call alloc_check(count,ndims) start=iodesc%start count=iodesc%count if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count,' ndims:',ndims end if else ! some compilers have problems passing ! unassociated pointers when they are intent in call alloc_check(start, 0) call alloc_check(count, 0) endif #ifdef TIMING call t_startf("pio_write_nf") #endif ierr = write_nf(File,IOBUF,varDesc,iodesc,start,count) #ifdef TIMING call t_stopf("pio_write_nf") #endif call dealloc_check(start) call dealloc_check(count) if(UseRearranger) call dealloc_check(IOBUF) ! call MPI_Barrier(File%iosystem%comp_comm,ierr) !-------------------------- ! set the error return code !-------------------------- iostat=ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_write_darray") #endif end subroutine write_darray_nf_int ! TYPE real,int,double !> !! @private !! @brief Write a 1D array of type double defined by varDesc using the decomposition described in iodesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 300 "piodarray.F90.in" subroutine write_darray_nf_double (File,varDesc,ioDesc,array, iostat, fillval) ! !DESCRIPTION: ! Writes a 2-d slab of TYPE to a netcdf file. ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor real(r8), target, intent(in) :: & array(:) ! array to be written real(r8), optional, intent(in) :: fillval ! rearrange receiver fill value integer(pio_offset), pointer :: start(:), count(:) integer(i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::write_darray_nf_double' real(r8), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if IO processor integer (i4) :: len, &! length of IO decomp segmap glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform ndims ! number of variable dimensions integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind) :: UseRearranger #if DEBUG_REARR real(r8), dimension(:), pointer :: array2 integer i #endif real(r8) :: rsum integer(i4) :: ierr #ifdef TIMING call t_startf("pio_write_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! --------------------------------------------------------- ! pull information out of the decomposition data structure ! --------------------------------------------------------- glen = iodesc%glen offset = iodesc%IOmap%start len = iodesc%IOmap%length if(Debug) & print *,__PIO_FILE__,__LINE__,' NAME : IAM: ',File%iosystem%comp_rank, & ' UseRearranger: ',UseRearranger, glen, offset, len #ifdef TIMING call t_startf("pio_rearrange_write") #endif if(UseRearranger) then if (IOproc) then if(Debug) print *, subName,': IAM: ',File%iosystem%comp_rank,'Before call to allocate(IOBUF): ',len call alloc_check(IOBUF,len,' TYPE :IOBUF') if (present(fillval)) then IOBUF=fillval else IOBUF= -1.0_r8 endif !------------------------------------------------ ! set the IO buffer to a particular test pattern !------------------------------------------------ !JMD IOBUF(:) = real(File%iosystem%io_rank,kind=r8) if(Debug) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len !,' IOBUF: ',IOBUF else call alloc_check(IOBUF,0) IOBUF= -1.0_r8 endif !------------------------------------ ! Rearrange data from comp->io decomp !------------------------------------ ! "array" is comp data call rearrange_comp2io(File%iosystem,iodesc, array, iobuf) #if DEBUG_REARR call alloc_check(array2,size(array),'array2') call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array2) do i=1,size(array) if (array(i) /= array2(i)) then print *, subName,': error: write ping-pong test failed on index',i end if end do print *, subName,': passed write ping-pong test' call dealloc_check(array2) !!!!!!! end debug #endif !-------------------------------------------- ! End data rearrange !-------------------------------------------- else iobuf=>array endif ! if(UseRearranger) #ifdef TIMING call t_stopf("pio_rearrange_write") #endif if (IOproc) then !---------------------------------------------- ! write the global 2-d slice from IO processors !---------------------------------------------- if (DebugIO.and.userearranger.and.len>1) then print *,__PIO_FILE__,__LINE__, & File%iosystem%comp_rank,': write IOBUF r8', & IOBUF(1:2),' ...',IOBUF(len-1:len), & glen,len,offset ! write an ascii version write(10+File%iosystem%comp_rank,*) IOBUF(1:len) close(10+File%iosystem%comp_rank) endif ! this is a time dependent multidimensional array if(vardesc%rec>=0 .and. iodesc%start(1)>=0) then ndims = size(iodesc%start)+1 call alloc_check(start,ndims) call alloc_check(count,ndims) start(1:ndims-1)=iodesc%start count(1:ndims-1)=iodesc%count start(ndims:ndims)=vardesc%rec count(ndims:ndims)=1 if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count,' ndims:',ndims, minval(iobuf),maxval(iobuf) ! this is a timedependent single value else if(vardesc%rec>=0) then call alloc_check(start,1) call alloc_check(count,1) start(1) = int(vardesc%rec,kind=PIO_Offset) count(1) = 1_PIO_Offset if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count ! this is a non-timedependent array else ndims = size(iodesc%start) call alloc_check(start,ndims) call alloc_check(count,ndims) start=iodesc%start count=iodesc%count if(Debug) print *, __PIO_FILE__,__LINE__,'start:',start,' count:',count,' ndims:',ndims end if else ! some compilers have problems passing ! unassociated pointers when they are intent in call alloc_check(start, 0) call alloc_check(count, 0) endif #ifdef TIMING call t_startf("pio_write_nf") #endif ierr = write_nf(File,IOBUF,varDesc,iodesc,start,count) #ifdef TIMING call t_stopf("pio_write_nf") #endif call dealloc_check(start) call dealloc_check(count) if(UseRearranger) call dealloc_check(IOBUF) ! call MPI_Barrier(File%iosystem%comp_comm,ierr) !-------------------------- ! set the error return code !-------------------------- iostat=ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_write_darray") #endif end subroutine write_darray_nf_double ! TYPE real,int,double !> !! @private !! @brief Write a 1D array of type real. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 529 "piodarray.F90.in" subroutine write_darray_bin_real(File,varDesc,ioDesc,array, iostat, fillval) ! !DESCRIPTION: ! Writes a 2-d slab of integers to a file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! varable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! IO decomp descriptor real(r4), dimension(:), intent(in), target :: & array ! array to be written real(r4), optional, intent(in) :: fillval integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::write_darray_bin_real' real(r4), dimension(:), pointer :: & IOBUF ! local IO buffer #if DEBUG_REARR real(r4), dimension(:), pointer :: array2 integer i #endif logical (log_kind) :: IOproc ! true if IO processor integer (i4) :: len, &! length of IO decomp segmap glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform varID ! variable ID integer (i4) :: ierr integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind) :: UseRearranger #ifdef TIMING call t_startf("pio_write_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! ------------------------------------------------- ! Pull information about the IO decomposition ! ------------------------------------------------- if(Debug) print *, __PIO_FILE__,__LINE__,iodesc%Write%elemTYPE, iodesc%Write%fileTYPE offset = iodesc%IOmap%start len = iodesc%IOmap%length glen = iodesc%glen ierr = pio_noerr if (IOproc) then if(userearranger) then call alloc_check(IOBUF,len,' TYPE :IOBUF') if (present(fillval)) then IOBUF=fillval else IOBUF=-1.0_r8 endif else iobuf=>array end if !------------------------------------------------ ! set the IO buffer to a particular test pattern !------------------------------------------------ !JMD IOBUF(:) = File%iosystem%io_rank if(Debug) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len ! ,' IOBUF: ',IOBUF else if(userearranger) call alloc_check(IOBUF,0,'write_darray_int:IOBUF') endif !----------------------------------------- !NEED HELP: ! ! Need a call to a data rearranger here ! ! call ESMF_rearrange() ! ! or ! ! call MCT_rearrange() !----------------------------------------- #ifdef TIMING call t_startf("pio_rearrange_write") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from comp->io decomp !------------------------------------ ! "array" is comp data call rearrange_comp2io(File%iosystem,iodesc,array,IOBUF) #if DEBUG_REARR call alloc_check(array2,size(array),'array2') call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array2) do i=1,size(array) if (array(i) /= array2(i)) then print *, subName,': error: int write ping-pong test failed on index',i end if end do print *, subName,': passed int write ping-pong test' call dealloc_check(array2) !!!!!!! end debug #endif !-------------------------------------------- ! End data rearrange !-------------------------------------------- endif #ifdef TIMING call t_stopf("pio_rearrange_write") #endif if(IOProc) then #ifdef TIMING call t_startf("pio_write_bin") #endif !---------------------------------------------- ! write the global 2-d slice from IO processors !---------------------------------------------- ierr = write_mpiio(File,IOBUF,varDesc,iodesc) #ifdef TIMING call t_stopf("pio_write_bin") #endif endif !-------------------------- ! deallocate the IO buffer !-------------------------- if(userearranger) call dealloc_check(IOBUF) ! call MPI_Barrier(File%iosystem%comp_comm,ierr) !-------------------------- ! set the error return code !-------------------------- iostat=ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_write_darray") #endif end subroutine write_darray_bin_real ! TYPE real,int,double !> !! @private !! @brief Write a 1D array of type int. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 529 "piodarray.F90.in" subroutine write_darray_bin_int(File,varDesc,ioDesc,array, iostat, fillval) ! !DESCRIPTION: ! Writes a 2-d slab of integers to a file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! varable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! IO decomp descriptor integer(i4), dimension(:), intent(in), target :: & array ! array to be written integer(i4), optional, intent(in) :: fillval integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::write_darray_bin_int' integer(i4), dimension(:), pointer :: & IOBUF ! local IO buffer #if DEBUG_REARR integer(i4), dimension(:), pointer :: array2 integer i #endif logical (log_kind) :: IOproc ! true if IO processor integer (i4) :: len, &! length of IO decomp segmap glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform varID ! variable ID integer (i4) :: ierr integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind) :: UseRearranger #ifdef TIMING call t_startf("pio_write_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! ------------------------------------------------- ! Pull information about the IO decomposition ! ------------------------------------------------- if(Debug) print *, __PIO_FILE__,__LINE__,iodesc%Write%elemTYPE, iodesc%Write%fileTYPE offset = iodesc%IOmap%start len = iodesc%IOmap%length glen = iodesc%glen ierr = pio_noerr if (IOproc) then if(userearranger) then call alloc_check(IOBUF,len,' TYPE :IOBUF') if (present(fillval)) then IOBUF=fillval else IOBUF=-1.0_r8 endif else iobuf=>array end if !------------------------------------------------ ! set the IO buffer to a particular test pattern !------------------------------------------------ !JMD IOBUF(:) = File%iosystem%io_rank if(Debug) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len ! ,' IOBUF: ',IOBUF else if(userearranger) call alloc_check(IOBUF,0,'write_darray_int:IOBUF') endif !----------------------------------------- !NEED HELP: ! ! Need a call to a data rearranger here ! ! call ESMF_rearrange() ! ! or ! ! call MCT_rearrange() !----------------------------------------- #ifdef TIMING call t_startf("pio_rearrange_write") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from comp->io decomp !------------------------------------ ! "array" is comp data call rearrange_comp2io(File%iosystem,iodesc,array,IOBUF) #if DEBUG_REARR call alloc_check(array2,size(array),'array2') call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array2) do i=1,size(array) if (array(i) /= array2(i)) then print *, subName,': error: int write ping-pong test failed on index',i end if end do print *, subName,': passed int write ping-pong test' call dealloc_check(array2) !!!!!!! end debug #endif !-------------------------------------------- ! End data rearrange !-------------------------------------------- endif #ifdef TIMING call t_stopf("pio_rearrange_write") #endif if(IOProc) then #ifdef TIMING call t_startf("pio_write_bin") #endif !---------------------------------------------- ! write the global 2-d slice from IO processors !---------------------------------------------- ierr = write_mpiio(File,IOBUF,varDesc,iodesc) #ifdef TIMING call t_stopf("pio_write_bin") #endif endif !-------------------------- ! deallocate the IO buffer !-------------------------- if(userearranger) call dealloc_check(IOBUF) ! call MPI_Barrier(File%iosystem%comp_comm,ierr) !-------------------------- ! set the error return code !-------------------------- iostat=ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_write_darray") #endif end subroutine write_darray_bin_int ! TYPE real,int,double !> !! @private !! @brief Write a 1D array of type double. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The data to be written !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !! @param fillval : An optional fill value to fill holes in the data written !< # 529 "piodarray.F90.in" subroutine write_darray_bin_double(File,varDesc,ioDesc,array, iostat, fillval) ! !DESCRIPTION: ! Writes a 2-d slab of integers to a file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! file information type (var_desc_t), intent(inout) :: & varDesc ! varable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! IO decomp descriptor real(r8), dimension(:), intent(in), target :: & array ! array to be written real(r8), optional, intent(in) :: fillval integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::write_darray_bin_double' real(r8), dimension(:), pointer :: & IOBUF ! local IO buffer #if DEBUG_REARR real(r8), dimension(:), pointer :: array2 integer i #endif logical (log_kind) :: IOproc ! true if IO processor integer (i4) :: len, &! length of IO decomp segmap glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform varID ! variable ID integer (i4) :: ierr integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind) :: UseRearranger #ifdef TIMING call t_startf("pio_write_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! ------------------------------------------------- ! Pull information about the IO decomposition ! ------------------------------------------------- if(Debug) print *, __PIO_FILE__,__LINE__,iodesc%Write%elemTYPE, iodesc%Write%fileTYPE offset = iodesc%IOmap%start len = iodesc%IOmap%length glen = iodesc%glen ierr = pio_noerr if (IOproc) then if(userearranger) then call alloc_check(IOBUF,len,' TYPE :IOBUF') if (present(fillval)) then IOBUF=fillval else IOBUF=-1.0_r8 endif else iobuf=>array end if !------------------------------------------------ ! set the IO buffer to a particular test pattern !------------------------------------------------ !JMD IOBUF(:) = File%iosystem%io_rank if(Debug) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len ! ,' IOBUF: ',IOBUF else if(userearranger) call alloc_check(IOBUF,0,'write_darray_int:IOBUF') endif !----------------------------------------- !NEED HELP: ! ! Need a call to a data rearranger here ! ! call ESMF_rearrange() ! ! or ! ! call MCT_rearrange() !----------------------------------------- #ifdef TIMING call t_startf("pio_rearrange_write") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from comp->io decomp !------------------------------------ ! "array" is comp data call rearrange_comp2io(File%iosystem,iodesc,array,IOBUF) #if DEBUG_REARR call alloc_check(array2,size(array),'array2') call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array2) do i=1,size(array) if (array(i) /= array2(i)) then print *, subName,': error: int write ping-pong test failed on index',i end if end do print *, subName,': passed int write ping-pong test' call dealloc_check(array2) !!!!!!! end debug #endif !-------------------------------------------- ! End data rearrange !-------------------------------------------- endif #ifdef TIMING call t_stopf("pio_rearrange_write") #endif if(IOProc) then #ifdef TIMING call t_startf("pio_write_bin") #endif !---------------------------------------------- ! write the global 2-d slice from IO processors !---------------------------------------------- ierr = write_mpiio(File,IOBUF,varDesc,iodesc) #ifdef TIMING call t_stopf("pio_write_bin") #endif endif !-------------------------- ! deallocate the IO buffer !-------------------------- if(userearranger) call dealloc_check(IOBUF) ! call MPI_Barrier(File%iosystem%comp_comm,ierr) !-------------------------- ! set the error return code !-------------------------- iostat=ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_write_darray") #endif end subroutine write_darray_bin_double ! TYPE real,int,double !> !! @private !! @brief Read a 1D array of type real defined by varDesc using the decomposition described in ioDesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 719 "piodarray.F90.in" subroutine read_darray_nf_real (File,varDesc,ioDesc,array, iostat) ! ! !DESCRIPTION: ! Reads a 2-d horizontal slice of integers from a binary file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! info about data file type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor ! !INPUT/OUTPUT PARAMETERS: real(r4), dimension(:), intent(out), target :: & array ! array to be read integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::read_darray_nf_real' real(r4), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if an IO processor integer (i4) :: len, &! local length of IO decomp glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform ndims integer(pio_offset), dimension(PIO_MAX_VAR_DIMS) :: start, count integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind), parameter :: Debug = .FALSE. logical(log_kind) :: UseRearranger integer(i4) :: ierr #if DEBUG_REARR real(r4), dimension(:), pointer :: iobuf2 integer i #endif #ifdef TIMING call t_startf("pio_read_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ierr = PIO_NOERR ! ----------------------------------------------------- ! Pull out information of DecompMap_t data structure ! ----------------------------------------------------- glen = iodesc%glen offset = iodesc%IOmap%start len = iodesc%IOmap%length start=0 count=0 if (IOproc) then !----------------------------- ! allocate temporary IO buffer !----------------------------- if(userearranger) then if(File%iotype == pio_iotype_netcdf .and. file%iosystem%io_rank==0) then call alloc_check(IOBUF,iodesc%maxiobuflen,' TYPE :IOBUF') else call alloc_check(IOBUF,len,' TYPE :IOBUF') end if else iobuf=>array end if !---------------------------------------------- ! read the global 2-d slice to IO processors !---------------------------------------------- ! this is a time dependent multidimensional array if(vardesc%rec>=0 .and. iodesc%start(1)>=0) then ndims = size(iodesc%start)+1 start(1:ndims-1)=iodesc%start count(1:ndims-1)=iodesc%count start(ndims:ndims)=vardesc%rec count(ndims:ndims)=1 ! this is a timedependent single value else if(vardesc%rec>=0) then ndims=1 start(1) = int(vardesc%rec,kind=PIO_Offset) count(1) = 1_PIO_Offset ! this is a non-timedependent array else ndims = size(iodesc%start) start(1:ndims)=iodesc%start count(1:ndims)=iodesc%count end if else ndims=1 if(userearranger) then call alloc_check(IOBUF,0,'IOBUF') end if endif #ifdef TIMING call t_startf("pio_read_nf") #endif ierr = read_nf(File,IOBUF,varDesc,iodesc,start(1:ndims),count(1:ndims)) #ifdef TIMING call t_stopf("pio_read_nf") #endif if(DebugIO) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len !,' IOBUF: ',IOBUF #ifdef TIMING call t_startf("pio_rearrange_read") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from io->comp decomp !------------------------------------ ! "array" is comp data call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array) #if DEBUG_REARR call alloc_check(iobuf2,size(IOBUF),'iobuf2') call rearrange_comp2io(File%iosystem,iodesc,array,iobuf2) do i=1,size(iobuf) if (iobuf(i) /= iobuf2(i)) then print *, subName,': error: int read ping-pong test failed on index',i end if end do print *, subName,': passed int read ping-pong test' call dealloc_check(iobuf2) !!!!!!! end debug #endif ! -------------------------- ! deallocate IO buffer ! -------------------------- call dealloc_check(IOBUF) endif #ifdef TIMING call t_stopf("pio_rearrange_read") #endif !---------------- ! set errror code !---------------- iostat = ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_read_darray") #endif end subroutine read_darray_nf_real ! TYPE real,int,double !> !! @private !! @brief Read a 1D array of type int defined by varDesc using the decomposition described in ioDesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 719 "piodarray.F90.in" subroutine read_darray_nf_int (File,varDesc,ioDesc,array, iostat) ! ! !DESCRIPTION: ! Reads a 2-d horizontal slice of integers from a binary file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! info about data file type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor ! !INPUT/OUTPUT PARAMETERS: integer(i4), dimension(:), intent(out), target :: & array ! array to be read integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::read_darray_nf_int' integer(i4), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if an IO processor integer (i4) :: len, &! local length of IO decomp glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform ndims integer(pio_offset), dimension(PIO_MAX_VAR_DIMS) :: start, count integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind), parameter :: Debug = .FALSE. logical(log_kind) :: UseRearranger integer(i4) :: ierr #if DEBUG_REARR integer(i4), dimension(:), pointer :: iobuf2 integer i #endif #ifdef TIMING call t_startf("pio_read_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ierr = PIO_NOERR ! ----------------------------------------------------- ! Pull out information of DecompMap_t data structure ! ----------------------------------------------------- glen = iodesc%glen offset = iodesc%IOmap%start len = iodesc%IOmap%length start=0 count=0 if (IOproc) then !----------------------------- ! allocate temporary IO buffer !----------------------------- if(userearranger) then if(File%iotype == pio_iotype_netcdf .and. file%iosystem%io_rank==0) then call alloc_check(IOBUF,iodesc%maxiobuflen,' TYPE :IOBUF') else call alloc_check(IOBUF,len,' TYPE :IOBUF') end if else iobuf=>array end if !---------------------------------------------- ! read the global 2-d slice to IO processors !---------------------------------------------- ! this is a time dependent multidimensional array if(vardesc%rec>=0 .and. iodesc%start(1)>=0) then ndims = size(iodesc%start)+1 start(1:ndims-1)=iodesc%start count(1:ndims-1)=iodesc%count start(ndims:ndims)=vardesc%rec count(ndims:ndims)=1 ! this is a timedependent single value else if(vardesc%rec>=0) then ndims=1 start(1) = int(vardesc%rec,kind=PIO_Offset) count(1) = 1_PIO_Offset ! this is a non-timedependent array else ndims = size(iodesc%start) start(1:ndims)=iodesc%start count(1:ndims)=iodesc%count end if else ndims=1 if(userearranger) then call alloc_check(IOBUF,0,'IOBUF') end if endif #ifdef TIMING call t_startf("pio_read_nf") #endif ierr = read_nf(File,IOBUF,varDesc,iodesc,start(1:ndims),count(1:ndims)) #ifdef TIMING call t_stopf("pio_read_nf") #endif if(DebugIO) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len !,' IOBUF: ',IOBUF #ifdef TIMING call t_startf("pio_rearrange_read") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from io->comp decomp !------------------------------------ ! "array" is comp data call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array) #if DEBUG_REARR call alloc_check(iobuf2,size(IOBUF),'iobuf2') call rearrange_comp2io(File%iosystem,iodesc,array,iobuf2) do i=1,size(iobuf) if (iobuf(i) /= iobuf2(i)) then print *, subName,': error: int read ping-pong test failed on index',i end if end do print *, subName,': passed int read ping-pong test' call dealloc_check(iobuf2) !!!!!!! end debug #endif ! -------------------------- ! deallocate IO buffer ! -------------------------- call dealloc_check(IOBUF) endif #ifdef TIMING call t_stopf("pio_rearrange_read") #endif !---------------- ! set errror code !---------------- iostat = ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_read_darray") #endif end subroutine read_darray_nf_int ! TYPE real,int,double !> !! @private !! @brief Read a 1D array of type double defined by varDesc using the decomposition described in ioDesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 719 "piodarray.F90.in" subroutine read_darray_nf_double (File,varDesc,ioDesc,array, iostat) ! ! !DESCRIPTION: ! Reads a 2-d horizontal slice of integers from a binary file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! info about data file type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor ! !INPUT/OUTPUT PARAMETERS: real(r8), dimension(:), intent(out), target :: & array ! array to be read integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::read_darray_nf_double' real(r8), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if an IO processor integer (i4) :: len, &! local length of IO decomp glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype, &! type of IO to perform ndims integer(pio_offset), dimension(PIO_MAX_VAR_DIMS) :: start, count integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind), parameter :: Debug = .FALSE. logical(log_kind) :: UseRearranger integer(i4) :: ierr #if DEBUG_REARR real(r8), dimension(:), pointer :: iobuf2 integer i #endif #ifdef TIMING call t_startf("pio_read_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ierr = PIO_NOERR ! ----------------------------------------------------- ! Pull out information of DecompMap_t data structure ! ----------------------------------------------------- glen = iodesc%glen offset = iodesc%IOmap%start len = iodesc%IOmap%length start=0 count=0 if (IOproc) then !----------------------------- ! allocate temporary IO buffer !----------------------------- if(userearranger) then if(File%iotype == pio_iotype_netcdf .and. file%iosystem%io_rank==0) then call alloc_check(IOBUF,iodesc%maxiobuflen,' TYPE :IOBUF') else call alloc_check(IOBUF,len,' TYPE :IOBUF') end if else iobuf=>array end if !---------------------------------------------- ! read the global 2-d slice to IO processors !---------------------------------------------- ! this is a time dependent multidimensional array if(vardesc%rec>=0 .and. iodesc%start(1)>=0) then ndims = size(iodesc%start)+1 start(1:ndims-1)=iodesc%start count(1:ndims-1)=iodesc%count start(ndims:ndims)=vardesc%rec count(ndims:ndims)=1 ! this is a timedependent single value else if(vardesc%rec>=0) then ndims=1 start(1) = int(vardesc%rec,kind=PIO_Offset) count(1) = 1_PIO_Offset ! this is a non-timedependent array else ndims = size(iodesc%start) start(1:ndims)=iodesc%start count(1:ndims)=iodesc%count end if else ndims=1 if(userearranger) then call alloc_check(IOBUF,0,'IOBUF') end if endif #ifdef TIMING call t_startf("pio_read_nf") #endif ierr = read_nf(File,IOBUF,varDesc,iodesc,start(1:ndims),count(1:ndims)) #ifdef TIMING call t_stopf("pio_read_nf") #endif if(DebugIO) print *, subName,': {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'offset: ',offset,'len: ',len !,' IOBUF: ',IOBUF #ifdef TIMING call t_startf("pio_rearrange_read") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from io->comp decomp !------------------------------------ ! "array" is comp data call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array) #if DEBUG_REARR call alloc_check(iobuf2,size(IOBUF),'iobuf2') call rearrange_comp2io(File%iosystem,iodesc,array,iobuf2) do i=1,size(iobuf) if (iobuf(i) /= iobuf2(i)) then print *, subName,': error: int read ping-pong test failed on index',i end if end do print *, subName,': passed int read ping-pong test' call dealloc_check(iobuf2) !!!!!!! end debug #endif ! -------------------------- ! deallocate IO buffer ! -------------------------- call dealloc_check(IOBUF) endif #ifdef TIMING call t_stopf("pio_rearrange_read") #endif !---------------- ! set errror code !---------------- iostat = ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_read_darray") #endif end subroutine read_darray_nf_double ! TYPE real,int,double !> !! @private !! @brief Read an array of type real defined by varDesc using the decomposition !! described in ioDesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 916 "piodarray.F90.in" subroutine read_darray_bin_real (File,varDesc,ioDesc,array, iostat) ! ! !DESCRIPTION: ! Reads a 2-d horizontal slice of integers from a binary file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! info about data file type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor ! !INPUT/OUTPUT PARAMETERS: real(r4), dimension(:), intent(out), target :: & array ! array to be read integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::read_darray_bin_real' real(r4), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if an IO processor integer (i4) :: len, &! local length of IO decomp glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype ! type of IO to perform integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind), parameter :: Debug = .FALSE. logical(log_kind) :: UseRearranger integer(i4) :: ierr #if DEBUG_REARR real(r4), dimension(:), pointer :: iobuf2 integer i #endif #ifdef TIMING call t_startf("pio_read_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! ----------------------------------------------------- ! Pull out information of DecompMap_t data structure ! ----------------------------------------------------- glen = iodesc%glen ! offset = iodesc%IOmap%start ! len = iodesc%IOmap%length len = iodesc%Read%n_words ierr = pio_noerr if (IOproc) then !----------------------------- ! allocate temporary IO buffer !----------------------------- if(userearranger) then call alloc_check(IOBUF,len,'read_darray_ :IOBUF') else iobuf=>array end if !---------------------------------------------- ! read the global 2-d slice to IO processors !---------------------------------------------- if(iotype.eq.pio_iotype_binary) then print *, subName,': TYPE : IO type not supported' iostat =-1 return end if #ifdef TIMING call t_startf("pio_read_bin") #endif ierr = read_mpiio(File,IOBUF,varDesc,iodesc) #ifdef TIMING call t_stopf("pio_read_bin") #endif if(DebugIO) print *, subName,': TYPE: {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'len: ',len,' IOBUF: ',IOBUF(1:4) else if(userearranger) then call alloc_check(IOBUF,0,'IOBUF') endif #ifdef TIMING call t_startf("pio_rearrange_read") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from io->comp decomp !------------------------------------ ! "array" is comp data call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array) #if DEBUG_REARR call alloc_check(iobuf2,size(IOBUF),'iobuf2') call rearrange_comp2io(File%iosystem,iodesc,array,iobuf2) do i=1,size(iobuf) if (iobuf(i) /= iobuf2(i)) then print *, subName,': error: int read ping-pong test failed on index',i end if end do print *, subName,': passed int read ping-pong test' call dealloc_check(iobuf2) !!!!!!! end debug #endif ! -------------------------- ! deallocate IO buffer ! -------------------------- call dealloc_check(IOBUF) endif #ifdef TIMING call t_stopf("pio_rearrange_read") #endif !---------------- ! set errror code !---------------- iostat = ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_read_darray") #endif end subroutine read_darray_bin_real ! TYPE real,int,double !> !! @private !! @brief Read an array of type int defined by varDesc using the decomposition !! described in ioDesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 916 "piodarray.F90.in" subroutine read_darray_bin_int (File,varDesc,ioDesc,array, iostat) ! ! !DESCRIPTION: ! Reads a 2-d horizontal slice of integers from a binary file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! info about data file type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor ! !INPUT/OUTPUT PARAMETERS: integer(i4), dimension(:), intent(out), target :: & array ! array to be read integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::read_darray_bin_int' integer(i4), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if an IO processor integer (i4) :: len, &! local length of IO decomp glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype ! type of IO to perform integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind), parameter :: Debug = .FALSE. logical(log_kind) :: UseRearranger integer(i4) :: ierr #if DEBUG_REARR integer(i4), dimension(:), pointer :: iobuf2 integer i #endif #ifdef TIMING call t_startf("pio_read_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! ----------------------------------------------------- ! Pull out information of DecompMap_t data structure ! ----------------------------------------------------- glen = iodesc%glen ! offset = iodesc%IOmap%start ! len = iodesc%IOmap%length len = iodesc%Read%n_words ierr = pio_noerr if (IOproc) then !----------------------------- ! allocate temporary IO buffer !----------------------------- if(userearranger) then call alloc_check(IOBUF,len,'read_darray_ :IOBUF') else iobuf=>array end if !---------------------------------------------- ! read the global 2-d slice to IO processors !---------------------------------------------- if(iotype.eq.pio_iotype_binary) then print *, subName,': TYPE : IO type not supported' iostat =-1 return end if #ifdef TIMING call t_startf("pio_read_bin") #endif ierr = read_mpiio(File,IOBUF,varDesc,iodesc) #ifdef TIMING call t_stopf("pio_read_bin") #endif if(DebugIO) print *, subName,': TYPE: {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'len: ',len,' IOBUF: ',IOBUF(1:4) else if(userearranger) then call alloc_check(IOBUF,0,'IOBUF') endif #ifdef TIMING call t_startf("pio_rearrange_read") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from io->comp decomp !------------------------------------ ! "array" is comp data call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array) #if DEBUG_REARR call alloc_check(iobuf2,size(IOBUF),'iobuf2') call rearrange_comp2io(File%iosystem,iodesc,array,iobuf2) do i=1,size(iobuf) if (iobuf(i) /= iobuf2(i)) then print *, subName,': error: int read ping-pong test failed on index',i end if end do print *, subName,': passed int read ping-pong test' call dealloc_check(iobuf2) !!!!!!! end debug #endif ! -------------------------- ! deallocate IO buffer ! -------------------------- call dealloc_check(IOBUF) endif #ifdef TIMING call t_stopf("pio_rearrange_read") #endif !---------------- ! set errror code !---------------- iostat = ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_read_darray") #endif end subroutine read_darray_bin_int ! TYPE real,int,double !> !! @private !! @brief Read an array of type double defined by varDesc using the decomposition !! described in ioDesc to the netcdf or pnetcdf file File. !! @details !! @param File @copydoc file_desc_t !! @param varDesc @copydoc var_desc_t !! @param ioDesc @copydoc io_desc_t !! @param array : The read data !! @param iostat : The status returned from this routine (see \ref PIO_seterrorhandling for details) !< # 916 "piodarray.F90.in" subroutine read_darray_bin_double (File,varDesc,ioDesc,array, iostat) ! ! !DESCRIPTION: ! Reads a 2-d horizontal slice of integers from a binary file ! ! !REVISION HISTORY: ! same as module ! !INPUT PARAMETERS: type (File_desc_t), intent(inout) :: & File ! info about data file type (var_desc_t), intent(inout) :: & varDesc ! variable descriptor type (io_desc_t), intent(inout) :: & ioDesc ! io decomp descriptor ! !INPUT/OUTPUT PARAMETERS: real(r8), dimension(:), intent(out), target :: & array ! array to be read integer (i4), intent(out) :: iostat !EOP !BOC !----------------------------------------------------------------------- ! ! local variables ! !----------------------------------------------------------------------- character(len=*), parameter :: subName=modName//'::read_darray_bin_double' real(r8), dimension(:), pointer :: & IOBUF ! local IO buffer logical (log_kind) :: IOproc ! true if an IO processor integer (i4) :: len, &! local length of IO decomp glen, &! global length of IO decomp in words reclen, &! global length of IO decomp in bytes iotype ! type of IO to perform integer(kind=PIO_OFFSET) :: offset ! local offset of IO decomp logical(log_kind), parameter :: Debug = .FALSE. logical(log_kind) :: UseRearranger integer(i4) :: ierr #if DEBUG_REARR real(r8), dimension(:), pointer :: iobuf2 integer i #endif #ifdef TIMING call t_startf("pio_read_darray") #endif ! ----------------------------------------------------- ! pull information from file_desc_t data structure ! ----------------------------------------------------- IOproc = File%iosystem%IOproc iotype = File%iotype UseRearranger = File%iosystem%UseRearranger ! ----------------------------------------------------- ! Pull out information of DecompMap_t data structure ! ----------------------------------------------------- glen = iodesc%glen ! offset = iodesc%IOmap%start ! len = iodesc%IOmap%length len = iodesc%Read%n_words ierr = pio_noerr if (IOproc) then !----------------------------- ! allocate temporary IO buffer !----------------------------- if(userearranger) then call alloc_check(IOBUF,len,'read_darray_ :IOBUF') else iobuf=>array end if !---------------------------------------------- ! read the global 2-d slice to IO processors !---------------------------------------------- if(iotype.eq.pio_iotype_binary) then print *, subName,': TYPE : IO type not supported' iostat =-1 return end if #ifdef TIMING call t_startf("pio_read_bin") #endif ierr = read_mpiio(File,IOBUF,varDesc,iodesc) #ifdef TIMING call t_stopf("pio_read_bin") #endif if(DebugIO) print *, subName,': TYPE: {comp,io}_rank: ',File%iosystem%comp_rank,File%iosystem%io_rank, & 'len: ',len,' IOBUF: ',IOBUF(1:4) else if(userearranger) then call alloc_check(IOBUF,0,'IOBUF') endif #ifdef TIMING call t_startf("pio_rearrange_read") #endif if(UseRearranger) then !------------------------------------ ! Rearrange data from io->comp decomp !------------------------------------ ! "array" is comp data call rearrange_io2comp(File%iosystem,iodesc,IOBUF,array) #if DEBUG_REARR call alloc_check(iobuf2,size(IOBUF),'iobuf2') call rearrange_comp2io(File%iosystem,iodesc,array,iobuf2) do i=1,size(iobuf) if (iobuf(i) /= iobuf2(i)) then print *, subName,': error: int read ping-pong test failed on index',i end if end do print *, subName,': passed int read ping-pong test' call dealloc_check(iobuf2) !!!!!!! end debug #endif ! -------------------------- ! deallocate IO buffer ! -------------------------- call dealloc_check(IOBUF) endif #ifdef TIMING call t_stopf("pio_rearrange_read") #endif !---------------- ! set errror code !---------------- iostat = ierr !----------------------------------------------------------------------- !EOC #ifdef TIMING call t_stopf("pio_read_darray") #endif end subroutine read_darray_bin_double end module piodarray