Commit dd3cbdf0 authored by dceresoli's avatar dceresoli

london_rvdw(1:nsx) to enter vdw radii


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/[email protected] c92efa57-630b-4861-b058-cf58834340f0
parent 167a87c8
......@@ -235,6 +235,7 @@ SUBROUTINE iosys()
B_field, fixed_magnetization, report, lspinorb,&
starting_spin_angle, assume_isolated,spline_ps,&
vdw_corr, london, london_s6, london_rcut, london_c6, &
london_rvdw, &
ts_vdw, ts_vdw_isolated, ts_vdw_econv_thr, &
xdm, xdm_a1, xdm_a2, lforcet, &
one_atom_occupations, &
......@@ -283,7 +284,7 @@ SUBROUTINE iosys()
!
USE constraints_module, ONLY : init_constraint
USE read_namelists_module, ONLY : read_namelists, sm_not_set
USE london_module, ONLY : init_london, lon_rcut, scal6, in_c6
USE london_module, ONLY : init_london, lon_rcut, scal6, in_c6, in_rvdw
USE xdm_module, ONLY : init_xdm, a1i, a2i
USE tsvdw_module, ONLY : vdw_isolated, vdw_econv_thr
USE us, ONLY : spline_ps_ => spline_ps
......@@ -1218,6 +1219,7 @@ SUBROUTINE iosys()
lon_rcut = london_rcut
scal6 = london_s6
in_c6(:) = london_c6(:)
in_rvdw(:) = london_rvdw(:)
END IF
IF ( lxdm ) THEN
a1i = xdm_a1
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment