--- /dev/null
+ SUBROUTINE ABRT
+ STOP 'IN ABRT'
+ END
+C*MODULE MTHLIB *DECK VCLR
+ SUBROUTINE VCLR(A,INCA,N)
+C
+ IMPLICIT DOUBLE PRECISION(A-H,O-Z)
+C
+ DIMENSION A(*)
+C
+ PARAMETER (ZERO=0.0D+00)
+C
+C ----- ZERO OUT VECTOR -A-, USING INCREMENT -INCA- -----
+C
+ IF (INCA .NE. 1) GO TO 200
+ DO 110 L=1,N
+ A(L) = ZERO
+ 110 CONTINUE
+ RETURN
+C
+ 200 CONTINUE
+ LA=1-INCA
+ DO 210 L=1,N
+ LA=LA+INCA
+ A(LA) = ZERO
+ 210 CONTINUE
+ RETURN
+ END