compare01/src/plot_field.F

C $Id: plot_field.F,v 1.1 1998/05/25 20:21:05 cnh Exp $
#include "CPP_OPTIONS.h"
#include "CPP_MACROS.h"
      SUBROUTINE PLOT_FIELD ( A, Nx, Ny )
      IMPLICIT NONE
C
      INTEGER Nx, Ny
      REAL    A(Nx,Ny)
C
      INTEGER NC
      PARAMETER ( NC = 28 )
      CHARACTER*(NC) CVAL
      PARAMETER ( CVAL =
     &'-abcdefghijklmnopqrstuvwxyz+' )
      REAL AMN, AMX, CINT
      CHARACTER BUFFER(1024)
      INTEGER I, J, IDX
C
      AMN =  1.D32
      AMX = -1.D32
      DO J=1,Ny
       DO I=1,Nx
        IF( A(I,J).NE.0.) THEN
         IF ( A(I,J) .LT. AMN ) AMN = A(I,J)
         IF ( A(I,J) .GT. AMX ) AMX = A(I,J)
        ENDIF
       ENDDO
      ENDDO
      CINT = (AMX-AMN)/FLOAT(NC-1)
      WRITE (0,*) ' Min = ', AMN,' Max = ',AMX, ' CINT = ', CINT
      DO J = Ny, 1 , -1
       DO I = 1, Nx
        IF ( CINT .NE. 0. ) THEN
         IDX = NINT(  FLOAT(NC-1)*(A(I,J)-AMN)/(AMX-AMN) )+1
         BUFFER(I) = CVAL(IDX:IDX)
        ELSE
         BUFFER(I) = '*'
        ENDIF
        IF ( A(I,J) .EQ. 0. ) BUFFER(I) = '.'
C       IF ( A(I,J) .LT. 0. ) BUFFER(I) = '-'
C       IF ( A(I,J) .GT. 0. ) BUFFER(I) = '+'
       ENDDO
       WRITE(0,'(I3,1X,1024A)') J,(BUFFER(I),I=1,Nx)
      ENDDO
C
      END
1	adcroft	1.2	C $Id: plot_field.F,v 1.1 1998/05/25 20:21:05 cnh Exp $
2	cnh	1.1	#include "CPP_OPTIONS.h"
3			#include "CPP_MACROS.h"
4			SUBROUTINE PLOT_FIELD ( A, Nx, Ny )
5			IMPLICIT NONE
6			C
7			INTEGER Nx, Ny
8			REAL A(Nx,Ny)
9			C
10			INTEGER NC
11			PARAMETER ( NC = 28 )
12			CHARACTER*(NC) CVAL
13			PARAMETER ( CVAL =
14			&'-abcdefghijklmnopqrstuvwxyz+' )
15			REAL AMN, AMX, CINT
16			CHARACTER BUFFER(1024)
17			INTEGER I, J, IDX
18			C
19			AMN = 1.D32
20			AMX = -1.D32
21			DO J=1,Ny
22			DO I=1,Nx
23			IF( A(I,J).NE.0.) THEN
24			IF ( A(I,J) .LT. AMN ) AMN = A(I,J)
25			IF ( A(I,J) .GT. AMX ) AMX = A(I,J)
26			ENDIF
27			ENDDO
28			ENDDO
29			CINT = (AMX-AMN)/FLOAT(NC-1)
30			WRITE (0,*) ' Min = ', AMN,' Max = ',AMX, ' CINT = ', CINT
31			DO J = Ny, 1 , -1
32			DO I = 1, Nx
33			IF ( CINT .NE. 0. ) THEN
34			IDX = NINT( FLOAT(NC-1)*(A(I,J)-AMN)/(AMX-AMN) )+1
35			BUFFER(I) = CVAL(IDX:IDX)
36			ELSE
37			BUFFER(I) = '*'
38			ENDIF
39			IF ( A(I,J) .EQ. 0. ) BUFFER(I) = '.'
40			C IF ( A(I,J) .LT. 0. ) BUFFER(I) = '-'
41			C IF ( A(I,J) .GT. 0. ) BUFFER(I) = '+'
42			ENDDO
43			WRITE(0,'(I3,1X,1024A)') J,(BUFFER(I),I=1,Nx)
44			ENDDO
45			C
46			END