$
C0M(I,J,4),C0M(I,J,5),C0M(I,J,6) )
9001
CONTINUE
9000
CONTINUE
DO 10000 I=1,NN-1
DO 10001 J=I+1,NN
TET(I,J)= PF(T,TETM(I,J,1),TETM(I,J,2),TETM(I,J,3),
$
TETM(I,J,4),TETM(I,J,5),TETM(I,J,6))
DO 10002 IA=1,NN
PSI(I,J,IA)=PF(T,PSIM(I,J,IA,1),PSIM(I,J,IA,2),
$
PSIM(I,J,IA,3),PSIM(I,J,IA,4),
$
PSIM(I,J,IA,5),PSIM(I,J,IA,6))
10002
CONTINUE
10001
CONTINUE
10000 CONTINUE
RETURN
2
CONTINUE
DO 11000 I = 1, 10
SUMCA(I) = 0.
SUMCAT(I) = 0.
SUMAN(I) = 0.
SUMZ(I)
= 0.
SUMAC(I) = 0.
SUMAA(I) = 0.
SUMCC(I) = 0.
SUMK(I)
= 0.
11000 CONTINUE
SMX=0.
ZZ=0.
UM=0.
DO 12000 I=IP+1,NMV
J=I-IP
XM(J)=EX(I)*55.50837
ZZ=ZZ+XM(J)*DABS(Z(J))
UM=UM+XM(J)*Z(J)**2
SMX=SMX+XM(J)
12000 CONTINUE
UM=UM/2.
SQ=SQRT(UM)
ALPHA=2.*SQ
ALPHA1=1.4*SQ
ALPHA2=12.*SQ
G1=2.*(1.-(1.+ALPHA1)*EXP(-ALPHA1))/ALPHA1**2
G2=2.*(1.-(1.+ALPHA2)*EXP(-ALPHA2))/ALPHA2**2
GPRI1=-2.*(1.-(1.+ALPHA1+ALPHA1**2/2.)*EXP(-ALPHA1))/ALPHA1**2
GPRI2=-2.*(1.-(1.+ALPHA2+ALPHA2**2/2.)*EXP(-ALPHA2))/ALPHA2**2
G=2.*(1-(1+ALPHA)*EXP(-ALPHA))/ALPHA**2
GPRIME=-2.*(1.-(1.+ALPHA+ALPHA**2/2.)*EXP(-ALPHA))/ALPHA**2
22