in_binary_512i_mod.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: IN_BIN_512I                                            C
!  Purpose: read an array in chunks of 512 bytes (INTEGER WORDS)       C
!                                                                      C
!  Author: P. Nicoletti                               Date: 02-JAN-92  C
!  Reviewer: P. Nicoletti, W. Rogers, M. Syamlal      Date: 24-JAN-92  C
!                                                                      C
!  Revision Number:                                                    C
!  Purpose:                                                            C
!  Author:                                            Date: dd-mmm-yy  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables: NWORDS, L, NSEG, NREM, LC, N1, N2                  C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
MODULE in_binary_512i

CONTAINS

      SUBROUTINE in_bin_512i(IUNIT, ARRAY, NN, NEXT_REC)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE machine
      IMPLICIT NONE
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      array to write out
      INTEGER          ARRAY(*)
!
!                      output unit number
      INTEGER          IUNIT
!
!                      number of elements in ARRAY
      INTEGER          NN
!
!                      next record number in direct access output file
      INTEGER          NEXT_REC
!
! local variables
!
!                      number of words for 512 bytes (nwords * 4 = 512)
      INTEGER          NWORDS
!
!                      loop counter
      INTEGER          L
!
!                      number of full 512 byte segments need to write N
!                      double precision words
      INTEGER          NSEG
!
!                      number of double precision words in the partially
!                      filled last record
      INTEGER          NREM
!
!                      loop counter
      INTEGER          LC
!
!                      write out array elements N1 to N2
      INTEGER          N1 , N2
!-----------------------------------------------
!
      nwords = nwords_i
      IF (nn <= nwords) THEN
         READ (iunit, rec=next_rec) (array(l),l=1,nn)
         next_rec = next_rec + 1
         RETURN
      ENDIF

      nseg = nn/nwords
      nrem = mod(nn,nwords)
      n1 = 1
      n2 = nwords
!
! read the full 512 byte segments
!
      DO lc = 1, nseg
         READ (iunit, rec=next_rec) (array(l),l=n1,n2)
         n1 = n1 + nwords
         n2 = n2 + nwords
         next_rec = next_rec + 1
      END DO
      IF (nrem /= 0) THEN
         READ (iunit, rec=next_rec) (array(l),l=n1,nn)
         next_rec = next_rec + 1
      ENDIF

      RETURN
      END SUBROUTINE in_bin_512i

      subroutine convert_from_io_i(arr_io,arr_internal,nn)

      use geometry
      use indices
      use compar
      use functions

      implicit none

      integer, intent(in) :: arr_io(:)
      integer, intent(out) :: arr_internal(:)
      integer   nn,i,j,k,ijk,ijk_io

!     write(*,*) 'C0:',C0
!     write(*,*) 'C1:',C1
!     write(*,*) 'C2:',C2

!     write(*,*) 'io:',size(arr_io)
!     write(*,*) 'int:',size(arr_internal)

     if(size(arr_io) == size(arr_internal)) then
         arr_internal = arr_io
     else
         do k = 1,kmax2
         do j = 1,jmax2
         do i = 1,imax2
            ijk    = funijk_gl(i,j,k)
            ijk_io = funijk_io(i,j,k)
            arr_internal(ijk) = arr_io(ijk_io)
           write(*,*)i,j,k,ijk, ijk_io
         end do
         end do
         end do
     endif

      return
    end subroutine convert_from_io_i

      subroutine convert_to_io_i(arr_internal,arr_io,nn)

      use geometry
      use indices
      use compar
      use functions

      implicit none

      integer   arr_io(*) , arr_internal(*)
      integer   nn,i,j,k,ijk,ijk_io

      do k = 1,kmax2
         do j = 1,jmax2
            do i = 1,imax2
               ijk    = funijk_gl(i,j,k)
               ijk_io = funijk_io(i,j,k)
               arr_io(ijk_io) = arr_internal(ijk)
            end do
         end do
      end do

      return
    end subroutine convert_to_io_i

      subroutine convert_to_io_c(arr_internal,arr_io,nn)

      use geometry
      use indices
      use compar
      use functions

      implicit none

      character(LEN=4)   arr_io(*) , arr_internal(*)
      integer       nn,i,j,k,ijk,ijk_io

      do k = 1,kmax2
         do j = 1,jmax2
            do i = 1,imax2
               ijk    = funijk_gl(i,j,k)
               ijk_io = funijk_io(i,j,k)
               arr_io(ijk_io) = arr_internal(ijk)
            end do
         end do
      end do
!
      return
    end subroutine convert_to_io_c

END MODULE in_binary_512i

!// Comments on the modifications for DMP version implementation
!// 001 Include header file and common declarations for parallelization