34 lines
870 B
Fortran
34 lines
870 B
Fortran
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
|