15
#include "mafdecls.fh"
15
17
#include "errquit.fh"
17
19
* **** local variables ****
20
character*255 full_filename
23
logical value,doit,oprint
25
character*255 full_filename,full_filename2
22
27
* **** external functions ****
23
28
logical control_Fei
24
29
external control_Fei
31
call Parallel_taskid(taskid)
32
oprint = (taskid.eq.MASTER)
34
if (.not.rtdb_get(rtdb,'nwpw:cpmd_properties',mt_log,1,doit))
39
write(luout,*) "================================"
40
write(luout,*) "==== CPMD Property Analysis ===="
41
write(luout,*) "================================"
44
write(luout,*) "CPMD property analysis is on. ",
45
> "To turn analysis off set cpmd_properties, i.e."
47
write(luout,*) " cpmd_properties off"
52
write(luout,*) "CPMD property analysis is off. ",
53
> "To turn analysis on set cpmd_properties, i.e."
55
write(luout,*) " cpmd_properties on"
26
63
if (.not.rtdb_cget(rtdb,'cpmd:ion_motion_filename',1,filename))
27
> call util_file_prefix('ion_motion',filename)
29
call util_file_name_noprefix(filename,.false.,
64
> call util_file_prefix('ion_motion',filename)
65
call util_file_name_noprefix(filename,.false.,
69
if (.not.rtdb_cget(rtdb,'cpmd:unfolded_xyz_filename',1,filename))
70
> call util_file_prefix('unfolded.xyz',filename)
71
call util_file_name_noprefix(filename,.false.,
34
75
call pair_dist(full_filename)
35
76
call pair_distab(full_filename)
36
77
call velocity_autocorrelation(full_filename)
78
call generate_unfolded_xyz(full_filename,full_filename2)
39
81
* **** perform shake force averaging if Fei set *****