Line data Source code
1 : module eos_support
2 :
3 : use eos_def
4 : use eos_lib
5 : use const_def
6 : use chem_def
7 : use math_lib
8 : use utils_lib, only: is_bad_num
9 :
10 : implicit none
11 :
12 : logical, parameter :: use_shared_data_dir = .true. ! MUST BE .true. FOR RELEASE
13 : !logical, parameter :: use_shared_data_dir = .false.
14 :
15 : real(dp) :: X, Z, Zinit, Y, dXC, dXO, XC, XO, abar, zbar, z2bar, z53bar, ye
16 : integer, parameter :: species = 7
17 : integer, parameter :: h1 = 1, he4 = 2, c12 = 3, n14 = 4, o16 = 5, ne20 = 6, mg24 = 7
18 : integer, target :: chem_id_array(species)
19 : integer, pointer, dimension(:) :: chem_id, net_iso
20 : real(dp) :: xa(species)
21 :
22 : real(dp), allocatable :: d_dxa(:, :) ! (num_d_dxa_basic_results,species)
23 :
24 : integer :: handle
25 : type(EoS_General_Info), pointer :: rq
26 :
27 : character(len=eos_name_length) :: eos_names(num_eos_basic_results)
28 :
29 : ! if false, then test using data from mesa/eos/data/eos_data
30 : ! if true, then test using data from mesa/data/eos_data
31 :
32 : contains
33 :
34 36 : subroutine Init_Composition(X_in, Zinit_in, XC_in, XO_in)
35 : use chem_lib
36 : real(dp), intent(IN) :: X_in, Zinit_in, XC_in, XO_in
37 :
38 : real(dp), parameter :: Zfrac_C = 0.173312d0
39 : real(dp), parameter :: Zfrac_N = 0.053177d0
40 : real(dp), parameter :: Zfrac_O = 0.482398d0
41 : real(dp), parameter :: Zfrac_Ne = 0.098675d0
42 :
43 576 : real(dp) :: Z, frac, dabar_dx(species), dzbar_dx(species), sumx, &
44 324 : mass_correction, dmc_dx(species)
45 :
46 36 : chem_id => chem_id_array
47 :
48 36 : allocate (net_iso(num_chem_isos))
49 282852 : net_iso(:) = 0
50 :
51 36 : chem_id(h1) = ih1; net_iso(ih1) = h1
52 36 : chem_id(he4) = ihe4; net_iso(ihe4) = he4
53 36 : chem_id(c12) = ic12; net_iso(ic12) = c12
54 36 : chem_id(n14) = in14; net_iso(in14) = n14
55 36 : chem_id(o16) = io16; net_iso(io16) = o16
56 36 : chem_id(ne20) = ine20; net_iso(ine20) = ne20
57 36 : chem_id(mg24) = img24; net_iso(img24) = mg24
58 :
59 36 : X = X_in
60 36 : Zinit = Zinit_in
61 36 : XC = XC_in; XO = XO_in
62 36 : Y = 1 - (X + Zinit + XC + XO)
63 36 : if (Y < 0) then ! adjust XC and XO
64 0 : if (XC + XO <= 0) then
65 0 : write (*, *) 'bad args to Init_Composition'
66 0 : call mesa_error(__FILE__, __LINE__)
67 : end if
68 0 : frac = (1 - X - Zinit)/(XC + XO)
69 0 : if (frac <= 0) call mesa_error(__FILE__, __LINE__, 'bad args to Init_Composition')
70 0 : XC = frac*XC; XO = frac*XO
71 0 : Y = 1 - (X + Zinit + XC + XO)
72 0 : if (Y < -1d-10) then
73 0 : write (*, *) 'screw up in Init_Composition'
74 0 : call mesa_error(__FILE__, __LINE__)
75 : end if
76 0 : if (Y < 0) Y = 0
77 : end if
78 :
79 36 : xa(h1) = X
80 36 : xa(he4) = Y
81 36 : xa(c12) = Zinit*Zfrac_C + XC
82 36 : xa(n14) = Zinit*Zfrac_N
83 36 : xa(o16) = Zinit*Zfrac_O + XO
84 36 : xa(ne20) = Zinit*Zfrac_Ne
85 252 : xa(species) = 1 - sum(xa(1:species - 1))
86 :
87 : call composition_info( &
88 : species, chem_id, xa, X, Y, Z, abar, zbar, z2bar, z53bar, &
89 36 : ye, mass_correction, sumx, dabar_dx, dzbar_dx, dmc_dx)
90 :
91 36 : end subroutine Init_Composition
92 :
93 2 : subroutine Setup_eos
94 36 : use chem_lib
95 : use const_lib, only: const_init
96 : !..allocate and load the eos tables
97 :
98 : character(len=256) :: my_mesa_dir
99 : integer :: info
100 : logical :: use_cache
101 :
102 2 : info = 0
103 2 : allocate (d_dxa(num_eos_d_dxa_results, species), stat=info)
104 2 : if (info /= 0) then
105 0 : write (*, *) 'allocate failed for Setup_eos'
106 0 : call mesa_error(__FILE__, __LINE__)
107 : end if
108 :
109 2 : my_mesa_dir = '../..'
110 2 : call const_init(my_mesa_dir, info)
111 2 : if (info /= 0) then
112 0 : write (*, *) 'const_init failed'
113 0 : call mesa_error(__FILE__, __LINE__)
114 : end if
115 :
116 2 : call math_init()
117 :
118 2 : call chem_init('isotopes.data', info)
119 2 : if (info /= 0) then
120 0 : write (*, *) 'chem_init failed'
121 0 : call mesa_error(__FILE__, __LINE__)
122 : end if
123 :
124 2 : use_cache = .true.
125 :
126 2 : call eos_init(' ', use_cache, info)
127 2 : if (info /= 0) then
128 0 : write (*, *) 'failed in eos_init'
129 0 : call mesa_error(__FILE__, __LINE__)
130 : end if
131 54 : eos_names = eosDT_result_names
132 :
133 2 : handle = alloc_eos_handle_using_inlist('inlist', info)
134 2 : if (info /= 0) then
135 0 : write (*, *) 'failed in alloc_eos_handle_using_inlist'
136 0 : call mesa_error(__FILE__, __LINE__)
137 : end if
138 :
139 2 : call eos_ptr(handle, rq, info)
140 2 : if (info /= 0) then
141 0 : write (*, *) 'failed in eos_ptr'
142 0 : call mesa_error(__FILE__, __LINE__)
143 : end if
144 :
145 10 : end subroutine Setup_eos
146 :
147 0 : subroutine Shutdown_eos
148 0 : call free_eos_handle(handle)
149 0 : call eos_shutdown
150 0 : deallocate (d_dxa)
151 2 : end subroutine Shutdown_eos
152 :
153 : end module eos_support
|
