SpectraRust/synspec/extracted/griem.f
2026-03-19 14:05:33 +08:00

19 lines
591 B
Fortran

SUBROUTINE GRIEM(ID,T,ANE,ION,FR,WGR,GAM)
C =========================================
C
C STARK DAMPING PARAMETER (GAM) CALCULATED FROM INPUT VALUES
C OF STARK WIDTHS FOR T=5000, 10000, 20000, 40000 K,
C AND FOR NE=1.E16 (FOR NEUTRALS) OR NE = 1.E17 (FOR IONS)
C
INCLUDE 'PARAMS.FOR'
INCLUDE 'MODELP.FOR'
DIMENSION WGR(4)
if(t.le.0.) return
J=JT(ID)
GAM=(TI0(ID)*WGR(J)+TI1(ID)*WGR(J-1)+TI2(ID)*WGR(J-2))
* *ANE*1.E-10*FR*1.E-10*FR*4.2E-14
IF(ION.GT.1) GAM=GAM*0.1
IF(GAM.LT.0.) GAM=0.
RETURN
END