C^^^^^1992 1/7 FOR STRESS ^^^^^^^^^^^
SUBROUTINE FORZFB
C^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
IMPLICIT REAL(A-H,O-Y)
IMPLICIT COMPLEX(Z)
INCLUDE 'PACVPP'
C VPP-PARALLEL
!XOCL SUBPROCESSOR PS(IPARA)=PQ(1:IPARA)
!XOCL INDEX PARTITION IP=(PS,INDEX=1:KNV3,PART=BAND)
DIMENSION SSS(KNG1,KNV3,KTYP,10)
& ,ZZZ(KNG1,KEG,KNV3),ZFC(KEG,KNV3,KATM,10)
!XOCL LOCAL SSS(:,/IP,:,:),ZZZ(:,:,/IP),ZFC(:,/IP,:,:)
EQUIVALENCE (SNL,SSS),(ZAJ,ZZZ),(ZFBB,ZFC)
C
cGEN
DIMENSION CWLNEW(10)
cGEN
C!XOCL PARALLEL REGION
!XOCL SPREAD DO /IP
DO 1000 IA=1,KATM
IF (NLSPD(KFTYPE(IA)).EQ.1) THEN
LNUM = 4
ELSE
LNUM = 9
END IF
DO 1003 L=1,LNUM
cGEN
*pdir pardo
DO 1001 IK=1,KV3
DO 1002 IBAN=NBD1,NBD2
ZFC(IBAN,IK,IA,L) =DCMPLX(0.0D0,0.0D0)
1002 CONTINUE
1001 CONTINUE
1003 CONTINUE
1000 CONTINUE
!XOCL END SPREAD
C!XOCL END PARALLEL
C VPP-PARALLEL START
C!XOCL PARALLEL REGION
!XOCL SPREAD DO /IP
cGEN
*pdir pardo
DO 2100 IK=1,KV3
DO 2110 IA=1,KATM
CS=1.0D0/(WS(KFTYPE(IA))*UNIVOL)
CP=1.0D0/(WP(KFTYPE(IA))*UNIVOL)
CWLNEW(1)=CS
CWLNEW(2)=CP
CWLNEW(3)=CP
CWLNEW(4)=CP
C
IF (NLSPD(KFTYPE(IA)).EQ.1) THEN
LNUM = 4
ELSE
LNUM = 9
CD=1.0D0/(WD(KFTYPE(IA))*UNIVOL)
CWLNEW(5)=CD
CWLNEW(6)=CD
CWLNEW(7)=CD
CWLNEW(8)=CD
CWLNEW(9)=CD
END IF
DO 3222 L=1,LNUM
DO 3200 IBAN=NBD1,NBD2
C 1991 11/28 I ---> I1
DO 3510 I=1,IBA(IK)
I1 = NBASE(I,IK)
L1 = IG1(I1)+KX1
L2 = IG2(I1)+KY1
L3 = IG3(I1)+KZ1
ZTMP=ZZZ(I,IBAN,IK)*DCONJG( ZFM2( I1 ,IA ) )
ZFC(IBAN,IK,IA,L)=ZFC(IBAN,IK,IA,L)+
& ZTMP*SSS(I,IK,KFTYPE(IA),L)
3510 CONTINUE
ZFC(IBAN,IK,IA,L)=CWLNEW(L)*ZFC(IBAN,IK,IA,L)
3200 CONTINUE
3222 CONTINUE
C
2110 CONTINUE
2100 CONTINUE
!XOCL END SPREAD
C!XOCL END PARALLEL
C
RETURN
END