SUBROUTINE INTLEM(PRFH,WL0,ILINE,ID)
C     ====================================
C
      INCLUDE 'IMPLIC.FOR'
      INCLUDE 'BASICS.FOR'
      INCLUDE 'MODELQ.FOR'
      PARAMETER (FOC1=1.25e-9,TTW=2./3.,VTBC=6.06e-9)
      DIMENSION PRFH(MHWL)
C
C     temperature is modified in order to account for the
C     effect of turbulent velocity on the Doppler width
C
      T=TEMP(ID)+VTBC*VTURBS(ID)*VTURBS(ID)
      ANE=ELEC(ID)
      TL=LOG10(T)
      ANEL=LOG10(ANE)
      F00=FOC1*EXP(TTW*LOG(ANE))
      XK=XK0(ILINE)
      FXK=F00*XK
      DOP=1.E8/WL0*SQRT(1.65E8*T)
      DBETA=WL0*WL0/2.997925E18/FXK
      BETAD=DBETA*DOP
C
C     interpolation to the actual values of temperature and electron
C     density
C
      NWL=NWLHYD(ILINE)
      DO IWL=1,NWL
         CALL INTHYD(PRFH0,TL,ANEL,IWL,ILINE)
         PRFH(IWL)=PRFH0
      END DO
      RETURN
      END