51 |
|
|
52 |
|
#include "applications/sequentialProps/SequentialPropsCmd.h" |
53 |
|
#include "applications/sequentialProps/SequentialAnalyzer.hpp" |
54 |
< |
#include "applications/sequentialProps/DensityAnalyzer.hpp" |
54 |
> |
#include "applications/sequentialProps/CenterOfMass.hpp" |
55 |
> |
#include "applications/sequentialProps/ContactAngle1.hpp" |
56 |
|
|
57 |
|
using namespace OpenMD; |
58 |
|
|
106 |
|
SimInfo* info = creator.createSim(dumpFileName, false); |
107 |
|
|
108 |
|
SequentialAnalyzer* analyzer; |
109 |
< |
if(args_info.density_given){ |
110 |
< |
analyzer = new DensityAnalyzer(info, dumpFileName, sele1); |
109 |
> |
if(args_info.com_given){ |
110 |
> |
analyzer = new CenterOfMass(info, dumpFileName, sele1); |
111 |
> |
} else if(args_info.ca1_given){ |
112 |
> |
RealType solidZ; |
113 |
> |
if (args_info.referenceZ_given) |
114 |
> |
solidZ = args_info.referenceZ_arg; |
115 |
> |
else { |
116 |
> |
sprintf( painCave.errMsg, |
117 |
> |
"--referenceZ must be set if --ca1 is used\n"); |
118 |
> |
painCave.severity = OPENMD_ERROR; |
119 |
> |
painCave.isFatal = 1; |
120 |
> |
simError(); |
121 |
> |
} |
122 |
> |
RealType dropletR; |
123 |
> |
if (args_info.dropletR_given) |
124 |
> |
dropletR = args_info.dropletR_arg; |
125 |
> |
else { |
126 |
> |
sprintf( painCave.errMsg, |
127 |
> |
"--dropletR must be set if --ca1 is used\n"); |
128 |
> |
painCave.severity = OPENMD_ERROR; |
129 |
> |
painCave.isFatal = 1; |
130 |
> |
simError(); |
131 |
> |
} |
132 |
> |
|
133 |
> |
analyzer = new ContactAngle1(info, dumpFileName, sele1, solidZ, dropletR); |
134 |
|
} |
135 |
|
|
136 |
|
if (args_info.output_given) { |