annotate Common/Src/decom.c @ 783:c31237d20491

Update digital O2 sensor information: In the previous version only one digital O2 sensor was supported and the additional data could be accessed via a separate menu item. Because of the intraduction of the multiplexer there are three datasets available not. Because of the existing menu structure it was not possible to add additional menu items. To solve the problem the sensor information has been moved beghind the sensor On/Off functionality. When selecting a sensor the sensor data will be shown and the de-/activation of a sensor may be done in scope of the sensor information dialog.
author Ideenmodellierer
date Mon, 29 May 2023 20:26:38 +0200
parents 9bc817e9e221
children 012f94ec2fe0
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>
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
40 #include "settings.h"
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
41 #include "calc_crush.h"
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
42
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
43 # define FRACTION_N2_AIR 0.7902
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
44
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
45 const float helium_time_constant[16] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
46 3.68695308808482E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
47 2.29518933960247E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
48 1.46853216220327E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
49 9.91626867753856E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
50 6.78890480470074E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
51 4.78692804254106E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
52 3.37626488338989E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
53 2.38113081607676E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
54 1.68239606932026E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
55 1.25592893741610E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
56 9.80544886914621E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
57 7.67264977374303E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
58 6.01220557342307E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
59 4.70185307665137E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
60 3.68225234041620E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
61 2.88775228329769E-003};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
62
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
63 const float nitrogen_time_constant[16] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
64 1.38629436111989E-001,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
65 8.66433975699932E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
66 5.54517744447956E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
67 3.74674151654024E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
68 2.56721177985165E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
69 1.80978376125312E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
70 1.27651414467762E-002,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
71 9.00191143584345E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
72 6.35914844550409E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
73 4.74758342849278E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
74 3.70666941475907E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
75 2.90019740820061E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
76 2.27261370675392E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
77 1.77730046297422E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
78 1.39186180835330E-003,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
79 1.09157036308653E-003};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
80
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
81
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
82 const float buehlmann_N2_a[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
83 1.1696f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
84 1.0000f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
85 0.8618f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
86 0.7562f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
87 0.6200f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
88 0.5043f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
89 0.4410f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
90 0.4000f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
91 0.3750f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
92 0.3500f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
93 0.3295f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
94 0.3065f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
95 0.2835f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
96 0.2610f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
97 0.2480f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
98 0.2327f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
99
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
100 const float buehlmann_N2_b[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
101 0.5578f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
102 0.6514f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
103 0.7222f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
104 0.7825f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
105 0.8126f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
106 0.8434f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
107 0.8693f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
108 0.8910f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
109 0.9092f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
110 0.9222f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
111 0.9319f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
112 0.9403f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
113 0.9477f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
114 0.9544f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
115 0.9602f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
116 0.9653f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
117
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
118 const float buehlmann_He_a[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
119 1.6189f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
120 1.3830f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
121 1.1919f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
122 1.0458f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
123 0.9220f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
124 0.8205f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
125 0.7305f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
126 0.6502f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
127 0.5950f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
128 0.5545f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
129 0.5333f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
130 0.5189f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
131 0.5181f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
132 0.5176f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
133 0.5172f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
134 0.5119f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
135
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
136 const float buehlmann_He_b[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
137 0.4770f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
138 0.5747f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
139 0.6527f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
140 0.7223f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
141 0.7582f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
142 0.7957f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
143 0.8279f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
144 0.8553f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
145 0.8757f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
146 0.8903f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
147 0.8997f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
148 0.9073f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
149 0.9122f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
150 0.9171f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
151 0.9217f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
152 0.9267f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
153
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
154 const float buehlmann_N2_t_halflife[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
155 5.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
156 8.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
157 12.5f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
158 18.5f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
159 27.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
160 38.3f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
161 54.3f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
162 77.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
163 109.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
164 146.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
165 187.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
166 239.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
167 305.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
168 390.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
169 498.0f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
170 635.0f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
171
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
172 const float buehlmann_He_t_halflife[] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
173 1.88f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
174 3.02f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
175 4.72f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
176 6.99f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
177 10.21f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
178 14.48f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
179 20.53f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
180 29.11f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
181 41.20f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
182 55.19f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
183 70.69f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
184 90.34f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
185 115.29f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
186 147.42f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
187 188.24f,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
188 240.03f};
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
189
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
190 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
191 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
192 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
193 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
194 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
195 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
196 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
197 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
198 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
199 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
200 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
201
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
202 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
203 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
204 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
205 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
206 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
207 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
208 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
209 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
210 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
211 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
212 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
213
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
214 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
215 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
216 float fraction_all_inertgases;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
217 float ppo2_fraction_setpoint;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
218 float diluent_divisor;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
219
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
220 *fraction_nitrogen = ((float)pGas->nitrogen_percentage) / 100.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
221 *fraction_helium = ((float)pGas->helium_percentage) / 100.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
222
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
223 if(pGas->AppliedDiveMode == DIVEMODE_CCR)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
224 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
225 // continue with CCR
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
226 fraction_all_inertgases = *fraction_nitrogen + *fraction_helium;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
227
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
228 ppo2_fraction_setpoint = (float)pGas->setPoint_cbar/ (100 * ambient_pressure_bar);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
229
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
230 diluent_divisor = (1.0f - ppo2_fraction_setpoint) / fraction_all_inertgases;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
231 if(diluent_divisor < 0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
232 diluent_divisor = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
233
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
234 *fraction_nitrogen *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
235 *fraction_helium *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
236 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
237 if(pGas->AppliedDiveMode == DIVEMODE_PSCR)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
238 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
239 fraction_all_inertgases = *fraction_nitrogen + *fraction_helium;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
240 ppo2_fraction_setpoint = decom_calc_SimppO2(ambient_pressure_bar, pGas) / ambient_pressure_bar;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
241 diluent_divisor = (1.0f - ppo2_fraction_setpoint) / fraction_all_inertgases;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
242 if(diluent_divisor < 0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
243 diluent_divisor = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
244
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
245 *fraction_nitrogen *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
246 *fraction_helium *= diluent_divisor;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
247 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
248 }
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 void decom_tissues_exposure(int period_in_seconds, SLifeData * pLifeData)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
252 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
253 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
254 }
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 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
258 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
259 int ci;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
260 float percent_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
261 float percent_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
262 float partial_pressure_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
263 float partial_pressure_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
264
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 int period_in_seconds_left;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
268
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
269 if(period_in_seconds > 0)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
270 {
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
271
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
272 decom_get_inert_gases(ambiant_pressure_bar, pActualGas, &percent_N2, &percent_He);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
273
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
274 partial_pressure_N2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
275 partial_pressure_He = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_He;
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
276 period_in_seconds_left = period_in_seconds;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
277
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
278 while(period_in_seconds_left)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
279 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
280 if(period_in_seconds_left >= 3600)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
281 period_in_seconds = 3600;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
282 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
283 if(period_in_seconds_left >= 800)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
284 period_in_seconds = 800;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
285 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
286 if(period_in_seconds_left >= 300)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
287 period_in_seconds = 300;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
288 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
289 if(period_in_seconds_left >= 100)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
290 period_in_seconds = 100;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
291 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
292 if(period_in_seconds_left >= 60)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
293 period_in_seconds = 60;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
294 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
295 if(period_in_seconds_left == 36)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
296 period_in_seconds = 18;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
297 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
298 if(period_in_seconds_left >= 20)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
299 period_in_seconds = 20;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
300 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
301 if(period_in_seconds_left >= 18)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
302 period_in_seconds = 18;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
303 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
304 if(period_in_seconds_left >= 10)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
305 period_in_seconds = 10;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
306 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
307 if(period_in_seconds_left >= 8)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
308 period_in_seconds = 8;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
309 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
310 if(period_in_seconds_left >= 3)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
311 period_in_seconds = 3;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
312 else
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
313 period_in_seconds = 1;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
314
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
315 period_in_seconds_left -= period_in_seconds;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
316
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
317 switch (period_in_seconds)
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
318 {
57
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
319 case 1:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
320 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
321 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
322 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
323 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
324 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
325 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
326 case 3:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
327 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
328 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
329 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
330 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
331 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
332 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
333 case 8:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
334 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
335 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
336 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
337 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
338 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
339 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
340 case 10:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
341 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
342 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
343 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
344 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
345 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
346 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
347 case 18:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
348 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
349 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
350 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
351 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
352 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
353 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
354 case 20:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
355 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
356 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
357 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
358 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
359 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
360 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
361 case 60:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
362 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
363 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
364 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
365 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
366 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
367 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
368 case 100:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
369 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
370 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
371 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
372 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
373 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
374 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
375 case 300:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
376 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
377 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
378 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
379 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
380 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
381 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
382 case 800:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
383 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
384 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
385 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
386 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
387 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
388 break;
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
389 case 3600:
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
390 for (ci=0;ci<16;ci++)
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
391 {
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
392 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
393 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
394 }
e941c9e49f73 Resolve format warnings
Ideenmodellierer
parents: 38
diff changeset
395 break;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
396 }
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 void decom_reset_with_1000mbar(SLifeData * pLifeData)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
402 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
403 double saturation = 1.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
404
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
405 saturation -= WATER_VAPOUR_PRESSURE;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
406 saturation *= FRACTION_N2_AIR;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
407
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
408 for(int i=0;i<16;i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
409 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
410 pLifeData->tissue_nitrogen_bar[i] = saturation;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
411 pLifeData->tissue_helium_bar[i] = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
412 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
413 pLifeData->otu = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
414 pLifeData->cns = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
415 pLifeData->desaturation_time_minutes = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
416 pLifeData->no_fly_time_minutes = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
417 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
418
129
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
419 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
420 {
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 float saturation = 1.0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
423 saturation = ambient;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
424 saturation -= WATER_VAPOUR_PRESSURE;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
425 saturation *= FRACTION_N2_AIR;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
426
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
427 for(int i=0;i<16;i++)
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
428 {
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
429 pLifeData->tissue_nitrogen_bar[i] = saturation;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
430 pLifeData->tissue_helium_bar[i] = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
431 }
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
432 pLifeData->otu = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
433 pLifeData->cns = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
434 pLifeData->desaturation_time_minutes = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
435 pLifeData->no_fly_time_minutes = 0;
239aa58b533d Added function to initialise structures with the actual ambient pressure
Ideenmodellierer
parents: 57
diff changeset
436 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
437
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
438 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
439 /* NOTE ABOUT PRESSURE UNITS USED IN CALCULATIONS: */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
440 /* It is the convention in decompression calculations to compute all gas */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
441 /* loadings, absolute pressures, partial pressures, etc., in the units of */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
442 /* depth pressure that you are diving - either feet of seawater (fsw) or */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
443 /* meters of seawater (msw). This program follows that convention with the */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
444 /* the exception that all VPM calculations are performed in SI units (by */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
445 /* necessity). Accordingly, there are several conversions back and forth */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
446 /* between the diving pressure units and the SI units. */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
447 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
448 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
449 /* FUNCTION SUBPROGRAM FOR GAS LOADING CALCULATIONS - ASCENT AND DESCENT */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
450 /* =============================================================================== */
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 float decom_schreiner_equation(float *initial_inspired_gas_pressure,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
454 float *rate_change_insp_gas_pressure,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
455 float *interval_time_minutes,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
456 const float *gas_time_constant,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
457 float *initial_gas_pressure)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
458 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
459 /* System generated locals */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
460 float ret_val;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
461 float time_null_pressure = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
462 float time_rest = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
463 float time = *interval_time_minutes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
464 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
465 /* Note: The Schreiner equation is applied when calculating the uptake or */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
466 /* elimination of compartment gases during linear ascents or descents at a */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
467 /* constant rate. For ascents, a negative number for rate must be used. */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
468 /* =============================================================================== */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
469 if( *rate_change_insp_gas_pressure < 0.0f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
470 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
471 time_null_pressure = -1.0f * *initial_inspired_gas_pressure / *rate_change_insp_gas_pressure;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
472 if(time > time_null_pressure )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
473 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
474 time_rest = time - time_null_pressure;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
475 time = time_null_pressure;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
476 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
477 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
478 ret_val =
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
479 *initial_inspired_gas_pressure +
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
480 *rate_change_insp_gas_pressure *
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
481 (time - 1.f / *gas_time_constant) -
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
482 (*initial_inspired_gas_pressure -
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
483 *initial_gas_pressure -
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
484 *rate_change_insp_gas_pressure / *gas_time_constant) *
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
485 expf(-(*gas_time_constant) * time);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
486
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
487 if(time_rest > 0.0f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
488 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
489 ret_val = ret_val * expf(-(*gas_time_constant) * time_rest);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
490 }
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 return ret_val;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
494 }; /* schreiner_equation__2 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
495
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
496 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
497 float* pTissue_nitrogen_bar, float* pTissue_helium_bar)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
498 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
499
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
500 float initial_pressure_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
501 float initial_pressure_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
502
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
503 float ending_pressure_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
504 float ending_pressure_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
505
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
506 float fraction_N2_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
507 float fraction_N2_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
508 float fraction_He_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
509 float fraction_He_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
510
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
511 float rate_N2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
512 float rate_He;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
513
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
514 float period_in_minutes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
515
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
516 int ci;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
517
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
518 if(period_in_seconds <= 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
519 return;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
520
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
521 decom_get_inert_gases(starting_ambient_pressure_bar, pGas, &fraction_N2_begin, &fraction_He_begin );
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
522 decom_get_inert_gases(ending_ambient_pressure_bar, pGas, &fraction_N2_end, &fraction_He_end );
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
523
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
524 initial_pressure_N2 = (starting_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_N2_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
525 initial_pressure_He = (starting_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_He_begin;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
526
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
527 ending_pressure_N2 = (ending_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_N2_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
528 ending_pressure_He = (ending_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_He_end;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
529
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
530 rate_N2 = (ending_pressure_N2 - initial_pressure_N2) / period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
531 rate_He = (ending_pressure_He - initial_pressure_He) / period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
532
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
533 period_in_minutes = ((float)period_in_seconds) / 60.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
534
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
535 for (ci=0;ci<16;ci++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
536 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
537 pTissue_nitrogen_bar[ci] =
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
538 decom_schreiner_equation(
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
539 &initial_pressure_N2,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
540 &rate_N2,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
541 &period_in_minutes,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
542 &nitrogen_time_constant[ci],
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
543 &pTissue_nitrogen_bar[ci]);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
544
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
545 pTissue_helium_bar[ci] =
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
546 decom_schreiner_equation(
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
547 &initial_pressure_He,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
548 &rate_He,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
549 &period_in_minutes,
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
550 &helium_time_constant[ci],
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
551 &pTissue_helium_bar[ci]);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
552 }
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 _Bool nextSetpointChange(SDiveSettings* pDiveSettings, uint8_t depth_meter, uint8_t* change_depth_meter, char* setpoint)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
556 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
557 uint8_t new_depth = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
558 char new_setpoint = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
559 for(int i = 1; i <= 5; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
560 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
561 if(pDiveSettings->setpoint[i].setpoint_cbar > 0 && pDiveSettings->setpoint[i].depth_meter > 0 )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
562 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
563 if( pDiveSettings->setpoint[i].depth_meter > new_depth && pDiveSettings->setpoint[i].depth_meter < depth_meter)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
564 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
565 new_depth = pDiveSettings->setpoint[i].depth_meter;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
566 new_setpoint = pDiveSettings->setpoint[i].setpoint_cbar;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
567 }
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 if(new_depth)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
571 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
572 * change_depth_meter = new_depth;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
573 * setpoint = new_setpoint;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
574 return 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
575 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
576 return 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
577 }
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
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
581 void decom_CreateGasChangeList(SDiveSettings* pInput, const SLifeData* pLifeData)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
582 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
583 int i=0, j = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
584 int count = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
585 for(i=0;i< 5;i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
586 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
587 //FirstGas
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
588
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
589 pInput->decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
590 pInput->decogaslist[i].GasIdInSettings = 255;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
591 pInput->decogaslist[i].setPoint_cbar = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
592 pInput->decogaslist[i].helium_percentage = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
593 pInput->decogaslist[i].nitrogen_percentage = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
594 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
595 //pInput->liveData.dive_time_seconds = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
596
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
597 /* FirstGas
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
598 * 0 = special gas, 1 to 5 ist OC gas, 6 to 10 is diluent
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
599 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
600
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
601
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
602
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
603 pInput->decogaslist[0] = pLifeData->actualGas;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
604
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
605 /* Add Deco Gases
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
606 * special (gasId == 0) is never a deco/travel gas but actual gas only
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
607 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
608 if(pInput->diveMode == DIVEMODE_OC)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
609 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
610
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
611 for(i=1;i<= 5;i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
612 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
613 if(pInput->gas[i].note.ub.active && pInput->gas[i].depth_meter
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
614 && (pLifeData->actualGas.GasIdInSettings != i)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
615 &&(pInput->gas[i].depth_meter < pLifeData->depth_meter ) )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
616 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
617 count = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
618 for(j=1;j<= 5;j++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
619 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
620 if( (pInput->gas[j].note.ub.active && pInput->gas[j].depth_meter > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
621 && (pLifeData->actualGas.GasIdInSettings != j) // new hw 160905
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
622 && (pInput->gas[j].depth_meter > pInput->gas[i].depth_meter))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
623 count++;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
624 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
625 pInput->decogaslist[count].change_during_ascent_depth_meter_otherwise_zero = pInput->gas[i].depth_meter;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
626 pInput->decogaslist[count].nitrogen_percentage = 100;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
627 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].oxygen_percentage;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
628 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].helium_percentage;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
629 pInput->decogaslist[count].helium_percentage = pInput->gas[i].helium_percentage;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
630 pInput->decogaslist[count].GasIdInSettings = i;
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
631 pInput->decogaslist[count].AppliedDiveMode = DIVEMODE_OC;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
632 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
633 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
634 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
635 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
636 {
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
637 //divmode CCR or PSCR
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
638 for(i=6; i <= 10; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
639 {
684
9bc817e9e221 Ignore automatic setpoint changes during ascend:
Ideenmodellierer
parents: 662
diff changeset
640 if((pInput->gas[i].note.ub.active) && (pInput->gas[i].depth_meter)
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
641 && (pLifeData->actualGas.GasIdInSettings != i)
684
9bc817e9e221 Ignore automatic setpoint changes during ascend:
Ideenmodellierer
parents: 662
diff changeset
642 && (pInput->gas[i].depth_meter < pLifeData->depth_meter ))
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
643 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
644 count = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
645 for(j=6;j<= 10;j++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
646 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
647 // if(pInput->gas[j].note.ub.active && pInput->gas[j].depth_meter > 0 &&pInput->gas[j].depth_meter > pInput->gas[i].depth_meter)
684
9bc817e9e221 Ignore automatic setpoint changes during ascend:
Ideenmodellierer
parents: 662
diff changeset
648 if(((pInput->gas[j].note.ub.active) && (pInput->gas[j].depth_meter > 0))
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
649 && (pLifeData->actualGas.GasIdInSettings != j) // new hw 160905
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
650 && (pInput->gas[j].depth_meter > pInput->gas[i].depth_meter))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
651 count++;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
652 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
653 pInput->decogaslist[count].change_during_ascent_depth_meter_otherwise_zero = pInput->gas[i].depth_meter;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
654 pInput->decogaslist[count].nitrogen_percentage = 100;
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
655 if(pInput->diveMode == DIVEMODE_PSCR)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
656 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
657 pInput->decogaslist[count].AppliedDiveMode = DIVEMODE_PSCR;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
658 pInput->decogaslist[count].setPoint_cbar = decom_calc_SimppO2_O2based((float)(pInput->gas[i].depth_meter / 10.0 + 1.0), pInput->gas[i].oxygen_percentage, pInput->decogaslist[count].pscr_factor ) * 100;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
659 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].oxygen_percentage;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
660 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
661 else
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
662 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
663 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].oxygen_percentage;
684
9bc817e9e221 Ignore automatic setpoint changes during ascend:
Ideenmodellierer
parents: 662
diff changeset
664 pInput->decogaslist[count].AppliedDiveMode = DIVEMODE_CCR;
9bc817e9e221 Ignore automatic setpoint changes during ascend:
Ideenmodellierer
parents: 662
diff changeset
665 pInput->decogaslist[count].setPoint_cbar = pInput->decogaslist[0].setPoint_cbar; /* assume that current setpoint is kept till end of the dive */
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
666 }
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
667 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].helium_percentage;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
668 pInput->decogaslist[count].helium_percentage = pInput->gas[i].helium_percentage;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
669 pInput->decogaslist[count].GasIdInSettings = i;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
670 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
671 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
672 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
673 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
674 void test_decom_CreateGasChangeList(void)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
675 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
676 SDiveSettings diveSetting;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
677 SLifeData lifeData;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
678 lifeData.depth_meter = 100;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
679 lifeData.actualGas.helium_percentage = 30;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
680 lifeData.actualGas.nitrogen_percentage = 60;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
681 lifeData.actualGas.setPoint_cbar = 18;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
682 lifeData.actualGas.GasIdInSettings = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
683 lifeData.actualGas.change_during_ascent_depth_meter_otherwise_zero = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
684 diveSetting.diveMode = DIVEMODE_CCR;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
685 diveSetting.gas[6].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
686 diveSetting.gas[6].helium_percentage = 30;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
687 diveSetting.gas[6].oxygen_percentage = 10;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
688 diveSetting.gas[6].note.ub.active = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
689
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
690 diveSetting.gas[7].depth_meter = 60;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
691 diveSetting.gas[7].helium_percentage = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
692 diveSetting.gas[7].oxygen_percentage = 10;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
693 diveSetting.gas[7].note.ub.active = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
694 diveSetting.gas[8].note.ub.active = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
695 diveSetting.gas[9].note.ub.active = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
696 diveSetting.gas[10].note.ub.active = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
697
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
698 diveSetting.setpoint[0].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
699 diveSetting.setpoint[1].depth_meter = 80;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
700 diveSetting.setpoint[1].setpoint_cbar = 20;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
701 diveSetting.setpoint[2].depth_meter = 60;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
702 diveSetting.setpoint[2].setpoint_cbar = 25;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
703 diveSetting.setpoint[3].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
704 diveSetting.setpoint[4].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
705 diveSetting.setpoint[5].depth_meter = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
706
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
707
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
708 decom_CreateGasChangeList(&diveSetting, &lifeData);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
709 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
710
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
711 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
712 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
713 float tissue_inertgas_saturation;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
714 float inertgas_a;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
715 float inertgas_b;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
716 float inertgas_tolerance;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
717 float gf_minus_1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
718
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
719 gf_minus_1 = GF_value - 1.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
720
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
721 for (int ci = 0; ci < 16; ci++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
722 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
723 if(Tissue_helium_bar[ci] == 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
724 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
725 tissue_inertgas_saturation = Tissue_nitrogen_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
726 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
727 inertgas_a = buehlmann_N2_a[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
728 inertgas_b = buehlmann_N2_b[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
729 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
730 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
731 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
732 tissue_inertgas_saturation = Tissue_nitrogen_bar[ci] + Tissue_helium_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
733 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
734 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
735 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
736 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
737 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
738 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
739 //
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
740 if(inertgas_tolerance < tissue_inertgas_saturation)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
741 return 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
742 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
743 return 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
744 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
745
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
746
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
747 void decom_tissues_desaturation_time(const SLifeData* pLifeData, SLifeData2* pOutput)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
748 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
749 float pressure_in_gas_for_complete;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
750 float pressure_in_gas_for_desat;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
751 float diff_to_complete;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
752 float diff_to_desatpoint;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
753 float necessary_halftimes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
754 float desattime;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
755
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
756 pressure_in_gas_for_complete = 0.7902f * ( pLifeData->pressure_surface_bar - 0.0627f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
757 pressure_in_gas_for_desat = 1.05f * pressure_in_gas_for_complete;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
758 for(int i=0; i<16; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
759 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
760 diff_to_complete = pressure_in_gas_for_complete - pLifeData->tissue_nitrogen_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
761 diff_to_desatpoint = pressure_in_gas_for_desat - pLifeData->tissue_nitrogen_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
762
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
763 if((diff_to_desatpoint >= 0) || (diff_to_complete >= 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
764 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
765 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
766 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
767 necessary_halftimes = (logf(1.0f - (diff_to_desatpoint/diff_to_complete)) / -0.6931f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
768 desattime = buehlmann_N2_t_halflife[i] * necessary_halftimes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
769 if(desattime <= (float)0xFFFF)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
770 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = desattime;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
771 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
772 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = 0xFFFF;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
773 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
774 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
775
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
776 for(int i=0; i<16; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
777 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
778 diff_to_desatpoint = 0.05f - pLifeData->tissue_helium_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
779 diff_to_complete = -1.0f * pLifeData->tissue_helium_bar[i];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
780
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
781 if((diff_to_desatpoint >= 0) || (diff_to_complete >= 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
782 pOutput->tissue_helium_desaturation_time_minutes[i] = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
783 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
784 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
785 necessary_halftimes = (logf(1.0f - (diff_to_desatpoint/diff_to_complete)) / -0.6931f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
786 desattime = buehlmann_He_t_halflife[i] * necessary_halftimes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
787 if(desattime <= (float)0xFFFF)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
788 pOutput->tissue_helium_desaturation_time_minutes[i] = desattime;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
789 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
790 pOutput->tissue_helium_desaturation_time_minutes[i] = 0xFFFF;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
791 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
792 }
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 #define MAX_DEGRADE_OTU_TIME_MINUTES (1440)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
796 //CNS&OTU:
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
797 #define OXY_TEN_MINUTES_IN_SECONDS (600)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
798 #define OXY_HALF_LIVE_OF_TEN_MINUTES__INVERSE_NINTH_ROOT_OF_TWO (0.92587471f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
799 #define OXY_NINE_DAYS_IN_TEN_MINUTES (1296)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
800 #define OXY_ONE_SIXTIETH_PART (0.0166667f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
801 #define OXY_NEGATIVE_FIVE_SIXTH_PARTS (-0.8333333f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
802 void decom_oxygen_calculate_otu(float* oxygen_otu, float pressure_oxygen_real)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
803 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
804 if(pressure_oxygen_real <= 0.5f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
805 return;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
806 *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
807 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
808
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
809 void decom_oxygen_calculate_otu_degrade(float* oxygen_otu, long seconds_since_last_dive)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
810 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
811 static long otu_time_ticker = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
812 static double otu_degrade_every_10_minutes = 999.9;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
813 long cycles_since_last_call;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
814
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
815 if((*oxygen_otu <= 0) || (seconds_since_last_dive == 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
816 *oxygen_otu = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
817 else if(seconds_since_last_dive < OXY_TEN_MINUTES_IN_SECONDS)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
818 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
819 otu_time_ticker = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
820 otu_degrade_every_10_minutes = *oxygen_otu / (MAX_DEGRADE_OTU_TIME_MINUTES / 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
821 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
822 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
823 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
824 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
825 *oxygen_otu -= ((double)cycles_since_last_call) * otu_degrade_every_10_minutes;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
826 otu_time_ticker += cycles_since_last_call;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
827 if((*oxygen_otu < 0) || (otu_time_ticker > (MAX_DEGRADE_OTU_TIME_MINUTES / 10)))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
828 *oxygen_otu = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
829 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
830 }
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 void decom_oxygen_calculate_cns_degrade(float* oxygen_cns, long seconds_since_last_dive)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
835 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
836 static long cns_time_ticker = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
837 int cns_max_cycles;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
838
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
839 if((*oxygen_cns <= 0.5f) || (seconds_since_last_dive == 0))
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
840 *oxygen_cns = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
841 else if(seconds_since_last_dive < OXY_TEN_MINUTES_IN_SECONDS)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
842 cns_time_ticker = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
843 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
844 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
845 cns_max_cycles = OXY_NINE_DAYS_IN_TEN_MINUTES;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
846 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
847 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
848 cns_time_ticker++;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
849 cns_max_cycles--;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
850 *oxygen_cns *= OXY_HALF_LIVE_OF_TEN_MINUTES__INVERSE_NINTH_ROOT_OF_TWO;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
851 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
852 }
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 // new hwOS style
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
857 void decom_oxygen_calculate_cns(float* oxygen_cns, float pressure_oxygen_real)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
858 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
859 uint8_t char_I_actual_ppO2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
860 float CNS_fraction = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
861 const float time_factor = 3000.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
862
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
863 if(pressure_oxygen_real < 0.15f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
864 char_I_actual_ppO2 = 15;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
865 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
866 if(pressure_oxygen_real >= 2.5f)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
867 char_I_actual_ppO2 = 255;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
868 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
869 char_I_actual_ppO2 = (uint8_t)(pressure_oxygen_real * 100);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
870
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
871 if (char_I_actual_ppO2 < 50)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
872 (void)0; // no changes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
873 //------------------------------------------------------------------------
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
874 // Below (and including) 1.60 bar
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
875 else if (char_I_actual_ppO2 < 61)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
876 CNS_fraction += time_factor/(-533.07f * char_I_actual_ppO2 + 54000.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
877 else if (char_I_actual_ppO2 < 71)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
878 CNS_fraction += time_factor/(-444.22f * char_I_actual_ppO2 + 48600.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
879 else if (char_I_actual_ppO2 < 81)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
880 CNS_fraction += time_factor/(-355.38f * char_I_actual_ppO2 + 42300.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
881 else if (char_I_actual_ppO2 < 91)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
882 CNS_fraction += time_factor/(-266.53f * char_I_actual_ppO2 + 35100.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
883 else if (char_I_actual_ppO2 < 111)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
884 CNS_fraction += time_factor/(-177.69f * char_I_actual_ppO2 + 27000.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
885 else if (char_I_actual_ppO2 < 152)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
886 CNS_fraction += time_factor/( -88.84f * char_I_actual_ppO2 + 17100.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
887 else if (char_I_actual_ppO2 < 167)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
888 CNS_fraction += time_factor/(-222.11f * char_I_actual_ppO2 + 37350.0f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
889 //------------------------------------------------------------------------
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
890 // Arieli et all.(2002): Modeling pulmonary and CNS O2 toxicity:
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
891 // J Appl Physiol 92: 248--256, 2002, doi:10.1152/japplphysiol.00434.2001
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
892 // Formula (A1) based on value for 1.55 and c=20
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
893 // example calculation: Sqrt((1.7/1.55)^20)*0.000404
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
894 else if (char_I_actual_ppO2 < 172)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
895 CNS_fraction += time_factor*0.00102f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
896 else if (char_I_actual_ppO2 < 177)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
897 CNS_fraction += time_factor*0.00136f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
898 else if (char_I_actual_ppO2 < 182)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
899 CNS_fraction += time_factor*0.00180f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
900 else if (char_I_actual_ppO2 < 187)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
901 CNS_fraction += time_factor*0.00237f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
902 else if (char_I_actual_ppO2 < 192)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
903 CNS_fraction += time_factor*0.00310f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
904 else if (char_I_actual_ppO2 < 198)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
905 CNS_fraction += time_factor*0.00401f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
906 else if (char_I_actual_ppO2 < 203)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
907 CNS_fraction += time_factor*0.00517f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
908 else if (char_I_actual_ppO2 < 233)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
909 CNS_fraction += time_factor*0.0209f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
910 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
911 CNS_fraction += time_factor*0.0482f; // value for 2.5
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
912
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
913 if( CNS_fraction > 999.0f) // Limit display to 999%
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
914 CNS_fraction = 999.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
915 if( CNS_fraction < 0.0f )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
916 CNS_fraction = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
917
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
918 //calculate cns for the actual ppo2 for 1 second
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
919 *oxygen_cns += OXY_ONE_SIXTIETH_PART * CNS_fraction;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
920
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
921 if( *oxygen_cns > 999.0f) // Limit display to 999%
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
922 *oxygen_cns = 999.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
923 if( *oxygen_cns < 0.0f )
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
924 *oxygen_cns = 0.0f;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
925 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
926
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
927 /* old DR5 style
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
928 void decom_oxygen_calculate_cns(float* oxygen_cns, float pressure_oxygen_real)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
929 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
930 int cns_no_range = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
931 _Bool not_found = 1;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
932 //for the cns calculation
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
933 const float cns_ppo2_ranges[60][2] = {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
934 {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
935 {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
936 {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
937 {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
938 {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
939 {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
940 {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
941 {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
942 {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
943 {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
944 //find the correct cns range for the corresponding ppo2
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
945 cns_no_range = 58;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
946 while (cns_no_range && not_found)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
947 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
948 if (pressure_oxygen_real > cns_ppo2_ranges[cns_no_range][0])
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
949 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
950 cns_no_range++;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
951 not_found = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
952 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
953 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
954 cns_no_range--;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
955 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
956
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
957 //calculate cns for the actual ppo2 for 1 second
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
958 *oxygen_cns += OXY_ONE_SIXTIETH_PART * cns_ppo2_ranges[cns_no_range][1];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
959 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
960 */
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
961
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
962 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
963 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
964 float pressure_oxygen_real;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
965 float one_second_cns;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
966
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
967 pressure_oxygen_real = decom_calc_ppO2(pressure_ambient_bar, pActualGas);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
968 one_second_cns = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
969 decom_oxygen_calculate_cns(&one_second_cns, pressure_oxygen_real);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
970 *oxygen_cns += one_second_cns * period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
971 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
972
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
973
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
974 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
975 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
976 if(ending_ambient_pressure_bar == starting_ambient_pressure_bar)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
977 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
978 decom_oxygen_calculate_cns_exposure(period_in_seconds, pGas, starting_ambient_pressure_bar, oxygen_cns);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
979 return;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
980 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
981
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
982 float pressure_oxygen_real;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
983 float initial_pressure_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
984 float ending_pressure_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
985 float rate_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
986
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
987 initial_pressure_oxygen = decom_calc_ppO2(starting_ambient_pressure_bar, pGas);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
988 ending_pressure_oxygen = decom_calc_ppO2(ending_ambient_pressure_bar, pGas);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
989
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
990 rate_oxygen = (ending_pressure_oxygen - initial_pressure_oxygen) / period_in_seconds;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
991
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
992 pressure_oxygen_real = initial_pressure_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
993 for(int i = 0; i < period_in_seconds; i++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
994 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
995 decom_oxygen_calculate_cns(oxygen_cns, pressure_oxygen_real);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
996 pressure_oxygen_real += rate_oxygen;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
997 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
998 }
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 float decom_calc_ppO2(const float ambiant_pressure_bar, const SGas* pGas)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1002 {
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1003 float percent_N2 = 0;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1004 float percent_He = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1005 float percent_O2 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1006
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1007 decom_get_inert_gases(ambiant_pressure_bar, pGas, &percent_N2, &percent_He);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1008 percent_O2 = 1 - percent_N2 - percent_He;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1009
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1010 return (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_O2;
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1011 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1012
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1013
662
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1014 float decom_calc_SimppO2(float ambiant_pressure_bar, const SGas* pGas)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1015 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1016 float o2Ratio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1017 float inertGasRatio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1018 float simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1019
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1020 o2Ratio = (100.0 - pGas->nitrogen_percentage - pGas->helium_percentage) / 100.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1021 inertGasRatio = 1.0 - o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1022 simulatedPSCRppo2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1023 simulatedPSCRppo2 -= (inertGasRatio * pGas->pscr_factor);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1024 if(simulatedPSCRppo2 < 0.0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1025 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1026 simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1027 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1028 return simulatedPSCRppo2;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1029 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1030
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1031 float decom_calc_SimppO2_O2based(float ambiant_pressure_bar, uint8_t O2PerCent, float factor)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1032 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1033 float o2Ratio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1034 float inertGasRatio = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1035 float simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1036
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1037 o2Ratio = O2PerCent / 100.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1038 inertGasRatio = 1.0 - o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1039 simulatedPSCRppo2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * o2Ratio;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1040 simulatedPSCRppo2 -= (inertGasRatio * factor);
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1041 if(simulatedPSCRppo2 < 0.0)
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1042 {
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1043 simulatedPSCRppo2 = 0.0;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1044 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1045 return simulatedPSCRppo2;
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1046 }
1b995079c045 PSCR Mode
heinrichs weikamp
parents: 129
diff changeset
1047
38
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1048 uint8_t decom_get_actual_deco_stop(SDiveState* pDiveState)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1049 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1050 SDecoinfo* pDecoinfo;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1051 uint8_t depthNext, depthLast, depthSecond, depthInc;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1052 if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1053 pDecoinfo = &pDiveState->decolistBuehlmann;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1054 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1055 pDecoinfo = &pDiveState->decolistVPM;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1056
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1057 depthLast = (uint8_t)(pDiveState->diveSettings.last_stop_depth_bar * 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1058 depthSecond = (uint8_t)(pDiveState->diveSettings.input_second_to_last_stop_depth_bar * 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1059 depthInc = (uint8_t)(pDiveState->diveSettings.input_next_stop_increment_depth_bar * 10);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1060 if(pDecoinfo->output_stop_length_seconds[0] > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1061 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1062 depthNext = depthLast;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1063 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1064 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1065 return 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1066 for(int i = DECOINFO_STRUCT_MAX_STOPS -1 ;i > 0; i--)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1067 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1068 if(pDecoinfo->output_stop_length_seconds[i] > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1069 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1070 depthNext = depthSecond + ( (i - 1) * depthInc);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1071 break;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1072 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1073 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1074 return depthNext;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1075 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1076
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 // decom_calc_desaturation_time
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1080 /// @brief This code is used to calculate desat, calculated by RTE and send to Firmware
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1081 /// similar but more technical in code than decom_tissues_desaturation_time()
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1082 /// the later has 0.05 for helium in contrast to this one.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1083 /// This one goes down to 70%, the oterh
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1084 ///
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1085 /// output is desat time in minutes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1086 // ===============================================================================
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1087 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
1088 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1089 const float N2_ratio = 0.7902; // FIXED sum as stated in b"uhlmann
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1090
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1091 float pres_surface;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1092 float temp_atem;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1093 float float_desaturation_multiplier;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1094 float temp1,temp2,temp3,temp4;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1095 int ci;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1096 int int_temp;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1097 int int_O_desaturation_time;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1098 pres_surface = ((float)surface_pressure_bar);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1099 temp_atem = N2_ratio * (pres_surface - 0.0627f);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1100
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1101 int_O_desaturation_time = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1102 float_desaturation_multiplier = 100 / 142.0f; // new in v.101 (70,42%/100.=142)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1103
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1104 for (ci=0;ci<16;ci++)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1105 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1106 // saturation_time (for flight) and N2_saturation in multiples of halftime
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1107 // 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
1108 // 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
1109 // 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
1110
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1111 // N2
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1112 temp1 = 1.05f * temp_atem;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1113 temp1 = temp1 - (float)Tissue_nitrogen_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1114 temp2 = temp_atem - (float)Tissue_nitrogen_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1115 if (temp2 >= 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1116 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1117 temp1 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1118 temp2 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1119 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1120 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1121 temp1 = temp1 / temp2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1122
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1123 if (temp1 > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1124 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1125 temp1 = logf(1.0f - temp1);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1126 temp1 = temp1 / -0.6931f; // temp1 is the multiples of half times necessary.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1127 // 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
1128 // minus because log is negative
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1129 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
1130 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1131 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1132 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1133 temp1 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1134 temp2 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1135 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1136
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1137 // He
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1138 temp3 = 0.1f - (float)Tissue_helium_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1139 if (temp3 >= 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1140 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1141 temp3 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1142 temp4 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1143 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1144 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1145 temp3 = -1.0f * temp3 / (float)Tissue_helium_bar[ci];
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1146 if (temp3 > 0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1147 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1148 temp3 = logf(1.0f - temp3);
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1149 temp3 = temp3 / -0.6931f; // temp1 is the multiples of half times necessary.
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1150 // 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
1151 // minus because log is negative
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1152 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
1153 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1154 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1155 {
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1156 temp3 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1157 temp4 = 0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1158 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1159
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1160 // saturation_time (for flight)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1161 if (temp4 > temp2)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1162 int_temp = (int)temp4;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1163 else
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1164 int_temp = (int)temp2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1165 if(int_temp > int_O_desaturation_time)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1166 int_O_desaturation_time = int_temp;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1167
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1168 /*// N2 saturation in multiples of halftime for display purposes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1169 temp2 = temp1 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1170 temp2 = temp2 + 80.0; // set center
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1171 if (temp2 < 0.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1172 temp2 = 0.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1173 if (temp2 > 255.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1174 temp2 = 255.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1175 U8_tissue_N2_saturation[ci] = (U8)temp2;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1176 // He saturation in multiples of halftime for display purposes
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1177 temp4 = temp3 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1178 temp4 = temp4 + 80.0; // set center
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1179 if (temp4 < 0.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1180 temp4 = 0.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1181 if (temp4 > 255.0)
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1182 temp4 = 255.0;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1183 U8_tissue_He_saturation[ci] = (char)temp4;*/
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1184 }
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1185
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1186 return int_O_desaturation_time;
5f11787b4f42 include in ostc4 repository
heinrichsweikamp
parents:
diff changeset
1187 }