File: DPRINT.FT of Tape: Various/ETH/prog1

(Source file text) 

---

      SUBROUTINE PRINT(IFUN)
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 /LINE/ A,B,DA,DB,SA,SB,SAB,TD,AMU,RR,RMSD,CHI
      COMMON /REGRES/ XL(10),YL(10),DY(10)
      COMMON/FRERUN/ I,JR,INIT,LS
      COMMON/ACCUR/ DDMU1,DDMU2,DCHI,ERRMU
C 
      LOGICAL ERRMU 
  500 FORMAT(1H1) 
  505 FORMAT(1X,T50,'VERSION  ',I3) 
  510 FORMAT(1X,T10,'RUNS:',T20,A6,T30,'TO',T40,A6,//)
  520 FORMAT(1X,T20,'ANGLE= ',F6.1,T40,'HI= ',F7.3,' T',T60 
     $,'DK= ',1PE10.2,' 1/T',//)
  525 FORMAT(1X,T20'(AMPLITUDE CALCULATED FOR I0/NP='I4' /'I4')'//)
  530 FORMAT(1X,T20,'FREQUENCY : ',F10.3,' T',T50'MU*P='
     $,T60,3PF7.3,T70,'1E-3*1/T'T85'MU='T92,0PF7.3) 
  540 FORMAT(//2X'RUN'T16'T[K]'T32'FREQ[T]'T51'PHASE'T71'AMP'T91, 
     $'ETA[T]'T107'A0'/)
  545 FORMAT(//2X'RUN'T16'T[K]'T33'TD[K]'T53'AMP'T72'ETA[T]'
     $T89'A0'/)
  550 FORMAT(1X,A6,2X,0PF6.3'+/-'F5.3,2X,F9.2'+/-'F5.2,3X,F6.1, 
     $'+/-'F4.1,3X,1PE9.2'+/-'E9.2,3X,0PF6.2'+/-'F5.2,3X,1PE9.2)
  555 FORMAT(1X,A6,2X,0PF6.3'+/-'F5.3,3X,F6.3'+/-'F5.3,
     $3X,1PE9.2'+/-'E9.2,3X,0PF6.2'+/-'F5.2,3X,1PE9.2)
  560 FORMAT(1X,//,T8,'MEAN FREQUENCY[T]='F9.2'+/-'F5.2,T95,
     $'MEAN A0='1PE9.2,/,T8,'PERIOD[1/T]=        '1PE9.2, 
     $T95'VAR.A0= '1PE9.2,////) 
  565 FORMAT(1X,//,T16,'MEAN TD ='2X,F6.3' +/-'F5.3,T77
     $'MEAN A0='1PE9.2,/,T77'VAR.A0 = 'E9.2,////) 
  570 FORMAT(1X,T6,'LOGAMP FIT:',T20,'MU =',T30,F7.3,T70, 
     $'MU*P =',T80,3PF7.3,T90'1E-3*1/T',/,T6,'...........',/,T33,
     $'+/-',T38,0PF7.3,'  (LINE FIT)'
     $,/,T38,F7.3,'  (VARIANCE)',//) 
  580 FORMAT(1X,T20'CHI='T30,F7.3,/,T20'RMSD='T30,F7.3,/,T20 
     $'RR=',T30,F7.3,//)
  585 FORMAT(1X,T20'CHI='T30,F7.3,/,T20'RMSD='T30,F7.3,//)
  590 FORMAT(1X,T6,'LOGAMP+ FIT:',T20,'MU =',T30,F7.3,T70, 
     $'MU*P =',T80,3PF7.3,T90'1E-3*1/T',//,T33,
     $'+/-',T38,0PF7.3,T83,'+/-',T88,F7.3,'  (LINE FIT)'
     $,/,T38,F7.3,T88,F7.3,'  (VARIANCE)',,//) 
  600 FORMAT(1X,T6,'LOGAMP0 FIT:',T20,'MU =',T30,F7.3,T70, 
     $'MU*P =',T80,3PF7.3,T90'1E-3*1/T',/,T6,'************',/,
     $T33'+/-'T38,0PF7.3,T83,'+/-',T88,F7.3,'  (LINE FIT)'
     $,/,T38,F7.3,T88,F7.3,'  (VARIANCE)',,//) 
C 
C       !!!!!!!!!!!!!!!!!!!!!!!!! 
      IVERS=8 
C       !!!!!!!!!!!!!!!!!!!!!!!!! 
      IF (IFUN.EQ.2) GO TO 2
      IF (IFUN.EQ.3) GO TO 3
      IF (IFUN.EQ.4) GO TO 4
      IF (IFUN.EQ.5) GO TO 5
      WRITE(4,500)
      WRITE(4,505) IVERS
      WRITE(4,510) ADENT(1),ADENT(MR) 
      WRITE(4,520) HANG,HI,DK 
      WRITE(4,525) I0,NP
      RETURN
    2 WRITE(4,500)
      AMP=AMU0/ANFM(I)
      WRITE(4,530) FREQ(I),AMP,AMU0
      WRITE(4,540)
      WRITE(4,550) (ADENT(J),T(J),DT(J),ANF(I,J),DANF(I,J), 
     1PH(I,J),DPH(I,J),AF0(J),DAF(I,J),ETA0(J),DETA(I,J),A0(J), 
     2J=1,MR) 
      PER=1./ANFM(I)
      WRITE(4,560) ANFM(I),DANFM(I),A0M,PER,DA0M
      RETURN
    5 DA0M=.0 
      DO 7 J=1,MR 
      A0(J)=A0M 
      TAU=ETA0(J)*DK*NP 
      EX=EXP(-TAU)
      AF0(J)=A0M*T(J)*(1.-EX)/(TAU*SQHI*EXP(ETA0(J)/HI))
    7 CONTINUE
      RETURN
    3 GOTO 4
    4 DMUL=DA
      DMUS=SQRT(SA)
      AMP=AMU/ANFM(I)
      IF (IFUN.EQ.3) WRITE(4,570) AMU,AMP,DMUL,DMUS
      IF (IFUN.EQ.4) WRITE(4,590) AMU,AMP,DMUL,DMUS
      IF (IFUN.EQ.5) WRITE(4,600) AMU,AMP,DMUL,DMUS
      WRITE(4,580) CHI,RMSD,RR
C     WRITE(4,530) FREQ(I),AMP,AMU0
      WRITE(4,545)
      WRITE(4,555) (ADENT(J),T(J),DT(J),TDJ(J),DTDJ(J),AF0(J),
     $DAF(I,J),ETA0(J),DETA(I,J),A0(J),J=1,MR)
      WRITE(4,565) TDM,DTDM,A0M,DA0M
      WRITE(4,585) CHI,RMSD
      RETURN
      END 

---