File: DLOGAM.FT of Tape: Various/ETH/f4

(Source file text) 

---

      SUBROUTINE LOGAMP
C
      COMMON/INI/ ANF(10,10),DANF(10,10),PH(10,10),DPH(10,10) 
     $,AF(10,10),DAF(10,10),ETA(10,10),DETA(10,10),EMU(10,10), 
     $T(10),DT(10),ANFM(10),DANFM(10)
CDC  $,F(10,512)
      COMMON/TESPAR/ A0M,DA0M,AMU0,AF0(10),ETA0(10),A0(10)
     $,TDJ(10),DTDJ(10),TDM,DTDM
      COMMON/PARAM/IDENT,HANG,HI,HF,DK,NP,PI,SQHI,DKHI,I0DH,H0AL
      COMMON/STEER/MR,ADENT(10),MF,FREQ(10),I0
      COMMON /REGRES/ XL(10),YL(10),DY(10)
      COMMON/FRERUN/ I,JR,INIT,LS
C 
      REAL KET,I0DH
C 450 FORMAT(1X,'YL(',I3,' ) = ',G10.3,'    DY(',I3,
C    $' ) = ',G10.3,/)
C 460 FORMAT(1X,'AF=',G10.3,'  TAUD=',G10.3,'  KET=',G10.3,
C    $'  I0DH=',G10.3,'  T=',G10.3)
C 
      DO 50 J=1,MR
      TAUD=ETA0(J)*DK*NP 
      DTAUD=DETA(I,J)*DK*NP 
      KET=(1.-EXP(-TAUD)) 
C     WRITE(4,460) AF0(J),TAUD,KET,I0DH,T(J)
      YL(J)=H0AL*(ALOG(ABS(AF0(J)*TAUD/T(J)/KET/I0DH))- 
     $I0*TAUD/NP)
      DY(J)=H0AL*SQRT((DAF(I,J)/AF0(J))**2+
     $((1./TAUD+EXP(-TAUD)/KET-I0/NP)*DTAUD)**2)
C     WRITE(4,450) J,YL(J),J,DY(J)
      XL(J)=T(J)
   50 CONTINUE
      RETURN
      END

---