Merge branch 'lipid' of mmka.chem.univ.gda.pl:unres into lipid

[unres.git] / source / unres / src_Eshel / SRC-SURPLUS / add.f

diff --git a/source/unres/src_Eshel/SRC-SURPLUS/add.f b/source/unres/src_Eshel/SRC-SURPLUS/add.f
new file mode 100644 (file)
index 0000000..fd91a70
--- /dev/null
+++ b/source/unres/src_Eshel/SRC-SURPLUS/add.f
@@ -0,0 +1,28 @@
+      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