2003-01-20 05:58:50 +08:00
|
|
|
!
|
|
|
|
|
! Copyright (C) 2002 CP90 group
|
|
|
|
|
! This file is distributed under the terms of the
|
|
|
|
|
! GNU General Public License. See the file `License'
|
|
|
|
|
! in the root directory of the present distribution,
|
|
|
|
|
! or http://www.gnu.org/copyleft/gpl.txt .
|
|
|
|
|
!
|
All namelists and cards moved to Modules/input_parameters.f90 .
From now on, all new input variables should be added
to this module, and then copied to the code internal
variables in the input.f90 subroutine
The namelists and cards parsers are in :
Modules/read_namelists.f90 and Modules/read_cards.f90
files input_parameters.f90 read_namelists.f90 read_cards.f90
are shared by all codes, while each code has its own version
of input.f90 ( used to copy input values into internals variables ).
EXAMPLE:
suppose you need to add a new input variable called "pippo"
to the namelist control, then:
1) add pippo to the input_parameters.f90 file containing the
namelist control
INTEGER :: pippo = 0
NAMELIST / control / ....., pippo
remember: always set an initialization value!
2) add pippo to the control_default subroutine
( cantained in module read_namelists.f90 )
subroutine control_default( prog )
...
IF( prog == 'PW' ) pippo = 10
...
end subroutine
this routine set the default value for pippo,
that could vary with the code
3) add pippo to the control_bcast subroutine
( cantained in module read_namelists.f90 )
subroutine control_bcast( )
...
call mp_bcast( pippo )
...
end subroutine
4) add pippo to the control_checkin subroutine
( cantained in module read_namelists.f90 )
subroutine control_checking( prog )
...
IF( pippo < 0 ) &
CALL error(' control_checkin ',' variable pippo less than 0 ', 1 )
...
end subroutine
5) Copy the value of pippo in the code internal variables
( file input.f90 )
subroutine iosys()
use input_parameters, only: ...., pippo
use pwcom, only: ....., myvar
...
call read_namelists( 'PW' )
...
myvar = pippo
...
end subroutine
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@282 c92efa57-630b-4861-b058-cf58834340f0
2003-07-31 21:24:20 +08:00
|
|
|
MODULE input_cp
|
2003-01-20 05:58:50 +08:00
|
|
|
|
All namelists and cards moved to Modules/input_parameters.f90 .
From now on, all new input variables should be added
to this module, and then copied to the code internal
variables in the input.f90 subroutine
The namelists and cards parsers are in :
Modules/read_namelists.f90 and Modules/read_cards.f90
files input_parameters.f90 read_namelists.f90 read_cards.f90
are shared by all codes, while each code has its own version
of input.f90 ( used to copy input values into internals variables ).
EXAMPLE:
suppose you need to add a new input variable called "pippo"
to the namelist control, then:
1) add pippo to the input_parameters.f90 file containing the
namelist control
INTEGER :: pippo = 0
NAMELIST / control / ....., pippo
remember: always set an initialization value!
2) add pippo to the control_default subroutine
( cantained in module read_namelists.f90 )
subroutine control_default( prog )
...
IF( prog == 'PW' ) pippo = 10
...
end subroutine
this routine set the default value for pippo,
that could vary with the code
3) add pippo to the control_bcast subroutine
( cantained in module read_namelists.f90 )
subroutine control_bcast( )
...
call mp_bcast( pippo )
...
end subroutine
4) add pippo to the control_checkin subroutine
( cantained in module read_namelists.f90 )
subroutine control_checking( prog )
...
IF( pippo < 0 ) &
CALL error(' control_checkin ',' variable pippo less than 0 ', 1 )
...
end subroutine
5) Copy the value of pippo in the code internal variables
( file input.f90 )
subroutine iosys()
use input_parameters, only: ...., pippo
use pwcom, only: ....., myvar
...
call read_namelists( 'PW' )
...
myvar = pippo
...
end subroutine
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@282 c92efa57-630b-4861-b058-cf58834340f0
2003-07-31 21:24:20 +08:00
|
|
|
IMPLICIT NONE
|
|
|
|
|
SAVE
|
2003-01-20 05:58:50 +08:00
|
|
|
|
All namelists and cards moved to Modules/input_parameters.f90 .
From now on, all new input variables should be added
to this module, and then copied to the code internal
variables in the input.f90 subroutine
The namelists and cards parsers are in :
Modules/read_namelists.f90 and Modules/read_cards.f90
files input_parameters.f90 read_namelists.f90 read_cards.f90
are shared by all codes, while each code has its own version
of input.f90 ( used to copy input values into internals variables ).
EXAMPLE:
suppose you need to add a new input variable called "pippo"
to the namelist control, then:
1) add pippo to the input_parameters.f90 file containing the
namelist control
INTEGER :: pippo = 0
NAMELIST / control / ....., pippo
remember: always set an initialization value!
2) add pippo to the control_default subroutine
( cantained in module read_namelists.f90 )
subroutine control_default( prog )
...
IF( prog == 'PW' ) pippo = 10
...
end subroutine
this routine set the default value for pippo,
that could vary with the code
3) add pippo to the control_bcast subroutine
( cantained in module read_namelists.f90 )
subroutine control_bcast( )
...
call mp_bcast( pippo )
...
end subroutine
4) add pippo to the control_checkin subroutine
( cantained in module read_namelists.f90 )
subroutine control_checking( prog )
...
IF( pippo < 0 ) &
CALL error(' control_checkin ',' variable pippo less than 0 ', 1 )
...
end subroutine
5) Copy the value of pippo in the code internal variables
( file input.f90 )
subroutine iosys()
use input_parameters, only: ...., pippo
use pwcom, only: ....., myvar
...
call read_namelists( 'PW' )
...
myvar = pippo
...
end subroutine
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@282 c92efa57-630b-4861-b058-cf58834340f0
2003-07-31 21:24:20 +08:00
|
|
|
PRIVATE
|
2003-01-20 05:58:50 +08:00
|
|
|
|
2003-08-25 06:25:53 +08:00
|
|
|
PUBLIC :: iosys
|
2004-02-17 17:53:14 +08:00
|
|
|
PUBLIC :: iosys_pseudo
|
2004-02-15 18:29:15 +08:00
|
|
|
PUBLIC :: read_input_file
|
2003-01-20 05:58:50 +08:00
|
|
|
|
All namelists and cards moved to Modules/input_parameters.f90 .
From now on, all new input variables should be added
to this module, and then copied to the code internal
variables in the input.f90 subroutine
The namelists and cards parsers are in :
Modules/read_namelists.f90 and Modules/read_cards.f90
files input_parameters.f90 read_namelists.f90 read_cards.f90
are shared by all codes, while each code has its own version
of input.f90 ( used to copy input values into internals variables ).
EXAMPLE:
suppose you need to add a new input variable called "pippo"
to the namelist control, then:
1) add pippo to the input_parameters.f90 file containing the
namelist control
INTEGER :: pippo = 0
NAMELIST / control / ....., pippo
remember: always set an initialization value!
2) add pippo to the control_default subroutine
( cantained in module read_namelists.f90 )
subroutine control_default( prog )
...
IF( prog == 'PW' ) pippo = 10
...
end subroutine
this routine set the default value for pippo,
that could vary with the code
3) add pippo to the control_bcast subroutine
( cantained in module read_namelists.f90 )
subroutine control_bcast( )
...
call mp_bcast( pippo )
...
end subroutine
4) add pippo to the control_checkin subroutine
( cantained in module read_namelists.f90 )
subroutine control_checking( prog )
...
IF( pippo < 0 ) &
CALL error(' control_checkin ',' variable pippo less than 0 ', 1 )
...
end subroutine
5) Copy the value of pippo in the code internal variables
( file input.f90 )
subroutine iosys()
use input_parameters, only: ...., pippo
use pwcom, only: ....., myvar
...
call read_namelists( 'PW' )
...
myvar = pippo
...
end subroutine
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@282 c92efa57-630b-4861-b058-cf58834340f0
2003-07-31 21:24:20 +08:00
|
|
|
CONTAINS
|
2003-01-20 05:58:50 +08:00
|
|
|
|
2004-02-15 18:29:15 +08:00
|
|
|
! subroutines
|
|
|
|
|
! ----------------------------------------------
|
|
|
|
|
! ----------------------------------------------
|
|
|
|
|
|
2004-07-06 23:21:25 +08:00
|
|
|
SUBROUTINE read_input_file( lneb, lsmd, lwf )
|
2004-02-15 18:29:15 +08:00
|
|
|
USE read_namelists_module, ONLY: read_namelists
|
|
|
|
|
USE read_cards_module, ONLY: read_cards
|
|
|
|
|
USE input_parameters, ONLY: calculation
|
|
|
|
|
IMPLICIT NONE
|
|
|
|
|
|
2004-07-06 23:21:25 +08:00
|
|
|
LOGICAL, INTENT(OUT) :: lneb, lsmd, lwf
|
2004-02-15 18:29:15 +08:00
|
|
|
|
|
|
|
|
! . Read NAMELISTS ..................................................!
|
|
|
|
|
|
|
|
|
|
CALL read_namelists( 'CP' )
|
|
|
|
|
|
|
|
|
|
! . Read CARDS ......................................................!
|
|
|
|
|
|
|
|
|
|
CALL read_cards( 'CP' )
|
|
|
|
|
|
|
|
|
|
lneb = ( TRIM( calculation ) == 'neb' )
|
2004-06-30 02:12:51 +08:00
|
|
|
lsmd = ( TRIM( calculation ) == 'smd' )
|
2004-07-06 23:21:25 +08:00
|
|
|
lwf = ( TRIM( calculation ) == 'cp-wf' )
|
2004-02-15 18:29:15 +08:00
|
|
|
|
|
|
|
|
RETURN
|
|
|
|
|
END SUBROUTINE
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2004-06-01 01:52:55 +08:00
|
|
|
subroutine iosys_pseudo( psfile_ , pseudo_dir_ , nsp_ )
|
|
|
|
|
use input_parameters, only: atom_pfile, pseudo_dir, ntyp
|
2004-02-17 17:53:14 +08:00
|
|
|
use parameters, only: nsx
|
|
|
|
|
implicit none
|
2004-10-26 17:32:48 +08:00
|
|
|
character(len=256) :: psfile_ ( nsx ) , pseudo_dir_
|
2004-06-01 01:52:55 +08:00
|
|
|
integer :: nsp_
|
2004-02-17 17:53:14 +08:00
|
|
|
nsp_ = ntyp
|
|
|
|
|
psfile_= ' '
|
|
|
|
|
psfile_ ( 1:nsp_ ) = atom_pfile( 1:nsp_ )
|
|
|
|
|
pseudo_dir_ = pseudo_dir
|
2004-02-19 01:00:32 +08:00
|
|
|
!
|
|
|
|
|
! read in pseudopotentials and wavefunctions files
|
|
|
|
|
!
|
|
|
|
|
call readpp()
|
2004-02-17 17:53:14 +08:00
|
|
|
return
|
|
|
|
|
end subroutine
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
!-----------------------------------------------------------------------
|
2004-08-27 18:20:42 +08:00
|
|
|
subroutine iosys( nbeg_ , ndr_ , ndw_ , nomore_ , iprint_ , isave_ &
|
2004-09-13 16:49:06 +08:00
|
|
|
& , delt_ , emass_ , emaec_ , tsde_ , frice_ , grease_ &
|
2004-06-30 02:12:51 +08:00
|
|
|
& , tortho_ , eps_ , max_ , trane_ , ampre_ , tranp_ , amprp_ &
|
|
|
|
|
& , tfor_ , tsdp_ , fricp_ , greasp_ , tcp_ , tcap_ , tolp_ , trhor_ , trhow_ , tvlocw_ &
|
|
|
|
|
& , tnosep_ , qnp_ , tempw_ , tnosee_ , qne_ , ekincw_ &
|
2004-08-27 18:20:42 +08:00
|
|
|
& , tpre_ , thdyn_ , thdiag_ , iforceh_ , wmass_ , frich_ , greash_ , press_ &
|
2004-10-14 16:38:05 +08:00
|
|
|
& , tnoseh_ , qnh_ , temph_ , celldm_ , ibrav_ , tau0_ , iforce_ &
|
2004-06-30 02:12:51 +08:00
|
|
|
& , nat_ , nsp_ , na_ , pmass_ , rcmax_ , f_ , nel_ , nspin_ , nupdwn_ &
|
|
|
|
|
& , iupdwn_ , n_ , nx_, nr1_ , nr2_ , nr3_ , omega_ , alat_ , a1_ , a2_ , a3_ &
|
2004-10-14 16:38:05 +08:00
|
|
|
& , nr1b_ , nr2b_ , nr3b_ , nr1s_ , nr2s_ , nr3s_ &
|
2004-06-30 16:49:41 +08:00
|
|
|
& , psfile_ , pseudo_dir_, iprsta_, ispin_ &
|
2004-06-30 02:12:51 +08:00
|
|
|
& , sm_p, smcp_, smlm_, smopt_, linr_, polm_, kwnp_, codfreq_, forfreq_, smwfreq_ &
|
|
|
|
|
& , tol_, lmfreq_, maxlm_ )
|
|
|
|
|
|
|
|
|
|
!-----------------------------------------------------------------------
|
|
|
|
|
! this subroutine reads control variables from standard input (unit 5)
|
|
|
|
|
! ------------------------------------------------------------------
|
|
|
|
|
|
|
|
|
|
use input_parameters, only: &
|
|
|
|
|
nr1, nr2, nr3, greash, press, nr2s, nr3s, nr1s, tolp, temph, grease, &
|
2004-09-13 16:49:06 +08:00
|
|
|
tempw, fnoseh, amprp, greasp, tranp, atomic_positions, nelec, &
|
2004-06-30 02:12:51 +08:00
|
|
|
if_pos, rd_ht, nelup, neldw, occupations, f_inp, pos, nr3b, pseudo_dir, &
|
2004-08-27 18:20:42 +08:00
|
|
|
nr1b, nr2b, sp_pos, atom_mass, atom_pfile, iprint, isave, orthogonalization, &
|
2004-06-30 02:12:51 +08:00
|
|
|
electron_velocities, startingwfc, ndr, ndw, ion_dynamics, ion_damping, &
|
|
|
|
|
cell_velocities, electron_dynamics, electron_damping, ion_velocities, &
|
|
|
|
|
celldm, nbnd, nspin, calculation, ntyp, ibrav, restart_mode, ion_positions, &
|
|
|
|
|
nstep, ecutwfc, ecutrho, ampre, ortho_eps, ortho_max, wmass, qcutz, q2sigma, &
|
|
|
|
|
ecfixed, ekincw, fnosep, nat, tstress, disk_io, fnosee, ion_temperature, &
|
|
|
|
|
cell_temperature, cell_dofree, cell_dynamics, cell_damping, electron_temperature, &
|
2004-09-23 22:11:33 +08:00
|
|
|
dt, emass, emass_cutoff, ion_radius, verbosity, tprnfor, &
|
2004-06-30 02:12:51 +08:00
|
|
|
ekin_conv_thr, etot_conv_thr, max_seconds, na_inp, rd_pos, atom_label, rd_vel, &
|
2004-06-30 16:49:41 +08:00
|
|
|
smd_polm, smd_kwnp, smd_linr, smd_stcd, smd_stcd1, smd_stcd2, smd_stcd3, smd_codf, &
|
2004-06-30 02:12:51 +08:00
|
|
|
smd_forf, smd_smwf, smd_lmfreq, smd_tol, smd_maxlm, smd_smcp, smd_smopt, smd_smlm, &
|
2004-09-23 22:11:33 +08:00
|
|
|
num_of_images, smd_ene_ini, smd_ene_fin
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
use read_namelists_module, only: read_namelists
|
|
|
|
|
use read_cards_module, only: read_cards
|
|
|
|
|
|
2004-08-27 18:20:42 +08:00
|
|
|
use constants, only: pi, scmass, factem, eps8, uma_au
|
2004-06-30 02:12:51 +08:00
|
|
|
use parameters, only: nsx, natx, nbndxx
|
|
|
|
|
use io_global, only: ionode, stdout
|
|
|
|
|
use control_flags, only: taurdr, tprnfor_ => tprnfor
|
2004-09-13 16:49:06 +08:00
|
|
|
use control_flags, only: tzerop, tzeroe, tzeroc
|
2004-06-30 02:12:51 +08:00
|
|
|
use mp, only: mp_bcast
|
|
|
|
|
USE control_flags, ONLY: tconvthrs, lneb, lsmd
|
|
|
|
|
USE check_stop, ONLY: check_stop_init
|
2004-08-27 18:20:42 +08:00
|
|
|
USE ions_base, ONLY: ions_base_init
|
2004-10-14 16:38:05 +08:00
|
|
|
|
|
|
|
|
use gvecw, only: agg => ecutz, sgg => ecsig, e0gg => ecfix
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
!
|
|
|
|
|
implicit none
|
|
|
|
|
!
|
|
|
|
|
!
|
|
|
|
|
real(kind=8) :: ampre_ , delt_ , ekincw_ , emass_ , emaec_ , eps_ , &
|
|
|
|
|
& frice_ , fricp_ , frich_ , grease_ , greasp_ , greash_ , &
|
|
|
|
|
& press_ , qnp_ , qne_ , qnh_ , tempw_ , temph_ , tolp_ , wmass_ , &
|
2004-10-14 16:38:05 +08:00
|
|
|
amprp_ ( nsx ), celldm_ ( 6 ), tau0_ ( 3, natx )
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-08-27 18:20:42 +08:00
|
|
|
integer :: nbeg_ , ndr_ , ndw_ , nomore_ , iprint_ , max_ , iforce_( 3, natx )
|
|
|
|
|
integer :: isave_
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-09-13 16:49:06 +08:00
|
|
|
logical :: trane_ , tsde_ , tortho_ , tnosee_ , tfor_ , tsdp_ , tcp_ , &
|
2004-06-30 02:12:51 +08:00
|
|
|
tcap_ , tnosep_ , trhor_ , trhow_ , tvlocw_ , tpre_ , thdyn_ , thdiag_ , &
|
2004-08-27 18:20:42 +08:00
|
|
|
tnoseh_ , tranp_ ( nsx )
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
integer :: nat_ , nsp_ , na_ ( nsx ), nel_ ( 2 ), nspin_ , &
|
|
|
|
|
& nupdwn_ ( 2 ), iupdwn_ ( 2 ), n_ , nx_ , nr1_ , nr2_ , nr3_ , &
|
2004-08-27 18:20:42 +08:00
|
|
|
& nr1b_ , nr2b_ , nr3b_ , nr1s_ , nr2s_ , nr3s_ , ibrav_, iprsta_, &
|
|
|
|
|
iforceh_( 3, 3 )
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
real(kind=8) :: pmass_ ( nsx ), rcmax_ ( nsx ), f_ ( nbndxx ), ispin_ ( nbndxx ), &
|
2004-10-14 16:38:05 +08:00
|
|
|
& omega_ , alat_ , a1_ ( 3 ), a2_ ( 3 ), a3_ ( 3 )
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
|
2004-10-26 17:32:48 +08:00
|
|
|
character(len=256) :: psfile_ ( nsx ) , pseudo_dir_
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! local variables
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
real(kind=8), parameter:: terahertz = 2.418D-5
|
2004-08-27 18:20:42 +08:00
|
|
|
real(kind=8) :: taus( 3, natx ), ocp, fsum
|
2004-06-30 02:12:51 +08:00
|
|
|
integer :: unit = 5, ionode_id = 0, i, ia, ios, is, iss, in, isa, smpm
|
|
|
|
|
|
|
|
|
|
integer, intent(out) :: sm_p, kwnp_, codfreq_, forfreq_, smwfreq_, &
|
|
|
|
|
& lmfreq_, maxlm_
|
|
|
|
|
logical, intent(out) :: smcp_, smlm_, smopt_, linr_, polm_
|
|
|
|
|
real(kind=8) :: tol_
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
IF( TRIM( calculation ) == 'nscf' ) trhor_ = .true.
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! translate from input to internals of SMCP,
|
2004-07-23 20:29:51 +08:00
|
|
|
!
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
! ... SM_P
|
|
|
|
|
|
|
|
|
|
sm_p = num_of_images -1
|
|
|
|
|
smpm = sm_p -1
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! ... what to do
|
|
|
|
|
|
|
|
|
|
smcp_ = smd_smcp
|
|
|
|
|
smopt_ = smd_smopt
|
|
|
|
|
smlm_ = smd_smlm
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! ... initial path info
|
|
|
|
|
|
|
|
|
|
linr_ = smd_linr
|
|
|
|
|
polm_ = smd_polm
|
|
|
|
|
kwnp_ = smd_kwnp
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! ... Frequencey of wiriting
|
|
|
|
|
|
|
|
|
|
codfreq_ = smd_codf
|
|
|
|
|
forfreq_ = smd_forf
|
|
|
|
|
smwfreq_ = smd_smwf
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! ... Lagrange multiplier info.
|
|
|
|
|
|
|
|
|
|
lmfreq_ = smd_lmfreq
|
|
|
|
|
tol_ = smd_tol
|
|
|
|
|
maxlm_ = smd_maxlm
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! ... if smlm
|
|
|
|
|
|
|
|
|
|
IF(smd_smlm) THEN
|
|
|
|
|
IF(smd_ene_ini >= 0.d0 .OR. smd_ene_fin >= 0.d0) &
|
|
|
|
|
& CALL errore(' start : ',' Check : ene_ini & ene_fin ', 1 )
|
|
|
|
|
ENDIF
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! translate from input to internals of CP, various checks
|
|
|
|
|
|
|
|
|
|
! ... Set the number of species
|
|
|
|
|
|
|
|
|
|
nsp_ = ntyp
|
|
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
if( .not. lneb ) then
|
|
|
|
|
CALL ions_base_init( ntyp , nat , na_inp , sp_pos , rd_pos , rd_vel, atom_mass, &
|
2004-06-30 02:12:51 +08:00
|
|
|
atom_label, if_pos, atomic_positions )
|
2004-07-23 20:29:51 +08:00
|
|
|
end if
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
! ... IBRAV and CELLDM
|
|
|
|
|
|
|
|
|
|
ibrav_ = ibrav
|
|
|
|
|
celldm_ = celldm
|
|
|
|
|
|
|
|
|
|
! ... Set Values for bands and spin
|
|
|
|
|
|
|
|
|
|
n_ = nbnd * nspin
|
|
|
|
|
nspin_ = nspin
|
|
|
|
|
|
|
|
|
|
! ... Set Values for the cutoff
|
|
|
|
|
|
2004-10-14 16:38:05 +08:00
|
|
|
CALL ecutoffs_setup( ecutwfc, ecutrho, ecfixed, qcutz, q2sigma )
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
ampre_ = ampre
|
|
|
|
|
SELECT CASE ( restart_mode )
|
|
|
|
|
CASE ('from_scratch')
|
|
|
|
|
nbeg_ = -2
|
|
|
|
|
nomore_ = nstep
|
|
|
|
|
trane_ = ( startingwfc == 'random' )
|
|
|
|
|
if ( ampre_ == 0.d0 ) ampre_ = 0.02
|
|
|
|
|
CASE ('reset_counters')
|
|
|
|
|
nbeg_ = 0
|
|
|
|
|
nomore_ = nstep
|
|
|
|
|
CASE ('restart')
|
|
|
|
|
nbeg_ = 1
|
|
|
|
|
nomore_ = nstep
|
|
|
|
|
if ( ion_positions == 'from_input' ) then
|
|
|
|
|
taurdr = .TRUE.
|
|
|
|
|
nbeg_ = -1
|
|
|
|
|
end if
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown restart_mode '//trim(restart_mode), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ndr_ = ndr
|
|
|
|
|
ndw_ = ndw
|
|
|
|
|
iprint_ = iprint
|
|
|
|
|
|
|
|
|
|
IF( .NOT. lneb ) THEN
|
|
|
|
|
CALL check_stop_init( max_seconds )
|
|
|
|
|
END IF
|
|
|
|
|
|
|
|
|
|
! ... TORTHO
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( orthogonalization )
|
|
|
|
|
CASE ('Gram-Schmidt')
|
|
|
|
|
tortho_ = .FALSE.
|
|
|
|
|
CASE ('ortho')
|
|
|
|
|
tortho_ = .TRUE.
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown orthogonalization '//&
|
|
|
|
|
trim(orthogonalization), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( electron_velocities )
|
|
|
|
|
CASE ('default')
|
|
|
|
|
continue
|
|
|
|
|
CASE ('zero')
|
|
|
|
|
print '("Warning: electron_velocities keyword has no effect")'
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' electron_velocities='// &
|
|
|
|
|
trim(electron_velocities)//' not implemented', 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
! ... TSDE
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( electron_dynamics )
|
|
|
|
|
CASE ('sd')
|
|
|
|
|
tsde_ = .TRUE.
|
|
|
|
|
frice_ = 0.d0
|
|
|
|
|
CASE ('verlet')
|
|
|
|
|
tsde_ = .FALSE.
|
|
|
|
|
frice_ = 0.d0
|
|
|
|
|
CASE ('damp')
|
|
|
|
|
tsde_ = .FALSE.
|
|
|
|
|
frice_ = electron_damping
|
|
|
|
|
CASE ('none')
|
|
|
|
|
tsde_ = .FALSE.
|
|
|
|
|
frice_ = 0.d0
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown electron_dynamics '//&
|
|
|
|
|
trim(electron_dynamics),1)
|
|
|
|
|
END SELECT
|
|
|
|
|
|
2004-09-13 16:49:06 +08:00
|
|
|
SELECT CASE ( electron_velocities )
|
|
|
|
|
CASE ('zero')
|
|
|
|
|
tzeroe = .TRUE.
|
|
|
|
|
CASE ('default')
|
|
|
|
|
tzeroe = .FALSE.
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown electron_velocities '//&
|
|
|
|
|
trim(electron_dynamics),1)
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
! Ion velocities
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( ion_velocities )
|
|
|
|
|
CASE ('default')
|
|
|
|
|
tcap_ = .false.
|
|
|
|
|
CASE ('random')
|
|
|
|
|
tcap_ = .true.
|
|
|
|
|
CASE ('zero')
|
2004-09-13 16:49:06 +08:00
|
|
|
tzerop = .TRUE.
|
2004-06-30 02:12:51 +08:00
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown ion_velocities '//trim(ion_velocities),1)
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
! ... TFOR TSDP
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( ion_dynamics )
|
|
|
|
|
CASE ('sd')
|
|
|
|
|
tsdp_ = .TRUE.
|
|
|
|
|
tfor_ = .TRUE.
|
|
|
|
|
fricp_= 0.d0
|
|
|
|
|
CASE ('verlet')
|
|
|
|
|
tsdp_ = .FALSE.
|
|
|
|
|
tfor_ = .TRUE.
|
|
|
|
|
fricp_= 0.d0
|
|
|
|
|
CASE ('damp')
|
|
|
|
|
tsdp_ = .FALSE.
|
|
|
|
|
tfor_ = .TRUE.
|
|
|
|
|
fricp_= ion_damping
|
|
|
|
|
CASE ('none')
|
|
|
|
|
tsdp_ = .FALSE.
|
|
|
|
|
tfor_ = .FALSE.
|
|
|
|
|
fricp_= 0.d0
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown ion_dynamics '//trim(ion_dynamics), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( cell_velocities )
|
|
|
|
|
CASE ('default')
|
2004-09-13 16:49:06 +08:00
|
|
|
tzeroc = .FALSE.
|
2004-06-30 02:12:51 +08:00
|
|
|
CASE ('zero')
|
2004-09-13 16:49:06 +08:00
|
|
|
tzeroc = .TRUE.
|
2004-06-30 02:12:51 +08:00
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown cell_velocities '//trim(cell_velocities),1)
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
! For SMD
|
|
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
IF( lsmd ) THEN
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( cell_dynamics )
|
|
|
|
|
CASE ('none')
|
|
|
|
|
tpre_ = .FALSE.
|
|
|
|
|
thdyn_= .FALSE.
|
|
|
|
|
frich_= 0.d0
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' smiosys ',' cell_dynamics not implemented : '//trim(cell_dynamics), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
ELSE
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( cell_dynamics )
|
|
|
|
|
CASE ('sd')
|
|
|
|
|
tpre_ = .TRUE.
|
|
|
|
|
thdyn_= .TRUE.
|
|
|
|
|
frich_= 0.d0
|
|
|
|
|
CASE ('pr')
|
|
|
|
|
tpre_ = .TRUE.
|
|
|
|
|
thdyn_= .TRUE.
|
|
|
|
|
frich_= 0.d0
|
|
|
|
|
CASE ('damp-pr')
|
|
|
|
|
tpre_ = .TRUE.
|
|
|
|
|
thdyn_= .TRUE.
|
|
|
|
|
frich_ = cell_damping
|
|
|
|
|
CASE ('none')
|
|
|
|
|
tpre_ = .FALSE.
|
|
|
|
|
thdyn_= .FALSE.
|
|
|
|
|
frich_= 0.d0
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown cell_dynamics '//trim(cell_dynamics), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
END IF
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( electron_temperature )
|
|
|
|
|
! temperature control of electrons via Nose' thermostat
|
|
|
|
|
! EKINW (REAL(DBL)) average kinetic energy (in atomic units)
|
|
|
|
|
! FNOSEE (REAL(DBL)) frequency (in terahertz)
|
|
|
|
|
CASE ('nose')
|
|
|
|
|
tnosee_ = .TRUE.
|
|
|
|
|
CASE ('not_controlled')
|
|
|
|
|
tnosee_ = .FALSE.
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown electron_temperature '//&
|
|
|
|
|
trim(electron_temperature), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( ion_temperature )
|
|
|
|
|
! temperature control of ions via Nose' thermostat
|
|
|
|
|
! TEMPW (REAL(DBL)) frequency (in which units?)
|
|
|
|
|
! FNOSEP (REAL(DBL)) temperature (in which units?)
|
|
|
|
|
CASE ('nose')
|
|
|
|
|
tnosep_ = .TRUE.
|
|
|
|
|
tcp_ = .false.
|
|
|
|
|
CASE ('not_controlled')
|
|
|
|
|
tnosep_ = .FALSE.
|
|
|
|
|
tcp_ = .false.
|
|
|
|
|
CASE ('rescaling' )
|
|
|
|
|
tnosep_ = .FALSE.
|
|
|
|
|
tcp_ = .true.
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown ion_temperature '//&
|
|
|
|
|
trim(ion_temperature), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( cell_temperature )
|
|
|
|
|
! cell temperature control of ions via Nose' thermostat
|
|
|
|
|
! FNOSEH (REAL(DBL)) frequency (in which units?)
|
|
|
|
|
! TEMPH (REAL(DBL)) temperature (in which units?)
|
|
|
|
|
CASE ('nose')
|
|
|
|
|
tnoseh_ = .TRUE.
|
|
|
|
|
CASE ('not_controlled')
|
|
|
|
|
tnoseh_ = .FALSE.
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown cell_temperature '//&
|
|
|
|
|
trim(cell_temperature), 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( cell_dofree )
|
|
|
|
|
CASE ('all')
|
|
|
|
|
thdiag_ =.false.
|
|
|
|
|
CASE ('xyz')
|
|
|
|
|
thdiag_ =.true.
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' unknown cell_dofree '//trim(cell_dofree), 1 )
|
|
|
|
|
END SELECT
|
2004-08-27 18:20:42 +08:00
|
|
|
if(thdyn_) then
|
|
|
|
|
if(thdiag_) then
|
|
|
|
|
iforceh_=0
|
|
|
|
|
do i=1,3
|
|
|
|
|
iforceh_(i,i)=1
|
|
|
|
|
enddo
|
|
|
|
|
else
|
|
|
|
|
iforceh_=1
|
|
|
|
|
endif
|
|
|
|
|
endif
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
|
|
|
|
|
! ... radii, masses
|
|
|
|
|
|
|
|
|
|
DO is = 1, nsp_
|
|
|
|
|
rcmax_ (is) = ion_radius(is)
|
|
|
|
|
IF( ion_radius(is) <= 0.d0 ) THEN
|
|
|
|
|
CALL errore(' iosys ',' invalid ion_radius ', is)
|
|
|
|
|
END IF
|
|
|
|
|
END DO
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! compatibility between FPMD and CP90
|
|
|
|
|
!
|
|
|
|
|
|
2004-09-13 16:49:06 +08:00
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
tconvthrs%active = .FALSE.
|
|
|
|
|
IF( ion_dynamics == 'none' .AND. cell_dynamics == 'none' ) THEN
|
|
|
|
|
tconvthrs%ekin = ekin_conv_thr
|
|
|
|
|
tconvthrs%derho = etot_conv_thr
|
|
|
|
|
tconvthrs%force = 10d+10
|
|
|
|
|
tconvthrs%active = .TRUE.
|
|
|
|
|
tconvthrs%nstep = 1
|
|
|
|
|
END IF
|
|
|
|
|
|
2004-08-27 18:20:42 +08:00
|
|
|
isave_ = isave
|
2004-06-30 02:12:51 +08:00
|
|
|
tprnfor_ = tprnfor
|
|
|
|
|
if ( trim( verbosity ) == 'high' ) then
|
|
|
|
|
iprsta_ = 3
|
|
|
|
|
else
|
|
|
|
|
iprsta_ = 1
|
|
|
|
|
end if
|
|
|
|
|
delt_ = dt
|
|
|
|
|
emass_ = emass
|
|
|
|
|
emaec_ = emass_cutoff
|
|
|
|
|
eps_ = ortho_eps
|
|
|
|
|
max_ = ortho_max
|
2004-08-27 18:20:42 +08:00
|
|
|
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
if ( tstress ) tpre_ = .true.
|
|
|
|
|
trhow_ = ( trim( disk_io ) == 'high' )
|
|
|
|
|
tvlocw_ = .false. ! temporaneo
|
|
|
|
|
!
|
|
|
|
|
qne_ = 0.0d0
|
|
|
|
|
qnp_ = 0.0d0
|
|
|
|
|
qnh_ = 0.0d0
|
|
|
|
|
if( fnosee > 0.0d0 ) qne_ = 4.d0*ekincw/(fnosee*(2.d0*pi)*terahertz)**2
|
|
|
|
|
if( fnosep > 0.0d0 ) qnp_ = 2.d0*(3*nat)*tempw/factem/(fnosep*(2.d0*pi)*terahertz)**2
|
|
|
|
|
if( fnoseh > 0.0d0 ) qnh_ = 2.d0*(3*3 )*temph/factem/(fnoseh*(2.d0*pi)*terahertz)**2
|
|
|
|
|
tempw_ = tempw
|
|
|
|
|
temph_ = temph
|
|
|
|
|
ekincw_ = ekincw
|
|
|
|
|
|
|
|
|
|
grease_ = grease
|
|
|
|
|
tranp_ ( 1 : nsp_ ) = tranp ( 1 : nsp_ )
|
|
|
|
|
amprp_ ( 1 : nsp_ ) = amprp ( 1 : nsp_ )
|
|
|
|
|
|
|
|
|
|
greasp_ = greasp
|
|
|
|
|
tolp_ = tolp
|
|
|
|
|
greash_ = greash
|
|
|
|
|
press_ = press
|
|
|
|
|
|
|
|
|
|
nr1_ = nr1
|
|
|
|
|
nr2_ = nr2
|
|
|
|
|
nr3_ = nr3
|
|
|
|
|
nr1s_ = nr1s
|
|
|
|
|
nr2s_ = nr2s
|
|
|
|
|
nr3s_ = nr3s
|
|
|
|
|
nr1b_ = nr1b
|
|
|
|
|
nr2b_ = nr2b
|
|
|
|
|
nr3b_ = nr3b
|
|
|
|
|
|
|
|
|
|
nat_ = nat
|
|
|
|
|
pseudo_dir_ = pseudo_dir
|
|
|
|
|
|
|
|
|
|
! cards parameters
|
|
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
IF( .NOT. lsmd ) THEN
|
|
|
|
|
tau0_ = 0.0
|
|
|
|
|
END IF
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
iforce_= 0
|
|
|
|
|
psfile_= ' '
|
|
|
|
|
a1_ = 0.0
|
|
|
|
|
a2_ = 0.0
|
|
|
|
|
a3_ = 0.0
|
|
|
|
|
|
2004-08-27 18:20:42 +08:00
|
|
|
! masses are brought to au
|
|
|
|
|
|
|
|
|
|
pmass_ ( 1:nsp_ ) = atom_mass( 1:nsp_ ) * scmass
|
2004-06-30 02:12:51 +08:00
|
|
|
psfile_ ( 1:nsp_ ) = atom_pfile( 1:nsp_ )
|
|
|
|
|
|
|
|
|
|
! ... Getting 'na
|
|
|
|
|
|
|
|
|
|
na_ = 0
|
|
|
|
|
isa = 0
|
|
|
|
|
do is = 1, nsp_
|
|
|
|
|
do ia = 1, nat_
|
|
|
|
|
if ( sp_pos(ia) == is) then
|
|
|
|
|
na_(is) = na_(is) + 1
|
2004-06-30 16:49:41 +08:00
|
|
|
isa = isa + 1
|
2004-06-30 02:12:51 +08:00
|
|
|
if( na_(is) > natx ) call errore(' cards',' na > natx', na_ (is) )
|
2004-07-23 20:29:51 +08:00
|
|
|
IF( .NOT. lsmd ) THEN
|
2004-08-27 18:20:42 +08:00
|
|
|
tau0_ (:, isa ) = rd_pos(:, ia)
|
2004-07-23 20:29:51 +08:00
|
|
|
END IF
|
2004-08-27 18:20:42 +08:00
|
|
|
iforce_(:, isa ) = if_pos(:, ia)
|
2004-06-30 02:12:51 +08:00
|
|
|
end if
|
|
|
|
|
end do
|
|
|
|
|
end do
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! set up atomic positions and crystal lattice
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
if ( ibrav_ == 0 ) then
|
|
|
|
|
a1_ = rd_ht( 1, 1:3 )
|
|
|
|
|
a2_ = rd_ht( 2, 1:3 )
|
|
|
|
|
a3_ = rd_ht( 3, 1:3 )
|
|
|
|
|
if ( celldm_ (1) == 0.d0 ) then
|
|
|
|
|
celldm_ (1) = sqrt( a1_ (1) ** 2 + a1_ (2) ** 2 + a1_ (3) ** 2 )
|
|
|
|
|
a1_(:) = a1_(:) / celldm_(1)
|
|
|
|
|
a2_(:) = a2_(:) / celldm_(1)
|
|
|
|
|
a3_(:) = a3_(:) / celldm_(1)
|
|
|
|
|
end if
|
|
|
|
|
else
|
|
|
|
|
call latgen( ibrav_ , celldm_ , a1_ , a2_ , a3_ , omega_ )
|
|
|
|
|
end if
|
|
|
|
|
alat_ = celldm_ (1)
|
|
|
|
|
|
|
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
IF( lsmd ) THEN
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! How to obtain the initial trial path.
|
|
|
|
|
!
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
IF(smd_smopt) THEN
|
|
|
|
|
|
|
|
|
|
CALL init_path(sm_p,kwnp_,smd_stcd,nsp_,nat_,alat_,nbeg_,1)
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
ELSEIF(smd_linr) THEN
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
CALL init_path(sm_p,kwnp_,smd_stcd,nsp_,nat_,alat_,nbeg_,2)
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
ELSEIF(smd_polm .AND. (smd_kwnp < num_of_images) ) THEN
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
CALL init_path(sm_p,kwnp_,smd_stcd,nsp_,nat_,alat_,nbeg_,3)
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
ELSEIF(smd_kwnp == num_of_images ) THEN
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
CALL init_path(sm_p,kwnp_,smd_stcd,nsp_,nat_,alat_,nbeg_,4)
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
ENDIF
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-07-23 20:29:51 +08:00
|
|
|
ELSE
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( atomic_positions )
|
|
|
|
|
!
|
|
|
|
|
! convert input atomic positions to internally used format:
|
|
|
|
|
! tau0 in atomic units
|
|
|
|
|
!
|
|
|
|
|
CASE ('alat')
|
|
|
|
|
!
|
|
|
|
|
! input atomic positions are divided by a0
|
|
|
|
|
!
|
|
|
|
|
tau0_ = tau0_ * alat_
|
|
|
|
|
CASE ('bohr')
|
|
|
|
|
!
|
|
|
|
|
! input atomic positions are in a.u.: do nothing
|
|
|
|
|
!
|
|
|
|
|
continue
|
|
|
|
|
CASE ('crystal')
|
|
|
|
|
!
|
|
|
|
|
! input atomic positions are in crystal axis ("scaled"):
|
|
|
|
|
!
|
|
|
|
|
taus = tau0_
|
2004-08-27 18:20:42 +08:00
|
|
|
isa = 0
|
2004-07-23 20:29:51 +08:00
|
|
|
do is = 1, nsp_
|
|
|
|
|
do ia = 1, na_(is)
|
2004-08-27 18:20:42 +08:00
|
|
|
isa = isa + 1
|
2004-07-23 20:29:51 +08:00
|
|
|
do i = 1, 3
|
2004-08-27 18:20:42 +08:00
|
|
|
tau0_ ( i, isa ) = a1_ (i) * taus( 1, isa ) &
|
|
|
|
|
+ a2_ (i) * taus( 2, isa ) &
|
|
|
|
|
+ a3_ (i) * taus( 3, isa )
|
2004-07-23 20:29:51 +08:00
|
|
|
end do
|
|
|
|
|
end do
|
|
|
|
|
end do
|
|
|
|
|
CASE ('angstrom')
|
|
|
|
|
!
|
|
|
|
|
! atomic positions in A
|
|
|
|
|
!
|
|
|
|
|
tau0_ = tau0_ / 0.529177
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' atomic_positions='//trim(atomic_positions)// &
|
|
|
|
|
' not implemented ', 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
END IF
|
2004-06-30 02:12:51 +08:00
|
|
|
|
2004-08-27 18:20:42 +08:00
|
|
|
! ... Set the default value for the cell mass
|
|
|
|
|
wmass_ = wmass
|
|
|
|
|
IF( wmass_ == 0.d0 ) THEN
|
|
|
|
|
wmass_ = 3.d0 / (4.d0 * pi**2 ) * SUM( atom_mass(1:ntyp)*na_inp(1:ntyp) )
|
|
|
|
|
wmass_ = wmass_ * UMA_AU
|
|
|
|
|
WRITE( stdout,999) wmass_
|
|
|
|
|
ELSE
|
|
|
|
|
WRITE( stdout,998) wmass_
|
|
|
|
|
END IF
|
|
|
|
|
998 format(' wmass (read from input) = ',f15.2,/)
|
|
|
|
|
999 format(' wmass (calculated) = ',f15.2,/)
|
|
|
|
|
|
2004-06-30 02:12:51 +08:00
|
|
|
|
|
|
|
|
!
|
|
|
|
|
!
|
|
|
|
|
! set occupancies
|
|
|
|
|
!
|
|
|
|
|
|
|
|
|
|
IF( nelec < 1 ) THEN
|
|
|
|
|
CALL errore(' iosys ',' nelec less than 1 ', int(nelec) )
|
|
|
|
|
END IF
|
|
|
|
|
IF( nint(nelec) - nelec > eps8 ) THEN
|
|
|
|
|
CALL errore(' iosys ',' nelec must be integer', int(nelec) )
|
|
|
|
|
END IF
|
|
|
|
|
|
|
|
|
|
if( mod( n_ , 2 ) .ne. 0 ) then
|
|
|
|
|
nx_ = n_ + 1
|
|
|
|
|
else
|
|
|
|
|
nx_= n_
|
|
|
|
|
end if
|
|
|
|
|
|
|
|
|
|
iupdwn_ ( 1 ) = 1
|
|
|
|
|
nel_ = 0
|
|
|
|
|
|
|
|
|
|
SELECT CASE ( TRIM(occupations) )
|
|
|
|
|
CASE ('bogus')
|
|
|
|
|
!
|
|
|
|
|
! empty-states calculation: occupancies have a (bogus) finite value
|
|
|
|
|
!
|
|
|
|
|
! bogus to ensure \sum_i f_i = Nelec (nelec is integer)
|
|
|
|
|
!
|
|
|
|
|
f_ ( : ) = nelec / n_
|
|
|
|
|
nel_ (1) = nint(nelec)
|
|
|
|
|
nupdwn_ (1) = n_
|
|
|
|
|
if ( nspin_ == 2 ) then
|
|
|
|
|
!
|
|
|
|
|
! bogus to ensure Nelec = Nup + Ndw
|
|
|
|
|
!
|
|
|
|
|
nel_ (1) = ( nint(nelec) + 1 ) / 2
|
|
|
|
|
nel_ (2) = nint(nelec) / 2
|
|
|
|
|
nupdwn_ (1)=nbnd
|
|
|
|
|
nupdwn_ (2)=nbnd
|
|
|
|
|
iupdwn_ (2)=nbnd+1
|
|
|
|
|
end if
|
|
|
|
|
CASE ('from_input')
|
|
|
|
|
!
|
|
|
|
|
! occupancies have been read from input
|
|
|
|
|
!
|
|
|
|
|
f_ ( 1:nbnd ) = f_inp( 1:nbnd, 1 )
|
|
|
|
|
if( nspin_ == 2 ) f_ ( nbnd+1 : 2*nbnd ) = f_inp( 1:nbnd, 2 )
|
|
|
|
|
if( nelec == 0.d0 ) nelec = SUM ( f_ ( 1:n_ ) )
|
|
|
|
|
if( nspin_ == 2 .and. nelup == 0) nelup = SUM ( f_ ( 1:nbnd ) )
|
|
|
|
|
if( nspin_ == 2 .and. neldw == 0) neldw = SUM ( f_ ( nbnd+1 : 2*nbnd ) )
|
|
|
|
|
|
|
|
|
|
if( nspin_ == 1 ) then
|
|
|
|
|
nel_ (1) = nint(nelec)
|
|
|
|
|
nupdwn_ (1) = n_
|
|
|
|
|
else
|
|
|
|
|
IF ( ABS (nelup + neldw - nelec) > eps8 ) THEN
|
|
|
|
|
CALL errore(' iosys ',' wrong # of up and down spin', 1 )
|
|
|
|
|
END IF
|
|
|
|
|
nel_ (1) = nint(nelup)
|
|
|
|
|
nel_ (2) = nint(neldw)
|
|
|
|
|
nupdwn_ (1)=nbnd
|
|
|
|
|
nupdwn_ (2)=nbnd
|
|
|
|
|
iupdwn_ (2)=nbnd+1
|
|
|
|
|
end if
|
|
|
|
|
|
|
|
|
|
CASE ('fixed')
|
|
|
|
|
|
|
|
|
|
if( nspin_ == 1 ) then
|
|
|
|
|
nel_ (1) = nint(nelec)
|
|
|
|
|
nupdwn_ (1) = n_
|
|
|
|
|
else
|
|
|
|
|
IF ( nelup + neldw /= nelec ) THEN
|
|
|
|
|
CALL errore(' iosys ',' wrong # of up and down spin', 1 )
|
|
|
|
|
END IF
|
|
|
|
|
nel_ (1) = nint(nelup)
|
|
|
|
|
nel_ (2) = nint(neldw)
|
|
|
|
|
nupdwn_ (1)=nbnd
|
|
|
|
|
nupdwn_ (2)=nbnd
|
|
|
|
|
iupdwn_ (2)=nbnd+1
|
|
|
|
|
end if
|
|
|
|
|
|
|
|
|
|
! ocp = 2 for spinless systems, ocp = 1 for spin-polarized systems
|
|
|
|
|
ocp = 2.d0 / nspin_
|
|
|
|
|
! default filling: attribute ocp electrons to each states
|
|
|
|
|
! until the good number of electrons is reached
|
|
|
|
|
do iss = 1, nspin_
|
|
|
|
|
fsum = 0.0d0
|
|
|
|
|
do in = iupdwn_ ( iss ), iupdwn_ ( iss ) - 1 + nupdwn_ ( iss )
|
|
|
|
|
if ( fsum + ocp < nel_ ( iss ) + 0.0001 ) then
|
|
|
|
|
f_ (in) = ocp
|
|
|
|
|
else
|
|
|
|
|
f_ (in) = max( nel_ ( iss ) - fsum, 0.d0 )
|
|
|
|
|
end if
|
|
|
|
|
fsum=fsum + f_(in)
|
|
|
|
|
end do
|
|
|
|
|
end do
|
|
|
|
|
|
|
|
|
|
CASE DEFAULT
|
|
|
|
|
CALL errore(' iosys ',' occupation method not implemented', 1 )
|
|
|
|
|
END SELECT
|
|
|
|
|
|
|
|
|
|
do iss = 1, nspin_
|
|
|
|
|
do in = iupdwn_(iss), iupdwn_(iss) - 1 + nupdwn_(iss)
|
|
|
|
|
ispin_(in) = iss
|
|
|
|
|
end do
|
|
|
|
|
end do
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
! --------------------------------------------------------
|
|
|
|
|
! print out heading
|
|
|
|
|
!
|
|
|
|
|
WRITE( stdout,500) nbeg_ , nomore_ , iprint_ , ndr_ , ndw_
|
|
|
|
|
WRITE( stdout,505) delt_
|
|
|
|
|
WRITE( stdout,510) emass_ , emaec_
|
|
|
|
|
!
|
|
|
|
|
if( tortho_ ) then
|
|
|
|
|
WRITE( stdout,511) eps_ , max_
|
|
|
|
|
else
|
|
|
|
|
WRITE( stdout,512)
|
|
|
|
|
endif
|
|
|
|
|
!
|
|
|
|
|
if( tsde_ ) then
|
|
|
|
|
WRITE( stdout,513)
|
|
|
|
|
else
|
|
|
|
|
if ( tnosee_ ) frice_ = 0.
|
|
|
|
|
WRITE( stdout,509)
|
|
|
|
|
WRITE( stdout,514) frice_ , grease_
|
|
|
|
|
endif
|
|
|
|
|
!
|
|
|
|
|
if ( trhor_ ) then
|
|
|
|
|
WRITE( stdout,720)
|
|
|
|
|
endif
|
|
|
|
|
!
|
|
|
|
|
if( .not. trhor_ .and. trhow_ )then
|
|
|
|
|
WRITE( stdout,721)
|
|
|
|
|
endif
|
|
|
|
|
!
|
|
|
|
|
if( tvlocw_ )then
|
|
|
|
|
WRITE( stdout,722)
|
|
|
|
|
endif
|
|
|
|
|
!
|
|
|
|
|
if( trane_ ) then
|
|
|
|
|
WRITE( stdout,515) ampre_
|
|
|
|
|
endif
|
|
|
|
|
WRITE( stdout,516)
|
|
|
|
|
do is =1, nsp_
|
|
|
|
|
if(tranp_(is)) WRITE( stdout,517) is, amprp_(is)
|
|
|
|
|
end do
|
|
|
|
|
!
|
|
|
|
|
if(tfor_) then
|
|
|
|
|
if(tnosep_) fricp_ = 0.
|
|
|
|
|
WRITE( stdout,520)
|
|
|
|
|
if(tsdp_)then
|
|
|
|
|
WRITE( stdout,521)
|
|
|
|
|
else
|
|
|
|
|
WRITE( stdout,522) fricp_ , greasp_
|
|
|
|
|
endif
|
|
|
|
|
else
|
|
|
|
|
WRITE( stdout,518)
|
|
|
|
|
endif
|
|
|
|
|
!
|
|
|
|
|
if( tfor_ ) then
|
|
|
|
|
if(( tcp_ .or. tcap_ .or. tnosep_ ) .and. tsdp_ ) then
|
|
|
|
|
call errore(' main',' t contr. for ions when tsdp=.t.',0)
|
|
|
|
|
endif
|
|
|
|
|
if(.not. tcp_ .and. .not. tcap_ .and. .not. tnosep_ ) then
|
|
|
|
|
WRITE( stdout,550)
|
|
|
|
|
else if(tcp_ .and. tcap_ ) then
|
|
|
|
|
call errore(' main',' tcp and tcap both true',0)
|
|
|
|
|
else if(tcp_ .and. tnosep_ ) then
|
|
|
|
|
call errore(' main',' tcp and tnosep both true',0)
|
|
|
|
|
else if(tcap_ .and. tnosep_ ) then
|
|
|
|
|
call errore(' main',' tcap and tnosep both true',0)
|
|
|
|
|
else if(tcp_ ) then
|
|
|
|
|
WRITE( stdout,555) tempw_ , tolp_
|
|
|
|
|
else if(tcap_) then
|
|
|
|
|
WRITE( stdout,560) tempw_ , tolp_
|
|
|
|
|
else if(tnosep_ ) then
|
|
|
|
|
WRITE( stdout,562) tempw_ , qnp_
|
|
|
|
|
end if
|
|
|
|
|
if(tnosee_) then
|
|
|
|
|
WRITE( stdout,566) ekincw_ , qne_
|
|
|
|
|
end if
|
|
|
|
|
end if
|
|
|
|
|
!
|
|
|
|
|
if(tpre_) then
|
|
|
|
|
WRITE( stdout,600)
|
|
|
|
|
if(thdyn_) then
|
|
|
|
|
if(thdiag_) WRITE( stdout,608)
|
|
|
|
|
if(tnoseh_) then
|
|
|
|
|
frich_=0.
|
|
|
|
|
WRITE( stdout,604) temph_,qnh_,press_
|
|
|
|
|
else
|
|
|
|
|
WRITE( stdout,602) frich_,greash_,press_
|
|
|
|
|
endif
|
|
|
|
|
else
|
|
|
|
|
WRITE( stdout,606)
|
|
|
|
|
endif
|
|
|
|
|
endif
|
2004-10-14 16:38:05 +08:00
|
|
|
if ( agg .ne. 0.d0) then
|
|
|
|
|
WRITE( stdout,650) agg, sgg, e0gg
|
2004-06-30 02:12:51 +08:00
|
|
|
end if
|
|
|
|
|
WRITE( stdout,700) iprsta_
|
|
|
|
|
|
|
|
|
|
!
|
|
|
|
|
500 format(// &
|
|
|
|
|
& ' nbeg=',i3,' nomore=',i7,3x,' iprint=',i4,/ &
|
|
|
|
|
& ' reads from',i3,' writes on',i3)
|
|
|
|
|
505 format(' time step = ',f9.4/)
|
|
|
|
|
510 format(' parameters for electron dynamics:'/ &
|
|
|
|
|
& ' emass= ',f10.2,2x,'emaec= ',f10.2,'ry')
|
|
|
|
|
511 format(' orthog. with lagrange multipliers: eps=',e10.2, &
|
|
|
|
|
& ' max=',i3)
|
|
|
|
|
512 format(' orthog. with gram-schmidt')
|
|
|
|
|
513 format(' electron dynamics with steepest descent')
|
|
|
|
|
509 format(' verlet algorithm for electron dynamics')
|
|
|
|
|
514 format(' with friction frice = ',f7.4,' , grease = ',f7.4)
|
|
|
|
|
720 format(' charge density is read from unit 47',/)
|
|
|
|
|
721 format(' charge density is written in unit 47',/)
|
|
|
|
|
722 format(' local potential is written in unit 46',/)
|
|
|
|
|
515 format(' initial random displacement of el. coordinates with ', &
|
|
|
|
|
& ' amplitude=',f10.6,/ &
|
|
|
|
|
& ' trane not to be used with mass preconditioning')
|
|
|
|
|
516 format(/)
|
|
|
|
|
517 format(' initial random displacement of ionic coord. for species ',&
|
|
|
|
|
& i4,' : amplitude=',f10.6)
|
|
|
|
|
518 format(' ions are not allowed to move'/)
|
|
|
|
|
520 format(' ions are allowed to move')
|
|
|
|
|
521 format(' ion dynamics with steepest descent')
|
|
|
|
|
522 format(' ion dynamics with fricp = ',f7.4,' and greasp = ',f7.4)
|
|
|
|
|
550 format(' ion dynamics: the temperature is not controlled'//)
|
|
|
|
|
555 format(' ion dynamics with rescaling of velocities:'/ &
|
|
|
|
|
& ' temperature required=',f10.5,'(kelvin)',' tolerance=', &
|
|
|
|
|
& f10.5//)
|
|
|
|
|
560 format(' ion dynamics with canonical temp. control:'/ &
|
|
|
|
|
& ' temperature required=',f10.5,'(kelvin)',' tolerance=', &
|
|
|
|
|
& f10.5//)
|
|
|
|
|
562 format(' ion dynamics with nose` temp. control:'/ &
|
|
|
|
|
& ' temperature required=',f10.5,'(kelvin)',' nose` mass = ',&
|
|
|
|
|
& f10.3//)
|
|
|
|
|
566 format(' electronic dynamics with nose` temp. control:'/ &
|
|
|
|
|
& ' elec. kin. en. required=',f10.5,'(hartree)', &
|
|
|
|
|
& ' nose` mass = ',f10.3//)
|
|
|
|
|
600 format(' internal stress tensor calculated')
|
|
|
|
|
602 format(' cell parameters dynamics with frich = ',f7.4, &
|
|
|
|
|
& ' and greash = ',f7.4,/ &
|
|
|
|
|
& ' external pressure = ',f11.7,'(gpa)'//)
|
|
|
|
|
604 format(' cell parameters dynamics with nose` temp. control:'/ &
|
|
|
|
|
& ' cell temperature required = ',f10.5,'(kelvin)', &
|
|
|
|
|
& ' nose` mass = ',f10.3,/ &
|
|
|
|
|
& ' external pressure = ',f11.7,'(gpa)'//)
|
|
|
|
|
606 format(' cell parameters are not allowed to move'//)
|
|
|
|
|
608 format(' frozen off-diagonal cell parameters'//)
|
|
|
|
|
650 format(' modified kinetic energy functional, with parameters:'/ &
|
|
|
|
|
& ' agg = ',f8.4,' sgg = ', f7.4,' e0gg = ',f6.2)
|
|
|
|
|
700 format(' iprsta = ',i2/)
|
|
|
|
|
return
|
|
|
|
|
end subroutine
|
|
|
|
|
!
|
|
|
|
|
|
All namelists and cards moved to Modules/input_parameters.f90 .
From now on, all new input variables should be added
to this module, and then copied to the code internal
variables in the input.f90 subroutine
The namelists and cards parsers are in :
Modules/read_namelists.f90 and Modules/read_cards.f90
files input_parameters.f90 read_namelists.f90 read_cards.f90
are shared by all codes, while each code has its own version
of input.f90 ( used to copy input values into internals variables ).
EXAMPLE:
suppose you need to add a new input variable called "pippo"
to the namelist control, then:
1) add pippo to the input_parameters.f90 file containing the
namelist control
INTEGER :: pippo = 0
NAMELIST / control / ....., pippo
remember: always set an initialization value!
2) add pippo to the control_default subroutine
( cantained in module read_namelists.f90 )
subroutine control_default( prog )
...
IF( prog == 'PW' ) pippo = 10
...
end subroutine
this routine set the default value for pippo,
that could vary with the code
3) add pippo to the control_bcast subroutine
( cantained in module read_namelists.f90 )
subroutine control_bcast( )
...
call mp_bcast( pippo )
...
end subroutine
4) add pippo to the control_checkin subroutine
( cantained in module read_namelists.f90 )
subroutine control_checking( prog )
...
IF( pippo < 0 ) &
CALL error(' control_checkin ',' variable pippo less than 0 ', 1 )
...
end subroutine
5) Copy the value of pippo in the code internal variables
( file input.f90 )
subroutine iosys()
use input_parameters, only: ...., pippo
use pwcom, only: ....., myvar
...
call read_namelists( 'PW' )
...
myvar = pippo
...
end subroutine
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@282 c92efa57-630b-4861-b058-cf58834340f0
2003-07-31 21:24:20 +08:00
|
|
|
END MODULE input_cp
|
2004-06-30 02:12:51 +08:00
|
|
|
|