57 |
|
#include "applications/staticProps/GofRAngle.hpp" |
58 |
|
#include "applications/staticProps/GofAngle2.hpp" |
59 |
|
#include "applications/staticProps/GofXyz.hpp" |
60 |
+ |
#include "applications/staticProps/TwoDGofR.hpp" |
61 |
|
#include "applications/staticProps/P2OrderParameter.hpp" |
62 |
|
#include "applications/staticProps/BondOrderParameter.hpp" |
63 |
|
#include "applications/staticProps/BOPofR.hpp" |
171 |
|
|
172 |
|
//parse md file and set up the system |
173 |
|
SimCreator creator; |
173 |
– |
std::cout << "dumpFile = " << dumpFileName << "\n"; |
174 |
|
SimInfo* info = creator.createSim(dumpFileName); |
175 |
|
|
176 |
|
RealType maxLen; |
216 |
|
painCave.isFatal = 1; |
217 |
|
simError(); |
218 |
|
} |
219 |
+ |
} else if (args_info.twodgofr_given){ |
220 |
+ |
if (args_info.dz_given) { |
221 |
+ |
analyser= new TwoDGofR(info, dumpFileName, sele1, sele2, maxLen, |
222 |
+ |
args_info.dz_arg, args_info.nbins_arg); |
223 |
+ |
} else { |
224 |
+ |
sprintf( painCave.errMsg, |
225 |
+ |
"A slab width (dz) must be specified when calculating TwoDGofR"); |
226 |
+ |
painCave.severity = OPENMD_ERROR; |
227 |
+ |
painCave.isFatal = 1; |
228 |
+ |
simError(); |
229 |
+ |
} |
230 |
|
} else if (args_info.p2_given) { |
231 |
|
analyser = new P2OrderParameter(info, dumpFileName, sele1, sele2); |
232 |
|
} else if (args_info.rp2_given){ |