Lines Matching refs:iwa

24       integer iwfr,ipww,number,iwm,iwpm,nax,iwpj,iwa,jwa,iax,iwmn,lwwptr
75       do 6 iwa=1,mwa
78       xi(nax+iwmn,1,iwa)=xw(iwfr+iwm,1,iwa)
79       xi(nax+iwmn,2,iwa)=xw(iwfr+iwm,2,iwa)
80       xi(nax+iwmn,3,iwa)=xw(iwfr+iwm,3,iwa)
81       xj(nax+iwmn,1,iwa)=xw(lwwptr,1,iwa)
82       xj(nax+iwmn,2,iwa)=xw(lwwptr,2,iwa)
83       xj(nax+iwmn,3,iwa)=xw(lwwptr,3,iwa)
84       pl(nax+iwmn,1,iwa)=pw(iwfr+iwm,1,iwa,2)
85       pl(nax+iwmn,2,iwa)=pw(iwfr+iwm,2,iwa,2)
86       pl(nax+iwmn,3,iwa)=pw(iwfr+iwm,3,iwa,2)
87       pj(nax+iwmn,1,iwa)=pw(lwwptr,1,iwa,2)
88       pj(nax+iwmn,2,iwa)=pw(lwwptr,2,iwa,2)
89       pj(nax+iwmn,3,iwa)=pw(lwwptr,3,iwa,2)
101       do 10 iwa=1,mwa
104       xi(nax+iwmn,1,iwa)=xw(iwfr+iwm,1,iwa)
105       xi(nax+iwmn,2,iwa)=xw(iwfr+iwm,2,iwa)
106       xi(nax+iwmn,3,iwa)=xw(iwfr+iwm,3,iwa)
107       xj(nax+iwmn,1,iwa)=xw(lwwptr,1,iwa)-rwx(iwmn,1)
108       xj(nax+iwmn,2,iwa)=xw(lwwptr,2,iwa)-rwx(iwmn,2)
109       xj(nax+iwmn,3,iwa)=xw(lwwptr,3,iwa)-rwx(iwmn,3)
110       pl(nax+iwmn,1,iwa)=pw(iwfr+iwm,1,iwa,2)
111       pl(nax+iwmn,2,iwa)=pw(iwfr+iwm,2,iwa,2)
112       pl(nax+iwmn,3,iwa)=pw(iwfr+iwm,3,iwa,2)
113       pj(nax+iwmn,1,iwa)=pw(lwwptr,1,iwa,2)
114       pj(nax+iwmn,2,iwa)=pw(lwwptr,2,iwa,2)
115       pj(nax+iwmn,3,iwa)=pw(lwwptr,3,iwa,2)
124       do 12 iwa=1,mwa
127       fi(iax,ix,iwa)=zero
128       fj(iax,ix,iwa)=zero
140       do 20 iwa=1,mwa
141       qai=qfaci*chg(iwq(iwa),1,iset)
142       pai=chg(iwq(iwa),2,iset)
147       rx=xj(iax,1,jwa)-xi(iax,1,iwa)
148       ry=xj(iax,2,jwa)-xi(iax,2,iwa)
149       rz=xj(iax,3,jwa)-xi(iax,3,iwa)
150       pix=pai*pl(iax,1,iwa)
151       piy=pai*pl(iax,2,iwa)
152       piz=pai*pl(iax,3,iwa)
160       fi(iax,1,iwa)=fi(iax,1,iwa)+((rmj-qaj)*rx-pjx)*ri3
161       fi(iax,2,iwa)=fi(iax,2,iwa)+((rmj-qaj)*ry-pjy)*ri3
162       fi(iax,3,iwa)=fi(iax,3,iwa)+((rmj-qaj)*rz-pjz)*ri3
175       do 24 iwa=1,mwa
178       pw(iwfr+iwm,1,iwa,1)=pw(iwfr+iwm,1,iwa,1)+fi(iax+iwmn,1,iwa)
179       pw(iwfr+iwm,2,iwa,1)=pw(iwfr+iwm,2,iwa,1)+fi(iax+iwmn,2,iwa)
180       pw(iwfr+iwm,3,iwa,1)=pw(iwfr+iwm,3,iwa,1)+fi(iax+iwmn,3,iwa)
181       pw(lwwptr,1,iwa,1)=pw(lwwptr,1,iwa,1)+fj(iax+iwmn,1,iwa)
182       pw(lwwptr,2,iwa,1)=pw(lwwptr,2,iwa,1)+fj(iax+iwmn,2,iwa)
183       pw(lwwptr,3,iwa,1)=pw(lwwptr,3,iwa,1)+fj(iax+iwmn,3,iwa)
198       do 32 iwa=1,mwa
201       pl(nax2+iwmn,1,iwa)=pwp(iwfr+iwm,1,iwa,ipset-1,2)
202       pl(nax2+iwmn,2,iwa)=pwp(iwfr+iwm,2,iwa,ipset-1,2)
203       pl(nax2+iwmn,3,iwa)=pwp(iwfr+iwm,3,iwa,ipset-1,2)
204       pj(nax2+iwmn,1,iwa)=pwp(lwwptr,1,iwa,ipset-1,2)
205       pj(nax2+iwmn,2,iwa)=pwp(lwwptr,2,iwa,ipset-1,2)
206       pj(nax2+iwmn,3,iwa)=pwp(lwwptr,3,iwa,ipset-1,2)
215       do 40 iwa=1,mwa
218       fi(iax,ix,iwa)=zero
219       fj(iax,ix,iwa)=zero
224       do 34 iwa=1,mwa
225       qai=qfaci*chg(iwq(iwa),1,ipset)
226       pai=chg(iwq(iwa),2,ipset)
231       rx=xj(iax,1,jwa)-xi(iax,1,iwa)
232       ry=xj(iax,2,jwa)-xi(iax,2,iwa)
233       rz=xj(iax,3,jwa)-xi(iax,3,iwa)
234       pix=pai*pl(iax,1,iwa)
235       piy=pai*pl(iax,2,iwa)
236       piz=pai*pl(iax,3,iwa)
244       fi(iax,1,iwa)=fi(iax,1,iwa)+((rmj-qaj)*rx-pjx)*ri3
245       fi(iax,2,iwa)=fi(iax,2,iwa)+((rmj-qaj)*ry-pjy)*ri3
246       fi(iax,3,iwa)=fi(iax,3,iwa)+((rmj-qaj)*rz-pjz)*ri3
259       do 38 iwa=1,mwa
262       pwp(iwfr+iwm,1,iwa,ipset-1,1)=pwp(iwfr+iwm,1,iwa,ipset-1,1)+
263      + fi(iax+iwmn,1,iwa)
264       pwp(iwfr+iwm,2,iwa,ipset-1,1)=pwp(iwfr+iwm,2,iwa,ipset-1,1)+
265      + fi(iax+iwmn,2,iwa)
266       pwp(iwfr+iwm,3,iwa,ipset-1,1)=pwp(iwfr+iwm,3,iwa,ipset-1,1)+
267      + fi(iax+iwmn,3,iwa)
268       pwp(lwwptr,1,iwa,ipset-1,1)=pwp(lwwptr,1,iwa,ipset-1,1)+
269      + fj(iax+iwmn,1,iwa)
270       pwp(lwwptr,2,iwa,ipset-1,1)=pwp(lwwptr,2,iwa,ipset-1,1)+
271      + fj(iax+iwmn,2,iwa)
272       pwp(lwwptr,3,iwa,ipset-1,1)=pwp(lwwptr,3,iwa,ipset-1,1)+
273      + fj(iax+iwmn,3,iwa)