22 lines
512 B
Fortran
22 lines
512 B
Fortran
FUNCTION GFREE1(ID,X)
|
|
C =====================
|
|
C
|
|
C Hydrogenic free-free Gaunt factor, for depth ID,
|
|
C frequency FR, and charge CH
|
|
C
|
|
INCLUDE 'IMPLIC.FOR'
|
|
INCLUDE 'BASICS.FOR'
|
|
INCLUDE 'MODELQ.FOR'
|
|
PARAMETER (XMIN=0.2D0,
|
|
* XMINI=UN/XMIN)
|
|
C
|
|
IF(X.LT.UN) THEN
|
|
GFREE1=((GF4(ID)*X-GF3(ID))*X+GF2(ID))*X+GF1(ID)
|
|
ELSE IF(X.LT.XMINI) THEN
|
|
GFREE1=GF0(ID)+GF5(ID)*X
|
|
ELSE
|
|
GFREE1=GF6(ID)
|
|
END IF
|
|
RETURN
|
|
END
|