78 |
|
#include "applications/staticProps/TetrahedralityParam.hpp" |
79 |
|
#include "applications/staticProps/TetrahedralityParamZ.hpp" |
80 |
|
#include "applications/staticProps/RNEMDStats.hpp" |
81 |
+ |
#include "applications/staticProps/NitrileFrequencyMap.hpp" |
82 |
+ |
#include "applications/staticProps/MultipoleSum.hpp" |
83 |
|
|
84 |
|
using namespace OpenMD; |
85 |
|
|
240 |
|
} else { |
241 |
|
sprintf( painCave.errMsg, |
242 |
|
"A cutoff radius (rcut) must be specified when calculating Bond Order Parameters"); |
243 |
+ |
painCave.severity = OPENMD_ERROR; |
244 |
+ |
painCave.isFatal = 1; |
245 |
+ |
simError(); |
246 |
+ |
} |
247 |
+ |
} else if (args_info.multipole_given){ |
248 |
+ |
if (args_info.rcut_given) { |
249 |
+ |
analyser = new MultipoleSum(info, dumpFileName, sele1, |
250 |
+ |
args_info.rcut_arg); |
251 |
+ |
} else { |
252 |
+ |
sprintf( painCave.errMsg, |
253 |
+ |
"A cutoff radius (rcut) must be specified when calculating Multipole Sums"); |
254 |
|
painCave.severity = OPENMD_ERROR; |
255 |
|
painCave.isFatal = 1; |
256 |
|
simError(); |
280 |
|
painCave.isFatal = 1; |
281 |
|
simError(); |
282 |
|
} |
283 |
< |
} else if (args_info.bor_given){ |
283 |
> |
} else if (args_info.ior_given){ |
284 |
|
if (args_info.rcut_given) { |
285 |
< |
analyser = new BOPofR(info, dumpFileName, sele1, args_info.rcut_arg, |
285 |
> |
analyser = new IcosahedralOfR(info, dumpFileName, sele1, |
286 |
> |
args_info.rcut_arg, |
287 |
> |
args_info.nbins_arg, maxLen); |
288 |
> |
} else { |
289 |
> |
sprintf( painCave.errMsg, |
290 |
> |
"A cutoff radius (rcut) must be specified when calculating Bond Order Parameters"); |
291 |
> |
painCave.severity = OPENMD_ERROR; |
292 |
> |
painCave.isFatal = 1; |
293 |
> |
simError(); |
294 |
> |
} |
295 |
> |
} else if (args_info.for_given){ |
296 |
> |
if (args_info.rcut_given) { |
297 |
> |
analyser = new FCCOfR(info, dumpFileName, sele1, args_info.rcut_arg, |
298 |
|
args_info.nbins_arg, maxLen); |
299 |
|
} else { |
300 |
|
sprintf( painCave.errMsg, |
339 |
|
} else if (args_info.rnemdrt_given) { |
340 |
|
analyser = new RNEMDRTheta(info, dumpFileName, sele1, |
341 |
|
args_info.nbins_arg, args_info.nanglebins_arg); |
342 |
+ |
} else if (args_info.nitrile_given) { |
343 |
+ |
analyser = new NitrileFrequencyMap(info, dumpFileName, sele1, |
344 |
+ |
args_info.nbins_arg); |
345 |
|
} else if (args_info.p_angle_given) { |
346 |
|
if (args_info.sele1_given) { |
347 |
|
if (args_info.sele2_given) |