out_array_k.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: OUT_ARRAY_K                                            C
!  Purpose: print out a 2D (constant k-plane) array to standard output C
!                                                                      C
!  Author: P.Nicoletti                                Date: 02-DEC-91  C
!  Reviewer: W. Rogers, M. Syamlal, S. Venkatesan     Date: 31-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: IMAX2, JMAX2                                  C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables: NCOL, NTAB, LL1, LL2, LL3, IFORM1, IFORM2, IJK, IJ2
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      SUBROUTINE out_array_k(ARRAY)
!...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 param
      USE param1
      USE geometry
      USE fldvar
      USE physprop
      USE indices
      USE funits
      USE compar
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
      DOUBLE PRECISION ARRAY(*)
!
!                      number of columns to print out across the page
      INTEGER          NCOL
!
!                      number of tables the 2D array must be split into
!                      for printing
      INTEGER          NTAB
!
!                      loop indices
      INTEGER          LL1, LL2, LL3
!
!                      start and end 'I' for current table
      INTEGER          IFORM1 , IFORM2
!
!                      start 'IJ' and end 'IJ' for a given 'J' to print out
      INTEGER          IJK , IJ2
!

!
!-----------------------------------------------
!
! NOTE:  IF NCOL IS CHANGED TO A NUMBER GREATER THAN 30, THEN THE "30"
!        IN FORMATS 5050 AND 5100 MUST BE CHANGED TO THAT NUMBER.
!
      ncol = 10
!// Adjust for cyclic in x direction
    IF(cyclic_x) then
      ntab = (imax2-1)/ncol + 1
    ELSE
      ntab = imax2/ncol + 1
    ENDIF

!// Adjust for cyclic in x direction
    IF(cyclic_x) then
      IF (mod(imax2-1,ncol) == 0) ntab = ntab - 1
    ELSE
      IF (mod(imax2,ncol) == 0) ntab = ntab - 1
    ENDIF
!
      DO ll1 = 1, ntab

!// Adjust for cyclic in x direction
       IF(cyclic_x.AND.ll1.eq.1) then
         iform1 = 2 + ncol*(ll1 - 1)
       ELSE
         iform1 = 1 + ncol*(ll1 - 1)
       ENDIF

         iform2 = ncol*ll1

!// Adjust for cyclic in x direction
       IF(cyclic_x) then
         iform2 = min(iform2,imax2-1)
       ELSE
         iform2 = min(iform2,imax2)
       ENDIF

         WRITE (unit_out, 5050) (ll3,ll3=iform1,iform2)
         DO ll2 = jmax2, 1, -1
            ijk = funijk_io(iform1,ll2,1)
            ij2 = funijk_io(iform2,ll2,1)
!efd
!            WRITE (UNIT_OUT, 5100) LL2, (ARRAY(LL3),LL3=IJK,IJ2)
            WRITE (unit_out, 5100) ll2,  &
                     (array(funijk_io(ll3,ll2,1)),ll3=iform1,iform2)
         END DO
      END DO
 5050 FORMAT(3x,'J',3x,'I=',3x,10(i3,9x))
 5100 FORMAT(1x,i3,3x,10(1pe12.4))
      RETURN
      END SUBROUTINE out_array_k