23 lines
616 B
Fortran
23 lines
616 B
Fortran
FUNCTION SBFHMI_old(FR)
|
|
C ===================
|
|
C
|
|
C Bound-free cross-section for H- (negative hydrogen ion)
|
|
C
|
|
INCLUDE 'IMPLIC.FOR'
|
|
PARAMETER (UN=1.D0)
|
|
SBFHMI_old=0.
|
|
FR0=1.8259D14
|
|
IF(FR.LT.FR0) RETURN
|
|
IF(FR.LT.2.111D14) GO TO 10
|
|
X=2.997925D15/FR
|
|
SBFHMI=(6.80133D-3+X*(1.78708D-1+X*(1.6479D-1+X*(-2.04842D-2+X*
|
|
* 5.95244D-4))))*1.D-17
|
|
SBFHMI_old=sbfhmi
|
|
RETURN
|
|
10 X=2.997925D15*(UN/FR0-UN/FR)
|
|
SBFHMI=(2.69818D-1+X*(2.2019D-1+X*(-4.11288D-2+X*2.73236D-3)))
|
|
1 *X*1.D-17
|
|
SBFHMI_old=sbfhmi
|
|
RETURN
|
|
END
|