File: DLOGAM.FT of Tape: Various/ETH/prog1
(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