add rigid 4-site water treatment to OPTIMIZE & CRITICAL

Author: jayponder

example/alatet.all

@@ -5,9 +5,7 @@
 ###################################################################
 
 parameters              NONE
-
-group-molecule
-polar-eps               0.000001
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              ARITHMETIC
@@ -16,6 +14,7 @@ radiussize              DIAMETER
 epsilonrule             GEOMETRIC
 dielectric              1.0
 polarization            MUTUAL
+polar-eps               0.000001
 mutual-11-scale         0.0
 
 ## Dang-Chang Parameters (J Chem Phys, 106, 8149-8159, 1997)

example/alatet.dvol

@@ -8,7 +8,7 @@ parameters              NONE
 verbose
 #debug
 
-group-molecule
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              GEOMETRIC

example/dang.key

@@ -5,8 +5,7 @@
 ################################################
 
 parameters              NONE
-
-group-molecule
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              ARITHMETIC

example/nada.key

@@ -5,8 +5,7 @@
 ################################################################
 
 parameters              NONE
-
-group-molecule
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              ARITHMETIC
@@ -19,6 +18,8 @@ polarization            MUTUAL
 polar-eps               0.000001
 mutual-11-scale         0.0
 
+## POL3 Parameters (J Chem Phys, 99, 6208-6219, 1995)
+
 atom             1    O    "POL3 Water O"     8     16.000     2
 atom             2    H    "POL3 Water H"     1      1.008     1
 

example/opc.key

@@ -5,8 +5,7 @@
 #########################################################
 
 parameters              NONE
-
-group-molecule
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              GEOMETRIC

example/pol3.key

@@ -5,8 +5,7 @@
 ###########################################################
 
 parameters              NONE
-
-group-molecule
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              GEOMETRIC

example/spc.key

@@ -5,8 +5,7 @@
 ##############################################################
 
 parameters              NONE
-
-group-molecule
+rattle                  WATER
 
 vdwtype                 LENNARD-JONES
 radiusrule              GEOMETRIC

example/tip3p.key

@@ -1699,7 +1699,7 @@ connolly.o: faces.o inform.o iounit.o math.o
 control.o: argue.o inform.o keys.o output.o
 correlate.o: ascii.o atomid.o atoms.o files.o inform.o iounit.o
 couple.o: sizes.o
-critical.o: atoms.o files.o inform.o iounit.o keys.o minima.o output.o sizes.o usage.o
+critical.o: atoms.o files.o freeze.o inform.o iounit.o keys.o minima.o output.o sizes.o usage.o
 crystal.o: atomid.o atoms.o bound.o boxes.o couple.o files.o iounit.o math.o molcul.o
 cspline.o: iounit.o
 ctrpot.o:
@@ -2052,7 +2052,7 @@ opbend.o:
 opdist.o:
 openend.o:
 openmp.o:
-optimize.o: atoms.o bound.o files.o inform.o iounit.o scales.o usage.o
+optimize.o: atoms.o bound.o files.o freeze.o inform.o iounit.o scales.o usage.o
 optinit.o: bound.o inform.o keys.o output.o potent.o
 optirot.o: files.o inform.o iounit.o math.o omega.o scales.o zcoord.o
 optrigid.o: files.o group.o inform.o iounit.o math.o output.o rigid.o