annotate Common/Src/decom.c @ 863:0c89c6fa949c Evo_2_23

Bugfix empty line in deco plan (VPM only): Floating numbers were used to calculate the target slot for the time entry of a deco stop. The float rounding caused a time to be written into one line above the intended one. In the next step the misplaced time was overwritten by the next shallower stop. To fix the problem the index calculation has been corrected and in addition digit numbers have generally been added to floating point operations to make the floating operation more visible.
author Ideenmodellierer
date Tue, 02 Jul 2024 20:05:08 +0200
parents 012f94ec2fe0
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1 ///////////////////////////////////////////////////////////////////////////////
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
2 /// -*- coding: UTF-8 -*-
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
3 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
4 /// \file Common/Src/decom.c
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
5 /// \brief This code is used to calculate desat, calculated by RTE and send to Firmware
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
6 /// \author heinrichs weikamp gmbh
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
7 /// \date 22-Feb-2016
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
8 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
9 /// $Id$
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
10 ///////////////////////////////////////////////////////////////////////////////
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
11 /// \par Copyright (c) 2014-2018 Heinrichs Weikamp gmbh
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
12 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
13 /// This program is free software: you can redistribute it and/or modify
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
14 /// it under the terms of the GNU General Public License as published by
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
15 /// the Free Software Foundation, either version 3 of the License, or
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
16 /// (at your option) any later version.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
17 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
18 /// This program is distributed in the hope that it will be useful,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
19 /// but WITHOUT ANY WARRANTY; without even the implied warranty of
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
20 /// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
21 /// GNU General Public License for more details.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
22 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
23 /// You should have received a copy of the GNU General Public License
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
24 /// along with this program. If not, see <http://www.gnu.org/licenses/>.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
25 //////////////////////////////////////////////////////////////////////////////
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
26 /**
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
27 @verbatim
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
28 ==============================================================================
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
29 ##### Changes #####
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
30 ==============================================================================
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
31 V1.0.2 1602220x decom_oxygen_calculate_cns() changed to hwOS version
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
32
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
33 @endverbatim
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
34 ******************************************************************************
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
35 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
36
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
37 #include "decom.h"
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
39 #include <math.h>
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
40 #include <string.h>
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
41 #include "settings.h"
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
42 #include "calc_crush.h"
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
43
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
44 #define FRACTION_N2_AIR 0.7902
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
45
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
46 const float helium_time_constant[16] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
47 3.68695308808482E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
48 2.29518933960247E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
49 1.46853216220327E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
50 9.91626867753856E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
51 6.78890480470074E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
52 4.78692804254106E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
53 3.37626488338989E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
54 2.38113081607676E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
55 1.68239606932026E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
56 1.25592893741610E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
57 9.80544886914621E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
58 7.67264977374303E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
59 6.01220557342307E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
60 4.70185307665137E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
61 3.68225234041620E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
62 2.88775228329769E-003};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
63
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
64 const float nitrogen_time_constant[16] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
65 1.38629436111989E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
66 8.66433975699932E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
67 5.54517744447956E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
68 3.74674151654024E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
69 2.56721177985165E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
70 1.80978376125312E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
71 1.27651414467762E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
72 9.00191143584345E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
73 6.35914844550409E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
74 4.74758342849278E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
75 3.70666941475907E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
76 2.90019740820061E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
77 2.27261370675392E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
78 1.77730046297422E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
79 1.39186180835330E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
80 1.09157036308653E-003};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
81
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
82
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
83 const float buehlmann_N2_a[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
84 1.1696f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
85 1.0000f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
86 0.8618f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
87 0.7562f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
88 0.6200f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
89 0.5043f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
90 0.4410f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
91 0.4000f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
92 0.3750f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
93 0.3500f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
94 0.3295f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
95 0.3065f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
96 0.2835f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
97 0.2610f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
98 0.2480f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
99 0.2327f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
100
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
101 const float buehlmann_N2_b[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
102 0.5578f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
103 0.6514f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
104 0.7222f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
105 0.7825f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
106 0.8126f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
107 0.8434f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
108 0.8693f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
109 0.8910f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
110 0.9092f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
111 0.9222f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
112 0.9319f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
113 0.9403f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
114 0.9477f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
115 0.9544f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
116 0.9602f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
117 0.9653f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
118
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
119 const float buehlmann_He_a[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
120 1.6189f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
121 1.3830f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
122 1.1919f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
123 1.0458f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
124 0.9220f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
125 0.8205f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
126 0.7305f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
127 0.6502f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
128 0.5950f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
129 0.5545f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
130 0.5333f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
131 0.5189f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
132 0.5181f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
133 0.5176f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
134 0.5172f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
135 0.5119f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
136
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
137 const float buehlmann_He_b[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
138 0.4770f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
139 0.5747f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
140 0.6527f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
141 0.7223f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
142 0.7582f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
143 0.7957f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
144 0.8279f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
145 0.8553f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
146 0.8757f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
147 0.8903f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
148 0.8997f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
149 0.9073f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
150 0.9122f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
151 0.9171f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
152 0.9217f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
153 0.9267f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
154
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
155 const float buehlmann_N2_t_halflife[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
156 5.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
157 8.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
158 12.5f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
159 18.5f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
160 27.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
161 38.3f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
162 54.3f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
163 77.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
164 109.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
165 146.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
166 187.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
167 239.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
168 305.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
169 390.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
170 498.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
171 635.0f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
172
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
173 const float buehlmann_He_t_halflife[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
174 1.88f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
175 3.02f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
176 4.72f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
177 6.99f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
178 10.21f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
179 14.48f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
180 20.53f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
181 29.11f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
182 41.20f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
183 55.19f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
184 70.69f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
185 90.34f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
186 115.29f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
187 147.42f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
188 188.24f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
189 240.03f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
190
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
191 const float float_buehlmann_N2_factor_expositon_one_second[] = { 2.30782347297664E-003f, 1.44301447809736E-003f, 9.23769302935806E-004f, 6.24261986779007E-004f, 4.27777107246730E-004f, 3.01585140931371E-004f, 2.12729727268379E-004f, 1.50020603047807E-004f, 1.05980191127841E-004f, 7.91232600646508E-005f, 6.17759153688224E-005f, 4.83354552742732E-005f, 3.78761777920511E-005f, 2.96212356654113E-005f, 2.31974277413727E-005f, 1.81926738960225E-005f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
192 const float float_buehlmann_N2_factor_expositon_003_second[] = { 6.90750456296407E-003f, 4.32279956671600E-003f, 2.76874864793053E-003f, 1.87161709452954E-003f, 1.28278242026003E-003f, 9.04482589432765E-004f, 6.38053429621421E-004f, 4.49994293975742E-004f, 3.17906879170993E-004f, 2.37350999218289E-004f, 1.85316297551252E-004f, 1.44999356986975E-004f, 1.13624229615916E-004f, 8.88610747694640E-005f, 6.95906688746861E-005f, 5.45770287740943E-005f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
193 const float float_buehlmann_N2_factor_expositon_008_second[] = { 1.83141447532454E-002f, 1.14859796471039E-002f, 7.36630472495203E-003f, 4.98319782231915E-003f, 3.41709742823104E-003f, 2.41013596224415E-003f, 1.70057124687550E-003f, 1.19953484034729E-003f, 8.47527105247492E-004f, 6.32810814525819E-004f, 4.94100480767923E-004f, 3.86618231662861E-004f, 3.02969256443353E-004f, 2.36945319086024E-004f, 1.85564355251966E-004f, 1.45532124251058E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
194 const float float_buehlmann_N2_factor_expositon_10_seconds[] = { 2.28400315657541E-002f, 1.43368013598124E-002f, 9.19938673477072E-003f, 6.22511239287027E-003f, 4.69545762670800E-003f, 3.01176178733265E-003f, 2.12526200031782E-003f, 1.49919365737827E-003f, 1.05929662305226E-03f, 7.909509380171760E-004f, 6.17587450108648E-004f, 4.83249432061905E-004f, 3.78697227222391E-004f, 2.61728759809380E-004f, 2.31950063482533E-004f, 1.81911845881011E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
195 const float float_buehlmann_N2_factor_expositon_18_seconds[] = { 4.07358806747357E-002f, 2.56581087982929E-002f, 1.64979259737517E-002f, 1.11772892486697E-002f, 7.67205373705648E-003f, 5.41463899418337E-003f, 3.82221908774349E-003f, 2.69693016270112E-003f, 1.90592594569927E-003f, 1.42326123023573E-003f, 1.11138278062062E-003f, 8.69680830683950E-004f, 6.81551750048359E-004f, 5.33048018290350E-004f, 4.17471377070378E-004f, 3.27417496114757E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
196 const float float_buehlmann_N2_factor_expositon_20_seconds[] = { 4.51583960895835E-002f, 2.84680588463941E-002f, 1.83141447532454E-002f, 1.24114727614367E-002f, 8.52086250432193E-003f, 6.01445286560154E-003f, 4.24600726206570E-003f, 2.99613973313428E-003f, 2.11747113676897E-003f, 1.58127627264804E-003f, 1.23479348595879E-003f, 9.66265334110261E-004f, 7.57251042854845E-004f, 5.92258033589421E-004f, 4.63846326133055E-004f, 3.63790599842373E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
197 const float float_buehlmann_N2_factor_expositon_one_minute[] = { 1.29449436703876E-001f, 8.29959567953288E-002f, 5.39423532744041E-002f, 3.67741962370398E-002f, 2.53453908775689E-002f, 1.79350552316596E-002f, 1.26840126026602E-002f, 8.96151553540825E-003f, 6.33897185233323E-003f, 4.73633146787078E-003f, 3.69980819572546E-003f, 2.89599589841472E-003f, 2.27003327536857E-003f, 1.77572199977927E-003f, 1.39089361795441E-003f, 1.09097481687104E-003f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
198 const float float_buehlmann_N2_factor_expositon_100_second[] = { 2.06299474015900E-001f, 1.34463438993857E-001f, 8.82775114417832E-002f, 6.05359181023788E-002f, 4.18844218114071E-002f, 2.97126970072147E-002f, 2.10505144045823E-002f, 1.48911986890571E-002f, 1.05426136839346E-002f, 7.88141652426455E-003f, 6.15873909572406E-003f, 4.82199900095137E-003f, 3.78052526350936E-003f, 2.95778454900952E-003f, 2.31708109427220E-003f, 1.81763004457269E-003f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
199 const float float_buehlmann_N2_factor_expositon_five_minutes[]= { 5.00000000000000E-001f, 3.51580222674495E-001f, 2.42141716744801E-001f, 1.70835801932547E-001f, 1.20463829104624E-001f, 8.65157896183918E-002f, 6.18314987350977E-002f, 4.40116547625051E-002f, 3.12955727186929E-002f, 2.34583889613009E-002f, 1.83626606868127E-002f, 1.43963540993090E-002f, 1.12987527093947E-002f, 8.84713405486026E-003f, 6.93514912851934E-003f, 5.44298480182925E-003f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
200 const float float_buehlmann_N2_factor_expositon_800_second[] = { 8.42509868763141E-001f, 6.85019737526282E-001f, 5.22579198044792E-001f, 3.93205767018569E-001f, 2.89861248917861E-001f, 2.14397627137602E-001f, 1.56505490290652E-001f, 1.13102166881646E-001f, 8.12935637814599E-002f, 6.13392112527207E-002f, 4.82208523469105E-002f, 3.79311861210304E-002f, 2.98470272862601E-002f, 2.34187624071612E-002f, 1.83870151711824E-002f, 1.44488700649190E-002f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
201 const float float_buehlmann_N2_factor_expositon_one_hour[]= { 9.99755859375000E-001f, 9.94475728271980E-001f, 9.64103176406343E-001f, 8.94394508891055E-001f, 7.85689004286732E-001f, 6.62392147498621E-001f, 5.35088626789486E-001f, 4.17318576947576E-001f, 3.17197008420226E-001f, 2.47876700002107E-001f, 1.99405069752929E-001f, 1.59713055172538E-001f, 1.27468761759271E-001f, 1.01149026804458E-001f, 8.01196838116008E-002f, 6.33955413542552E-002f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
202
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
203 const float float_buehlmann_He_factor_expositon_one_second[] = { 6.12608039419837E-003f, 3.81800836683133E-003f, 2.44456078654209E-003f, 1.65134647076792E-003f, 1.13084424730725E-003f, 7.97503165599123E-004f, 5.62552521860549E-004f, 3.96776399429366E-004f, 2.80360036664540E-004f, 2.09299583354805E-004f, 1.63410794820518E-004f, 1.27869320250551E-004f, 1.00198406028040E-004f, 7.83611475491108E-005f, 6.13689891868496E-005f, 4.81280465299827E-005f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
204 const float float_buehlmann_He_factor_expositon_003_second[] = { 1.82658845044263E-002f, 1.14103491926518E-002f, 7.31576933570466E-003f, 4.94586307993539E-003f, 3.38869776192019E-003f, 2.39060197012086E-003f, 1.68670834759044E-003f, 1.18985696621965E-003f, 8.40844326779777E-004f, 6.27767340286467E-004f, 4.90152279561396E-004f, 3.83558911153159E-004f, 3.00565099928485E-004f, 2.35065021719993E-004f, 1.84095669333084E-004f, 1.44377190774980E-004f}; // 3 He
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
205 const float float_buehlmann_He_factor_expositon_008_second[] = { 4.79706116082057E-002f, 3.01390075707096E-002f, 1.93899772993034E-002f, 1.31346689569831E-002f, 9.01102820363486E-003f, 6.36224538449637E-003f, 4.49156910795023E-003f, 3.16980660943422E-003f, 2.24068067793926E-003f, 1.67317060331207E-003f, 1.30653891641375E-003f, 1.02249686330114E-003f, 8.01306192375617E-004f, 6.26717274191169E-004f, 4.90846474157092E-004f, 3.84959521834594E-004f}; // 8 He
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
206 const float float_buehlmann_He_factor_expositon_10_seconds[] = { 5.95993001714799E-002f, 3.75307444923134E-002f, 2.41784389107607E-002f, 1.63912909924208E-002f, 1.25106927410620E-002f, 7.94647192918641E-003f, 5.61130562069978E-003f, 3.96068706690245E-003f, 2.80006593100546E-003f, 2.09102564918129E-003f, 1.63290683272987E-003f, 1.27795767799976E-003f, 1.00153239354972E-003f, 7.33352120986130E-004f, 6.13520442722559E-004f, 4.81176244777948E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
207 const float float_buehlmann_He_factor_expositon_18_seconds[] = { 1.04710896899039E-001f, 6.65386126706349E-002f, 4.30995968284519E-002f, 2.93106657684409E-002f, 2.01607137751910E-002f, 1.42581599093282E-002f, 1.00776711616688E-002f, 7.11793906429403E-003f, 5.03447255531631E-003f, 3.76069760984632E-003f, 2.93731229281968E-003f, 2.29914783358365E-003f, 1.80203605181650E-003f, 1.40956155658090E-003f, 1.10406577253352E-003f, 8.65950533235460E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
208 const float float_buehlmann_He_factor_expositon_20_seconds[] = { 1.15646523762030E-001f, 7.36529322024796E-002f, 4.77722809133601E-002f, 3.25139075644434E-002f, 2.23755519884017E-002f, 1.58297974422514E-002f, 1.11911244906306E-002f, 7.90568709176287E-003f, 5.59229149279306E-003f, 4.17767891009702E-003f, 3.26314728073529E-003f, 2.55428218017273E-003f, 2.00206171996409E-003f, 1.56605681014277E-003f, 1.22666447811148E-003f, 9.62120958977297E-004f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
209 const float float_buehlmann_He_factor_expositon_one_minute[] = { 3.08363886219441E-001f, 2.05084082411030E-001f, 1.36579295730211E-001f, 9.44046323514587E-002f, 6.56358626478964E-002f, 4.67416115355790E-002f, 3.31990512604121E-002f, 2.35300557146709E-002f, 1.66832281977395E-002f, 1.24807506400979E-002f, 9.75753219809561E-003f, 7.64329013320042E-003f, 5.99416843126677E-003f, 4.69081666943783E-003f, 3.67548116287808E-003f, 2.88358673732592E-003f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
210 const float float_buehlmann_He_factor_expositon_100_second[] = { 4.59084487437744E-001f, 3.17867635141657E-001f, 2.17103957783539E-001f, 1.52336166567559E-001f, 1.06981885584572E-001f, 7.66825160768219E-002f, 5.47171474343117E-002f, 3.89083581201959E-002f, 2.76504642556165E-002f, 2.07145921483078E-002f, 1.62096019995457E-002f, 1.27063337640768E-002f, 9.97030625587825E-003f, 7.80579708939710E-003f, 6.11829377951190E-003f, 4.80135692933603E-003f}; // 100 He
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
211 const float float_buehlmann_He_factor_expositon_five_minutes[]= { 8.41733751018722E-001f, 6.82600697933713E-001f, 5.20142493735619E-001f, 3.90924736715930E-001f, 2.87834706153524E-001f, 2.12857832580192E-001f, 1.55333364924147E-001f, 1.12242395185686E-001f, 8.06788883581406E-002f, 6.08653819753062E-002f, 4.78448115000141E-002f, 3.76366999883051E-002f, 2.96136888654287E-002f, 2.32350754744602E-002f, 1.82428098114835E-002f, 1.43350223887367E-002f}; // thre
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
212 const float float_buehlmann_He_factor_expositon_800_second[] = { 9.92671155759686E-001f, 9.53124140216102E-001f, 8.58865632718416E-001f, 7.33443528431762E-001f, 5.95533881446524E-001f, 4.71787742036413E-001f, 3.62479376011699E-001f, 2.72021750877104E-001f, 2.00940186773410E-001f, 1.54187175639359E-001f, 1.22553521140786E-001f, 9.72431193565182E-002f, 7.70338702477497E-002f, 6.07666995543268E-002f, 4.79109397391700E-002f, 3.77715319879068E-002f}; // 800 He
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
213 const float float_buehlmann_He_factor_expositon_one_hour[]= { 9.99999999753021E-001f, 9.99998954626205E-001f, 9.99850944669188E-001f, 9.97393537149572E-001f, 9.82979603888650E-001f, 9.43423231328217E-001f, 8.68106292901111E-001f, 7.60374619482322E-001f, 6.35576141220644E-001f, 5.29310840978539E-001f, 4.44744511849213E-001f, 3.68942936079581E-001f, 3.02834419265355E-001f, 2.45810174422126E-001f, 1.98231319020275E-001f, 1.59085372294989E-001f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
214
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
215 void decom_get_inert_gases(const float ambient_pressure_bar,const SGas* pGas, float* fraction_nitrogen, float* fraction_helium )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
216 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
217 float fraction_all_inertgases;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
218 float ppo2_fraction_setpoint;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
219 float diluent_divisor;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
220
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
221 *fraction_nitrogen = ((float)pGas->nitrogen_percentage) / 100.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
222 *fraction_helium = ((float)pGas->helium_percentage) / 100.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
223
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
224 if(pGas->AppliedDiveMode == DIVEMODE_CCR)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
225 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
226 // continue with CCR
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
227 fraction_all_inertgases = *fraction_nitrogen + *fraction_helium;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
228
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
229 ppo2_fraction_setpoint = (float)pGas->setPoint_cbar/ (100 * ambient_pressure_bar);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
230
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
231 diluent_divisor = (1.0f - ppo2_fraction_setpoint) / fraction_all_inertgases;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
232 if(diluent_divisor < 0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
233 diluent_divisor = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
234
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
235 *fraction_nitrogen *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
236 *fraction_helium *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
237 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
238 if(pGas->AppliedDiveMode == DIVEMODE_PSCR)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
239 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
240 fraction_all_inertgases = *fraction_nitrogen + *fraction_helium;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
241 ppo2_fraction_setpoint = decom_calc_SimppO2(ambient_pressure_bar, pGas) / ambient_pressure_bar;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
242 diluent_divisor = (1.0f - ppo2_fraction_setpoint) / fraction_all_inertgases;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
243 if(diluent_divisor < 0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
244 diluent_divisor = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
245
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
246 *fraction_nitrogen *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
247 *fraction_helium *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
248 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
249 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
250
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
251
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
252 void decom_tissues_exposure(int period_in_seconds, SLifeData * pLifeData)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
253 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
254 decom_tissues_exposure2(period_in_seconds, &pLifeData->actualGas, pLifeData->pressure_ambient_bar, pLifeData->tissue_nitrogen_bar, pLifeData->tissue_helium_bar);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
255 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
256
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
257
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
258 void decom_tissues_exposure2(int period_in_seconds, SGas* pActualGas, float ambiant_pressure_bar, float *tissue_N2_selected_stage, float *tissue_He_selected_stage)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
259 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
260 int ci;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
261 float percent_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
262 float percent_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
263 float partial_pressure_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
264 float partial_pressure_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
265
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
266
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
267
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
268 int period_in_seconds_left;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
269
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
270 if(period_in_seconds > 0)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
271 {
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
272
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
273 decom_get_inert_gases(ambiant_pressure_bar, pActualGas, &percent_N2, &percent_He);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
274
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
275 partial_pressure_N2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
276 partial_pressure_He = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_He;
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
277 period_in_seconds_left = period_in_seconds;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
278
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
279 while(period_in_seconds_left)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
280 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
281 if(period_in_seconds_left >= 3600)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
282 period_in_seconds = 3600;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
283 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
284 if(period_in_seconds_left >= 800)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
285 period_in_seconds = 800;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
286 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
287 if(period_in_seconds_left >= 300)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
288 period_in_seconds = 300;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
289 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
290 if(period_in_seconds_left >= 100)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
291 period_in_seconds = 100;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
292 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
293 if(period_in_seconds_left >= 60)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
294 period_in_seconds = 60;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
295 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
296 if(period_in_seconds_left == 36)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
297 period_in_seconds = 18;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
298 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
299 if(period_in_seconds_left >= 20)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
300 period_in_seconds = 20;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
301 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
302 if(period_in_seconds_left >= 18)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
303 period_in_seconds = 18;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
304 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
305 if(period_in_seconds_left >= 10)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
306 period_in_seconds = 10;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
307 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
308 if(period_in_seconds_left >= 8)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
309 period_in_seconds = 8;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
310 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
311 if(period_in_seconds_left >= 3)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
312 period_in_seconds = 3;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
313 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
314 period_in_seconds = 1;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
315
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
316 period_in_seconds_left -= period_in_seconds;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
317
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
318 switch (period_in_seconds)
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
319 {
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
320 case 1:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
321 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
322 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
323 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_one_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
324 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_one_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
325 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
326 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
327 case 3:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
328 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
329 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
330 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_003_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
331 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_003_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
332 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
333 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
334 case 8:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
335 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
336 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
337 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_008_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
338 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_008_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
339 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
340 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
341 case 10:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
342 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
343 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
344 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_10_seconds[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
345 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_10_seconds[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
346 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
347 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
348 case 18:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
349 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
350 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
351 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_18_seconds[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
352 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_18_seconds[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
353 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
354 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
355 case 20:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
356 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
357 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
358 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_20_seconds[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
359 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_20_seconds[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
360 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
361 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
362 case 60:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
363 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
364 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
365 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_one_minute[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
366 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_one_minute[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
367 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
368 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
369 case 100:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
370 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
371 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
372 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_100_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
373 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_100_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
374 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
375 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
376 case 300:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
377 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
378 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
379 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_five_minutes[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
380 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_five_minutes[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
381 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
382 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
383 case 800:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
384 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
385 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
386 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_800_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
387 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_800_second[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
388 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
389 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
390 case 3600:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
391 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
392 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
393 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_one_hour[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
394 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_one_hour[ci];
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
395 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
396 break;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
397 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
398 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
399 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
400 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
401
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
402 void decom_reset_with_1000mbar(SLifeData * pLifeData)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
403 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
404 double saturation = 1.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
405
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
406 saturation -= WATER_VAPOUR_PRESSURE;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
407 saturation *= FRACTION_N2_AIR;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
408
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
409 for(int i=0;i<16;i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
410 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
411 pLifeData->tissue_nitrogen_bar[i] = saturation;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
412 pLifeData->tissue_helium_bar[i] = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
413 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
414 pLifeData->otu = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
415 pLifeData->cns = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
416 pLifeData->desaturation_time_minutes = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
417 pLifeData->no_fly_time_minutes = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
418 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
419
129
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
420 void decom_reset_with_ambientmbar(float ambient, SLifeData * pLifeData)
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
421 {
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
422
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
423 float saturation = 1.0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
424 saturation = ambient;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
425 saturation -= WATER_VAPOUR_PRESSURE;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
426 saturation *= FRACTION_N2_AIR;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
427
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
428 for(int i=0;i<16;i++)
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
429 {
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
430 pLifeData->tissue_nitrogen_bar[i] = saturation;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
431 pLifeData->tissue_helium_bar[i] = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
432 }
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
433 pLifeData->otu = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
434 pLifeData->cns = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
435 pLifeData->desaturation_time_minutes = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
436 pLifeData->no_fly_time_minutes = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
437 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
438
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
439 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
440 /* NOTE ABOUT PRESSURE UNITS USED IN CALCULATIONS: */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
441 /* It is the convention in decompression calculations to compute all gas */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
442 /* loadings, absolute pressures, partial pressures, etc., in the units of */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
443 /* depth pressure that you are diving - either feet of seawater (fsw) or */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
444 /* meters of seawater (msw). This program follows that convention with the */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
445 /* the exception that all VPM calculations are performed in SI units (by */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
446 /* necessity). Accordingly, there are several conversions back and forth */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
447 /* between the diving pressure units and the SI units. */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
448 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
449 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
450 /* FUNCTION SUBPROGRAM FOR GAS LOADING CALCULATIONS - ASCENT AND DESCENT */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
451 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
452
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
453
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
454 float decom_schreiner_equation(float *initial_inspired_gas_pressure,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
455 float *rate_change_insp_gas_pressure,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
456 float *interval_time_minutes,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
457 const float *gas_time_constant,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
458 float *initial_gas_pressure)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
459 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
460 /* System generated locals */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
461 float ret_val;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
462 float time_null_pressure = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
463 float time_rest = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
464 float time = *interval_time_minutes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
465 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
466 /* Note: The Schreiner equation is applied when calculating the uptake or */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
467 /* elimination of compartment gases during linear ascents or descents at a */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
468 /* constant rate. For ascents, a negative number for rate must be used. */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
469 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
470 if( *rate_change_insp_gas_pressure < 0.0f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
471 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
472 time_null_pressure = -1.0f * *initial_inspired_gas_pressure / *rate_change_insp_gas_pressure;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
473 if(time > time_null_pressure )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
474 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
475 time_rest = time - time_null_pressure;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
476 time = time_null_pressure;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
477 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
478 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
479 ret_val =
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
480 *initial_inspired_gas_pressure +
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
481 *rate_change_insp_gas_pressure *
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
482 (time - 1.f / *gas_time_constant) -
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
483 (*initial_inspired_gas_pressure -
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
484 *initial_gas_pressure -
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
485 *rate_change_insp_gas_pressure / *gas_time_constant) *
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
486 expf(-(*gas_time_constant) * time);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
487
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
488 if(time_rest > 0.0f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
489 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
490 ret_val = ret_val * expf(-(*gas_time_constant) * time_rest);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
491 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
492
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
493
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
494 return ret_val;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
495 }; /* schreiner_equation__2 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
496
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
497 void decom_tissues_exposure_stage_schreiner(int period_in_seconds, SGas* pGas, float starting_ambient_pressure_bar, float ending_ambient_pressure_bar,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
498 float* pTissue_nitrogen_bar, float* pTissue_helium_bar)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
499 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
500
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
501 float initial_pressure_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
502 float initial_pressure_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
503
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
504 float ending_pressure_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
505 float ending_pressure_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
506
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
507 float fraction_N2_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
508 float fraction_N2_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
509 float fraction_He_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
510 float fraction_He_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
511
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
512 float rate_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
513 float rate_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
514
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
515 float period_in_minutes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
516
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
517 int ci;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
518
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
519 if(period_in_seconds <= 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
520 return;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
521
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
522 decom_get_inert_gases(starting_ambient_pressure_bar, pGas, &fraction_N2_begin, &fraction_He_begin );
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
523 decom_get_inert_gases(ending_ambient_pressure_bar, pGas, &fraction_N2_end, &fraction_He_end );
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
524
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
525 initial_pressure_N2 = (starting_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_N2_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
526 initial_pressure_He = (starting_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_He_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
527
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
528 ending_pressure_N2 = (ending_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_N2_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
529 ending_pressure_He = (ending_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_He_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
530
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
531 rate_N2 = (ending_pressure_N2 - initial_pressure_N2) / period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
532 rate_He = (ending_pressure_He - initial_pressure_He) / period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
533
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
534 period_in_minutes = ((float)period_in_seconds) / 60.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
535
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
536 for (ci=0;ci<16;ci++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
537 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
538 pTissue_nitrogen_bar[ci] =
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
539 decom_schreiner_equation(
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
540 &initial_pressure_N2,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
541 &rate_N2,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
542 &period_in_minutes,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
543 &nitrogen_time_constant[ci],
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
544 &pTissue_nitrogen_bar[ci]);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
545
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
546 pTissue_helium_bar[ci] =
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
547 decom_schreiner_equation(
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
548 &initial_pressure_He,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
549 &rate_He,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
550 &period_in_minutes,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
551 &helium_time_constant[ci],
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
552 &pTissue_helium_bar[ci]);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
553 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
554 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
555
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
556 _Bool nextSetpointChange(SDiveSettings* pDiveSettings, uint8_t depth_meter, uint8_t* change_depth_meter, char* setpoint)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
557 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
558 uint8_t new_depth = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
559 char new_setpoint = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
560 for(int i = 1; i <= 5; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
561 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
562 if(pDiveSettings->setpoint[i].setpoint_cbar > 0 && pDiveSettings->setpoint[i].depth_meter > 0 )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
563 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
564 if( pDiveSettings->setpoint[i].depth_meter > new_depth && pDiveSettings->setpoint[i].depth_meter < depth_meter)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
565 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
566 new_depth = pDiveSettings->setpoint[i].depth_meter;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
567 new_setpoint = pDiveSettings->setpoint[i].setpoint_cbar;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
568 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
569 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
570 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
571 if(new_depth)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
572 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
573 * change_depth_meter = new_depth;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
574 * setpoint = new_setpoint;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
575 return 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
576 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
577 return 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
578 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
579
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
580
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
581 void insertGasIntoList(SGasLine* pGas, SGasLine** pGasList, uint8_t gasInSettings, uint8_t* pGasInSettingsList, uint8_t GasListLength)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
582 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
583 uint8_t localGasIndex = GasListLength;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
584 if(pGas != 0)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
585 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
586 while(localGasIndex != 0) /* first entry */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
587 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
588 if(pGasList[localGasIndex-1]->depth_meter > pGas->depth_meter) /* switch depth of existing gas is deeper then new one => move down */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
589 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
590 pGasList[localGasIndex] = pGasList[localGasIndex-1];
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
591 pGasInSettingsList[localGasIndex] = pGasInSettingsList[localGasIndex - 1];
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
592 localGasIndex--;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
593 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
594 else
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
595 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
596 break;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
597 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
598 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
599 pGasList[localGasIndex] = pGas;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
600 pGasInSettingsList[localGasIndex] = gasInSettings;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
601 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
602 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
603
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
604 void decom_CreateGasChangeList(SDiveSettings* pInput, const SLifeData* pLifeData)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
605 {
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
606 SGasLine localPSCRFirst;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
607 SGasLine *pLocalGasList[5] = {0,0,0,0,0};
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
608 uint8_t localGasInSettingsList[5] = {0,0,0,0,0};
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
609
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
610 uint8_t gasStart = 1;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
611 uint8_t gasIndex = 0;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
612 uint8_t gasEntryCnt = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
613
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
614 int i=0;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
615 for(i=0;i< 5;i++) /* reset list */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
616 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
617 pInput->decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
618 pInput->decogaslist[i].GasIdInSettings = 255;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
619 pInput->decogaslist[i].setPoint_cbar = 0;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
620 pInput->decogaslist[i].helium_percentage = 0;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
621 pInput->decogaslist[i].nitrogen_percentage = 0;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
622 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
623 /* FirstGas
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
624 * 0 = special gas, 1 to 5 is OC gas, 6 to 10 is diluent
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
625 */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
626 pInput->decogaslist[0] = pLifeData->actualGas;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
627 /* Add Deco Gases
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
628 * special (gasId == 0) is never a deco/travel gas but actual gas only
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
629 */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
630
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
631 if(pInput->diveMode != DIVEMODE_OC)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
632 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
633 gasStart = 6; /* CCR or PSCR => CC gaslist */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
634 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
635
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
636 if(pInput->diveMode == DIVEMODE_PSCR) /* Handle first gas as deco gas */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
637 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
638 for(gasIndex = gasStart; gasIndex < gasStart + 5; gasIndex++)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
639 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
640 if(pInput->gas[gasIndex].note.ub.first)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
641 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
642 if (pLifeData->actualGas.GasIdInSettings != gasIndex)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
643 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
644 memcpy(&localPSCRFirst, &pInput->gas[gasIndex], sizeof(SGasLine));
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
645 localPSCRFirst.depth_meter = calc_MOD(gasIndex);
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
646 insertGasIntoList(&localPSCRFirst, pLocalGasList, gasIndex, localGasInSettingsList, gasEntryCnt);
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
647 gasEntryCnt++;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
648 break;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
649 }
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
650 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
651 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
652 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
653
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
654
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
655 for(gasIndex = gasStart; gasIndex < gasStart + 5; gasIndex++)
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
656 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
657 if(((pInput->gas[gasIndex].note.ub.active) && (pInput->gas[gasIndex].depth_meter)) /* ready for deco calculation */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
658 && (pLifeData->actualGas.GasIdInSettings != gasIndex) /* not the actual gas */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
659 && (pInput->gas[gasIndex].depth_meter < pLifeData->depth_meter )) /* a gas which is on the way to surface */
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
660 {
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
661 insertGasIntoList(&pInput->gas[gasIndex], pLocalGasList, gasIndex, localGasInSettingsList, gasEntryCnt);
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
662 gasEntryCnt++;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
663 }
855
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
664 }
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
665 for(gasIndex = 1; gasIndex < gasEntryCnt+1; gasIndex++) /* move gasLine Information into deco List */
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
666 {
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
667 pInput->decogaslist[gasIndex].change_during_ascent_depth_meter_otherwise_zero = pLocalGasList[gasIndex-1]->depth_meter;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
668 pInput->decogaslist[gasIndex].nitrogen_percentage = 100;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
669 pInput->decogaslist[gasIndex].nitrogen_percentage -= pLocalGasList[gasIndex-1]->oxygen_percentage;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
670 pInput->decogaslist[gasIndex].nitrogen_percentage -= pLocalGasList[gasIndex-1]->helium_percentage;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
671 pInput->decogaslist[gasIndex].helium_percentage = pLocalGasList[gasIndex-1]->helium_percentage;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
672 pInput->decogaslist[gasIndex].GasIdInSettings = localGasInSettingsList[gasIndex-1];
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
673 pInput->decogaslist[gasIndex].AppliedDiveMode = pInput->diveMode;
012f94ec2fe0 Redefine Deco gas creation function:
Ideenmodellierer
parents: 684
diff changeset
674 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
675 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
676 void test_decom_CreateGasChangeList(void)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
677 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
678 SDiveSettings diveSetting;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
679 SLifeData lifeData;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
680 lifeData.depth_meter = 100;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
681 lifeData.actualGas.helium_percentage = 30;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
682 lifeData.actualGas.nitrogen_percentage = 60;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
683 lifeData.actualGas.setPoint_cbar = 18;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
684 lifeData.actualGas.GasIdInSettings = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
685 lifeData.actualGas.change_during_ascent_depth_meter_otherwise_zero = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
686 diveSetting.diveMode = DIVEMODE_CCR;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
687 diveSetting.gas[6].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
688 diveSetting.gas[6].helium_percentage = 30;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
689 diveSetting.gas[6].oxygen_percentage = 10;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
690 diveSetting.gas[6].note.ub.active = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
691
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
692 diveSetting.gas[7].depth_meter = 60;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
693 diveSetting.gas[7].helium_percentage = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
694 diveSetting.gas[7].oxygen_percentage = 10;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
695 diveSetting.gas[7].note.ub.active = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
696 diveSetting.gas[8].note.ub.active = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
697 diveSetting.gas[9].note.ub.active = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
698 diveSetting.gas[10].note.ub.active = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
699
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
700 diveSetting.setpoint[0].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
701 diveSetting.setpoint[1].depth_meter = 80;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
702 diveSetting.setpoint[1].setpoint_cbar = 20;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
703 diveSetting.setpoint[2].depth_meter = 60;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
704 diveSetting.setpoint[2].setpoint_cbar = 25;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
705 diveSetting.setpoint[3].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
706 diveSetting.setpoint[4].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
707 diveSetting.setpoint[5].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
708
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
709
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
710 decom_CreateGasChangeList(&diveSetting, &lifeData);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
711 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
712
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
713 uint8_t decom_tissue_test_tolerance(float* Tissue_nitrogen_bar, float* Tissue_helium_bar, float GF_value, float depth_in_bar_absolute)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
714 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
715 float tissue_inertgas_saturation;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
716 float inertgas_a;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
717 float inertgas_b;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
718 float inertgas_tolerance;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
719 float gf_minus_1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
720
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
721 gf_minus_1 = GF_value - 1.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
722
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
723 for (int ci = 0; ci < 16; ci++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
724 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
725 if(Tissue_helium_bar[ci] == 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
726 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
727 tissue_inertgas_saturation = Tissue_nitrogen_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
728 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
729 inertgas_a = buehlmann_N2_a[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
730 inertgas_b = buehlmann_N2_b[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
731 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
732 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
733 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
734 tissue_inertgas_saturation = Tissue_nitrogen_bar[ci] + Tissue_helium_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
735 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
736 inertgas_a = ( ( buehlmann_N2_a[ci] * Tissue_nitrogen_bar[ci]) + ( buehlmann_He_a[ci] * Tissue_helium_bar[ci]) ) / tissue_inertgas_saturation;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
737 inertgas_b = ( ( buehlmann_N2_b[ci] * Tissue_nitrogen_bar[ci]) + ( buehlmann_He_b[ci] * Tissue_helium_bar[ci]) ) / tissue_inertgas_saturation;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
738 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
739 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
740 inertgas_tolerance = ( (GF_value / inertgas_b - gf_minus_1) * depth_in_bar_absolute ) + ( GF_value * inertgas_a );
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
741 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
742 if(inertgas_tolerance < tissue_inertgas_saturation)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
743 return 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
744 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
745 return 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
746 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
747
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
748
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
749 void decom_tissues_desaturation_time(const SLifeData* pLifeData, SLifeData2* pOutput)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
750 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
751 float pressure_in_gas_for_complete;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
752 float pressure_in_gas_for_desat;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
753 float diff_to_complete;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
754 float diff_to_desatpoint;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
755 float necessary_halftimes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
756 float desattime;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
757
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
758 pressure_in_gas_for_complete = 0.7902f * ( pLifeData->pressure_surface_bar - 0.0627f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
759 pressure_in_gas_for_desat = 1.05f * pressure_in_gas_for_complete;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
760 for(int i=0; i<16; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
761 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
762 diff_to_complete = pressure_in_gas_for_complete - pLifeData->tissue_nitrogen_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
763 diff_to_desatpoint = pressure_in_gas_for_desat - pLifeData->tissue_nitrogen_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
764
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
765 if((diff_to_desatpoint >= 0) || (diff_to_complete >= 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
766 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
767 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
768 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
769 necessary_halftimes = (logf(1.0f - (diff_to_desatpoint/diff_to_complete)) / -0.6931f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
770 desattime = buehlmann_N2_t_halflife[i] * necessary_halftimes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
771 if(desattime <= (float)0xFFFF)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
772 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = desattime;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
773 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
774 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = 0xFFFF;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
775 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
776 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
777
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
778 for(int i=0; i<16; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
779 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
780 diff_to_desatpoint = 0.05f - pLifeData->tissue_helium_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
781 diff_to_complete = -1.0f * pLifeData->tissue_helium_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
782
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
783 if((diff_to_desatpoint >= 0) || (diff_to_complete >= 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
784 pOutput->tissue_helium_desaturation_time_minutes[i] = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
785 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
786 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
787 necessary_halftimes = (logf(1.0f - (diff_to_desatpoint/diff_to_complete)) / -0.6931f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
788 desattime = buehlmann_He_t_halflife[i] * necessary_halftimes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
789 if(desattime <= (float)0xFFFF)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
790 pOutput->tissue_helium_desaturation_time_minutes[i] = desattime;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
791 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
792 pOutput->tissue_helium_desaturation_time_minutes[i] = 0xFFFF;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
793 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
794 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
795 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
796
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
797 #define MAX_DEGRADE_OTU_TIME_MINUTES (1440)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
798 //CNS&OTU:
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
799 #define OXY_TEN_MINUTES_IN_SECONDS (600)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
800 #define OXY_HALF_LIVE_OF_TEN_MINUTES__INVERSE_NINTH_ROOT_OF_TWO (0.92587471f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
801 #define OXY_NINE_DAYS_IN_TEN_MINUTES (1296)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
802 #define OXY_ONE_SIXTIETH_PART (0.0166667f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
803 #define OXY_NEGATIVE_FIVE_SIXTH_PARTS (-0.8333333f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
804 void decom_oxygen_calculate_otu(float* oxygen_otu, float pressure_oxygen_real)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
805 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
806 if(pressure_oxygen_real <= 0.5f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
807 return;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
808 *oxygen_otu += (pow((double)(0.5f / (pressure_oxygen_real - 0.5f)),OXY_NEGATIVE_FIVE_SIXTH_PARTS)) * OXY_ONE_SIXTIETH_PART;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
809 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
810
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
811 void decom_oxygen_calculate_otu_degrade(float* oxygen_otu, long seconds_since_last_dive)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
812 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
813 static long otu_time_ticker = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
814 static double otu_degrade_every_10_minutes = 999.9;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
815 long cycles_since_last_call;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
816
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
817 if((*oxygen_otu <= 0) || (seconds_since_last_dive == 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
818 *oxygen_otu = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
819 else if(seconds_since_last_dive < OXY_TEN_MINUTES_IN_SECONDS)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
820 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
821 otu_time_ticker = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
822 otu_degrade_every_10_minutes = *oxygen_otu / (MAX_DEGRADE_OTU_TIME_MINUTES / 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
823 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
824 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
825 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
826 cycles_since_last_call = seconds_since_last_dive / (otu_time_ticker * OXY_TEN_MINUTES_IN_SECONDS);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
827 *oxygen_otu -= ((double)cycles_since_last_call) * otu_degrade_every_10_minutes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
828 otu_time_ticker += cycles_since_last_call;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
829 if((*oxygen_otu < 0) || (otu_time_ticker > (MAX_DEGRADE_OTU_TIME_MINUTES / 10)))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
830 *oxygen_otu = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
831 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
832 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
833
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
834
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
835
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
836 void decom_oxygen_calculate_cns_degrade(float* oxygen_cns, long seconds_since_last_dive)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
837 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
838 static long cns_time_ticker = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
839 int cns_max_cycles;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
840
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
841 if((*oxygen_cns <= 0.5f) || (seconds_since_last_dive == 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
842 *oxygen_cns = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
843 else if(seconds_since_last_dive < OXY_TEN_MINUTES_IN_SECONDS)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
844 cns_time_ticker = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
845 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
846 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
847 cns_max_cycles = OXY_NINE_DAYS_IN_TEN_MINUTES;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
848 while((*oxygen_cns >= 0.5f) && ((cns_time_ticker * OXY_TEN_MINUTES_IN_SECONDS) < seconds_since_last_dive) && cns_max_cycles)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
849 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
850 cns_time_ticker++;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
851 cns_max_cycles--;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
852 *oxygen_cns *= OXY_HALF_LIVE_OF_TEN_MINUTES__INVERSE_NINTH_ROOT_OF_TWO;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
853 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
854 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
855 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
856
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
857
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
858 // new hwOS style
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
859 void decom_oxygen_calculate_cns(float* oxygen_cns, float pressure_oxygen_real)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
860 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
861 uint8_t char_I_actual_ppO2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
862 float CNS_fraction = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
863 const float time_factor = 3000.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
864
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
865 if(pressure_oxygen_real < 0.15f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
866 char_I_actual_ppO2 = 15;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
867 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
868 if(pressure_oxygen_real >= 2.5f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
869 char_I_actual_ppO2 = 255;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
870 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
871 char_I_actual_ppO2 = (uint8_t)(pressure_oxygen_real * 100);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
872
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
873 if (char_I_actual_ppO2 < 50)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
874 (void)0; // no changes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
875 //------------------------------------------------------------------------
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
876 // Below (and including) 1.60 bar
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
877 else if (char_I_actual_ppO2 < 61)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
878 CNS_fraction += time_factor/(-533.07f * char_I_actual_ppO2 + 54000.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
879 else if (char_I_actual_ppO2 < 71)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
880 CNS_fraction += time_factor/(-444.22f * char_I_actual_ppO2 + 48600.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
881 else if (char_I_actual_ppO2 < 81)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
882 CNS_fraction += time_factor/(-355.38f * char_I_actual_ppO2 + 42300.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
883 else if (char_I_actual_ppO2 < 91)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
884 CNS_fraction += time_factor/(-266.53f * char_I_actual_ppO2 + 35100.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
885 else if (char_I_actual_ppO2 < 111)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
886 CNS_fraction += time_factor/(-177.69f * char_I_actual_ppO2 + 27000.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
887 else if (char_I_actual_ppO2 < 152)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
888 CNS_fraction += time_factor/( -88.84f * char_I_actual_ppO2 + 17100.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
889 else if (char_I_actual_ppO2 < 167)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
890 CNS_fraction += time_factor/(-222.11f * char_I_actual_ppO2 + 37350.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
891 //------------------------------------------------------------------------
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
892 // Arieli et all.(2002): Modeling pulmonary and CNS O2 toxicity:
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
893 // J Appl Physiol 92: 248--256, 2002, doi:10.1152/japplphysiol.00434.2001
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
894 // Formula (A1) based on value for 1.55 and c=20
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
895 // example calculation: Sqrt((1.7/1.55)^20)*0.000404
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
896 else if (char_I_actual_ppO2 < 172)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
897 CNS_fraction += time_factor*0.00102f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
898 else if (char_I_actual_ppO2 < 177)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
899 CNS_fraction += time_factor*0.00136f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
900 else if (char_I_actual_ppO2 < 182)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
901 CNS_fraction += time_factor*0.00180f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
902 else if (char_I_actual_ppO2 < 187)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
903 CNS_fraction += time_factor*0.00237f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
904 else if (char_I_actual_ppO2 < 192)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
905 CNS_fraction += time_factor*0.00310f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
906 else if (char_I_actual_ppO2 < 198)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
907 CNS_fraction += time_factor*0.00401f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
908 else if (char_I_actual_ppO2 < 203)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
909 CNS_fraction += time_factor*0.00517f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
910 else if (char_I_actual_ppO2 < 233)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
911 CNS_fraction += time_factor*0.0209f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
912 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
913 CNS_fraction += time_factor*0.0482f; // value for 2.5
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
914
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
915 if( CNS_fraction > 999.0f) // Limit display to 999%
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
916 CNS_fraction = 999.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
917 if( CNS_fraction < 0.0f )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
918 CNS_fraction = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
919
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
920 //calculate cns for the actual ppo2 for 1 second
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
921 *oxygen_cns += OXY_ONE_SIXTIETH_PART * CNS_fraction;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
922
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
923 if( *oxygen_cns > 999.0f) // Limit display to 999%
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
924 *oxygen_cns = 999.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
925 if( *oxygen_cns < 0.0f )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
926 *oxygen_cns = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
927 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
928
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
929 /* old DR5 style
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
930 void decom_oxygen_calculate_cns(float* oxygen_cns, float pressure_oxygen_real)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
931 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
932 int cns_no_range = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
933 _Bool not_found = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
934 //for the cns calculation
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
935 const float cns_ppo2_ranges[60][2] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
936 {0.50f, 0.00f}, {0.60f, 0.14f}, {0.64f, 0.15f}, {0.66f, 0.16f}, {0.68f, 0.17f}, {0.70f, 0.18f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
937 {0.74f, 0.19f}, {0.76f, 0.20f}, {0.78f, 0.21f}, {0.80f, 0.22f}, {0.82f, 0.23f}, {0.84f, 0.24f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
938 {0.86f, 0.25f}, {0.88f, 0.26f}, {0.90f, 0.28f}, {0.92f, 0.29f}, {0.94f, 0.30f}, {0.96f, 0.31f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
939 {0.98f, 0.32f}, {1.00f, 0.33f}, {1.02f, 0.35f}, {1.04f, 0.36f}, {1.06f, 0.38f}, {1.08f, 0.40f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
940 {1.10f, 0.42f}, {1.12f, 0.43f}, {1.14f, 0.43f}, {1.16f, 0.44f}, {1.18f, 0.46f}, {1.20f, 0.47f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
941 {1.22f, 0.48f}, {1.24f, 0.51f}, {1.26f, 0.52f}, {1.28f, 0.54f}, {1.30f, 0.56f}, {1.32f, 0.57f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
942 {1.34f, 0.60f}, {1.36f, 0.62f}, {1.38f, 0.63f}, {1.40f, 0.65f}, {1.42f, 0.68f}, {1.44f, 0.71f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
943 {1.46f, 0.74f}, {1.48f, 0.78f}, {1.50f, 0.83f}, {1.52f, 0.93f}, {1.54f, 1.04f}, {1.56f, 1.19f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
944 {1.58f, 1.47f}, {1.60f, 2.22f}, {1.62f, 5.00f}, {1.65f, 6.25f}, {1.67f, 7.69f}, {1.70f, 10.0f},
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
945 {1.72f,12.50f}, {1.74f,20.00f}, {1.77f,25.00f}, {1.79f,31.25f}, {1.80f,50.00f}, {1.82f,100.0f}};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
946 //find the correct cns range for the corresponding ppo2
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
947 cns_no_range = 58;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
948 while (cns_no_range && not_found)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
949 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
950 if (pressure_oxygen_real > cns_ppo2_ranges[cns_no_range][0])
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
951 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
952 cns_no_range++;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
953 not_found = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
954 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
955 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
956 cns_no_range--;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
957 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
958
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
959 //calculate cns for the actual ppo2 for 1 second
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
960 *oxygen_cns += OXY_ONE_SIXTIETH_PART * cns_ppo2_ranges[cns_no_range][1];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
961 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
962 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
963
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
964 void decom_oxygen_calculate_cns_exposure(int period_in_seconds, SGas* pActualGas, float pressure_ambient_bar, float* oxygen_cns)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
965 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
966 float pressure_oxygen_real;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
967 float one_second_cns;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
968
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
969 pressure_oxygen_real = decom_calc_ppO2(pressure_ambient_bar, pActualGas);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
970 one_second_cns = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
971 decom_oxygen_calculate_cns(&one_second_cns, pressure_oxygen_real);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
972 *oxygen_cns += one_second_cns * period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
973 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
974
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
975
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
976 void decom_oxygen_calculate_cns_stage_SchreinerStyle(int period_in_seconds, SGas* pGas, float starting_ambient_pressure_bar, float ending_ambient_pressure_bar, float* oxygen_cns)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
977 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
978 if(ending_ambient_pressure_bar == starting_ambient_pressure_bar)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
979 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
980 decom_oxygen_calculate_cns_exposure(period_in_seconds, pGas, starting_ambient_pressure_bar, oxygen_cns);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
981 return;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
982 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
983
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
984 float pressure_oxygen_real;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
985 float initial_pressure_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
986 float ending_pressure_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
987 float rate_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
988
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
989 initial_pressure_oxygen = decom_calc_ppO2(starting_ambient_pressure_bar, pGas);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
990 ending_pressure_oxygen = decom_calc_ppO2(ending_ambient_pressure_bar, pGas);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
991
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
992 rate_oxygen = (ending_pressure_oxygen - initial_pressure_oxygen) / period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
993
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
994 pressure_oxygen_real = initial_pressure_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
995 for(int i = 0; i < period_in_seconds; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
996 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
997 decom_oxygen_calculate_cns(oxygen_cns, pressure_oxygen_real);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
998 pressure_oxygen_real += rate_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
999 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1000 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1001
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1002
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1003 float decom_calc_ppO2(const float ambiant_pressure_bar, const SGas* pGas)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1004 {
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1005 float percent_N2 = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1006 float percent_He = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1007 float percent_O2 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1008
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1009 decom_get_inert_gases(ambiant_pressure_bar, pGas, &percent_N2, &percent_He);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1010 percent_O2 = 1 - percent_N2 - percent_He;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1011
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1012 return (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_O2;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1013 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1014
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1015
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1016 float decom_calc_SimppO2(float ambiant_pressure_bar, const SGas* pGas)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1017 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1018 float o2Ratio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1019 float inertGasRatio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1020 float simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1021
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1022 o2Ratio = (100.0 - pGas->nitrogen_percentage - pGas->helium_percentage) / 100.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1023 inertGasRatio = 1.0 - o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1024 simulatedPSCRppo2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1025 simulatedPSCRppo2 -= (inertGasRatio * pGas->pscr_factor);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1026 if(simulatedPSCRppo2 < 0.0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1027 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1028 simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1029 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1030 return simulatedPSCRppo2;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1031 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1032
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1033 float decom_calc_SimppO2_O2based(float ambiant_pressure_bar, uint8_t O2PerCent, float factor)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1034 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1035 float o2Ratio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1036 float inertGasRatio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1037 float simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1038
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1039 o2Ratio = O2PerCent / 100.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1040 inertGasRatio = 1.0 - o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1041 simulatedPSCRppo2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1042 simulatedPSCRppo2 -= (inertGasRatio * factor);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1043 if(simulatedPSCRppo2 < 0.0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1044 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1045 simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1046 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1047 return simulatedPSCRppo2;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1048 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1049
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1050 uint8_t decom_get_actual_deco_stop(SDiveState* pDiveState)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1051 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1052 SDecoinfo* pDecoinfo;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1053 uint8_t depthNext, depthLast, depthSecond, depthInc;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1054 if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1055 pDecoinfo = &pDiveState->decolistBuehlmann;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1056 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1057 pDecoinfo = &pDiveState->decolistVPM;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1058
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1059 depthLast = (uint8_t)(pDiveState->diveSettings.last_stop_depth_bar * 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1060 depthSecond = (uint8_t)(pDiveState->diveSettings.input_second_to_last_stop_depth_bar * 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1061 depthInc = (uint8_t)(pDiveState->diveSettings.input_next_stop_increment_depth_bar * 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1062 if(pDecoinfo->output_stop_length_seconds[0] > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1063 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1064 depthNext = depthLast;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1065 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1066 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1067 return 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1068 for(int i = DECOINFO_STRUCT_MAX_STOPS -1 ;i > 0; i--)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1069 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1070 if(pDecoinfo->output_stop_length_seconds[i] > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1071 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1072 depthNext = depthSecond + ( (i - 1) * depthInc);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1073 break;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1074 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1075 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1076 return depthNext;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1077 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1078
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1079
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1080 // ===============================================================================
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1081 // decom_calc_desaturation_time
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1082 /// @brief This code is used to calculate desat, calculated by RTE and send to Firmware
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1083 /// similar but more technical in code than decom_tissues_desaturation_time()
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1084 /// the later has 0.05 for helium in contrast to this one.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1085 /// This one goes down to 70%, the oterh
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1086 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1087 /// output is desat time in minutes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1088 // ===============================================================================
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1089 int decom_calc_desaturation_time(float* Tissue_nitrogen_bar, float* Tissue_helium_bar, float surface_pressure_bar)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1090 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1091 const float N2_ratio = 0.7902; // FIXED sum as stated in b"uhlmann
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1092
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1093 float pres_surface;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1094 float temp_atem;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1095 float float_desaturation_multiplier;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1096 float temp1,temp2,temp3,temp4;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1097 int ci;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1098 int int_temp;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1099 int int_O_desaturation_time;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1100 pres_surface = ((float)surface_pressure_bar);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1101 temp_atem = N2_ratio * (pres_surface - 0.0627f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1102
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1103 int_O_desaturation_time = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1104 float_desaturation_multiplier = 100 / 142.0f; // new in v.101 (70,42%/100.=142)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1105
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1106 for (ci=0;ci<16;ci++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1107 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1108 // saturation_time (for flight) and N2_saturation in multiples of halftime
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1109 // version v.100: 1.1 = 10 percent distance to totally clean (totally clean is not possible, would take infinite time )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1110 // new in version v.101: 1.07 = 7 percent distance to totally clean (totally clean is not possible, would take infinite time )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1111 // changes in v.101: 1.05 = 5 percent dist to totally clean is new desaturation point for display and noFly calculations
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1112
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1113 // N2
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1114 temp1 = 1.05f * temp_atem;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1115 temp1 = temp1 - (float)Tissue_nitrogen_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1116 temp2 = temp_atem - (float)Tissue_nitrogen_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1117 if (temp2 >= 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1118 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1119 temp1 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1120 temp2 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1121 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1122 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1123 temp1 = temp1 / temp2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1124
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1125 if (temp1 > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1126 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1127 temp1 = logf(1.0f - temp1);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1128 temp1 = temp1 / -0.6931f; // temp1 is the multiples of half times necessary.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1129 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1130 // minus because log is negative
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1131 temp2 = buehlmann_N2_t_halflife[ci] * temp1 / float_desaturation_multiplier; // time necessary (in minutes ) for complete desaturation (see comment about 10 percent) , new in v.101: float_desaturation_multiplier
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1132 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1133 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1134 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1135 temp1 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1136 temp2 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1137 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1138
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1139 // He
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1140 temp3 = 0.1f - (float)Tissue_helium_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1141 if (temp3 >= 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1142 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1143 temp3 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1144 temp4 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1145 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1146 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1147 temp3 = -1.0f * temp3 / (float)Tissue_helium_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1148 if (temp3 > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1149 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1150 temp3 = logf(1.0f - temp3);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1151 temp3 = temp3 / -0.6931f; // temp1 is the multiples of half times necessary.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1152 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1153 // minus because log is negative
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1154 temp4 = buehlmann_He_t_halflife[ci] * temp3 / float_desaturation_multiplier; // time necessary (in minutes ) for "complete" desaturation, new in v.101 float_desaturation_multiplier
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1155 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1156 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1157 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1158 temp3 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1159 temp4 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1160 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1161
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1162 // saturation_time (for flight)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1163 if (temp4 > temp2)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1164 int_temp = (int)temp4;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1165 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1166 int_temp = (int)temp2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1167 if(int_temp > int_O_desaturation_time)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1168 int_O_desaturation_time = int_temp;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1169
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1170 /*// N2 saturation in multiples of halftime for display purposes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1171 temp2 = temp1 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1172 temp2 = temp2 + 80.0; // set center
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1173 if (temp2 < 0.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1174 temp2 = 0.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1175 if (temp2 > 255.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1176 temp2 = 255.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1177 U8_tissue_N2_saturation[ci] = (U8)temp2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1178 // He saturation in multiples of halftime for display purposes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1179 temp4 = temp3 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1180 temp4 = temp4 + 80.0; // set center
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1181 if (temp4 < 0.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1182 temp4 = 0.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1183 if (temp4 > 255.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1184 temp4 = 255.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1185 U8_tissue_He_saturation[ci] = (char)temp4;*/
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1186 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1187
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1188 return int_O_desaturation_time;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1189 }