source/unres/src_MD/add.f

   1       SUBROUTINE ABRT
   2       STOP 'IN ABRT'
   3       END
   4 C*MODULE MTHLIB  *DECK VCLR
   5       SUBROUTINE VCLR(A,INCA,N)
   6 C
   7       IMPLICIT DOUBLE PRECISION(A-H,O-Z)
   8 C
   9       DIMENSION A(*)
  10 C
  11       PARAMETER (ZERO=0.0D+00)
  12 C
  13 C     ----- ZERO OUT VECTOR -A-, USING INCREMENT -INCA- -----
  14 C
  15       IF (INCA .NE. 1) GO TO 200
  16       DO 110 L=1,N
  17          A(L) = ZERO
  18   110 CONTINUE
  19       RETURN
  20 C
  21   200 CONTINUE
  22       LA=1-INCA
  23       DO 210 L=1,N
  24          LA=LA+INCA
  25          A(LA) = ZERO
  26   210 CONTINUE
  27       RETURN
  28       END