File: DINGLE.FT of Tape: Various/ETH/f2

(Source file text) 

---

C        PROGRAM DINGLE(OUTPUT=240B,TAPE4=OUTPUT,TAPE5,TAPE6)
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)
        COMMON /LINE/ DY(6),A,B,DA,DB,SA,SB,SAB,RR,RMSD,CHISQ 
C 
        DIMENSION X(400),Y(400),NB(400) 
	DIMENSION XL(10),YL(10)
C 
100     FORMAT(I3)
110     FORMAT((F8.1))
130     FORMAT((A6))
C 
        READ(5,100) MF
        READ(5,110) (FREQ(I),I=1,MF)
        READ(5,100) MR
        READ(5,130) (ADENT(I),I=1,MR) 
        CALL INPUT
        CALL AZERO
        CALL PRINT(1,0) 
	N=MR
        DO 50 I=1,MF
        DO 51 J=1,MR
        YL(J)=ETA(I,J)
        DY(J)=DETA(I,J) 
        XL(J)=T(J) 
51	CONTINUE
        CALL LINREG(XL,YL,N)
	CALL PRINT(2,I)
	TMAX=0.
	DO 10 K=1,MR
10	TMAX=AMAX1(TMAX,T(K))
	LMIN=1
	LMAX=(400-4*MR)/2
	DO 61 L=1,LMAX
	X(L)=1.+(TMAX-1.)*(L-LMIN)/FLOAT(LMAX-LMIN)
	Y(L)=A*X(L)+B
	NB(L)=4
61	CONTINUE
        DO 52 J=1,MR
        YL(J)=ETA0(I,J)
        DY(J)=DETA(I,J) 
        XL(J)=T(J) 
52	CONTINUE
        CALL LINREG(XL,YL,N)
        CALL PRINT(3,I) 
	LMIN=LMAX+1
	LMAX=LMAX*2
	DO 62 L=LMIN,LMAX
	X(L)=1.+(TMAX-1.)*(L-LMIN)/FLOAT(LMAX-LMIN)
	Y(L)=A*X(L)+B
	NB(L)=10
62	CONTINUE
	DO 60 J=1,MR
	K=4*(J-1)+LMAX+1
	X(K)=T(J)
	Y(K)=ETA(I,J)+DETA(I,J)
	NB(K)=2
	X(K+1)=T(J)
	Y(K+1)=ETA(I,J)-DETA(I,J)
	NB(K+1)=2
	X(K+2)=T(J)
	Y(K+2)=ETA(I,J)
	NB(K+2)=1
	X(K+3)=T(J)
	Y(K+3)=ETA0(I,J)
	NB(K+3)=3
60	CONTINUE
	NPOINT=400
	CALL GRAPH(X,Y,NB,NPOINT)
50      CONTINUE
        END 

---