File: OUTFOU.FT of Tape: Various/ETH/f2
(Source file text)
SUBROUTINE OUTFOU(LOGOUT)
C
COMMON /TITLE/ ADENT,GLOBAL,RUN,ICREAT
COMMON /FITPAR/ JIT,FISTOP,AMARQI,UPMARQ,DNMARQ,VARI,
$ UPVAR,DNVAR,VARMIN,MLOOP,LOOPLW,ISTART,ISTOP,ITEST,IDEFIX
COMMON /PARAM/ DKHI,SHFACT,CHISQ,ITORQ,IDILA,DK,DL,HANG
$ ,HUP,HDOWN,AMI,PS,PT,SOLLT,TEMP,TV,ES,AKL,AKV,CAP,EICH
COMMON /PEAKS/ TASWI,KBLOW,KBLOW1,KMAX,K1,K3,CHI,
$ Q(11,4),HIGHT(21),INDEX(21)
COMMON /VECT/ C(44),D1(11,4),DIAGEL(44),D(44),E(44)
COMMON /DATIN/ F(512),NP
COMMON /CONST/ PI,TWOPI
COMMON /PLOTC/ PLTBUF(400)
COMMON /FAF/ A(2049),B(2049),N
DIMENSION AR(44,44),Z(44,44)
EQUIVALENCE (A,AR),(B,Z)
INTEGER RUN
REAL K3
LOGICAL ITORQ,IDILA,ITEST,TASWI
C
C LOOKOUT WITH DIM!!
DIMENSION QN(11,4),DN(11,4)
EQUIVALENCE (QN,D),(DN,E)
C
599 FORMAT(1H1,T10,'JOSS =',1L1,' DILA =',1L1,' RUN =',
$ I2,' ANGLE =',F6.1)
600 FORMAT(//T10,'P0',T30,'P1',T50,'P2',T70,'P3'/)
601 FORMAT(T11,F8.3,T30,F8.3,T50,F8.3,T70,F8.3)
602 FORMAT(/T10,'NF',T30,'PH',T50,'AF',T70,'TF'/)
603 FORMAT(T10,F9.3,T30,F8.3,T50,1PE12.3,T70,0PF8.3)
604 FORMAT(//T10,'STANDARD DEVIATIONS'/)
605 FORMAT(1H0,'X',T10,'Q(X+1)',T30,'Q(X+2)',T50,'Q(X+3)',
$ T70,'Q(X+4)')
606 FORMAT(1H ,I2,T10,F8.3,T30,F8.3,T50,F8.3,T70,F8.3)
610 FORMAT('$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$')
620 FORMAT(A1,A2,A1,I2)
630 FORMAT(I2)
640 FORMAT(8(1PE10.3))
C
ICHAN=3
IF (LOGOUT.NE.0) GO TO 90
WRITE(ICHAN,605)
DO 50 K=1,K1
KX=(K-1)*4
50 WRITE(ICHAN,606) KX,(Q(K,J),J=1,4)
RETURN
90 DO 120 J=1,4
QN(1,J)=Q(1,J)
120 DN(1,J)=SQRT(AR(J,J))
DO 100 K=2,K1
KK=4*(K-1)
QN(K,1)=Q(K,1)*NP/(DL*TWOPI)
QN(K,2)=Q(K,2)*180./PI
QN(K,3)=Q(K,3)*EICH
QN(K,4)=Q(K,4)/DL
DN(K,1)=SQRT(AR(KK+1,KK+1))*NP/(DL*TWOPI)
DN(K,2)=SQRT(AR(KK+2,KK+2))*180./PI
DN(K,3)=SQRT(AR(KK+3,KK+3))*EICH
100 DN(K,4)=SQRT(AR(KK+4,KK+4))/DL
IF (LOGOUT.EQ.-1) GO TO 200
WRITE(ICHAN,599) ITORQ,IDILA,RUN,HANG
WRITE(ICHAN,600)
WRITE(ICHAN,601) (QN(1,J),J=1,4)
WRITE(ICHAN,602)
WRITE(ICHAN,603) ((QN(I,J),J=1,4),I=2,K1)
WRITE(ICHAN,604)
WRITE(ICHAN,600)
WRITE(ICHAN,601) (DN(1,J),J=1,4)
WRITE(ICHAN,602)
WRITE(ICHAN,603) ((DN(I,J),J=1,4),I=2,K1)
130 CONTINUE
RETURN
200 ICHAN=7
ICREAT=1HM
IF (ITORQ.AND.IDILA) ICREAT=1HD
IF (ITORQ.AND..NOT.IDILA) ICREAT=1HT
WRITE(ICHAN,610)
IMIST=1
WRITE(ICHAN,630) IMIST
WRITE(ICHAN,620) ADENT,GLOBAL,ICREAT,RUN
IMIST=8
WRITE(ICHAN,630) IMIST
FNP=FLOAT(NP)
WRITE(ICHAN,640) DK,HANG,HUP,PS,TEMP,TV,EICH,FNP
WRITE(ICHAN,610)
IMIST=2
WRITE(ICHAN,630) IMIST
WRITE(ICHAN,630) K1
WRITE(ICHAN,640) ((QN(I,J),J=1,4),(DN(I,J),J=1,4),I=1,K1)
RETURN
END