File: AZERO.FT of Tape: Various/ETH/f2
(Source file text)
SUBROUTINE AZERO C COMMON /IN/ ANF(6,6),AF(6,6),ETA(6,6),DAF(6,6),DETA(6,6) $ ,DK,HANG,HI,NP,T(6),DT(6),MR,ADENT(6),MF,FREQ(6) COMMON /MOD/ AM(6,6),A0(6,6),A0M(6),DA0M(6),AF0(6,6) $ ,ANFM(6),ETA0(6,6) C DO 10 J=1,MR DO 11 I=1,MF TAU=ETA(I,J)*DK*NP AM(I,J)=AF(I,J)*TAU/(1-EXP(-TAU)) 11 A0(I,J)=AM(I,J)*SQRT(HI)*EXP(ETA(I,J)/HI)/T(J) 10 CONTINUE DO 12 I=1,MF SP=0. SA=0. SB=0. SC=0. DO 13 J=1,MR A0P=A0(I,J)**2*((DAF(I,J)/AF(I,J))**2+(DETA(I,J)/ETA(I,J))**2) A0P=1./A0P SA=SA+A0(I,J)*A0P SB=SB+A0(I,J)**2*A0P SC=SC+ANF(I,J) SP=SP+A0P 13 CONTINUE A0M(I)=SA/SP DA0M(I)=SQRT((SB/SP-A0M(I)**2)/(MR-1)) ANFM(I)=SC/MR 12 CONTINUE DO 14 I=1,MF DO 15 J=1,MR ETA0(I,J)=HI*ALOG(A0M(I)*T(J)/SQRT(HI)/AM(I,J)) TAU0=ETA0(I,J)*DK*NP AF0(I,J)=AM(I,J)*(1.-EXP(-TAU0))/TAU0 15 CONTINUE 14 CONTINUE RETURN END