11
11
integer zv(2),amass(2),rc(2)
12
12
integer lmmax(2),lmax(2),locp(2),nmax(2)
13
13
integer nprj(2),n_projector(2),l_projector(2),m_projector(2)
14
integer b_projector(2)
14
15
integer npsp,nprj_max,nprj_mult
15
16
integer psp_type(2)
19
19
integer nkatmx,lmmax_max,lmax_max,nmax_max
20
20
parameter (nkatmx=50,lmax_max=4,lmmax_max=16,nmax_max=3)
23
23
common / psp_block / prjtmp,vl,vnl,Gijl,zv,amass,rc,
24
24
> lmmax,lmax,locp,nmax,
25
25
> nprj,n_projector,l_projector,m_projector,
26
> psp_type,core_kin,npsp,nprj_max,nprj_mult,
26
> b_projector,r3_matrix,
27
> psp_type,npsp,nprj_max,nprj_mult,
30
30
* **** common block for psp stress ****
33
33
common / psp2_block / dvl,dvnl
36
* ***** PAW variables *****
37
integer log_amesh(2),r1(2),rmax(2),sigma(2),zion(2)
38
integer n1dgrid(2),n1dbasis(2),nae(2),nps(2),lps(2),icut(2),eig(2)
39
integer phi_ae(2),dphi_ae(2),phi_ps(2),dphi_ps(2)
40
integer core_ae(2),core_ps(2),core_ae_prime(2),core_ps_prime(2)
41
integer rgrid(2),core_kin(2),core_ion(2)
42
integer hartree_matrix(2),comp_charge_matrix(2),comp_pot_matrix(2)
44
common /psp_paw_block/ log_amesh,r1,rmax,sigma,zion,
45
> n1dgrid,n1dbasis,nae,nps,lps,icut,eig,
46
> phi_ae,dphi_ae,phi_ps,dphi_ps,
47
> core_ae,core_ps,core_ae_prime,core_ps_prime,
48
> rgrid,core_kin,core_ion,
50
> comp_charge_matrix,comp_pot_matrix,
53
* **** common block for psp_xc_data ****
57
integer paw_xc_dylm_theta(2)
58
integer paw_xc_dylm_phi(2)
59
integer paw_xc_angle_phi(2)
60
integer paw_xc_cos_theta(2)
61
integer paw_xc_w_phi(2)
62
integer paw_xc_w_theta(2)
64
integer i_paw_xc_pot(2)
67
integer paw_dvxc_ae(2)
68
integer paw_dvxc_ps(2)
69
integer paw_xc_pot_size
71
integer paw_xc_matr(2)
72
integer paw_xc_dmatr(2)
73
common /psp_xc_data/ paw_xc_nphi,paw_xc_ntheta,paw_xc_e,
74
> i_paw_xc_pot,paw_xc_angle_phi,
75
> paw_xc_cos_theta, paw_xc_ylm,
76
> paw_xc_dylm_theta,paw_xc_dylm_phi,
77
> paw_xc_w_phi,paw_xc_w_theta,paw_vxc_ae,
78
> paw_vxc_ps,paw_dvxc_ae,paw_dvxc_ps,
79
> paw_xc_pot,paw_xc_pot_size,paw_xc_matr,
85
* **** common block for pspspin ****
87
integer pspspin_upl,pspspin_downl
88
integer pspspin_upions(2),pspspin_downions(2)
89
real*8 pspspin_upscale,pspspin_downscale
90
common / pspspin_block / pspspin_upscale,pspspin_downscale,
91
> pspspin_upions, pspspin_downions,
92
> pspspin_upl, pspspin_downl,