Commit 1d2ceae5 authored by dalcorso's avatar dalcorso

Minor cleanup of the phonon code (part 2). USE phcom and USE pwcom

substituted by multiple USE statements.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@5368 c92efa57-630b-4861-b058-cf58834340f0
parent 8deedecc
......@@ -12,12 +12,21 @@ subroutine addnlcc (imode0, drhoscf, npe)
!
#include "f_defs.h"
USE kinds, only : DP
USE ions_base, ONLY : nat
use funct, only : dft_is_gradient
use pwcom
USE cell_base, ONLY : omega, alat
use scf, only : rho, rho_core
USE kinds, only : DP
use phcom
USE gvect, ONLY : nrxx, g, ngm, nl, nrx1, nrx2, nrx3, nr1, nr2, nr3
USE lsda_mod, ONLY : nspin
USE spin_orb, ONLY : domag
USE dynmat, ONLY : dyn, dyn_rec
USE modes, ONLY : nirr, npert, npertx
USE gc_ph, ONLY: grho, dvxc_rr, dvxc_sr, dvxc_ss, dvxc_s
USE eqv, ONLY : dmuxc
USE nlcc_ph, ONLY : nlcc_any
USE qpoint, ONLY : xq
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -11,11 +11,21 @@
SUBROUTINE addnlcc_zstar_eu_us( drhoscf )
!----------===================-------------------
USE kinds, ONLY : DP
USE funct, only : dft_is_gradient
USE pwcom
USE scf, only : rho, rho_core
USE kinds, ONLY : DP
USE phcom
USE cell_base, ONLY : omega, alat
USE lsda_mod, ONLY : nspin
USE gvect, ONLY : nrxx, ngm, nl, g, nrx1, nrx2, nrx3, nr1, nr2, nr3
USE spin_orb, ONLY : domag
USE efield_mod, ONLY : zstareu0
USE qpoint, ONLY : xq
USE nlcc_ph, ONLY : nlcc_any
USE modes, ONLY : npert, nirr
USE eqv, ONLY : dmuxc
USE gc_ph, ONLY: grho, dvxc_rr, dvxc_sr, dvxc_ss, dvxc_s
USE mp_global, ONLY : my_pool_id
......
......@@ -22,7 +22,9 @@ subroutine addusdbec (ik, wgt, psi, dbecsum)
USE wvfct, only: npw, npwx, nbnd
USE uspp, only: nkb, vkb, okvan
USE uspp_param, only: upf, nh, nhm
use phcom
USE phus, ONLY : becp1
USE qpoint, ONLY : npwq
USE control_ph, ONLY : nbnd_occ, lgamma
implicit none
!
! the dummy variables
......
......@@ -24,7 +24,10 @@ subroutine addusdbec_nc (ik, wgt, psi, dbecsum_nc)
USE uspp, only: nkb, vkb, okvan
USE noncollin_module, ONLY : noncolin, npol
USE uspp_param, only: upf, nh, nhm
use phcom
USE qpoint, ONLY : npwq
USE phus, ONLY : becp1, becp1_nc
USE control_ph, ONLY : nbnd_occ, lgamma
implicit none
!
! the dummy variables
......
......@@ -27,10 +27,16 @@ subroutine addusddens (drhoscf, dbecsum, irr, mode0, npe, iflag)
USE cell_base, ONLY : tpiba
USE ions_base, ONLY : nat, ityp, ntyp => nsp
USE wavefunctions_module, ONLY: psic
use phcom
USE uspp_param, ONLY: upf, lmaxq, nh, nhm
USE paw_variables, ONLY : okpaw
USE modes, ONLY : u, npert, npertx
USE qpoint, ONLY : xq, eigqts
USE phus, ONLY : becsumort, alphasum
USE units_ph, ONLY : iudrhous, lrdrhous
USE control_ph, ONLY : lgamma
implicit none
!
! the dummy variables
......
......@@ -17,12 +17,16 @@ subroutine addusddense (drhoscf, dbecsum)
!
#include "f_defs.h"
USE ions_base, ONLY : nat, ityp, ntyp => nsp
use pwcom
use phcom
USE kinds, only : DP
USE ions_base, ONLY : nat, ityp, ntyp => nsp
USE gvect, ONLY : nrxx, nrx1, nrx2, nrx3, nr1, nr2, nr3, nl, g, gg, &
ngm, eigts1, eigts2, eigts3, ig1, ig2, ig3
USE uspp, ONLY: okvan
USE uspp_param, ONLY: upf, lmaxq, nh, nhm
USE qpoint, ONLY : eigqts
USE lsda_mod, ONLY : nspin
USE spin_orb, ONLY : domag
implicit none
!
! the dummy variables
......
......@@ -17,13 +17,18 @@ subroutine addusdynmat (dynwrk)
!
#include "f_defs.h"
USE kinds, only : DP
USE ions_base, ONLY : nat, ityp
use pwcom
USE noncollin_module, ONLY : noncolin, npol
USE kinds, only : DP
USE uspp, ONLY: okvan, becsum
USE uspp_param, only: upf, nh
use phcom
USE lsda_mod, ONLY : nspin
USE spin_orb, ONLY : lspinorb
USE phus, ONLY : int1, int1_nc, int2, int2_so, int4, int4_nc, &
int5, int5_so, alphasum, alphasum_nc, becsum, becsum_nc
USE modes, ONLY : nmodes
implicit none
complex(DP) :: dynwrk (3 * nat, 3 * nat)
......
......@@ -13,12 +13,11 @@ subroutine allocate_part
! dynamical allocation of arrays for the control of partial computation
! of the dynamical matrix
!
#include "f_defs.h"
USE ions_base, ONLY : nat
use pwcom
USE kinds, only : DP
use phcom
USE ions_base, ONLY : nat
USE partial, ONLY : comp_irr, done_irr, list, atomo
implicit none
!
......
......@@ -15,20 +15,35 @@ subroutine allocate_phq
!
#include "f_defs.h"
USE kinds, only : DP
USE ions_base, ONLY : nat, ntyp => nsp
use pwcom
USE klist, only : nks
USE wvfct, ONLY : nbnd, igk, npwx
USE gvect, ONLY : nrxx, ngm
USE lsda_mod, ONLY : nspin
USE noncollin_module, ONLY : noncolin, npol
USE wavefunctions_module, ONLY: evc
USE kinds, only : DP
USE spin_orb, ONLY : lspinorb
use phcom
use el_phon
USE becmod, ONLY: becp, becp_nc
USE uspp, ONLY: okvan, nkb
USE paw_variables, ONLY : okpaw
USE uspp_param, ONLY: nhm
USE ramanm, ONLY: ramtns, lraman
USE paw_variables, ONLY : okpaw
USE qpoint, ONLY : nksq, eigqts, igkq
USE phus, ONLY : int1, int1_nc, int2, int2_so, int3, int3_nc, int3_paw, &
int4, int4_nc, int5, int5_so, becsumort, dpqq, &
dpqq_so, alphasum, alphasum_nc, becsum_nc, &
becp1, becp1_nc, alphap, alphap_nc
USE efield_mod, ONLY : zstareu, zstareu0, zstarue0, zstarue0_rec, zstarue
USE eqv, ONLY : dpsi, evq, vlocq, dmuxc, dvpsi
USE units_ph, ONLY : this_pcxpsi_is_on_file, this_dvkb3_is_on_file
USE dynmat, ONLY : dyn00, dyn, dyn_rec, w2
USE modes, ONLY : u, ubar, rtau, max_irr_dim, npert, t, tmq
USE control_ph, ONLY : elph, lgamma
USE el_phon, ONLY : el_ph_mat
implicit none
!
! allocate space for the quantities needed in the phonon program
......
......@@ -15,11 +15,14 @@ subroutine cch_psi_all (n, h, ah, e, ik, m)
!
#include "f_defs.h"
use pwcom
use becmod
USE uspp, ONLY: nkb, vkb
USE kinds, only : DP
use phcom
USE becmod, ONLY : becp, calbec
USE uspp, ONLY: nkb, vkb
USE wvfct, ONLY : npwx, nbnd
USE control_ph, ONLY : alpha_pv, nbnd_occ, lgamma
USE eqv, ONLY : evq
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -15,12 +15,15 @@ subroutine ch_psi_all (n, h, ah, e, ik, m)
!
#include "f_defs.h"
use pwcom
use becmod
USE kinds, only : DP
USE wvfct, ONLY : npwx, nbnd
USE becmod, ONLY : becp, becp_nc, calbec
USE uspp, ONLY: nkb, vkb
USE noncollin_module, ONLY : noncolin, npol
USE kinds, only : DP
use phcom
USE control_ph, ONLY : alpha_pv, nbnd_occ, lgamma
USE eqv, ONLY : evq
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -18,13 +18,19 @@ subroutine compute_alphasum
!
#include "f_defs.h"
USE kinds, only : DP
USE ions_base, ONLY : nat, ityp, ntyp => nsp
use pwcom
USE lsda_mod, ONLY : current_spin, isk, lsda
USE wvfct, ONLY : nbnd, wg
USE noncollin_module, ONLY : noncolin, npol
USE kinds, only : DP
USE uspp, ONLY: okvan
USE uspp_param, ONLY: upf, nh
use phcom
USE phus, ONLY : alphasum, alphasum_nc, becp1, becp1_nc, alphap, &
alphap_nc
USE qpoint, ONLY : nksq
USE control_ph, ONLY : nbnd_occ, lgamma
implicit none
integer :: ik, ikk, ikq, ijkb0, ijh, ikb, jkb, ih, jh, na, nt, &
......
......@@ -18,13 +18,19 @@ subroutine compute_becsum_ph
!
#include "f_defs.h"
!
USE kinds, only : DP
USE ions_base, ONLY : nat, ityp, ntyp => nsp
use pwcom
USE lsda_mod, ONLY : current_spin, isk, lsda
USE wvfct, ONLY : nbnd, wg
USE noncollin_module, ONLY : noncolin, npol
USE kinds, only : DP
USE uspp, ONLY: okvan, becsum
USE uspp_param, ONLY: upf, nh
use phcom
USE phus, ONLY : alphasum, alphasum_nc, becp1, becp1_nc, alphap, &
alphap_nc, becsum_nc
USE qpoint, ONLY : nksq
USE control_ph, ONLY : nbnd_occ, lgamma
implicit none
integer :: ik, ikk, ikq, ijkb0, ijh, ikb, jkb, ih, jh, na, nt, ibnd
......
......@@ -26,7 +26,12 @@ subroutine compute_drhous (drhous, dbecsum, wgg, becq, alpq)
USE gvect, ONLY : nrxx
USE gsmooth, ONLY : nrxxs, nr1s,nr2s,nr3s, nrx1s,nrx2s,nrx3s, nls
USE wvfct, ONLY : npw, nbnd, igk
use phcom
USE qpoint, ONLY : nksq, igkq, npwq
USE eqv, ONLY : evq
USE units_ph, ONLY : iuwfc, lrwfc
USE control_ph, ONLY : lgamma
implicit none
!
! the dummy variables
......
......@@ -28,7 +28,12 @@ subroutine compute_drhous_nc (drhous, dbecsum, wgg, becq, alpq)
USE io_files, ONLY: iunigk
USE uspp, ONLY: okvan, nkb, vkb
USE uspp_param, ONLY: nhm
use phcom
USE qpoint, ONLY : nksq, igkq, npwq
USE eqv, ONLY : evq
USE units_ph, ONLY : lrwfc, iuwfc
USE control_ph, ONLY : lgamma
implicit none
!
! the dummy variables
......
......@@ -26,7 +26,13 @@ subroutine compute_nldyn (wdyn, wgg, becq, alpq)
USE uspp_param,ONLY : nh
USE spin_orb, ONLY : lspinorb
USE wvfct, ONLY : nbnd, et
use phcom
USE qpoint, ONLY : nksq
USE modes, ONLY : u
USE phus, ONLY : becp1, becp1_nc, alphap, alphap_nc, int1, int2, &
int2_so, int1_nc
USE control_ph, ONLY : nbnd_occ, lgamma
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -18,7 +18,8 @@ subroutine compute_weight (wgg)
USE klist, ONLY : wk, lgauss, degauss, ngauss
USE ener, ONLY : ef
USE wvfct, ONLY : nbnd, wg, et
use phcom
USE qpoint, ONLY : nksq
USE control_ph, ONLY : lgamma
implicit none
real(DP) :: wgg (nbnd, nbnd, nksq)
......
......@@ -11,12 +11,31 @@ subroutine deallocate_phq
!
! deallocates the variables allocated by allocate_phq
!
use phcom
USE noncollin_module, ONLY : m_loc
USE becmod, ONLY: becp, becp_nc
USE wavefunctions_module, ONLY: evc
USE ramanm, ONLY: ramtns
use el_phon
USE modes, ONLY : tmq, t, npert, u, ubar, rtau
USE qpoint, ONLY : eigqts, igkq
USE efield_mod, ONLY : zstareu, zstarue, zstarue0, zstareu0, &
zstarue0_rec
USE phus, ONLY : int1, int1_nc, int2, int2_so, int3, int3_nc, int3_paw, &
int4, int4_nc, int5, int5_so, becsum_nc, &
becsumort, alphasum, alphasum_nc, dpqq, dpqq_so, &
becp1, becp1_nc, alphap, alphap_nc
USE gc_ph, ONLY: grho, gmag, dvxc_rr, dvxc_sr, dvxc_ss, dvxc_s, &
vsgga, segni
USE gamma_gamma, ONLY : with_symmetry, has_equivalent, equiv_atoms, &
n_equiv_atoms
USE eqv, ONLY : dmuxc, vlocq, dpsi, dvpsi, evq
USE nlcc_ph, ONLY : drc
USE units_ph, ONLY : this_dvkb3_is_on_file, this_pcxpsi_is_on_file
USE dynmat, ONLY : dyn00, dyn_rec, dyn, w2
USE control_ph, ONLY : lgamma
USE el_phon, ONLY : el_ph_mat
IMPLICIT NONE
if(allocated(ramtns)) deallocate (ramtns)
if (lgamma) then
......
......@@ -30,7 +30,14 @@ subroutine drho
USE wvfct, ONLY : nbnd
USE spin_orb, ONLY : domag
USE paw_variables, ONLY : okpaw
use phcom
USE control_ph, ONLY : ldisp, all_done, lgamma, recover
USE dynmat, ONLY : dyn00
USE qpoint, ONLY : nksq
USE modes, ONLY : npertx, npert, nirr
USE phus, ONLY : becsumort, alphap, becp1, alphap_nc, becp1_nc
USE units_ph, ONLY : lrdrhous, iudrhous
USE mp_global, ONLY : inter_pool_comm, intra_pool_comm
USE mp, ONLY : mp_sum
......
......@@ -26,7 +26,14 @@ subroutine drhodv (nu_i0, nper, drhoscf)
USE io_global, ONLY : stdout
USE noncollin_module, ONLY : noncolin, npol
USE io_files, ONLY: iunigk
use phcom
USE dynmat, ONLY : dyn, dyn_rec
USE modes, ONLY : u, npertx
USE qpoint, ONLY : nksq, npwq, igkq
USE eqv, ONLY : dpsi
USE units_ph, ONLY : lrdwf, iuwfc, iudwf
USE control_ph, ONLY : lgamma
USE mp_global, ONLY : inter_pool_comm, intra_pool_comm
USE mp, ONLY : mp_sum
......
......@@ -21,7 +21,8 @@ subroutine drhodvloc (nu_i0, nper, drhoscf, wdyn)
USE gsmooth, ONLY : nrxxs, nr1s, nr2s, nr3s
USE cell_base, ONLY : omega
USE lsda_mod, ONLY : nspin
use phcom
USE modes, ONLY : npertx
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -29,7 +29,12 @@ subroutine drhodvus (irr, imode0, dvscfin, npe)
USE io_global, ONLY : stdout
USE uspp_param, ONLY : upf, nh
USE paw_variables, ONLY : okpaw
use phcom
USE modes, ONLY : npert, npertx, nirr
USE dynmat, ONLY : dyn, dyn_rec
USE phus, ONLY : becsumort, int3_paw
USE units_ph, ONLY : iudrhous, lrdrhous
USE mp_global, ONLY : inter_pool_comm, intra_pool_comm
USE mp, ONLY : mp_sum
implicit none
......
......@@ -23,7 +23,13 @@ subroutine dv_of_drho (mode, dvscf, flag)
USE spin_orb, ONLY : domag
USE funct, ONLY : dft_is_gradient
USE scf, ONLY : rho, rho_core
use phcom
USE eqv, ONLY : dmuxc
USE nlcc_ph, ONLY : nlcc_any
USE qpoint, ONLY : xq
USE gc_ph, ONLY : grho, dvxc_rr, dvxc_sr, dvxc_ss, dvxc_s
USE control_ph, ONLY : lrpa
implicit none
integer :: mode
......
......@@ -27,9 +27,14 @@ subroutine dvanqq
use spin_orb, only : lspinorb, domag
use scf, only : v, vltot
use noncollin_module, ONLY : noncolin, npol
use phcom
USE uspp, ONLY: okvan
USE uspp_param, ONLY: upf, lmaxq, nh
USE phus, ONLY : int1, int2, int4, int4_nc, int5, int5_so
USE eqv, ONLY : vlocq
USE qpoint, ONLY : eigqts, xq
USE control_ph, ONLY : recover, ldisp, lgamma
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -29,7 +29,11 @@ subroutine dvqpsi_us (ik, mode, uact, addnlcc)
use uspp_param,ONLY : upf
USE wvfct, ONLY : nbnd, npw, npwx, igk
USE wavefunctions_module, ONLY: evc
use phcom
USE nlcc_ph, ONLY : nlcc_any, drc
USE eqv, ONLY : dvpsi, dmuxc, vlocq
USE qpoint, ONLY : npwq, igkq, xq, eigqts
USE control_ph, ONLY : lgamma
implicit none
!
! The dummy variables
......
......@@ -28,7 +28,12 @@ subroutine dvqpsi_us_only (ik, mode, uact)
USE noncollin_module, ONLY : noncolin, npol
USE uspp, ONLY: okvan, nkb, vkb, qq, qq_so, deeq, deeq_nc
USE uspp_param, ONLY: nh
use phcom
USE qpoint, ONLY : igkq, npwq
USE phus, ONLY : int1, int1_nc, int2, int2_so, alphap, alphap_nc, &
becp1, becp1_nc
USE eqv, ONLY : dvpsi
USE control_ph, ONLY : lgamma
implicit none
!
! The dummy variables
......
......@@ -31,9 +31,16 @@ SUBROUTINE dynmat_us()
USE uspp_param, ONLY : nh
USE noncollin_module, ONLY : noncolin, npol
USE spin_orb, ONLY : lspinorb
USE phcom
USE becmod, ONLY : calbec
USE io_global, ONLY : stdout
USE qpoint, ONLY : npwq, nksq, igkq
USE modes, ONLY : u
USE dynmat, ONLY : dyn
USE phus, ONLY : becp1, becp1_nc, alphap, alphap_nc
USE control_ph, ONLY : nbnd_occ, lgamma
USE units_ph, ONLY : iuwfc, lrwfc
USE mp_global, ONLY : my_pool_id, inter_pool_comm, intra_pool_comm
USE mp, ONLY : mp_sum
......
......@@ -22,7 +22,10 @@ subroutine incdrhoscf (drhoscf, weight, ik, dbecsum)
USE wvfct, ONLY : npw, igk
USE uspp_param,ONLY: nhm
USE wavefunctions_module, ONLY: evc
use phcom
USE eqv, ONLY : dpsi
USE qpoint, ONLY : npwq, igkq
USE control_ph, ONLY : nbnd_occ, lgamma
implicit none
integer :: ik
......
......@@ -26,7 +26,10 @@ subroutine incdrhoscf_nc (drhoscf, weight, ik, dbecsum)
USE uspp_param,ONLY : nhm
USE wvfct, ONLY : npw, npwx, igk
USE wavefunctions_module, ONLY: evc
use phcom
USE qpoint, ONLY : npwq, igkq
USE eqv, ONLY : dpsi
USE control_ph, ONLY : nbnd_occ, lgamma
implicit none
integer :: ik
......
......@@ -28,7 +28,11 @@ subroutine newdq (dvscf, npe)
USE uspp_param, ONLY : upf, lmaxq, nh, nhm
USE spin_orb, ONLY : domag
USE paw_variables, ONLY : okpaw
use phcom
USE phus, ONLY : int3, int3_paw
USE qpoint, ONLY : xq, eigqts
USE control_ph, ONLY : lgamma
USE mp_global, ONLY: intra_pool_comm
USE mp, ONLY: mp_sum
......
......@@ -365,7 +365,8 @@ MODULE ph_restart
! ... this routine reads the format version of the current xml datafile
!
USE parser, ONLY : version_compare
USE xml_io_base
USE xml_io_base, ONLY : attr
IMPLICIT NONE
!
CHARACTER(LEN=*), INTENT(IN) :: dirname
......
......@@ -50,7 +50,7 @@ SUBROUTINE phq_readin()
USE control_flags, ONLY : twfcollect
USE paw_variables, ONLY : okpaw
USE ramanm, ONLY : eth_rps, eth_ns, lraman, elop, dek
USE freq_ph
USE freq_ph, ONLY : fpol, fiu, nfs
!
IMPLICIT NONE
!
......
......@@ -71,14 +71,31 @@ subroutine phq_setup
USE rap_point_group, ONLY : code_group, nclass, nelem, elem, which_irr,&
char_mat, name_rap, gname, name_class, ir_ram
USE rap_point_group_is, ONLY : code_group_is, gname_is
use phcom
USE nlcc_ph, ONLY : drc, nlcc_any
USE eqv, ONLY : dmuxc
USE control_ph, ONLY : rec_code, lgamma_gamma, search_sym, start_irr, &
last_irr, niter_ph, alpha_mix, all_done, &
epsil, lgamma, recover, where_rec, alpha_pv, &
nbnd_occ
USE output, ONLY : fildrho
USE modes, ONLY : u, ubar, npertx, npert, gi, gimq, nirr, &
max_irr_dim, t, tmq, irotmq, irgq, minus_q, &
nsymq, nmodes, rtau, invs
USE dynmat, ONLY : dyn, dyn_rec, dyn00
USE efield_mod, ONLY : epsilon, zstareu
USE qpoint, ONLY : xq
USE partial, ONLY : comp_irr, atomo, nat_todo, list, nrapp, all_comp, &
npert, done_irr
USE gamma_gamma, ONLY : has_equivalent, asr, nasr, n_diff_sites, &
equiv_atoms, n_equiv_atoms, with_symmetry
USE control_ph, ONLY : rec_code
USE ph_restart, ONLY : ph_writefile
USE control_flags, ONLY : iverbosity, modenum, noinv
USE funct, ONLY : dmxc, dmxc_spin, dmxc_nc, dft_is_gradient
USE ramanm, ONLY : lraman, elop, ramtns, eloptns
USE mp, ONLY : mp_max, mp_min
USE mp_global, ONLY : inter_pool_comm
USE ramanm, ONLY : lraman, elop, ramtns, eloptns
implicit none
......
......@@ -19,7 +19,7 @@ subroutine symdvscf (nper, irr, dvtosym)
USE cell_base, ONLY : at
USE symme, ONLY : s, ftau
USE lsda_mod, ONLY: nspin
use phcom
USE modes, ONLY : minus_q, irotmq, nsymq, irgq, gi, t, tmq, gimq
implicit none
integer :: nper, irr
......
......@@ -22,7 +22,6 @@ subroutine syme (dvsym)
USE symme, only : nsym, s, ftau
USE lsda_mod, only : nspin
USE kinds, only : DP
use phcom
implicit none
complex(DP) :: dvsym (nrx1, nrx2, nrx3, nspin, 3)
......
......@@ -31,8 +31,17 @@ subroutine zstar_eu_us
USE wavefunctions_module, ONLY : evc
USE uspp_param, ONLY : upf, nhm, nh
USE noncollin_module, ONLY : noncolin, npol
USE qpoint, ONLY : nksq, npwq
USE control_ph, ONLY : nbnd_occ
USE efield_mod, ONLY : zstareu0
USE eqv, ONLY : dvpsi, dpsi
USE phus, ONLY : becsumort, int3, int3_paw
USE modes, ONLY : u, npert, nirr
USE nlcc_ph, ONLY : nlcc_any
USE units_ph, ONLY : lrdwf, iucom, lrcom, lrebar, iuebar, lrdrhous, &
iudrhous, iudwf, lrwfc, iuwfc
USE mp_global, ONLY : nproc_pool
use phcom
!
implicit none
integer :: ibnd, jbnd, ipol, jpol, imode0, irr, imode, nrec, mode
......
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