1 files changed, 118 insertions, 0 deletions
diff --git a/sys/gio/ncarutil/tests/conrect.f b/sys/gio/ncarutil/tests/conrect.f
new file mode 100644
index 00000000..401aad9b
--- /dev/null
+++ b/sys/gio/ncarutil/tests/conrect.f
@@ -0,0 +1,118 @@
+      SUBROUTINE TCONRE	(nplot, IERROR)
+C
+C LATEST REVISION	 JUNE 1984
+C
+C PURPOSE		 TO PROVIDE A SIMPLE DEMONSTRATION OF
+C			 CONREC	AND TO TEST CONREC ON A	SINGLE
+C			 PROBLEM
+C
+C USAGE			 CALL TCONRE (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			   CONREC 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 LANGUAGE		 FORTRAN
+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			 TCONRE	CALL SUBROUTINES EZCNTR, CONREC, AND
+C			 PWRIT TO DRAW TWO LABELLED CONTOUR PLOTS OF THE
+C			 ARRAY Z.
+C
+C PORTABILITY		 FORTRAN77
+C
+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
+C     DATA TX/.3955/, TY/.9765/
+	data tx/.4267/, ty/.97/
+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 NUMBER 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 CONREC
+C
+c +noao: flag added to plot either EZCNTR or CONREC
+      if (nplot .eq. 1) then
+      CALL WTSTR ( TX, TY,
+     1		 'DEMONSTRATION PLOT FOR EZCNTR ENTRY OF CONREC',2,0,0 )
+      CALL EZCNTR (Z,21,25)
+      endif
+c -noao
+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 CONREC
+C
+c +noao:  flag added to plot either EZCNTR of CONREC
+      if (nplot .eq. 2) then
+      CALL WTSTR ( TX ,TY,
+     1		 'DEMONSTRATION PLOT FOR CONREC ENTRY OF CONREC',2,0,0 )
+      CALL CONREC (Z,21,21,25,-4.5,4.5,.3,0,0,0)
+      endif
+c -noao
+c     CALL NEWFM
+C
+C     WRITE (6,1001)
+      RETURN
+C
+C1001 FORMAT ('	    CONREC TEST	SUCCESSFUL',24X,
+C    1	      'SEE PLOT	TO VERIFY PERFORMANCE')
+C
+      END