3 FUNCTION cooling_dalgarno (self, temp)
10 class(cooling_t):: self
11 real(kind=dp) cooling_dalgarno, temp, qq
12 real(kind=dp),
parameter:: &
38 if (temp .le. t_mc)
then 40 else if (temp .le. te2)
then 42 else if (temp .le. te3)
then 44 else if (temp .le. te4)
then 46 else if (temp .le. te5)
then 48 else if (temp .le. te6)
then 50 else if (temp .le. te7)
then 52 else if (temp .le. te8)
then 57 cooling_dalgarno=qq*4e-28