Mercurial > public > ostc4
annotate Common/Src/decom.c @ 677:1953cd9145b3 Betatest
Move include file to common to ...
... get things compiled. The configuration.h file is also used from the RTE
code, so it should be in the Common tree.
Signed-off-by: Jan Mulder <jan@jlmulder.nl>
| author | Jan Mulder <jan@jlmulder.nl> |
|---|---|
| date | Wed, 20 Apr 2022 20:25:04 +0200 |
| parents | 1b995079c045 |
| children | 9bc817e9e221 |
| rev | line source |
|---|---|
| 38 | 1 /////////////////////////////////////////////////////////////////////////////// |
| 2 /// -*- coding: UTF-8 -*- | |
| 3 /// | |
| 4 /// \file Common/Src/decom.c | |
| 5 /// \brief This code is used to calculate desat, calculated by RTE and send to Firmware | |
| 6 /// \author heinrichs weikamp gmbh | |
| 7 /// \date 22-Feb-2016 | |
| 8 /// | |
| 9 /// $Id$ | |
| 10 /////////////////////////////////////////////////////////////////////////////// | |
| 11 /// \par Copyright (c) 2014-2018 Heinrichs Weikamp gmbh | |
| 12 /// | |
| 13 /// This program is free software: you can redistribute it and/or modify | |
| 14 /// it under the terms of the GNU General Public License as published by | |
| 15 /// the Free Software Foundation, either version 3 of the License, or | |
| 16 /// (at your option) any later version. | |
| 17 /// | |
| 18 /// This program is distributed in the hope that it will be useful, | |
| 19 /// but WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 20 /// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
| 21 /// GNU General Public License for more details. | |
| 22 /// | |
| 23 /// You should have received a copy of the GNU General Public License | |
| 24 /// along with this program. If not, see <http://www.gnu.org/licenses/>. | |
| 25 ////////////////////////////////////////////////////////////////////////////// | |
| 26 /** | |
| 27 @verbatim | |
| 28 ============================================================================== | |
| 29 ##### Changes ##### | |
| 30 ============================================================================== | |
| 31 V1.0.2 1602220x decom_oxygen_calculate_cns() changed to hwOS version | |
| 32 | |
| 33 @endverbatim | |
| 34 ****************************************************************************** | |
| 35 */ | |
| 36 | |
| 37 #include "decom.h" | |
| 38 | |
| 39 #include <math.h> | |
| 40 #include "settings.h" | |
| 41 #include "calc_crush.h" | |
| 42 | |
| 43 # define FRACTION_N2_AIR 0.7902 | |
| 44 | |
| 45 const float helium_time_constant[16] = { | |
| 46 3.68695308808482E-001, | |
| 47 2.29518933960247E-001, | |
| 48 1.46853216220327E-001, | |
| 49 9.91626867753856E-002, | |
| 50 6.78890480470074E-002, | |
| 51 4.78692804254106E-002, | |
| 52 3.37626488338989E-002, | |
| 53 2.38113081607676E-002, | |
| 54 1.68239606932026E-002, | |
| 55 1.25592893741610E-002, | |
| 56 9.80544886914621E-003, | |
| 57 7.67264977374303E-003, | |
| 58 6.01220557342307E-003, | |
| 59 4.70185307665137E-003, | |
| 60 3.68225234041620E-003, | |
| 61 2.88775228329769E-003}; | |
| 62 | |
| 63 const float nitrogen_time_constant[16] = { | |
| 64 1.38629436111989E-001, | |
| 65 8.66433975699932E-002, | |
| 66 5.54517744447956E-002, | |
| 67 3.74674151654024E-002, | |
| 68 2.56721177985165E-002, | |
| 69 1.80978376125312E-002, | |
| 70 1.27651414467762E-002, | |
| 71 9.00191143584345E-003, | |
| 72 6.35914844550409E-003, | |
| 73 4.74758342849278E-003, | |
| 74 3.70666941475907E-003, | |
| 75 2.90019740820061E-003, | |
| 76 2.27261370675392E-003, | |
| 77 1.77730046297422E-003, | |
| 78 1.39186180835330E-003, | |
| 79 1.09157036308653E-003}; | |
| 80 | |
| 81 | |
| 82 const float buehlmann_N2_a[] = { | |
| 83 1.1696f, | |
| 84 1.0000f, | |
| 85 0.8618f, | |
| 86 0.7562f, | |
| 87 0.6200f, | |
| 88 0.5043f, | |
| 89 0.4410f, | |
| 90 0.4000f, | |
| 91 0.3750f, | |
| 92 0.3500f, | |
| 93 0.3295f, | |
| 94 0.3065f, | |
| 95 0.2835f, | |
| 96 0.2610f, | |
| 97 0.2480f, | |
| 98 0.2327f}; | |
| 99 | |
| 100 const float buehlmann_N2_b[] = { | |
| 101 0.5578f, | |
| 102 0.6514f, | |
| 103 0.7222f, | |
| 104 0.7825f, | |
| 105 0.8126f, | |
| 106 0.8434f, | |
| 107 0.8693f, | |
| 108 0.8910f, | |
| 109 0.9092f, | |
| 110 0.9222f, | |
| 111 0.9319f, | |
| 112 0.9403f, | |
| 113 0.9477f, | |
| 114 0.9544f, | |
| 115 0.9602f, | |
| 116 0.9653f}; | |
| 117 | |
| 118 const float buehlmann_He_a[] = { | |
| 119 1.6189f, | |
| 120 1.3830f, | |
| 121 1.1919f, | |
| 122 1.0458f, | |
| 123 0.9220f, | |
| 124 0.8205f, | |
| 125 0.7305f, | |
| 126 0.6502f, | |
| 127 0.5950f, | |
| 128 0.5545f, | |
| 129 0.5333f, | |
| 130 0.5189f, | |
| 131 0.5181f, | |
| 132 0.5176f, | |
| 133 0.5172f, | |
| 134 0.5119f}; | |
| 135 | |
| 136 const float buehlmann_He_b[] = { | |
| 137 0.4770f, | |
| 138 0.5747f, | |
| 139 0.6527f, | |
| 140 0.7223f, | |
| 141 0.7582f, | |
| 142 0.7957f, | |
| 143 0.8279f, | |
| 144 0.8553f, | |
| 145 0.8757f, | |
| 146 0.8903f, | |
| 147 0.8997f, | |
| 148 0.9073f, | |
| 149 0.9122f, | |
| 150 0.9171f, | |
| 151 0.9217f, | |
| 152 0.9267f}; | |
| 153 | |
| 154 const float buehlmann_N2_t_halflife[] = { | |
| 155 5.0f, | |
| 156 8.0f, | |
| 157 12.5f, | |
| 158 18.5f, | |
| 159 27.0f, | |
| 160 38.3f, | |
| 161 54.3f, | |
| 162 77.0f, | |
| 163 109.0f, | |
| 164 146.0f, | |
| 165 187.0f, | |
| 166 239.0f, | |
| 167 305.0f, | |
| 168 390.0f, | |
| 169 498.0f, | |
| 170 635.0f}; | |
| 171 | |
| 172 const float buehlmann_He_t_halflife[] = { | |
| 173 1.88f, | |
| 174 3.02f, | |
| 175 4.72f, | |
| 176 6.99f, | |
| 177 10.21f, | |
| 178 14.48f, | |
| 179 20.53f, | |
| 180 29.11f, | |
| 181 41.20f, | |
| 182 55.19f, | |
| 183 70.69f, | |
| 184 90.34f, | |
| 185 115.29f, | |
| 186 147.42f, | |
| 187 188.24f, | |
| 188 240.03f}; | |
| 189 | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 201 | |
| 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}; | |
| 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 | |
| 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 | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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}; | |
| 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 | |
| 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 | |
| 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 | |
| 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}; | |
| 213 | |
| 214 void decom_get_inert_gases(const float ambient_pressure_bar,const SGas* pGas, float* fraction_nitrogen, float* fraction_helium ) | |
| 215 { | |
| 216 float fraction_all_inertgases; | |
| 217 float ppo2_fraction_setpoint; | |
| 218 float diluent_divisor; | |
| 219 | |
| 220 *fraction_nitrogen = ((float)pGas->nitrogen_percentage) / 100.0f; | |
| 221 *fraction_helium = ((float)pGas->helium_percentage) / 100.0f; | |
| 222 | |
| 662 | 223 if(pGas->AppliedDiveMode == DIVEMODE_CCR) |
| 224 { | |
| 225 // continue with CCR | |
| 226 fraction_all_inertgases = *fraction_nitrogen + *fraction_helium; | |
| 38 | 227 |
| 662 | 228 ppo2_fraction_setpoint = (float)pGas->setPoint_cbar/ (100 * ambient_pressure_bar); |
| 229 | |
| 230 diluent_divisor = (1.0f - ppo2_fraction_setpoint) / fraction_all_inertgases; | |
| 231 if(diluent_divisor < 0) | |
| 232 diluent_divisor = 0; | |
| 38 | 233 |
| 662 | 234 *fraction_nitrogen *= diluent_divisor; |
| 235 *fraction_helium *= diluent_divisor; | |
| 236 } | |
| 237 if(pGas->AppliedDiveMode == DIVEMODE_PSCR) | |
| 238 { | |
| 239 fraction_all_inertgases = *fraction_nitrogen + *fraction_helium; | |
| 240 ppo2_fraction_setpoint = decom_calc_SimppO2(ambient_pressure_bar, pGas) / ambient_pressure_bar; | |
| 241 diluent_divisor = (1.0f - ppo2_fraction_setpoint) / fraction_all_inertgases; | |
| 242 if(diluent_divisor < 0) | |
| 243 diluent_divisor = 0; | |
| 38 | 244 |
| 662 | 245 *fraction_nitrogen *= diluent_divisor; |
| 246 *fraction_helium *= diluent_divisor; | |
| 247 } | |
| 38 | 248 } |
| 249 | |
| 250 | |
| 251 void decom_tissues_exposure(int period_in_seconds, SLifeData * pLifeData) | |
| 252 { | |
| 253 decom_tissues_exposure2(period_in_seconds, &pLifeData->actualGas, pLifeData->pressure_ambient_bar, pLifeData->tissue_nitrogen_bar, pLifeData->tissue_helium_bar); | |
| 254 } | |
| 255 | |
| 256 | |
| 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) | |
| 258 { | |
| 259 int ci; | |
| 260 float percent_N2; | |
| 261 float percent_He; | |
| 262 float partial_pressure_N2; | |
| 263 float partial_pressure_He; | |
| 264 | |
| 265 | |
| 266 | |
| 267 int period_in_seconds_left; | |
| 268 | |
| 57 | 269 if(period_in_seconds > 0) |
| 270 { | |
| 38 | 271 |
| 272 decom_get_inert_gases(ambiant_pressure_bar, pActualGas, &percent_N2, &percent_He); | |
| 273 | |
| 274 partial_pressure_N2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_N2; | |
| 275 partial_pressure_He = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_He; | |
| 57 | 276 period_in_seconds_left = period_in_seconds; |
| 38 | 277 |
| 57 | 278 while(period_in_seconds_left) |
| 279 { | |
| 280 if(period_in_seconds_left >= 3600) | |
| 281 period_in_seconds = 3600; | |
| 282 else | |
| 283 if(period_in_seconds_left >= 800) | |
| 284 period_in_seconds = 800; | |
| 285 else | |
| 286 if(period_in_seconds_left >= 300) | |
| 287 period_in_seconds = 300; | |
| 288 else | |
| 289 if(period_in_seconds_left >= 100) | |
| 290 period_in_seconds = 100; | |
| 291 else | |
| 292 if(period_in_seconds_left >= 60) | |
| 293 period_in_seconds = 60; | |
| 294 else | |
| 295 if(period_in_seconds_left == 36) | |
| 296 period_in_seconds = 18; | |
| 297 else | |
| 298 if(period_in_seconds_left >= 20) | |
| 299 period_in_seconds = 20; | |
| 300 else | |
| 301 if(period_in_seconds_left >= 18) | |
| 302 period_in_seconds = 18; | |
| 303 else | |
| 304 if(period_in_seconds_left >= 10) | |
| 305 period_in_seconds = 10; | |
| 306 else | |
| 307 if(period_in_seconds_left >= 8) | |
| 308 period_in_seconds = 8; | |
| 309 else | |
| 310 if(period_in_seconds_left >= 3) | |
| 311 period_in_seconds = 3; | |
| 312 else | |
| 313 period_in_seconds = 1; | |
| 38 | 314 |
| 57 | 315 period_in_seconds_left -= period_in_seconds; |
| 38 | 316 |
| 57 | 317 switch (period_in_seconds) |
| 38 | 318 { |
| 57 | 319 case 1: |
| 320 for (ci=0;ci<16;ci++) | |
| 321 { | |
| 322 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_one_second[ci]; | |
| 323 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_one_second[ci]; | |
| 324 } | |
| 325 break; | |
| 326 case 3: | |
| 327 for (ci=0;ci<16;ci++) | |
| 328 { | |
| 329 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_003_second[ci]; | |
| 330 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_003_second[ci]; | |
| 331 } | |
| 332 break; | |
| 333 case 8: | |
| 334 for (ci=0;ci<16;ci++) | |
| 335 { | |
| 336 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_008_second[ci]; | |
| 337 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_008_second[ci]; | |
| 338 } | |
| 339 break; | |
| 340 case 10: | |
| 341 for (ci=0;ci<16;ci++) | |
| 342 { | |
| 343 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_10_seconds[ci]; | |
| 344 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_10_seconds[ci]; | |
| 345 } | |
| 346 break; | |
| 347 case 18: | |
| 348 for (ci=0;ci<16;ci++) | |
| 349 { | |
| 350 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_18_seconds[ci]; | |
| 351 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_18_seconds[ci]; | |
| 352 } | |
| 353 break; | |
| 354 case 20: | |
| 355 for (ci=0;ci<16;ci++) | |
| 356 { | |
| 357 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_20_seconds[ci]; | |
| 358 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_20_seconds[ci]; | |
| 359 } | |
| 360 break; | |
| 361 case 60: | |
| 362 for (ci=0;ci<16;ci++) | |
| 363 { | |
| 364 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_one_minute[ci]; | |
| 365 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_one_minute[ci]; | |
| 366 } | |
| 367 break; | |
| 368 case 100: | |
| 369 for (ci=0;ci<16;ci++) | |
| 370 { | |
| 371 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_100_second[ci]; | |
| 372 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_100_second[ci]; | |
| 373 } | |
| 374 break; | |
| 375 case 300: | |
| 376 for (ci=0;ci<16;ci++) | |
| 377 { | |
| 378 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_five_minutes[ci]; | |
| 379 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_five_minutes[ci]; | |
| 380 } | |
| 381 break; | |
| 382 case 800: | |
| 383 for (ci=0;ci<16;ci++) | |
| 384 { | |
| 385 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_800_second[ci]; | |
| 386 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_800_second[ci]; | |
| 387 } | |
| 388 break; | |
| 389 case 3600: | |
| 390 for (ci=0;ci<16;ci++) | |
| 391 { | |
| 392 tissue_N2_selected_stage[ci] += (partial_pressure_N2 - tissue_N2_selected_stage[ci]) * float_buehlmann_N2_factor_expositon_one_hour[ci]; | |
| 393 tissue_He_selected_stage[ci] += (partial_pressure_He - tissue_He_selected_stage[ci]) * float_buehlmann_He_factor_expositon_one_hour[ci]; | |
| 394 } | |
| 395 break; | |
| 38 | 396 } |
| 397 } | |
| 398 } | |
| 399 } | |
| 400 | |
| 401 void decom_reset_with_1000mbar(SLifeData * pLifeData) | |
| 402 { | |
| 403 double saturation = 1.0; | |
| 404 | |
| 405 saturation -= WATER_VAPOUR_PRESSURE; | |
| 406 saturation *= FRACTION_N2_AIR; | |
| 407 | |
| 408 for(int i=0;i<16;i++) | |
| 409 { | |
| 410 pLifeData->tissue_nitrogen_bar[i] = saturation; | |
| 411 pLifeData->tissue_helium_bar[i] = 0; | |
| 412 } | |
| 413 pLifeData->otu = 0; | |
| 414 pLifeData->cns = 0; | |
| 415 pLifeData->desaturation_time_minutes = 0; | |
| 416 pLifeData->no_fly_time_minutes = 0; | |
| 417 } | |
| 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 | 437 |
| 438 /* =============================================================================== */ | |
| 439 /* NOTE ABOUT PRESSURE UNITS USED IN CALCULATIONS: */ | |
| 440 /* It is the convention in decompression calculations to compute all gas */ | |
| 441 /* loadings, absolute pressures, partial pressures, etc., in the units of */ | |
| 442 /* depth pressure that you are diving - either feet of seawater (fsw) or */ | |
| 443 /* meters of seawater (msw). This program follows that convention with the */ | |
| 444 /* the exception that all VPM calculations are performed in SI units (by */ | |
| 445 /* necessity). Accordingly, there are several conversions back and forth */ | |
| 446 /* between the diving pressure units and the SI units. */ | |
| 447 /* =============================================================================== */ | |
| 448 /* =============================================================================== */ | |
| 449 /* FUNCTION SUBPROGRAM FOR GAS LOADING CALCULATIONS - ASCENT AND DESCENT */ | |
| 450 /* =============================================================================== */ | |
| 451 | |
| 452 | |
| 453 float decom_schreiner_equation(float *initial_inspired_gas_pressure, | |
| 454 float *rate_change_insp_gas_pressure, | |
| 455 float *interval_time_minutes, | |
| 456 const float *gas_time_constant, | |
| 457 float *initial_gas_pressure) | |
| 458 { | |
| 459 /* System generated locals */ | |
| 460 float ret_val; | |
| 461 float time_null_pressure = 0.0f; | |
| 462 float time_rest = 0.0f; | |
| 463 float time = *interval_time_minutes; | |
| 464 /* =============================================================================== */ | |
| 465 /* Note: The Schreiner equation is applied when calculating the uptake or */ | |
| 466 /* elimination of compartment gases during linear ascents or descents at a */ | |
| 467 /* constant rate. For ascents, a negative number for rate must be used. */ | |
| 468 /* =============================================================================== */ | |
| 469 if( *rate_change_insp_gas_pressure < 0.0f) | |
| 470 { | |
| 471 time_null_pressure = -1.0f * *initial_inspired_gas_pressure / *rate_change_insp_gas_pressure; | |
| 472 if(time > time_null_pressure ) | |
| 473 { | |
| 474 time_rest = time - time_null_pressure; | |
| 475 time = time_null_pressure; | |
| 476 } | |
| 477 } | |
| 478 ret_val = | |
| 479 *initial_inspired_gas_pressure + | |
| 480 *rate_change_insp_gas_pressure * | |
| 481 (time - 1.f / *gas_time_constant) - | |
| 482 (*initial_inspired_gas_pressure - | |
| 483 *initial_gas_pressure - | |
| 484 *rate_change_insp_gas_pressure / *gas_time_constant) * | |
| 485 expf(-(*gas_time_constant) * time); | |
| 486 | |
| 487 if(time_rest > 0.0f) | |
| 488 { | |
| 489 ret_val = ret_val * expf(-(*gas_time_constant) * time_rest); | |
| 490 } | |
| 491 | |
| 492 | |
| 493 return ret_val; | |
| 494 }; /* schreiner_equation__2 */ | |
| 495 | |
| 496 void decom_tissues_exposure_stage_schreiner(int period_in_seconds, SGas* pGas, float starting_ambient_pressure_bar, float ending_ambient_pressure_bar, | |
| 497 float* pTissue_nitrogen_bar, float* pTissue_helium_bar) | |
| 498 { | |
| 499 | |
| 500 float initial_pressure_N2; | |
| 501 float initial_pressure_He; | |
| 502 | |
| 503 float ending_pressure_N2; | |
| 504 float ending_pressure_He; | |
| 505 | |
| 506 float fraction_N2_begin; | |
| 507 float fraction_N2_end; | |
| 508 float fraction_He_begin; | |
| 509 float fraction_He_end; | |
| 510 | |
| 511 float rate_N2; | |
| 512 float rate_He; | |
| 513 | |
| 514 float period_in_minutes; | |
| 515 | |
| 516 int ci; | |
| 517 | |
| 518 if(period_in_seconds <= 0) | |
| 519 return; | |
| 520 | |
| 521 decom_get_inert_gases(starting_ambient_pressure_bar, pGas, &fraction_N2_begin, &fraction_He_begin ); | |
| 522 decom_get_inert_gases(ending_ambient_pressure_bar, pGas, &fraction_N2_end, &fraction_He_end ); | |
| 523 | |
| 524 initial_pressure_N2 = (starting_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_N2_begin; | |
| 525 initial_pressure_He = (starting_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_He_begin; | |
| 526 | |
| 527 ending_pressure_N2 = (ending_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_N2_end; | |
| 528 ending_pressure_He = (ending_ambient_pressure_bar - WATER_VAPOUR_PRESSURE) * fraction_He_end; | |
| 529 | |
| 530 rate_N2 = (ending_pressure_N2 - initial_pressure_N2) / period_in_seconds; | |
| 531 rate_He = (ending_pressure_He - initial_pressure_He) / period_in_seconds; | |
| 532 | |
| 533 period_in_minutes = ((float)period_in_seconds) / 60.0f; | |
| 534 | |
| 535 for (ci=0;ci<16;ci++) | |
| 536 { | |
| 537 pTissue_nitrogen_bar[ci] = | |
| 538 decom_schreiner_equation( | |
| 539 &initial_pressure_N2, | |
| 540 &rate_N2, | |
| 541 &period_in_minutes, | |
| 542 &nitrogen_time_constant[ci], | |
| 543 &pTissue_nitrogen_bar[ci]); | |
| 544 | |
| 545 pTissue_helium_bar[ci] = | |
| 546 decom_schreiner_equation( | |
| 547 &initial_pressure_He, | |
| 548 &rate_He, | |
| 549 &period_in_minutes, | |
| 550 &helium_time_constant[ci], | |
| 551 &pTissue_helium_bar[ci]); | |
| 552 } | |
| 553 } | |
| 554 | |
| 555 _Bool nextSetpointChange(SDiveSettings* pDiveSettings, uint8_t depth_meter, uint8_t* change_depth_meter, char* setpoint) | |
| 556 { | |
| 557 uint8_t new_depth = 0; | |
| 558 char new_setpoint = 0; | |
| 559 for(int i = 1; i <= 5; i++) | |
| 560 { | |
| 561 if(pDiveSettings->setpoint[i].setpoint_cbar > 0 && pDiveSettings->setpoint[i].depth_meter > 0 ) | |
| 562 { | |
| 563 if( pDiveSettings->setpoint[i].depth_meter > new_depth && pDiveSettings->setpoint[i].depth_meter < depth_meter) | |
| 564 { | |
| 565 new_depth = pDiveSettings->setpoint[i].depth_meter; | |
| 566 new_setpoint = pDiveSettings->setpoint[i].setpoint_cbar; | |
| 567 } | |
| 568 } | |
| 569 } | |
| 570 if(new_depth) | |
| 571 { | |
| 572 * change_depth_meter = new_depth; | |
| 573 * setpoint = new_setpoint; | |
| 574 return 1; | |
| 575 } | |
| 576 return 0; | |
| 577 } | |
| 578 | |
| 579 | |
| 580 | |
| 581 void decom_CreateGasChangeList(SDiveSettings* pInput, const SLifeData* pLifeData) | |
| 582 { | |
| 583 int i=0, j = 0; | |
| 584 int count = 0; | |
| 585 for(i=0;i< 5;i++) | |
| 586 { | |
| 587 //FirstGas | |
| 588 | |
| 589 pInput->decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0; | |
| 590 pInput->decogaslist[i].GasIdInSettings = 255; | |
| 591 pInput->decogaslist[i].setPoint_cbar = 0; | |
| 592 pInput->decogaslist[i].helium_percentage = 0; | |
| 593 pInput->decogaslist[i].nitrogen_percentage = 0; | |
| 594 } | |
| 595 //pInput->liveData.dive_time_seconds = 0; | |
| 596 | |
| 597 /* FirstGas | |
| 598 * 0 = special gas, 1 to 5 ist OC gas, 6 to 10 is diluent | |
| 599 */ | |
| 600 | |
| 601 | |
| 602 | |
| 603 pInput->decogaslist[0] = pLifeData->actualGas; | |
| 604 | |
| 605 /* Add Deco Gases | |
| 606 * special (gasId == 0) is never a deco/travel gas but actual gas only | |
| 607 */ | |
| 608 if(pInput->diveMode == DIVEMODE_OC) | |
| 609 { | |
| 610 | |
| 611 for(i=1;i<= 5;i++) | |
| 612 { | |
| 613 if(pInput->gas[i].note.ub.active && pInput->gas[i].depth_meter | |
| 614 && (pLifeData->actualGas.GasIdInSettings != i) | |
| 615 &&(pInput->gas[i].depth_meter < pLifeData->depth_meter ) ) | |
| 616 { | |
| 617 count = 1; | |
| 618 for(j=1;j<= 5;j++) | |
| 619 { | |
| 620 if( (pInput->gas[j].note.ub.active && pInput->gas[j].depth_meter > 0) | |
| 621 && (pLifeData->actualGas.GasIdInSettings != j) // new hw 160905 | |
| 622 && (pInput->gas[j].depth_meter > pInput->gas[i].depth_meter)) | |
| 623 count++; | |
| 624 } | |
| 625 pInput->decogaslist[count].change_during_ascent_depth_meter_otherwise_zero = pInput->gas[i].depth_meter; | |
| 626 pInput->decogaslist[count].nitrogen_percentage = 100; | |
| 627 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].oxygen_percentage; | |
| 628 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].helium_percentage; | |
| 629 pInput->decogaslist[count].helium_percentage = pInput->gas[i].helium_percentage; | |
| 630 pInput->decogaslist[count].GasIdInSettings = i; | |
| 662 | 631 pInput->decogaslist[count].AppliedDiveMode = DIVEMODE_OC; |
| 38 | 632 } |
| 633 } | |
| 634 } | |
| 635 else | |
| 636 { | |
| 662 | 637 //divmode CCR or PSCR |
| 38 | 638 for(i=6; i <= 10; i++) |
| 639 { | |
| 640 if(pInput->gas[i].note.ub.active && pInput->gas[i].depth_meter | |
| 641 && (pLifeData->actualGas.GasIdInSettings != i) | |
| 642 &&(pInput->gas[i].depth_meter < pLifeData->depth_meter ) ) | |
| 643 { | |
| 644 count = 1; | |
| 645 for(j=6;j<= 10;j++) | |
| 646 { | |
| 647 // if(pInput->gas[j].note.ub.active && pInput->gas[j].depth_meter > 0 &&pInput->gas[j].depth_meter > pInput->gas[i].depth_meter) | |
| 648 if( (pInput->gas[j].note.ub.active && pInput->gas[j].depth_meter > 0) | |
| 649 && (pLifeData->actualGas.GasIdInSettings != j) // new hw 160905 | |
| 650 && (pInput->gas[j].depth_meter > pInput->gas[i].depth_meter)) | |
| 651 count++; | |
| 652 } | |
| 653 pInput->decogaslist[count].change_during_ascent_depth_meter_otherwise_zero = pInput->gas[i].depth_meter; | |
| 654 pInput->decogaslist[count].nitrogen_percentage = 100; | |
| 662 | 655 if(pInput->diveMode == DIVEMODE_PSCR) |
| 656 { | |
| 657 pInput->decogaslist[count].AppliedDiveMode = DIVEMODE_PSCR; | |
| 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; | |
| 659 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].oxygen_percentage; | |
| 660 } | |
| 661 else | |
| 662 { | |
| 663 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].oxygen_percentage; | |
| 664 } | |
| 38 | 665 pInput->decogaslist[count].nitrogen_percentage -= pInput->gas[i].helium_percentage; |
| 666 pInput->decogaslist[count].helium_percentage = pInput->gas[i].helium_percentage; | |
| 667 pInput->decogaslist[count].GasIdInSettings = i; | |
| 668 | |
| 669 } | |
| 670 } | |
| 662 | 671 if(pInput->diveMode == DIVEMODE_CCR) |
| 38 | 672 { |
| 662 | 673 /* Include Setpoint Changes */ |
| 674 for(j=0; j <= count; j++) | |
| 38 | 675 { |
| 662 | 676 uint8_t depth = 0; |
| 677 uint8_t changedepth = 0; | |
| 678 char newSetpoint; | |
| 679 | |
| 680 pInput->decogaslist[j].AppliedDiveMode = DIVEMODE_CCR; | |
| 681 if(j == 0) | |
| 682 { | |
| 683 depth = pLifeData->depth_meter; | |
| 684 } | |
| 685 else | |
| 686 { | |
| 687 //no setpointchange ? | |
| 688 pInput->decogaslist[j].setPoint_cbar = pInput->decogaslist[j - 1].setPoint_cbar; | |
| 689 depth = pInput->decogaslist[j].change_during_ascent_depth_meter_otherwise_zero + 0.1f; | |
| 690 } | |
| 691 /* Setpoint change at the same depth as gas changes */ | |
| 692 if(nextSetpointChange(pInput,depth + 1, &changedepth,&newSetpoint) && changedepth == depth) | |
| 38 | 693 { |
| 662 | 694 pInput->decogaslist[j].setPoint_cbar = newSetpoint; |
| 38 | 695 } |
| 662 | 696 /* Setpoint changes inbetween gas changes */ |
| 697 while(nextSetpointChange(pInput, depth, &changedepth,&newSetpoint) | |
| 698 && ( | |
| 699 ( (j < count) && (changedepth > pInput->decogaslist[j + 1].change_during_ascent_depth_meter_otherwise_zero)) | |
| 700 || ((j == count) && (changedepth > 0)) | |
| 701 )) | |
| 702 { | |
| 703 //Include new entry with setpoint change in decogaslist | |
| 704 for(int k = count; k > j; k--) | |
| 705 { | |
| 706 pInput->decogaslist[k+1] = pInput->decogaslist[k]; | |
| 707 } | |
| 708 pInput->decogaslist[j + 1] = pInput->decogaslist[j]; | |
| 709 pInput->decogaslist[j + 1].setPoint_cbar = newSetpoint; | |
| 710 j++; | |
| 711 count++; | |
| 712 depth = changedepth; | |
| 713 } | |
| 714 | |
| 38 | 715 } |
| 716 | |
| 717 } | |
| 718 } | |
| 719 } | |
| 720 void test_decom_CreateGasChangeList(void) | |
| 721 { | |
| 722 SDiveSettings diveSetting; | |
| 723 SLifeData lifeData; | |
| 724 lifeData.depth_meter = 100; | |
| 725 lifeData.actualGas.helium_percentage = 30; | |
| 726 lifeData.actualGas.nitrogen_percentage = 60; | |
| 727 lifeData.actualGas.setPoint_cbar = 18; | |
| 728 lifeData.actualGas.GasIdInSettings = 0; | |
| 729 lifeData.actualGas.change_during_ascent_depth_meter_otherwise_zero = 0; | |
| 730 diveSetting.diveMode = DIVEMODE_CCR; | |
| 731 diveSetting.gas[6].depth_meter = 0; | |
| 732 diveSetting.gas[6].helium_percentage = 30; | |
| 733 diveSetting.gas[6].oxygen_percentage = 10; | |
| 734 diveSetting.gas[6].note.ub.active = 1; | |
| 735 | |
| 736 diveSetting.gas[7].depth_meter = 60; | |
| 737 diveSetting.gas[7].helium_percentage = 0; | |
| 738 diveSetting.gas[7].oxygen_percentage = 10; | |
| 739 diveSetting.gas[7].note.ub.active = 1; | |
| 740 diveSetting.gas[8].note.ub.active = 0; | |
| 741 diveSetting.gas[9].note.ub.active = 0; | |
| 742 diveSetting.gas[10].note.ub.active = 0; | |
| 743 | |
| 744 diveSetting.setpoint[0].depth_meter = 0; | |
| 745 diveSetting.setpoint[1].depth_meter = 80; | |
| 746 diveSetting.setpoint[1].setpoint_cbar = 20; | |
| 747 diveSetting.setpoint[2].depth_meter = 60; | |
| 748 diveSetting.setpoint[2].setpoint_cbar = 25; | |
| 749 diveSetting.setpoint[3].depth_meter = 0; | |
| 750 diveSetting.setpoint[4].depth_meter = 0; | |
| 751 diveSetting.setpoint[5].depth_meter = 0; | |
| 752 | |
| 753 | |
| 754 decom_CreateGasChangeList(&diveSetting, &lifeData); | |
| 755 } | |
| 756 | |
| 757 uint8_t decom_tissue_test_tolerance(float* Tissue_nitrogen_bar, float* Tissue_helium_bar, float GF_value, float depth_in_bar_absolute) | |
| 758 { | |
| 759 float tissue_inertgas_saturation; | |
| 760 float inertgas_a; | |
| 761 float inertgas_b; | |
| 762 float inertgas_tolerance; | |
| 763 float gf_minus_1; | |
| 764 | |
| 765 gf_minus_1 = GF_value - 1.0f; | |
| 766 | |
| 767 for (int ci = 0; ci < 16; ci++) | |
| 768 { | |
| 769 if(Tissue_helium_bar[ci] == 0) | |
| 770 { | |
| 771 tissue_inertgas_saturation = Tissue_nitrogen_bar[ci]; | |
| 772 // | |
| 773 inertgas_a = buehlmann_N2_a[ci]; | |
| 774 inertgas_b = buehlmann_N2_b[ci]; | |
| 775 } | |
| 776 else | |
| 777 { | |
| 778 tissue_inertgas_saturation = Tissue_nitrogen_bar[ci] + Tissue_helium_bar[ci]; | |
| 779 // | |
| 780 inertgas_a = ( ( buehlmann_N2_a[ci] * Tissue_nitrogen_bar[ci]) + ( buehlmann_He_a[ci] * Tissue_helium_bar[ci]) ) / tissue_inertgas_saturation; | |
| 781 inertgas_b = ( ( buehlmann_N2_b[ci] * Tissue_nitrogen_bar[ci]) + ( buehlmann_He_b[ci] * Tissue_helium_bar[ci]) ) / tissue_inertgas_saturation; | |
| 782 } | |
| 783 // | |
| 784 inertgas_tolerance = ( (GF_value / inertgas_b - gf_minus_1) * depth_in_bar_absolute ) + ( GF_value * inertgas_a ); | |
| 785 // | |
| 786 if(inertgas_tolerance < tissue_inertgas_saturation) | |
| 787 return 0; | |
| 788 } | |
| 789 return 1; | |
| 790 } | |
| 791 | |
| 792 | |
| 793 void decom_tissues_desaturation_time(const SLifeData* pLifeData, SLifeData2* pOutput) | |
| 794 { | |
| 795 float pressure_in_gas_for_complete; | |
| 796 float pressure_in_gas_for_desat; | |
| 797 float diff_to_complete; | |
| 798 float diff_to_desatpoint; | |
| 799 float necessary_halftimes; | |
| 800 float desattime; | |
| 801 | |
| 802 pressure_in_gas_for_complete = 0.7902f * ( pLifeData->pressure_surface_bar - 0.0627f); | |
| 803 pressure_in_gas_for_desat = 1.05f * pressure_in_gas_for_complete; | |
| 804 for(int i=0; i<16; i++) | |
| 805 { | |
| 806 diff_to_complete = pressure_in_gas_for_complete - pLifeData->tissue_nitrogen_bar[i]; | |
| 807 diff_to_desatpoint = pressure_in_gas_for_desat - pLifeData->tissue_nitrogen_bar[i]; | |
| 808 | |
| 809 if((diff_to_desatpoint >= 0) || (diff_to_complete >= 0)) | |
| 810 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = 0; | |
| 811 else | |
| 812 { | |
| 813 necessary_halftimes = (logf(1.0f - (diff_to_desatpoint/diff_to_complete)) / -0.6931f); | |
| 814 desattime = buehlmann_N2_t_halflife[i] * necessary_halftimes; | |
| 815 if(desattime <= (float)0xFFFF) | |
| 816 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = desattime; | |
| 817 else | |
| 818 pOutput->tissue_nitrogen_desaturation_time_minutes[i] = 0xFFFF; | |
| 819 } | |
| 820 } | |
| 821 | |
| 822 for(int i=0; i<16; i++) | |
| 823 { | |
| 824 diff_to_desatpoint = 0.05f - pLifeData->tissue_helium_bar[i]; | |
| 825 diff_to_complete = -1.0f * pLifeData->tissue_helium_bar[i]; | |
| 826 | |
| 827 if((diff_to_desatpoint >= 0) || (diff_to_complete >= 0)) | |
| 828 pOutput->tissue_helium_desaturation_time_minutes[i] = 0; | |
| 829 else | |
| 830 { | |
| 831 necessary_halftimes = (logf(1.0f - (diff_to_desatpoint/diff_to_complete)) / -0.6931f); | |
| 832 desattime = buehlmann_He_t_halflife[i] * necessary_halftimes; | |
| 833 if(desattime <= (float)0xFFFF) | |
| 834 pOutput->tissue_helium_desaturation_time_minutes[i] = desattime; | |
| 835 else | |
| 836 pOutput->tissue_helium_desaturation_time_minutes[i] = 0xFFFF; | |
| 837 } | |
| 838 } | |
| 839 } | |
| 840 | |
| 841 #define MAX_DEGRADE_OTU_TIME_MINUTES (1440) | |
| 842 //CNS&OTU: | |
| 843 #define OXY_TEN_MINUTES_IN_SECONDS (600) | |
| 844 #define OXY_HALF_LIVE_OF_TEN_MINUTES__INVERSE_NINTH_ROOT_OF_TWO (0.92587471f) | |
| 845 #define OXY_NINE_DAYS_IN_TEN_MINUTES (1296) | |
| 846 #define OXY_ONE_SIXTIETH_PART (0.0166667f) | |
| 847 #define OXY_NEGATIVE_FIVE_SIXTH_PARTS (-0.8333333f) | |
| 848 void decom_oxygen_calculate_otu(float* oxygen_otu, float pressure_oxygen_real) | |
| 849 { | |
| 850 if(pressure_oxygen_real <= 0.5f) | |
| 851 return; | |
| 852 *oxygen_otu += (pow((double)(0.5f / (pressure_oxygen_real - 0.5f)),OXY_NEGATIVE_FIVE_SIXTH_PARTS)) * OXY_ONE_SIXTIETH_PART; | |
| 853 } | |
| 854 | |
| 855 void decom_oxygen_calculate_otu_degrade(float* oxygen_otu, long seconds_since_last_dive) | |
| 856 { | |
| 857 static long otu_time_ticker = 0; | |
| 858 static double otu_degrade_every_10_minutes = 999.9; | |
| 859 long cycles_since_last_call; | |
| 860 | |
| 861 if((*oxygen_otu <= 0) || (seconds_since_last_dive == 0)) | |
| 862 *oxygen_otu = 0; | |
| 863 else if(seconds_since_last_dive < OXY_TEN_MINUTES_IN_SECONDS) | |
| 864 { | |
| 865 otu_time_ticker = 1; | |
| 866 otu_degrade_every_10_minutes = *oxygen_otu / (MAX_DEGRADE_OTU_TIME_MINUTES / 10); | |
| 867 } | |
| 868 else | |
| 869 { | |
| 870 cycles_since_last_call = seconds_since_last_dive / (otu_time_ticker * OXY_TEN_MINUTES_IN_SECONDS); | |
| 871 *oxygen_otu -= ((double)cycles_since_last_call) * otu_degrade_every_10_minutes; | |
| 872 otu_time_ticker += cycles_since_last_call; | |
| 873 if((*oxygen_otu < 0) || (otu_time_ticker > (MAX_DEGRADE_OTU_TIME_MINUTES / 10))) | |
| 874 *oxygen_otu = 0; | |
| 875 } | |
| 876 } | |
| 877 | |
| 878 | |
| 879 | |
| 880 void decom_oxygen_calculate_cns_degrade(float* oxygen_cns, long seconds_since_last_dive) | |
| 881 { | |
| 882 static long cns_time_ticker = 0; | |
| 883 int cns_max_cycles; | |
| 884 | |
| 885 if((*oxygen_cns <= 0.5f) || (seconds_since_last_dive == 0)) | |
| 886 *oxygen_cns = 0; | |
| 887 else if(seconds_since_last_dive < OXY_TEN_MINUTES_IN_SECONDS) | |
| 888 cns_time_ticker = 1; | |
| 889 else | |
| 890 { | |
| 891 cns_max_cycles = OXY_NINE_DAYS_IN_TEN_MINUTES; | |
| 892 while((*oxygen_cns >= 0.5f) && ((cns_time_ticker * OXY_TEN_MINUTES_IN_SECONDS) < seconds_since_last_dive) && cns_max_cycles) | |
| 893 { | |
| 894 cns_time_ticker++; | |
| 895 cns_max_cycles--; | |
| 896 *oxygen_cns *= OXY_HALF_LIVE_OF_TEN_MINUTES__INVERSE_NINTH_ROOT_OF_TWO; | |
| 897 } | |
| 898 } | |
| 899 } | |
| 900 | |
| 901 | |
| 902 // new hwOS style | |
| 903 void decom_oxygen_calculate_cns(float* oxygen_cns, float pressure_oxygen_real) | |
| 904 { | |
| 905 uint8_t char_I_actual_ppO2; | |
| 906 float CNS_fraction = 0; | |
| 907 const float time_factor = 3000.0f; | |
| 908 | |
| 909 if(pressure_oxygen_real < 0.15f) | |
| 910 char_I_actual_ppO2 = 15; | |
| 911 else | |
| 912 if(pressure_oxygen_real >= 2.5f) | |
| 913 char_I_actual_ppO2 = 255; | |
| 914 else | |
| 915 char_I_actual_ppO2 = (uint8_t)(pressure_oxygen_real * 100); | |
| 916 | |
| 917 if (char_I_actual_ppO2 < 50) | |
| 918 (void)0; // no changes | |
| 919 //------------------------------------------------------------------------ | |
| 920 // Below (and including) 1.60 bar | |
| 921 else if (char_I_actual_ppO2 < 61) | |
| 922 CNS_fraction += time_factor/(-533.07f * char_I_actual_ppO2 + 54000.0f); | |
| 923 else if (char_I_actual_ppO2 < 71) | |
| 924 CNS_fraction += time_factor/(-444.22f * char_I_actual_ppO2 + 48600.0f); | |
| 925 else if (char_I_actual_ppO2 < 81) | |
| 926 CNS_fraction += time_factor/(-355.38f * char_I_actual_ppO2 + 42300.0f); | |
| 927 else if (char_I_actual_ppO2 < 91) | |
| 928 CNS_fraction += time_factor/(-266.53f * char_I_actual_ppO2 + 35100.0f); | |
| 929 else if (char_I_actual_ppO2 < 111) | |
| 930 CNS_fraction += time_factor/(-177.69f * char_I_actual_ppO2 + 27000.0f); | |
| 931 else if (char_I_actual_ppO2 < 152) | |
| 932 CNS_fraction += time_factor/( -88.84f * char_I_actual_ppO2 + 17100.0f); | |
| 933 else if (char_I_actual_ppO2 < 167) | |
| 934 CNS_fraction += time_factor/(-222.11f * char_I_actual_ppO2 + 37350.0f); | |
| 935 //------------------------------------------------------------------------ | |
| 936 // Arieli et all.(2002): Modeling pulmonary and CNS O2 toxicity: | |
| 937 // J Appl Physiol 92: 248--256, 2002, doi:10.1152/japplphysiol.00434.2001 | |
| 938 // Formula (A1) based on value for 1.55 and c=20 | |
| 939 // example calculation: Sqrt((1.7/1.55)^20)*0.000404 | |
| 940 else if (char_I_actual_ppO2 < 172) | |
| 941 CNS_fraction += time_factor*0.00102f; | |
| 942 else if (char_I_actual_ppO2 < 177) | |
| 943 CNS_fraction += time_factor*0.00136f; | |
| 944 else if (char_I_actual_ppO2 < 182) | |
| 945 CNS_fraction += time_factor*0.00180f; | |
| 946 else if (char_I_actual_ppO2 < 187) | |
| 947 CNS_fraction += time_factor*0.00237f; | |
| 948 else if (char_I_actual_ppO2 < 192) | |
| 949 CNS_fraction += time_factor*0.00310f; | |
| 950 else if (char_I_actual_ppO2 < 198) | |
| 951 CNS_fraction += time_factor*0.00401f; | |
| 952 else if (char_I_actual_ppO2 < 203) | |
| 953 CNS_fraction += time_factor*0.00517f; | |
| 954 else if (char_I_actual_ppO2 < 233) | |
| 955 CNS_fraction += time_factor*0.0209f; | |
| 956 else | |
| 957 CNS_fraction += time_factor*0.0482f; // value for 2.5 | |
| 958 | |
| 959 if( CNS_fraction > 999.0f) // Limit display to 999% | |
| 960 CNS_fraction = 999.0f; | |
| 961 if( CNS_fraction < 0.0f ) | |
| 962 CNS_fraction = 0.0f; | |
| 963 | |
| 964 //calculate cns for the actual ppo2 for 1 second | |
| 965 *oxygen_cns += OXY_ONE_SIXTIETH_PART * CNS_fraction; | |
| 966 | |
| 967 if( *oxygen_cns > 999.0f) // Limit display to 999% | |
| 968 *oxygen_cns = 999.0f; | |
| 969 if( *oxygen_cns < 0.0f ) | |
| 970 *oxygen_cns = 0.0f; | |
| 971 } | |
| 972 | |
| 973 /* old DR5 style | |
| 974 void decom_oxygen_calculate_cns(float* oxygen_cns, float pressure_oxygen_real) | |
| 975 { | |
| 976 int cns_no_range = 0; | |
| 977 _Bool not_found = 1; | |
| 978 //for the cns calculation | |
| 979 const float cns_ppo2_ranges[60][2] = { | |
| 980 {0.50f, 0.00f}, {0.60f, 0.14f}, {0.64f, 0.15f}, {0.66f, 0.16f}, {0.68f, 0.17f}, {0.70f, 0.18f}, | |
| 981 {0.74f, 0.19f}, {0.76f, 0.20f}, {0.78f, 0.21f}, {0.80f, 0.22f}, {0.82f, 0.23f}, {0.84f, 0.24f}, | |
| 982 {0.86f, 0.25f}, {0.88f, 0.26f}, {0.90f, 0.28f}, {0.92f, 0.29f}, {0.94f, 0.30f}, {0.96f, 0.31f}, | |
| 983 {0.98f, 0.32f}, {1.00f, 0.33f}, {1.02f, 0.35f}, {1.04f, 0.36f}, {1.06f, 0.38f}, {1.08f, 0.40f}, | |
| 984 {1.10f, 0.42f}, {1.12f, 0.43f}, {1.14f, 0.43f}, {1.16f, 0.44f}, {1.18f, 0.46f}, {1.20f, 0.47f}, | |
| 985 {1.22f, 0.48f}, {1.24f, 0.51f}, {1.26f, 0.52f}, {1.28f, 0.54f}, {1.30f, 0.56f}, {1.32f, 0.57f}, | |
| 986 {1.34f, 0.60f}, {1.36f, 0.62f}, {1.38f, 0.63f}, {1.40f, 0.65f}, {1.42f, 0.68f}, {1.44f, 0.71f}, | |
| 987 {1.46f, 0.74f}, {1.48f, 0.78f}, {1.50f, 0.83f}, {1.52f, 0.93f}, {1.54f, 1.04f}, {1.56f, 1.19f}, | |
| 988 {1.58f, 1.47f}, {1.60f, 2.22f}, {1.62f, 5.00f}, {1.65f, 6.25f}, {1.67f, 7.69f}, {1.70f, 10.0f}, | |
| 989 {1.72f,12.50f}, {1.74f,20.00f}, {1.77f,25.00f}, {1.79f,31.25f}, {1.80f,50.00f}, {1.82f,100.0f}}; | |
| 990 //find the correct cns range for the corresponding ppo2 | |
| 991 cns_no_range = 58; | |
| 992 while (cns_no_range && not_found) | |
| 993 { | |
| 994 if (pressure_oxygen_real > cns_ppo2_ranges[cns_no_range][0]) | |
| 995 { | |
| 996 cns_no_range++; | |
| 997 not_found = 0; | |
| 998 } | |
| 999 else | |
| 1000 cns_no_range--; | |
| 1001 } | |
| 1002 | |
| 1003 //calculate cns for the actual ppo2 for 1 second | |
| 1004 *oxygen_cns += OXY_ONE_SIXTIETH_PART * cns_ppo2_ranges[cns_no_range][1]; | |
| 1005 } | |
| 1006 */ | |
| 1007 | |
| 1008 void decom_oxygen_calculate_cns_exposure(int period_in_seconds, SGas* pActualGas, float pressure_ambient_bar, float* oxygen_cns) | |
| 1009 { | |
| 1010 float pressure_oxygen_real; | |
| 1011 float one_second_cns; | |
| 1012 | |
| 1013 pressure_oxygen_real = decom_calc_ppO2(pressure_ambient_bar, pActualGas); | |
| 1014 one_second_cns = 0; | |
| 1015 decom_oxygen_calculate_cns(&one_second_cns, pressure_oxygen_real); | |
| 1016 *oxygen_cns += one_second_cns * period_in_seconds; | |
| 1017 } | |
| 1018 | |
| 1019 | |
| 1020 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) | |
| 1021 { | |
| 1022 if(ending_ambient_pressure_bar == starting_ambient_pressure_bar) | |
| 1023 { | |
| 1024 decom_oxygen_calculate_cns_exposure(period_in_seconds, pGas, starting_ambient_pressure_bar, oxygen_cns); | |
| 1025 return; | |
| 1026 } | |
| 1027 | |
| 1028 float pressure_oxygen_real; | |
| 1029 float initial_pressure_oxygen; | |
| 1030 float ending_pressure_oxygen; | |
| 1031 float rate_oxygen; | |
| 1032 | |
| 1033 initial_pressure_oxygen = decom_calc_ppO2(starting_ambient_pressure_bar, pGas); | |
| 1034 ending_pressure_oxygen = decom_calc_ppO2(ending_ambient_pressure_bar, pGas); | |
| 1035 | |
| 1036 rate_oxygen = (ending_pressure_oxygen - initial_pressure_oxygen) / period_in_seconds; | |
| 1037 | |
| 1038 pressure_oxygen_real = initial_pressure_oxygen; | |
| 1039 for(int i = 0; i < period_in_seconds; i++) | |
| 1040 { | |
| 1041 decom_oxygen_calculate_cns(oxygen_cns, pressure_oxygen_real); | |
| 1042 pressure_oxygen_real += rate_oxygen; | |
| 1043 } | |
| 1044 } | |
| 1045 | |
| 1046 | |
| 1047 float decom_calc_ppO2(const float ambiant_pressure_bar, const SGas* pGas) | |
| 1048 { | |
| 662 | 1049 float percent_N2 = 0; |
| 38 | 1050 float percent_He = 0; |
| 1051 float percent_O2 = 0; | |
| 1052 | |
| 662 | 1053 decom_get_inert_gases(ambiant_pressure_bar, pGas, &percent_N2, &percent_He); |
| 1054 percent_O2 = 1 - percent_N2 - percent_He; | |
| 1055 | |
| 1056 return (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * percent_O2; | |
| 38 | 1057 } |
| 1058 | |
| 1059 | |
| 662 | 1060 float decom_calc_SimppO2(float ambiant_pressure_bar, const SGas* pGas) |
| 1061 { | |
| 1062 float o2Ratio = 0.0; | |
| 1063 float inertGasRatio = 0.0; | |
| 1064 float simulatedPSCRppo2 = 0.0; | |
| 1065 | |
| 1066 o2Ratio = (100.0 - pGas->nitrogen_percentage - pGas->helium_percentage) / 100.0; | |
| 1067 inertGasRatio = 1.0 - o2Ratio; | |
| 1068 simulatedPSCRppo2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * o2Ratio; | |
| 1069 simulatedPSCRppo2 -= (inertGasRatio * pGas->pscr_factor); | |
| 1070 if(simulatedPSCRppo2 < 0.0) | |
| 1071 { | |
| 1072 simulatedPSCRppo2 = 0.0; | |
| 1073 } | |
| 1074 return simulatedPSCRppo2; | |
| 1075 } | |
| 1076 | |
| 1077 float decom_calc_SimppO2_O2based(float ambiant_pressure_bar, uint8_t O2PerCent, float factor) | |
| 1078 { | |
| 1079 float o2Ratio = 0.0; | |
| 1080 float inertGasRatio = 0.0; | |
| 1081 float simulatedPSCRppo2 = 0.0; | |
| 1082 | |
| 1083 o2Ratio = O2PerCent / 100.0; | |
| 1084 inertGasRatio = 1.0 - o2Ratio; | |
| 1085 simulatedPSCRppo2 = (ambiant_pressure_bar - WATER_VAPOUR_PRESSURE) * o2Ratio; | |
| 1086 simulatedPSCRppo2 -= (inertGasRatio * factor); | |
| 1087 if(simulatedPSCRppo2 < 0.0) | |
| 1088 { | |
| 1089 simulatedPSCRppo2 = 0.0; | |
| 1090 } | |
| 1091 return simulatedPSCRppo2; | |
| 1092 } | |
| 1093 | |
| 38 | 1094 uint8_t decom_get_actual_deco_stop(SDiveState* pDiveState) |
| 1095 { | |
| 1096 SDecoinfo* pDecoinfo; | |
| 1097 uint8_t depthNext, depthLast, depthSecond, depthInc; | |
| 1098 if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE) | |
| 1099 pDecoinfo = &pDiveState->decolistBuehlmann; | |
| 1100 else | |
| 1101 pDecoinfo = &pDiveState->decolistVPM; | |
| 1102 | |
| 1103 depthLast = (uint8_t)(pDiveState->diveSettings.last_stop_depth_bar * 10); | |
| 1104 depthSecond = (uint8_t)(pDiveState->diveSettings.input_second_to_last_stop_depth_bar * 10); | |
| 1105 depthInc = (uint8_t)(pDiveState->diveSettings.input_next_stop_increment_depth_bar * 10); | |
| 1106 if(pDecoinfo->output_stop_length_seconds[0] > 0) | |
| 1107 { | |
| 1108 depthNext = depthLast; | |
| 1109 } | |
| 1110 else | |
| 1111 return 0; | |
| 1112 for(int i = DECOINFO_STRUCT_MAX_STOPS -1 ;i > 0; i--) | |
| 1113 { | |
| 1114 if(pDecoinfo->output_stop_length_seconds[i] > 0) | |
| 1115 { | |
| 1116 depthNext = depthSecond + ( (i - 1) * depthInc); | |
| 1117 break; | |
| 1118 } | |
| 1119 } | |
| 1120 return depthNext; | |
| 1121 } | |
| 1122 | |
| 1123 | |
| 1124 // =============================================================================== | |
| 1125 // decom_calc_desaturation_time | |
| 1126 /// @brief This code is used to calculate desat, calculated by RTE and send to Firmware | |
| 1127 /// similar but more technical in code than decom_tissues_desaturation_time() | |
| 1128 /// the later has 0.05 for helium in contrast to this one. | |
| 1129 /// This one goes down to 70%, the oterh | |
| 1130 /// | |
| 1131 /// output is desat time in minutes | |
| 1132 // =============================================================================== | |
| 1133 int decom_calc_desaturation_time(float* Tissue_nitrogen_bar, float* Tissue_helium_bar, float surface_pressure_bar) | |
| 1134 { | |
| 1135 const float N2_ratio = 0.7902; // FIXED sum as stated in b"uhlmann | |
| 1136 | |
| 1137 float pres_surface; | |
| 1138 float temp_atem; | |
| 1139 float float_desaturation_multiplier; | |
| 1140 float temp1,temp2,temp3,temp4; | |
| 1141 int ci; | |
| 1142 int int_temp; | |
| 1143 int int_O_desaturation_time; | |
| 1144 pres_surface = ((float)surface_pressure_bar); | |
| 1145 temp_atem = N2_ratio * (pres_surface - 0.0627f); | |
| 1146 | |
| 1147 int_O_desaturation_time = 0; | |
| 1148 float_desaturation_multiplier = 100 / 142.0f; // new in v.101 (70,42%/100.=142) | |
| 1149 | |
| 1150 for (ci=0;ci<16;ci++) | |
| 1151 { | |
| 1152 // saturation_time (for flight) and N2_saturation in multiples of halftime | |
| 1153 // version v.100: 1.1 = 10 percent distance to totally clean (totally clean is not possible, would take infinite time ) | |
| 1154 // new in version v.101: 1.07 = 7 percent distance to totally clean (totally clean is not possible, would take infinite time ) | |
| 1155 // changes in v.101: 1.05 = 5 percent dist to totally clean is new desaturation point for display and noFly calculations | |
| 1156 | |
| 1157 // N2 | |
| 1158 temp1 = 1.05f * temp_atem; | |
| 1159 temp1 = temp1 - (float)Tissue_nitrogen_bar[ci]; | |
| 1160 temp2 = temp_atem - (float)Tissue_nitrogen_bar[ci]; | |
| 1161 if (temp2 >= 0) | |
| 1162 { | |
| 1163 temp1 = 0; | |
| 1164 temp2 = 0; | |
| 1165 } | |
| 1166 else | |
| 1167 temp1 = temp1 / temp2; | |
| 1168 | |
| 1169 if (temp1 > 0) | |
| 1170 { | |
| 1171 temp1 = logf(1.0f - temp1); | |
| 1172 temp1 = temp1 / -0.6931f; // temp1 is the multiples of half times necessary. | |
| 1173 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested. | |
| 1174 // minus because log is negative | |
| 1175 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 | |
| 1176 } | |
| 1177 else | |
| 1178 { | |
| 1179 temp1 = 0; | |
| 1180 temp2 = 0; | |
| 1181 } | |
| 1182 | |
| 1183 // He | |
| 1184 temp3 = 0.1f - (float)Tissue_helium_bar[ci]; | |
| 1185 if (temp3 >= 0) | |
| 1186 { | |
| 1187 temp3 = 0; | |
| 1188 temp4 = 0; | |
| 1189 } | |
| 1190 else | |
| 1191 temp3 = -1.0f * temp3 / (float)Tissue_helium_bar[ci]; | |
| 1192 if (temp3 > 0) | |
| 1193 { | |
| 1194 temp3 = logf(1.0f - temp3); | |
| 1195 temp3 = temp3 / -0.6931f; // temp1 is the multiples of half times necessary. | |
| 1196 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested. | |
| 1197 // minus because log is negative | |
| 1198 temp4 = buehlmann_He_t_halflife[ci] * temp3 / float_desaturation_multiplier; // time necessary (in minutes ) for "complete" desaturation, new in v.101 float_desaturation_multiplier | |
| 1199 } | |
| 1200 else | |
| 1201 { | |
| 1202 temp3 = 0; | |
| 1203 temp4 = 0; | |
| 1204 } | |
| 1205 | |
| 1206 // saturation_time (for flight) | |
| 1207 if (temp4 > temp2) | |
| 1208 int_temp = (int)temp4; | |
| 1209 else | |
| 1210 int_temp = (int)temp2; | |
| 1211 if(int_temp > int_O_desaturation_time) | |
| 1212 int_O_desaturation_time = int_temp; | |
| 1213 | |
| 1214 /*// N2 saturation in multiples of halftime for display purposes | |
| 1215 temp2 = temp1 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8 | |
| 1216 temp2 = temp2 + 80.0; // set center | |
| 1217 if (temp2 < 0.0) | |
| 1218 temp2 = 0.0; | |
| 1219 if (temp2 > 255.0) | |
| 1220 temp2 = 255.0; | |
| 1221 U8_tissue_N2_saturation[ci] = (U8)temp2; | |
| 1222 // He saturation in multiples of halftime for display purposes | |
| 1223 temp4 = temp3 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8 | |
| 1224 temp4 = temp4 + 80.0; // set center | |
| 1225 if (temp4 < 0.0) | |
| 1226 temp4 = 0.0; | |
| 1227 if (temp4 > 255.0) | |
| 1228 temp4 = 255.0; | |
| 1229 U8_tissue_He_saturation[ci] = (char)temp4;*/ | |
| 1230 } | |
| 1231 | |
| 1232 return int_O_desaturation_time; | |
| 1233 } |