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