Initial commit

1 files changed, 101 insertions, 0 deletions

diff --git a/sys/gio/ncarutil/tests/strmlnt.f b/sys/gio/ncarutil/tests/strmlnt.f
new file mode 100644
index 00000000..f2b40c69
--- /dev/null
+++ b/sys/gio/ncarutil/tests/strmlnt.f
@@ -0,0 +1,101 @@
+      SUBROUTINE TSTRML	(IERROR)
+C
+C LATEST REVISION	 JUNE 1984
+C
+C PURPOSE		 TO PROVIDE A SIMPLE DEMONSTRATION OF
+C			 ROUTINE STRMLN.
+C
+C USAGE			 CALL TSTRML (IERROR)
+C
+C ARGUMENTS
+C
+C ON OUTPUT		 IERROR
+C			   AN INTEGER VARIABLE
+C			   =0 IF THERE IS A NORMAL EXIT	FROM THE
+C			      ROUTINE STRMLN.
+C			   =1 OTHERWISE
+C
+C I/O			 IF THERE IS A NORMAL EXIT FROM	THE ROUTINE
+C			 STRMLN	THE MESSAGE
+C			   STRMLN TEST SUCCESSFUL . . .	SEE PLOT TO
+C			   VERIFY PERFORMANCE
+C			 IS PRINTED.
+C
+C PRECISION		 SINGLE
+C
+C
+C LANGUAGE		 FORTRAN
+C
+C ALGORITHM		 ROUTINE TSTRML	CALLS ROUTINE STRMLN TO
+C			 PRODUCE A PLOT	REPRESENTING THE  FLOW AND
+C			 MAGNITUDE OF A	VECTOR FIELD.
+C
+C PORTABILITY		 FORTRAN77
+C
+C
+C
+      REAL	      U(21,25)	 ,V(21,25)   ,WRK(1050)
+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/.5/,TY/.9765/
+C
+C     SET DIMENSIONS
+C
+      DATA NH,NV/21,25/
+C
+C INITIALIZE ERROR PARAMETER
+C
+      IERROR = 1
+C
+C SPECIFY HORIZONTAL AND VERTICAL VECTOR COMPONENTS U AND V ON
+C THE RECTANGULAR GRID
+C
+      TPIMX = 2.*3.14/FLOAT(NH)
+      TPJMX = 2.*3.14/FLOAT(NV)
+      DO  20 J=1,NV
+	 DO  10	I=1,NH
+	 U(I,J)	= SIN(TPIMX*(FLOAT(I)-1.))
+	 V(I,J)	= SIN(TPJMX*(FLOAT(J)-1.))
+   10	 CONTINUE
+   20 CONTINUE
+C
+C SELECT NORMALIZATION TRANSFORMATION 0
+C
+      CALL GSELNT (0)
+C
+C CALL WTSTR FOR STRMLN	PLOT TITLE
+C
+      CALL WTSTR (TX,TY,'DEMONSTRATION PLOT FOR ROUTINE STRMLN',2,
+     1		  0,0)
+C
+C DEFINE NORMALIZATION TRANSFORMATION 1, AND SET UP LOG	SCALING
+C
+      CALL GSVP	( 1, 0.1, 0.9, 0.1, 0.9	)
+      CALL GSWN	( 1, 1.0, 21., 1.0, 25.	)
+      CALL SETUSV ( 'LS' , 1 )
+C
+C SELECT NORMALIZATION TRANSFORMATION 1
+C
+      CALL GSELNT (1)
+C
+C     DRAW PERIMETER
+C
+c     CALL PERIM(1,0,1,0)
+C
+C CALL STRMLN FOR VECTOR FIELD STREAMLINES PLOT
+C
+      CALL STRMLN (U,V,WRK,NH,NH,NV,0,IER)
+C
+c     CALL NEWFM
+C
+      IERROR = 0
+c     WRITE (6,1001)
+      RETURN
+C
+c1001 FORMAT ('	    STRMLN TEST	SUCCESSFUL',24X,
+c    1	      'SEE PLOT	TO VERIFY PERFORMANCE')
+C
+      END