C  FIT TO ORG 40 PAR FOR 4.2/ CGW
      SUBROUTINE HD42
      COMMON NQI,QT(6),TWL(6),D(10),DL(10)
      COMMON TIM,DXT,HT,HML,SXS,SYS,SZS,TIVCH,UT,BR,SF,TMP,ALFA,SY100,
     1BETA,SZ100,Z,RC,V,QS
      COMMON HDM(33),IPR(1),ND,IPO,I2MC,IMA,OPC,IMM,IDD,IHS,NOV,INP,MRL,
     1NMU,ID2,IDEP,IMTCH,IM,IR,IL,IRL,ISM,IVD
      DIMENSION FWT(10),FST(3),ALBT4(3),SZR4(3)
      DATA FWT/1.6,.89,.64,.52,.43,.37,.33,.3,.27,.25/
      DATA FST/.7,1.,1.25/
      DATA SYR4/6.82/
      DATA ALBT4/3.33,1.4,1.02/
      DATA SZR4/287.,12.84,6.97/
      TF=1.8*TMP+32.
      IF (TF.LT.50.) GO TO 4
      IUT=UT+.5
      IF (IUT.GT.10) IUT=10
      IMT=(IM+1)/2
      IF (TF.LE.80.) FT=EXP(3.6889-.046052*TF)
      IF (TF.GT.80.) FT=EXP(3.4239-.042799*TF)
C     22. TIME AFTER FUNCTIONING
 1    CALL READA (22,IRT,3,TIM)
      IF (IRT.LT.0) GO TO 1
      FEL=ALOG(FST(IMT)*FWT(IUT)*FT*120./TIM)
      R=1.
      IF (FEL.LT.-1.2) GO TO 3
      R=EXP(.36464-.86189*FEL)
      IF (FEL.GT..4) R=R-EXP(-.05129-1.6767*FEL)
      IF (FEL.LE..4) R=R-EXP(-.24846-1.1373*FEL)
 3    QT(1)=QS*R*NMU
      SY100=SYR4
      SZ100=SZR4(IMT)
      ALFA=ALBT4(IMT)
      BETA=ALBT4(IMT)
      SYS=3.8
      SZS=.2
      SXS=SYS
      RETURN
 4    WRITE(*,5)
 5    FORMAT(' TEMPERATURE FUNCTION NOT DEFINED BELOW 10 DEG C')
      QT(1)=0.
      RETURN
      END
