1
block data initial_smd_cutoff_data
3
#include "smd_cutoff_data.fh"
7
data rcut_verlet /0.0d0/
12
subroutine smd_cutoff_get_rcut(arcut)
14
#include "smd_cutoff_data.fh"
16
double precision arcut
22
subroutine smd_cutoff_set_rcut(arcut)
24
#include "smd_cutoff_data.fh"
26
double precision arcut
32
subroutine smd_cutoff_get_rcut_verlet(arcut_verlet)
34
#include "smd_cutoff_data.fh"
36
double precision arcut_verlet
38
arcut_verlet = rcut_verlet
42
subroutine smd_cutoff_set_rcut_verlet(arcut_verlet)
44
#include "smd_cutoff_data.fh"
46
double precision arcut_verlet
48
rcut_verlet = arcut_verlet
52
subroutine smd_cutoff_print(un)
54
#include "smd_cutoff_data.fh"
58
write(un,*) "cutoff data"
59
write(un,*) "cutoff radius",rcut
60
write(un,*) "verlet cutoff radius",rcut_verlet
64
subroutine smd_cutoff_init()
67
#include "mafdecls.fh"
69
#include "smd_cutoff_data.fh"
76
pname = "smd_cutoff_init"
78
call smd_rtdb_get_handle(rtdb)
80
if (.not.rtdb_get(rtdb,tag,mt_dbl,1,rcut))
81
> call errquit(pname//'failed to get'//tag,0,
85
if (.not.rtdb_get(rtdb,tag,mt_dbl,1,rcut_verlet))
86
> call errquit(pname//'failed to get'//tag,0,