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

(Source file text) 

---

      SUBROUTINE CHIRMSD
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/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
C 
      DIMENSION V(10),P(10) 
C 
  150 FORMAT(1X'SUBR. CHIRMSD: J,DY(J),P(J)  'I2,2E11.4)
C     RMSD=ROOT MEAN SQUARE DEVIATION
      SP=0. 
      RMSD=0
      CHI=0
      Q=0
      QA=0
      N=MR
      DO 60 J=1,N 
      XL(J)=T(J)
      YL(J)=ETA0(J)
      DY(J)=DETA(I,J)
      P(J)=1. 
      IF (DY(J).EQ.0.) GO TO 72 
      P(J)=1./DY(J)**2
      WRITE(4,150) J,DY(J),P(J)
   72 SP=SP+P(J)
      A=14.69*AMU0
      B=A*TDM
      V(J)=YL(J)-XL(J)*A-B
      Q=Q+V(J)*V(J)*P(J)
      QA=QA+(V(J))**2 
   60 CONTINUE
      CHI=SQRT(Q/(N-2)) 
      RMSD=SQRT(QA/(N-2)) 
      RETURN
      END 

---