Initial commit

1 files changed, 114 insertions, 0 deletions

diff --git a/sys/gio/ncarutil/tests/conrcqckt.f b/sys/gio/ncarutil/tests/conrcqckt.f
new file mode 100644
index 00000000..d9d2f827
--- /dev/null
+++ b/sys/gio/ncarutil/tests/conrcqckt.f
@@ -0,0 +1,114 @@
+      SUBROUTINE TCNQCK	(IERROR)
+C
+C LATEST REVISION	 JUNE 1984
+C
+C PURPOSE		 TO PROVIDE A SIMPLE DEMONSTRATION OF
+C			 CONRECQCK AND TO TEST CONRECQCK ON A SINGLE
+C			 PROBLEM
+C
+C USAGE			 CALL TCNQCK (IERROR)
+C
+C ARGUMENTS
+C
+C ON OUTPUT		 IERROR
+C			   AN INTEGER VARIABLE
+C			   = 0,	IF THE TEST WAS	SUCCESSFUL,
+C			   = 1,	OTHERWISE
+C
+C I/O			 IF THE	TEST IS	SUCCESSFUL, THE	MESSAGE
+C
+C			   CONRECQCK TEST SUCCESSFUL  .	. . SEE	PLOT TO
+C			   VERIFY PERFORMANCE
+C
+C			 IS PRINTED ON UNIT 6.
+C			 IN ADDITION, TWO FRAMES CONTAINING THE	CONTOUR
+C			 PLOT ARE PRODUCED ON THE MACHINE GRAPHICS
+C			 DEVICE.  IN ORDER TO DETERMINE	IF THE TEST
+C			 WAS SUCCESSFUL, IT IS NECESSARY TO EXAMINE
+C			 THESE PLOTS.
+C
+C PRECISION		 SINGLE
+C
+C ALGORITHM		 THE FUNCTION
+C			   Z(X,Y) = X +	Y + 1./((X-.1)**2+Y**2+.09)
+C				    -1./((X+.1)**2+Y**2+.09)
+C			 FOR X = -1. TO	+1. IN INCREMENTS OF .1	AND
+C			     Y = -1.2 TO +1.2 IN INCREMENTS OF .1
+C			 IS COMPUTED.
+C			 TCNQCK	CALLS SUBROUTINES EZCNTR, CONREC, AND
+C			 PWRIT TO DRAW TWO LABELLED CONTOUR PLOTS OF THE
+C			 ARRAY Z.
+C
+C PORTABILITY		 ANSI FORTRAN77
+C
+C Z CONTAINS THE VALUES	TO BE PLOTTED.
+C
+      REAL	      Z(21,25)
+C
+C SPECIFY COORDINATES FOR PLOT TITLES.	ON AN ABSTRACT GRID WHERE
+C THE INTEGER COORDINATES RANGE	FROM 0.0 TO 1.0, THE VALUES TX AND TY
+C DEFINE THE CENTER OF THE TITLE STRING.
+C
+      DATA TX/.4267/, TY/.9765/
+C
+C
+C INITIALIZE ERROR PARAMETER
+C
+      IERROR = 0
+C
+C FILL TWO DIMENSIONAL ARRAY TO	BE PLOTTED
+C
+      DO  20 I=1,21
+	 X = .1*FLOAT(I-11)
+	 DO  10	J=1,25
+	    Y =	.1*FLOAT(J-13)
+	    Z(I,J) = X+Y+1./((X-.10)**2+Y**2+.09)-
+     1		     1./((X+.10)**2+Y**2+.09)
+   10	 CONTINUE
+   20 CONTINUE
+C
+C SELECT NORMALIZATION TRANSFORMATION 0
+C
+      CALL GSELNT (0)
+C
+C ENTRY	EZCNTR REQUIRES	ONLY THE ARRAY NAME AND	ITS DIMENSIONS
+C
+C THE TITLE FOR	THIS PLOT IS
+C
+C  DEMONSTRATION PLOT FOR EZCNTR ENTRY OF CONRECQCK
+C
+      CALL WTSTR (TX,TY,
+     1		 'DEMONSTRATION	PLOT FOR EZCNTR	ENTRY OF CONRECQCK',
+     2		 2,0,0)
+      CALL EZCNTR (Z,21,25)
+C
+C
+C ENTRY	CONREC ALLOWS USER SPECIFICATION OF PLOT PARAMETERS, IF	DESIRED
+C
+C IN THIS EXAMPLE, THE LOWEST CONTOUR LEVEL (-4.5), THE	HIGHEST	CONTOUR
+C LEVEL	(4.5), AND THE INCREMENT BETWEEN CONTOUR LEVELS	(0.3) ARE
+C SPECIFIED.
+C
+C THE TITLE FOR	THIS PLOT IS
+C
+C  DEMONSTRATION PLOT FOR CONREC ENTRY OF CONRECQCK
+C
+      CALL WTSTR (TX,TY,
+     1		 'DEMONSTRATION	PLOT FOR CONREC	ENTRY OF CONRECQCK',
+     2		 2,0,0)
+      CALL CONREC (Z,21,21,25,-4.5,4.5,.3,0,0,0)
+c     CALL NEWFM
+C
+c     WRITE (6,1001)
+      RETURN
+C
+c1001 FORMAT ('	 CONRECQCK TEST	SUCCESSFUL',24X,
+c    1	      'SEE PLOT	TO VERIFY PERFORMANCE')
+C
+C---------------------------------------------------------------------
+C REVISION HISTORY
+C
+C JUNE 1984			CONVERTED TO FORTRAN 77	AND GKS
+C
+C---------------------------------------------------------------------
+      END