24 lines
604 B
Fortran
24 lines
604 B
Fortran
FUNCTION EPS(T,ANE,ALAM,ION,N)
|
|
C ==============================
|
|
C
|
|
C NLTE PARAMETER EPSILON (COLLISIONAL/SPONTANEOUS DEEXCITATION)
|
|
C AFTER KASTNER, 1981, J.Q.S.R.T. 26, 377
|
|
C
|
|
INCLUDE 'PARAMS.FOR'
|
|
DATA CK0,CK1 /7.75E-8, 2.58E-8/
|
|
X=1.438E8/ALAM/T
|
|
XKT=12390./ALAM
|
|
TT=0.75*X
|
|
T1=TT+1.
|
|
A=4.36E7*XKT*XKT/(1.-EXP(-X))
|
|
IF(ION.EQ.1) GO TO 10
|
|
B=1.1+LOG(T1/TT)-0.4/T1/T1
|
|
C=X*B*SQRT(T)/XKT/XKT*ANE
|
|
IF(N.EQ.0) C=CK0*C
|
|
IF(N.NE.0) C=CK1*C
|
|
GO TO 20
|
|
10 C=2.16/T/SQRT(T)/X**1.68*ANE
|
|
20 EPS=C/(C+A)
|
|
RETURN
|
|
END
|