-
Notifications
You must be signed in to change notification settings - Fork 1
/
abundances.f90
49 lines (40 loc) · 1.17 KB
/
abundances.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
!>
!! \brief This module sets elemental abundances
!!
!! Module for Capreole / C2-Ray (f90)
!!
!! \b Author: Garrelt Mellema
!!
!! \b Date: 2008-05-27
!!
!! This module is also accepted by the F compiler (Dec 9, 2003)
!<
module abundances
! This module sets elemental abundances
use precision, only: dp
implicit none
private
!> Helium abundance (by number)
real(kind=dp),public,parameter :: abu_he=0.074
!real(kind=dp),parameter :: abu_he=0.08 ! for EoR KP sims
!> Carbon abundance (by number)
real(kind=dp),public,parameter :: abu_c=7.1e-7
!> Mean molecular weight
real(kind=dp),public,parameter :: mu=(1.0-abu_he)+4.0*abu_he
!-----------------------------------------------------------------------
! set heavy element abundances
! from Osterbrock (1989), Table 5.13
!-----------------------------------------------------------------------
!parameter(abu_he=0.1)
!abu_he=0.074 !cosmic abundance of He (24.2% by mass)
!parameter(abu_c=7.1e-7)
!abu_n=1.3e-7
!abu_o=4.2e-7
!abu_ne=1.0e-7
!abu_s=1.0e-8
!abu_c=3.3e-4
! abu_n=0.0
! abu_o=6.6e-4
! abu_ne=8.3e-5
! abu_s=1.6e-5
end module abundances