Mercurial > public > ostc4
annotate Discovery/Src/vpm.c @ 863:0c89c6fa949c Evo_2_23 tip
Bugfix empty line in deco plan (VPM only):
Floating numbers were used to calculate the target slot for the time entry of a deco stop. The float rounding caused a time to be written into one line above the intended one. In the next step the misplaced time was overwritten by the next shallower stop. To fix the problem the index calculation has been corrected and in addition digit numbers have generally been added to floating point operations to make the floating operation more visible.
| author | Ideenmodellierer |
|---|---|
| date | Tue, 02 Jul 2024 20:05:08 +0200 |
| parents | b7d93ff6b3b2 |
| children |
| rev | line source |
|---|---|
| 38 | 1 /////////////////////////////////////////////////////////////////////////////// |
| 2 /// -*- coding: UTF-8 -*- | |
| 3 /// | |
| 4 /// \file Discovery/Src/vpm.c | |
| 5 /// \brief critical_volume comment by hw | |
| 6 /// \author Heinrichs Weikamp, Erik C. Baker | |
| 7 /// \date 19-April-2014 | |
| 8 /// | |
| 9 /// \details | |
| 10 /// | |
| 11 /// $Id$ | |
| 12 /////////////////////////////////////////////////////////////////////////////// | |
| 13 /// \par Copyright (c) 2014-2018 Heinrichs Weikamp gmbh | |
| 14 /// | |
| 15 /// This program is free software: you can redistribute it and/or modify | |
| 16 /// it under the terms of the GNU General Public License as published by | |
| 17 /// the Free Software Foundation, either version 3 of the License, or | |
| 18 /// (at your option) any later version. | |
| 19 /// | |
| 20 /// This program is distributed in the hope that it will be useful, | |
| 21 /// but WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 22 /// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
| 23 /// GNU General Public License for more details. | |
| 24 /// | |
| 25 /// You should have received a copy of the GNU General Public License | |
| 26 /// along with this program. If not, see <http://www.gnu.org/licenses/>. | |
| 27 ////////////////////////////////////////////////////////////////////////////// | |
| 28 /// \par Varying Permeability Model (VPM) Decompression Program in c (converted from FORTRAN) | |
| 29 /// | |
| 30 /// Author: Erik C. Baker | |
| 31 /// | |
| 32 /// "DISTRIBUTE FREELY - CREDIT THE AUTHORS" | |
| 33 /// | |
| 34 /// This program extends the 1986 VPM algorithm (Yount & Hoffman) to include | |
| 35 /// mixed gas, repetitive, and altitude diving. Developments to the algorithm | |
| 36 /// were made by David E. Yount, Eric B. Maiken, and Erik C. Baker over a | |
| 37 /// period from 1999 to 2001. This work is dedicated in remembrance of | |
| 38 /// Professor David E. Yount who passed away on April 27, 2000. | |
| 39 /// | |
| 40 /// Notes: | |
| 41 /// 1. This program uses the sixteen (16) half-time compartments of the | |
| 42 /// Buhlmann ZH-L16 model. The optional Compartment 1b is used here with | |
| 43 /// half-times of 1.88 minutes for helium and 5.0 minutes for nitrogen. | |
| 44 /// | |
| 45 /// 2. This program uses various DEC, IBM, and Microsoft extensions which | |
| 46 /// may not be supported by all FORTRAN compilers. Comments are made with | |
| 47 /// a capital "C" in the first column or an exclamation point "!" placed | |
| 48 /// in a line after code. An asterisk "*" in column 6 is a continuation | |
| 49 /// of the previous line. All code, except for line numbers, starts in | |
| 50 /// column 7. | |
| 51 /// | |
| 52 /// 3. Comments and suggestions for improvements are welcome. Please | |
| 53 /// respond by e-mail to: EBaker@se.aeieng.com | |
| 54 /// | |
| 55 /// Acknowledgment: Thanks to Kurt Spaugh for recommendations on how to clean | |
| 56 /// up the code. | |
| 57 /// =============================================================================== | |
| 58 /// Converted to vpmdeco.c using f2c; R.McGinnis (CABER Swe) 5/01 | |
| 59 /// =============================================================================== | |
| 60 /// | |
| 61 /// ************************ Heirichs Weipkamp ************************************** | |
| 62 /// | |
| 63 /// The original Yount & Baker code has been adjusted for real life calculation. | |
| 64 /// | |
| 65 /// 1) The original main function has been split in several functions | |
| 66 /// | |
| 67 /// 2) When the deco zone is reached (while ascending) the gradient factors are kept fix | |
| 68 /// and critical volume algorithm is switched of. maxfirststopdepth is kept fix | |
| 69 /// to make shure Boeyls Law algorithm works correctly | |
| 70 /// | |
| 71 /// 4) gas_loadings_ascent_descend heeds all gaschanges and CCR support has been added | |
| 72 /// | |
| 73 | |
| 74 #include <stdio.h> | |
| 75 #include <stdlib.h> | |
| 76 #include <string.h> | |
| 77 #include <math.h> | |
| 78 #include <time.h> | |
| 79 | |
| 80 #include "vpm.h" | |
| 81 #include "decom.h" | |
| 82 | |
| 83 #define GAS_N2 0 | |
| 84 #define GAS_HE 1 | |
| 85 | |
| 290 | 86 static const _Bool buehlmannSafety = true; |
| 38 | 87 /* Common Block Declarations */ |
| 88 | |
| 89 extern const float SURFACE_TENSION_GAMMA; //!Adj. Range: 0.015 to 0.065 N/m | |
| 90 extern const float SKIN_COMPRESSION_GAMMAC; //!Adj. Range: 0.160 to 0.290 N/m | |
| 91 extern const float UNITS_FACTOR; | |
| 92 extern const float WATER_VAPOR_PRESSURE; // (Schreiner value) based on respiratory quotien | |
| 93 extern const float CRIT_VOLUME_PARAMETER_LAMBDA; //!Adj. Range: 6500 to 8300 fsw-min | |
| 290 | 94 //extern const float GRADIENT_ONSET_OF_IMPERM_ATM; //!Adj. Range: 5.0 to 10.0 atm |
| 38 | 95 extern const float REGENERATION_TIME_CONSTANT; //!Adj. Range: 10080 to 51840 min |
| 290 | 96 //extern const float PRESSURE_OTHER_GASES_MMHG; //!Constant value for PO2 up to 2 atm |
| 38 | 97 extern const float CONSTANT_PRESSURE_OTHER_GASES; // PRESSURE_OTHER_GASES_MMHG / 760. * UNITS_FACTOR; |
| 98 | |
| 99 extern const float HELIUM_TIME_CONSTANT[]; | |
| 100 extern const float NITROGEN_TIME_CONSTANT[]; | |
| 101 | |
| 290 | 102 static float minimum_deco_stop_time; |
| 103 static float run_time, run_time_first_stop; | |
| 104 static float segment_time; | |
| 105 static short mix_number; | |
| 106 static float barometric_pressure; | |
| 107 static _Bool altitude_dive_algorithm_off; | |
| 108 static _Bool units_equal_fsw, units_equal_msw; | |
| 38 | 109 |
| 110 /* by hw 11.06.2015 to allow */ | |
| 290 | 111 static float gCNS_VPM; |
| 38 | 112 |
| 290 | 113 static float helium_pressure[16], nitrogen_pressure[16]; |
| 114 static float surface_phase_volume_time[16]; | |
| 115 static float regenerated_radius_he[16], regenerated_radius_n2[16]; | |
| 116 static float allowable_gradient_he[16], allowable_gradient_n2[16]; | |
| 38 | 117 |
| 118 //_Bool deco_zone_reached; | |
| 290 | 119 static _Bool critical_volume_algorithm_off; |
| 120 static float max_first_stop_depth; | |
| 121 static float max_deco_ceiling_depth; | |
| 38 | 122 //Boylslaw compensation |
| 290 | 123 static float deco_gradient_he[16]; |
| 124 static float deco_gradient_n2[16]; | |
| 125 static int vpm_calc_what; | |
| 126 static int count_critical_volume_iteration; | |
| 127 static short number_of_changes; | |
| 128 static float depth_change[11]; | |
| 129 static float step_size_change[11]; | |
| 130 static float rate_change[11]; | |
| 131 static short mix_change[11]; | |
| 38 | 132 |
| 290 | 133 static const _Bool vpm_b = true; |
| 38 | 134 |
| 135 extern const float float_buehlmann_N2_factor_expositon_20_seconds[]; | |
| 136 extern const float float_buehlmann_He_factor_expositon_20_seconds[]; | |
| 137 extern const float float_buehlmann_N2_factor_expositon_one_minute[]; | |
| 138 extern const float float_buehlmann_He_factor_expositon_one_minute[]; | |
| 139 extern const float float_buehlmann_N2_factor_expositon_five_minutes[]; | |
| 140 extern const float float_buehlmann_He_factor_expositon_five_minutes[]; | |
| 141 extern const float float_buehlmann_N2_factor_expositon_one_hour[]; | |
| 142 extern const float float_buehlmann_He_factor_expositon_one_hour[]; | |
| 143 | |
| 290 | 144 static float depth_start_of_deco_calc; |
| 145 static float depth_start_of_deco_zone; | |
| 146 static float first_stop_depth; | |
| 147 static float run_time_start_of_deco_zone; | |
| 38 | 148 |
| 290 | 149 static float r_nint(float *x); |
| 150 static float r_int(float *x); | |
| 151 static _Bool nullzeit_unter60; | |
| 152 static int vpm_calc_status; | |
| 153 static _Bool buehlmann_wait_exceeded = false; | |
| 38 | 154 |
| 290 | 155 static SLifeData* pInput = NULL; |
| 156 static SVpm* pVpm = NULL; | |
| 157 static SDecoinfo* pDecoInfo = NULL; | |
| 158 static SDiveSettings* pDiveSettings = NULL; | |
| 159 | |
| 160 static float r_nint(float *x) | |
| 38 | 161 { |
| 162 return( (*x)>=0 ? | |
| 163 floorf(*x + 0.5f) : -floorf(0.5f - *x) ); | |
| 164 } | |
| 165 | |
| 290 | 166 static float r_int(float *x) |
| 38 | 167 { |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
168 return( (*x>0.0) ? floorf(*x) : -floorf(- *x) ); |
| 38 | 169 } |
| 170 | |
| 171 /** private functions | |
| 172 */ | |
| 290 | 173 extern int radius_root_finder (float *a, float *b, float *c,float *low_bound, float *high_bound, float *ending_radius); |
| 38 | 174 |
| 290 | 175 static int nuclear_regeneration(float *dive_time);// clock_(); |
| 176 static int calc_deco_ceiling(float *deco_ceiling_depth,_Bool fallowablw); | |
| 177 static int critical_volume(float *deco_phase_volume_time); ; | |
| 178 static int calc_start_of_deco_zone(float *starting_depth, float *rate, float *depth_start_of_deco_zone); | |
| 179 static int calc_initial_allowable_gradient(void); | |
| 180 static void decompression_stop(float *deco_stop_depth, float *step_size, _Bool final_deco_calculation); | |
| 181 static int gas_loadings_ascent_descen(float* helium_pressure, float* nitrogen_pressure, float starting_depth,float ending_depth, float rate,_Bool check_gas_change); | |
| 182 static int calc_surface_phase_volume_time(void); | |
| 183 static int calc_max_actual_gradient(float *deco_stop_depth); | |
| 184 static int projected_ascent(float *starting_depth, float *rate, float *deco_stop_depth, float *step_size); | |
| 185 static void vpm_calc_deco(void); | |
| 186 static int vpm_calc_critcal_volume(_Bool begin,_Bool calc_nulltime); | |
| 187 static int vpm_check_converged(_Bool calc_nulltime); | |
| 188 static int vpm_calc_final_deco(_Bool begin); | |
| 189 static void BOYLES_LAW_COMPENSATION (float* First_Stop_Depth,float * Deco_Stop_Depth,float* Step_Size); | |
| 292 | 190 static int vpm_calc_ndl(void); |
| 290 | 191 static void vpm_init_1(void); |
| 192 static void vpm_calc_deco_ceiling(void); | |
| 38 | 193 |
| 290 | 194 static void vpm_init_1(void) |
| 38 | 195 { |
| 196 units_equal_msw = true; | |
| 197 units_equal_fsw = false; | |
| 198 altitude_dive_algorithm_off= true; //!Options: ON or OFF | |
| 199 minimum_deco_stop_time=1.0; //!Options: float positive number | |
| 200 critical_volume_algorithm_off= false; //!Options: ON or OFF | |
| 201 run_time = 0.; | |
| 202 //barometric_pressure = dive_data.surface * 10; | |
| 203 | |
| 204 //mix_number = dive_data.selected_gas + 1; | |
| 205 | |
| 206 max_first_stop_depth = 0; | |
| 207 max_deco_ceiling_depth = 0; | |
| 208 //deco_zone_reached = false; | |
| 209 depth_start_of_deco_calc = 0; | |
| 210 depth_start_of_deco_zone = 0; | |
| 211 first_stop_depth = 0; | |
| 212 run_time_start_of_deco_zone = 0; | |
| 213 | |
| 214 gCNS_VPM = 0; | |
| 215 } | |
| 216 | |
| 217 float vpm_get_CNS(void) | |
| 218 { | |
| 219 return gCNS_VPM; | |
| 220 } | |
| 221 | |
| 222 int vpm_calc(SLifeData* pINPUT, | |
| 223 SDiveSettings* pSettings, | |
| 224 SVpm* pVPM, | |
| 225 SDecoinfo* | |
| 226 pDECOINFO, | |
| 227 int calc_what) | |
| 228 { | |
| 229 vpm_init_1(); | |
| 230 //decom_CreateGasChangeList(pSettings, pINPUT); | |
| 231 vpm_calc_what = calc_what; | |
| 232 /**clear decoInfo*/ | |
| 233 pDECOINFO->output_time_to_surface_seconds = 0; | |
| 234 pDECOINFO->output_ndl_seconds = 0; | |
| 235 pDECOINFO->output_ceiling_meter = 0; | |
|
247
3949781096d4
feature: Relative GF to Saturation renames
Jan Mulder <jlmulder@xs4all.nl>
parents:
149
diff
changeset
|
236 pDECOINFO->super_saturation = 0; |
| 38 | 237 uint8_t tmp_calc_status; |
| 238 for(int i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 239 { | |
| 240 pDECOINFO->output_stop_length_seconds[i] = 0; | |
| 241 } | |
| 242 | |
|
305
305f251cc981
bugfix, consistency: show deco/NDL really after 1 minute divetime
Jan Mulder <jlmulder@xs4all.nl>
parents:
293
diff
changeset
|
243 if(pINPUT->dive_time_seconds_without_surface_time < 60) |
| 38 | 244 { |
| 292 | 245 vpm_calc_status = CALC_NDL; |
| 38 | 246 return vpm_calc_status; |
| 247 } | |
| 248 pVpm = pVPM; | |
| 249 pInput = pINPUT; | |
| 250 pDecoInfo = pDECOINFO; | |
| 251 pDiveSettings = pSettings; | |
| 252 | |
| 292 | 253 if(vpm_calc_status == CALC_NDL) |
| 38 | 254 { |
| 292 | 255 tmp_calc_status = vpm_calc_ndl(); |
| 38 | 256 } |
| 257 else | |
| 258 { | |
| 259 tmp_calc_status = CALC_BEGIN; | |
| 260 } | |
| 261 //Normal Deco calculation | |
| 292 | 262 if(tmp_calc_status != CALC_NDL) |
| 38 | 263 { |
| 264 max_first_stop_depth = pVpm->max_first_stop_depth_save; | |
| 265 run_time_start_of_deco_zone = pVpm->run_time_start_of_deco_zone_save; | |
| 266 depth_start_of_deco_zone = pVpm->depth_start_of_deco_zone_save; | |
| 267 for (int i = 0; i < 16; ++i) { | |
| 268 helium_pressure[i] = pInput->tissue_helium_bar[i] * 10; | |
| 269 nitrogen_pressure[i] = pInput->tissue_nitrogen_bar[i] * 10; | |
| 270 } | |
| 271 vpm_calc_deco(); | |
| 272 tmp_calc_status = vpm_calc_critcal_volume(true,false); | |
| 273 if(vpm_calc_what == DECOSTOPS) | |
| 274 { | |
| 275 pVpm->max_first_stop_depth_save = max_first_stop_depth; | |
| 276 pVpm->run_time_start_of_deco_zone_save = run_time_start_of_deco_zone; | |
| 277 pVpm->depth_start_of_deco_zone_save = depth_start_of_deco_zone; | |
| 278 } | |
| 279 } | |
| 280 | |
| 281 //Only Decostops not futute stops | |
| 282 if(vpm_calc_what == DECOSTOPS) | |
| 283 vpm_calc_status = tmp_calc_status; | |
| 284 return vpm_calc_status; | |
| 285 } | |
| 286 | |
| 287 void vpm_saturation_after_ascent(SLifeData* input) | |
| 288 { | |
| 289 int i = 0; | |
| 290 for (i = 0; i < 16; ++i) { | |
| 291 pInput->tissue_helium_bar[i] = helium_pressure[i] / 10; | |
| 292 pInput->tissue_nitrogen_bar[i] = nitrogen_pressure[i] / 10; | |
| 293 } | |
| 294 pInput->pressure_ambient_bar = pInput->pressure_surface_bar; | |
| 295 } | |
| 296 /* =============================================================================== */ | |
| 297 /* NOTE ABOUT PRESSURE UNITS USED IN CALCULATIONS: */ | |
| 298 /* It is the convention in decompression calculations to compute all gas */ | |
| 299 /* loadings, absolute pressures, partial pressures, etc., in the units of */ | |
| 300 /* depth pressure that you are diving - either feet of seawater (fsw) or */ | |
| 301 /* meters of seawater (msw). This program follows that convention with the */ | |
| 302 /* the exception that all VPM calculations are performed in SI units (by */ | |
| 303 /* necessity). Accordingly, there are several conversions back and forth */ | |
| 304 /* between the diving pressure units and the SI units. */ | |
| 305 /* =============================================================================== */ | |
| 306 /* =============================================================================== */ | |
| 307 /* FUNCTION SUBPROGRAM FOR GAS LOADING CALCULATIONS - ASCENT AND DESCENT */ | |
| 308 /* =============================================================================== */ | |
| 309 | |
| 310 | |
| 311 | |
| 312 /* =============================================================================== */ | |
| 313 /* SUBROUTINE GAS_LOADINGS_ASCENT_DESCENT */ | |
| 314 /* Purpose: This subprogram applies the Schreiner equation to update the */ | |
| 315 /* gas loadings (partial pressures of helium and nitrogen) in the half-time */ | |
| 316 /* compartments due to a linear ascent or descent segment at a constant rate. */ | |
| 317 /* =============================================================================== */ | |
| 318 | |
| 290 | 319 static int gas_loadings_ascent_descen(float* helium_pressure, |
| 38 | 320 float* nitrogen_pressure, |
| 321 float starting_depth, | |
| 322 float ending_depth, | |
| 323 float rate,_Bool check_gas_change) | |
| 324 { | |
| 325 short i; | |
| 326 float initial_inspired_n2_pressure, | |
| 327 initial_inspired_he_pressure, nitrogen_rate, | |
| 328 last_run_time, | |
| 329 starting_ambient_pressure, | |
| 330 ending_ambient_pressure; | |
| 331 float initial_helium_pressure[16]; | |
| 332 float initial_nitrogen_pressure[16]; | |
| 333 float helium_rate; | |
| 334 float fraction_helium_begin; | |
| 335 float fraction_helium_end; | |
| 336 float fraction_nitrogen_begin; | |
| 337 float fraction_nitrogen_end; | |
| 338 float ending_depth_tmp = ending_depth; | |
| 339 float segment_time_tmp = 0; | |
| 340 /* loop */ | |
| 341 /* =============================================================================== */ | |
| 342 /* CALCULATIONS */ | |
| 343 /* =============================================================================== */ | |
| 344 segment_time = (ending_depth_tmp - starting_depth) / rate; | |
| 345 last_run_time = run_time; | |
| 346 run_time = last_run_time + segment_time; | |
| 347 do { | |
| 348 ending_depth_tmp = ending_depth; | |
| 349 if (starting_depth > ending_depth && check_gas_change && number_of_changes > 1) | |
| 350 { | |
| 351 for (i = 1; i < number_of_changes; ++i) | |
| 352 { | |
| 353 if (depth_change[i] < starting_depth && depth_change[i] > ending_depth) | |
| 354 { | |
| 355 ending_depth_tmp = depth_change[i]; | |
| 356 break; | |
| 357 } | |
| 358 } | |
| 359 for (i = 1; i < number_of_changes; ++i) | |
| 360 { | |
| 361 if (depth_change[i] >= starting_depth) | |
| 362 { | |
| 363 mix_number = mix_change[i]; | |
| 364 } | |
| 365 } | |
| 366 } | |
| 367 segment_time_tmp = (ending_depth_tmp - starting_depth) / rate; | |
| 368 ending_ambient_pressure = ending_depth_tmp + barometric_pressure; | |
| 369 starting_ambient_pressure = starting_depth + barometric_pressure; | |
| 370 decom_get_inert_gases( starting_ambient_pressure / 10, (&pDiveSettings->decogaslist[mix_number]), &fraction_nitrogen_begin, &fraction_helium_begin ); | |
| 371 decom_get_inert_gases( ending_ambient_pressure / 10, (&pDiveSettings->decogaslist[mix_number]), &fraction_nitrogen_end, &fraction_helium_end ); | |
| 372 | |
| 373 initial_inspired_he_pressure = (starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_helium_begin; | |
| 374 initial_inspired_n2_pressure = (starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_nitrogen_begin; | |
| 375 //helium_rate = *rate * fraction_helium[mix_number - 1]; | |
| 376 helium_rate = ((ending_ambient_pressure - WATER_VAPOR_PRESSURE)* fraction_helium_end - initial_inspired_he_pressure)/segment_time_tmp; | |
| 377 //nitrogen_rate2 = *rate * fraction_nitrogen[mix_number - 1]; | |
| 378 nitrogen_rate = ((ending_ambient_pressure - WATER_VAPOR_PRESSURE)* fraction_nitrogen_end - initial_inspired_n2_pressure)/segment_time_tmp; | |
| 379 | |
| 380 | |
| 381 decom_oxygen_calculate_cns_stage_SchreinerStyle(segment_time_tmp,&pDiveSettings->decogaslist[mix_number],starting_ambient_pressure/10,ending_ambient_pressure/10,&gCNS_VPM); | |
| 382 //if(fabs(nitrogen_rate - nitrogen_rate2) > 0.000001) | |
| 383 //return -2; | |
| 384 for (i = 1; i <= 16; ++i) | |
| 385 { | |
| 386 initial_helium_pressure[i - 1] = helium_pressure[i - 1]; | |
| 387 initial_nitrogen_pressure[i - 1] = nitrogen_pressure[i - 1]; | |
| 388 helium_pressure[i - 1] = | |
| 389 schreiner_equation__2(&initial_inspired_he_pressure, | |
| 390 &helium_rate, | |
| 391 &segment_time, | |
| 392 &HELIUM_TIME_CONSTANT[i - 1], | |
| 393 &initial_helium_pressure[i - 1]); | |
| 394 nitrogen_pressure[i - 1] = | |
| 395 schreiner_equation__2(&initial_inspired_n2_pressure, | |
| 396 &nitrogen_rate, | |
| 397 &segment_time, | |
| 398 &NITROGEN_TIME_CONSTANT[i - 1], | |
| 399 &initial_nitrogen_pressure[i - 1]); | |
| 400 | |
| 401 //nextround??? | |
| 402 | |
| 403 } | |
| 404 starting_depth = ending_depth_tmp; | |
| 405 } while(ending_depth_tmp > ending_depth); | |
| 406 | |
| 407 return 0; | |
| 408 } /* gas_loadings_ascent_descen */ | |
| 409 | |
| 290 | 410 static float last_phase_volume_time[16]; |
| 411 static float n2_pressure_start_of_deco_zone[16]; | |
| 412 static float he_pressure_start_of_deco_zone[16]; | |
| 413 static float phase_volume_time[16]; | |
| 414 static float n2_pressure_start_of_ascent[16]; | |
| 415 static float he_pressure_start_of_ascent[16]; | |
| 416 static float run_time_start_of_deco_calc; | |
| 417 static float starting_depth; | |
| 418 static float last_run_time; | |
| 419 static float deco_phase_volume_time; | |
| 420 static float run_time_start_of_ascent; | |
| 421 static float rate; | |
| 422 static float step_size; | |
| 423 static _Bool vpm_violates_buehlmann; | |
| 38 | 424 |
| 290 | 425 static void vpm_calc_deco(void) |
| 38 | 426 { |
| 427 /* System generated locals */ | |
| 428 | |
| 429 //float deepest_possible_stop_depth; | |
| 430 // altitude_of_dive, | |
| 431 short i; | |
| 432 int j = 0; | |
| 433 | |
| 434 // float rounding_operation; | |
| 435 | |
| 436 /* =============================================================================== */ | |
| 437 /* INPUT PARAMETERS TO BE USED FOR STAGED DECOMPRESSION AND SAVE IN ARRAYS. */ | |
| 438 /* ASSIGN INITAL PARAMETERS TO BE USED AT START OF ASCENT */ | |
| 439 /* The user has the ability to change mix, ascent rate, and step size in any */ | |
| 440 /* combination at any depth during the ascent. */ | |
| 441 /* =============================================================================== */ | |
| 442 | |
| 443 run_time = ((float)pInput->dive_time_seconds )/ 60; | |
| 444 count_critical_volume_iteration = 0; | |
| 445 number_of_changes = 1; | |
| 446 | |
| 447 barometric_pressure = pInput->pressure_surface_bar * 10; | |
| 448 depth_change[0] =(pInput->pressure_ambient_bar - pInput->pressure_surface_bar)* 10; | |
| 449 mix_change[0] = 0; | |
| 450 rate_change[0 ] = -10;// neu 160215 hw, zuvor: -12; | |
| 451 step_size_change[0] = 3; | |
| 452 vpm_violates_buehlmann = false; | |
| 453 | |
| 454 for (i = 1; i < BUEHLMANN_STRUCT_MAX_GASES; i++) | |
| 455 { | |
| 456 depth_change[i] = 0; | |
| 457 mix_change[i] = 0; | |
| 458 } | |
| 459 j = 0; | |
| 460 | |
| 461 for (i = 1; i < BUEHLMANN_STRUCT_MAX_GASES; i++) | |
| 462 { | |
|
830
b7d93ff6b3b2
Added selection if an active gas shall be used for deco calculation or not:
Ideenmodellierer
parents:
305
diff
changeset
|
463 if((pDiveSettings->decogaslist[i].change_during_ascent_depth_meter_otherwise_zero >= depth_change[0] + 1) |
|
b7d93ff6b3b2
Added selection if an active gas shall be used for deco calculation or not:
Ideenmodellierer
parents:
305
diff
changeset
|
464 && (pDiveSettings->gas[pDiveSettings->decogaslist[i].GasIdInSettings].note.ub.decocalc)) |
| 38 | 465 continue; |
| 466 | |
|
830
b7d93ff6b3b2
Added selection if an active gas shall be used for deco calculation or not:
Ideenmodellierer
parents:
305
diff
changeset
|
467 if((pDiveSettings->decogaslist[i].change_during_ascent_depth_meter_otherwise_zero <= 0) |
|
b7d93ff6b3b2
Added selection if an active gas shall be used for deco calculation or not:
Ideenmodellierer
parents:
305
diff
changeset
|
468 || (pDiveSettings->gas[pDiveSettings->decogaslist[i].GasIdInSettings].note.ub.decocalc == 0)) |
| 38 | 469 break; |
| 470 | |
| 471 j++; | |
| 472 number_of_changes ++; | |
| 473 depth_change[j] = pDiveSettings->decogaslist[i].change_during_ascent_depth_meter_otherwise_zero ; | |
| 474 mix_change[j] = i; | |
| 475 rate_change[j] = -10;// neu 160215 hw, zuvor: -12; | |
| 476 step_size_change[j] = 3; | |
| 477 } | |
| 478 | |
| 479 starting_depth = depth_change[0] ; | |
| 480 mix_number = mix_change[0] ; | |
| 481 rate = rate_change[0]; | |
| 482 step_size = step_size_change[0]; | |
| 483 | |
| 484 for (i = 0; i < 16; ++i) { | |
| 485 he_pressure_start_of_ascent[i ] = helium_pressure[i]; | |
| 486 n2_pressure_start_of_ascent[i] = nitrogen_pressure[i]; | |
| 487 } | |
| 488 run_time_start_of_ascent = run_time; | |
| 489 if(starting_depth <= depth_start_of_deco_zone && vpm_calc_what == DECOSTOPS) | |
| 490 { | |
| 491 pVpm->deco_zone_reached = true; | |
| 492 depth_start_of_deco_calc = starting_depth; | |
| 493 critical_volume_algorithm_off = true; | |
| 494 } | |
| 495 else | |
| 496 { | |
| 497 //if(deco_zone_reached) | |
| 498 //{ | |
| 499 pVpm->deco_zone_reached = false; | |
| 500 critical_volume_algorithm_off = false; | |
| 501 //max_first_stop_depth = 0; | |
| 502 //max_first_stop_depth_save = 0; | |
| 503 //} | |
| 504 /* =============================================================================== */ | |
| 505 /* BEGIN PROCESS OF ASCENT AND DECOMPRESSION */ | |
| 506 /* First, calculate the regeneration of critical radii that takes place over */ | |
| 507 /* the dive time. The regeneration time constant has a time scale of weeks */ | |
| 508 /* so this will have very little impact on dives of normal length, but will */ | |
| 509 /* have major impact for saturation dives. */ | |
| 510 /* =============================================================================== */ | |
| 511 | |
| 512 nuclear_regeneration(&run_time); | |
| 513 | |
| 514 /* =============================================================================== */ | |
| 515 /* CALCULATE INITIAL ALLOWABLE GRADIENTS FOR ASCENT */ | |
| 516 /* This is based on the maximum effective crushing pressure on critical radii */ | |
| 517 /* in each compartment achieved during the dive profile. */ | |
| 518 /* =============================================================================== */ | |
| 519 | |
| 520 calc_initial_allowable_gradient(); | |
| 521 | |
| 522 /* =============================================================================== */ | |
| 523 /* SAVE VARIABLES AT START OF ASCENT (END OF BOTTOM TIME) SINCE THESE WILL */ | |
| 524 /* BE USED LATER TO COMPUTE THE FINAL ASCENT PROFILE THAT IS WRITTEN TO THE */ | |
| 525 /* OUTPUT FILE. */ | |
| 526 /* The VPM uses an iterative process to compute decompression schedules so */ | |
| 527 /* there will be more than one pass through the decompression loop. */ | |
| 528 /* =============================================================================== */ | |
| 529 | |
| 530 /* =============================================================================== */ | |
| 531 /* CALCULATE THE DEPTH WHERE THE DECOMPRESSION ZONE BEGINS FOR THIS PROFILE */ | |
| 532 /* BASED ON THE INITIAL ASCENT PARAMETERS AND WRITE THE DEEPEST POSSIBLE */ | |
| 533 /* DECOMPRESSION STOP DEPTH TO THE OUTPUT FILE */ | |
| 534 /* Knowing where the decompression zone starts is very important. Below */ | |
| 535 /* that depth there is no possibility for bubble formation because there */ | |
| 536 /* will be no supersaturation gradients. Deco stops should never start */ | |
| 537 /* below the deco zone. The deepest possible stop deco stop depth is */ | |
| 538 /* defined as the next "standard" stop depth above the point where the */ | |
| 539 /* leading compartment enters the deco zone. Thus, the program will not */ | |
| 540 /* base this calculation on step sizes larger than 10 fsw or 3 msw. The */ | |
| 541 /* deepest possible stop depth is not used in the program, per se, rather */ | |
| 542 /* it is information to tell the diver where to start putting on the brakes */ | |
| 543 /* during ascent. This should be prominently displayed by any deco program. */ | |
| 544 /* =============================================================================== */ | |
| 545 | |
| 546 calc_start_of_deco_zone(&starting_depth, &rate, &depth_start_of_deco_zone); | |
| 547 /* =============================================================================== */ | |
| 548 /* TEMPORARILY ASCEND PROFILE TO THE START OF THE DECOMPRESSION ZONE, SAVE */ | |
| 549 /* VARIABLES AT THIS POINT, AND INITIALIZE VARIABLES FOR CRITICAL VOLUME LOOP */ | |
| 550 /* The iterative process of the VPM Critical Volume Algorithm will operate */ | |
| 551 /* only in the decompression zone since it deals with excess gas volume */ | |
| 552 /* released as a result of supersaturation gradients (not possible below the */ | |
| 553 /* decompression zone). */ | |
| 554 /* =============================================================================== */ | |
| 555 gas_loadings_ascent_descen(helium_pressure,nitrogen_pressure, starting_depth, depth_start_of_deco_zone, rate, true); | |
| 556 | |
| 557 run_time_start_of_deco_zone = run_time; | |
| 558 depth_start_of_deco_calc = depth_start_of_deco_zone; | |
| 559 | |
| 560 for (i = 0; i < 16; ++i) | |
| 561 { | |
| 562 pVpm->max_actual_gradient[i] = 0.; | |
| 563 } | |
| 564 } | |
| 565 | |
| 566 for (i = 0; i < 16; ++i) | |
| 567 { | |
| 568 surface_phase_volume_time[i] = 0.; | |
| 569 last_phase_volume_time[i] = 0.; | |
| 570 he_pressure_start_of_deco_zone[i] = helium_pressure[i]; | |
| 571 n2_pressure_start_of_deco_zone[i] = nitrogen_pressure[i]; | |
| 572 //pVpm->max_actual_gradient[i] = 0.; | |
| 573 } | |
| 574 run_time_start_of_deco_calc = run_time; | |
| 575 } | |
| 576 /* =============================================================================== */ | |
| 577 /* START OF CRITICAL VOLUME LOOP */ | |
| 578 /* This loop operates between Lines 50 and 100. If the Critical Volume */ | |
| 579 /* Algorithm is toggled "off" in the program settings, there will only be */ | |
| 580 /* one pass through this loop. Otherwise, there will be two or more passes */ | |
| 581 /* through this loop until the deco schedule is "converged" - that is when a */ | |
| 582 /* comparison between the phase volume time of the present iteration and the */ | |
| 583 /* last iteration is less than or equal to one minute. This implies that */ | |
| 584 /* the volume of released gas in the most recent iteration differs from the */ | |
| 585 /* "critical" volume limit by an acceptably small amount. The critical */ | |
| 586 /* volume limit is set by the Critical Volume Parameter Lambda in the program */ | |
| 587 /* settings (default setting is 7500 fsw-min with adjustability range from */ | |
| 588 /* from 6500 to 8300 fsw-min according to Bruce Wienke). */ | |
| 589 /* =============================================================================== */ | |
| 590 /* L50: */ | |
| 591 | |
| 290 | 592 static float deco_stop_depth; |
| 593 static int vpm_calc_critcal_volume(_Bool begin, | |
| 38 | 594 _Bool calc_nulltime) |
| 595 { /* loop will run continuous there is an exit stateme */ | |
| 596 | |
| 597 short i; | |
| 598 | |
| 599 float rounding_operation2; | |
| 600 //float ending_depth; | |
| 601 float deco_ceiling_depth; | |
| 602 | |
| 603 //float deco_time; | |
| 604 int count = 0; | |
| 605 _Bool first_stop; | |
| 606 int dp = 0; | |
| 607 float tissue_He_saturation[16]; | |
| 608 float tissue_N2_saturation[16]; | |
| 609 float vpm_buehlmann_safety_gradient = 1.0f - (((float)pDiveSettings->vpm_conservatism) / 40); | |
| 610 /* =============================================================================== */ | |
| 611 /* CALCULATE CURRENT DECO CEILING BASED ON ALLOWABLE SUPERSATURATION */ | |
| 612 /* GRADIENTS AND SET FIRST DECO STOP. CHECK TO MAKE SURE THAT SELECTED STEP */ | |
| 613 /* SIZE WILL NOT ROUND UP FIRST STOP TO A DEPTH THAT IS BELOW THE DECO ZONE. */ | |
| 614 /* =============================================================================== */ | |
| 615 if(begin) | |
| 616 { | |
| 617 if(depth_start_of_deco_calc < max_first_stop_depth ) | |
| 618 { | |
| 619 if(vpm_b) | |
| 620 { | |
| 621 BOYLES_LAW_COMPENSATION(&max_first_stop_depth, &depth_start_of_deco_calc, &step_size); | |
| 622 } | |
| 623 calc_deco_ceiling(&deco_ceiling_depth, false); | |
| 624 } | |
| 625 else | |
| 626 calc_deco_ceiling(&deco_ceiling_depth, true); | |
| 627 | |
| 628 | |
| 629 if (deco_ceiling_depth <= 0.0f) { | |
| 630 deco_stop_depth = 0.0f; | |
| 631 } else { | |
| 632 rounding_operation2 = deco_ceiling_depth / step_size + ( float)0.5f; | |
| 633 deco_stop_depth = r_nint(&rounding_operation2) * step_size; | |
| 634 } | |
| 635 | |
| 636 // buehlmann safety | |
| 637 if(buehlmannSafety) | |
| 638 { | |
| 639 for (i = 0; i < 16; i++) | |
| 640 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
641 tissue_He_saturation[i] = helium_pressure[i] / 10.0; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
642 tissue_N2_saturation[i] = nitrogen_pressure[i] / 10.0; |
| 38 | 643 } |
| 644 | |
| 645 if(!decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (deco_stop_depth / 10.0f) + pInput->pressure_surface_bar)) | |
| 646 { | |
| 647 | |
| 648 vpm_violates_buehlmann = true; | |
| 649 do { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
650 deco_stop_depth += 3.0; |
| 38 | 651 } while (!decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (deco_stop_depth / 10.0f) + pInput->pressure_surface_bar)); |
| 652 } | |
| 653 } | |
| 654 | |
| 655 /* =============================================================================== */ | |
| 656 /* PERFORM A SEPARATE "PROJECTED ASCENT" OUTSIDE OF THE MAIN PROGRAM TO MAKE */ | |
| 657 /* SURE THAT AN INCREASE IN GAS LOADINGS DURING ASCENT TO THE FIRST STOP WILL */ | |
| 658 /* NOT CAUSE A VIOLATION OF THE DECO CEILING. IF SO, ADJUST THE FIRST STOP */ | |
| 659 /* DEEPER BASED ON STEP SIZE UNTIL A SAFE ASCENT CAN BE MADE. */ | |
| 660 /* Note: this situation is a possibility when ascending from extremely deep */ | |
| 661 /* dives or due to an unusual gas mix selection. */ | |
| 662 /* CHECK AGAIN TO MAKE SURE THAT ADJUSTED FIRST STOP WILL NOT BE BELOW THE */ | |
| 663 /* DECO ZONE. */ | |
| 664 /* =============================================================================== */ | |
| 665 if (deco_stop_depth < depth_start_of_deco_calc) | |
| 666 { | |
| 667 projected_ascent(&depth_start_of_deco_calc, &rate, &deco_stop_depth, &step_size); | |
| 668 } | |
| 669 | |
| 670 /*if (deco_stop_depth > depth_start_of_deco_zone) { | |
| 671 printf("\t\n"); | |
| 672 printf(fmt_905); | |
| 673 printf(fmt_900); | |
| 674 printf("\nPROGRAM TERMINATED\n"); | |
| 675 exit(1); | |
| 676 }*/ | |
| 677 | |
| 678 /* =============================================================================== */ | |
| 679 /* HANDLE THE SPECIAL CASE WHEN NO DECO STOPS ARE REQUIRED - ASCENT CAN BE */ | |
| 680 /* MADE DIRECTLY TO THE SURFACE */ | |
| 681 /* Write ascent data to output file and exit the Critical Volume Loop. */ | |
| 682 /* =============================================================================== */ | |
| 683 | |
| 684 if (deco_stop_depth == 0.0f) | |
| 685 { | |
| 686 if(calc_nulltime) | |
| 687 { | |
| 688 return CALC_END; | |
| 689 } | |
| 690 if(pVpm->deco_zone_reached) | |
| 691 { | |
| 692 for(dp = 0;dp < DECOINFO_STRUCT_MAX_STOPS;dp++) | |
| 693 { | |
| 694 pDecoInfo->output_stop_length_seconds[dp] = 0; | |
| 695 } | |
| 696 pDecoInfo->output_ndl_seconds = 0; | |
| 697 } | |
| 698 | |
| 292 | 699 return CALC_NDL; |
| 38 | 700 /* exit the critical volume l */ |
| 701 } | |
| 702 | |
| 703 /* =============================================================================== */ | |
| 704 /* ASSIGN VARIABLES FOR ASCENT FROM START OF DECO ZONE TO FIRST STOP. SAVE */ | |
| 705 /* FIRST STOP DEPTH FOR LATER USE WHEN COMPUTING THE FINAL ASCENT PROFILE */ | |
| 706 /* =============================================================================== */ | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
707 deco_stop_depth = fmaxf(deco_stop_depth,(float)pDiveSettings->last_stop_depth_bar * 10.0); |
| 38 | 708 starting_depth = depth_start_of_deco_calc; |
| 709 first_stop_depth = deco_stop_depth; | |
| 710 first_stop = true; | |
| 711 } | |
| 712 /* =============================================================================== */ | |
| 713 /* DECO STOP LOOP BLOCK WITHIN CRITICAL VOLUME LOOP */ | |
| 714 /* This loop computes a decompression schedule to the surface during each */ | |
| 715 /* iteration of the critical volume loop. No output is written from this */ | |
| 716 /* loop, rather it computes a schedule from which the in-water portion of the */ | |
| 717 /* total phase volume time (Deco_Phase_Volume_Time) can be extracted. Also, */ | |
| 718 /* the gas loadings computed at the end of this loop are used the subroutine */ | |
| 719 /* which computes the out-of-water portion of the total phase volume time */ | |
| 720 /* (Surface_Phase_Volume_Time) for that schedule. */ | |
| 721 | |
| 722 /* Note that exit is made from the loop after last ascent is made to a deco */ | |
| 723 /* stop depth that is less than or equal to zero. A final deco stop less */ | |
| 724 /* than zero can happen when the user makes an odd step size change during */ | |
| 725 /* ascent - such as specifying a 5 msw step size change at the 3 msw stop! */ | |
| 726 /* =============================================================================== */ | |
| 727 | |
| 728 while(true) /* loop will run continuous there is an break statement */ | |
| 729 { | |
| 730 if(starting_depth > deco_stop_depth ) | |
| 731 gas_loadings_ascent_descen(helium_pressure, nitrogen_pressure, starting_depth, deco_stop_depth, rate,first_stop); | |
| 732 | |
| 733 first_stop = false; | |
| 734 if (deco_stop_depth <= 0.0f) | |
| 735 { | |
| 736 break; | |
| 737 } | |
| 738 if (number_of_changes > 1) | |
| 739 { | |
| 740 int i1 = number_of_changes; | |
| 741 for (i = 2; i <= i1; ++i) { | |
| 742 if (depth_change[i - 1] >= deco_stop_depth) | |
| 743 { | |
| 744 mix_number = mix_change[i - 1]; | |
| 745 rate = rate_change[i - 1]; | |
| 746 step_size = step_size_change[i - 1]; | |
| 747 } | |
| 748 } | |
| 749 } | |
| 750 if(vpm_b) | |
| 751 { | |
| 752 float fist_stop_depth2 = fmaxf(first_stop_depth,max_first_stop_depth); | |
| 753 BOYLES_LAW_COMPENSATION(&fist_stop_depth2, &deco_stop_depth, &step_size); | |
| 754 } | |
| 755 decompression_stop(&deco_stop_depth, &step_size, false); | |
| 756 starting_depth = deco_stop_depth; | |
| 757 | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
758 if(deco_stop_depth == (float)pDiveSettings->last_stop_depth_bar * 10.0) |
| 38 | 759 deco_stop_depth = 0; |
| 760 else | |
| 761 { | |
| 762 deco_stop_depth = deco_stop_depth - step_size; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
763 deco_stop_depth = fmaxf(deco_stop_depth,(float)pDiveSettings->last_stop_depth_bar * 10.0); |
| 38 | 764 } |
| 765 | |
| 766 count++; | |
| 767 //if(count > 14) | |
| 768 //return CALC_CRITICAL2; | |
| 769 /* L60: */ | |
| 770 } | |
| 771 | |
| 772 return vpm_check_converged(calc_nulltime); | |
| 773 } | |
| 774 /* =============================================================================== */ | |
| 775 /* COMPUTE TOTAL PHASE VOLUME TIME AND MAKE CRITICAL VOLUME COMPARISON */ | |
| 776 /* The deco phase volume time is computed from the run time. The surface */ | |
| 777 /* phase volume time is computed in a subroutine based on the surfacing gas */ | |
| 778 /* loadings from previous deco loop block. Next the total phase volume time */ | |
| 779 /* (in-water + surface) for each compartment is compared against the previous */ | |
| 780 /* total phase volume time. The schedule is converged when the difference is */ | |
| 781 /* less than or equal to 1 minute in any one of the 16 compartments. */ | |
| 782 | |
| 783 /* Note: the "phase volume time" is somewhat of a mathematical concept. */ | |
| 784 /* It is the time divided out of a total integration of supersaturation */ | |
| 785 /* gradient x time (in-water and surface). This integration is multiplied */ | |
| 786 /* by the excess bubble number to represent the amount of free-gas released */ | |
| 787 /* as a result of allowing a certain number of excess bubbles to form. */ | |
| 788 /* =============================================================================== */ | |
| 789 /* end of deco stop loop */ | |
| 790 | |
| 290 | 791 static int vpm_check_converged(_Bool calc_nulltime) |
| 38 | 792 { |
| 793 | |
| 794 short i; | |
| 795 float critical_volume_comparison; | |
| 796 float r1; | |
| 797 _Bool schedule_converged = false; | |
| 798 | |
| 799 | |
| 800 deco_phase_volume_time = run_time - run_time_start_of_deco_zone; | |
| 801 calc_surface_phase_volume_time(); | |
| 802 | |
| 803 for (i = 1; i <= 16; ++i) | |
| 804 { | |
| 805 phase_volume_time[i - 1] = | |
| 806 deco_phase_volume_time + surface_phase_volume_time[i - 1]; | |
| 807 critical_volume_comparison = (r1 = phase_volume_time[i - 1] - last_phase_volume_time[i - 1], fabs(r1)); | |
| 808 | |
| 809 if (critical_volume_comparison <= 1.0f) | |
| 810 { | |
| 811 schedule_converged = true; | |
| 812 } | |
| 813 } | |
| 814 | |
| 815 /* =============================================================================== */ | |
| 816 /* CRITICAL VOLUME DECISION TREE BETWEEN LINES 70 AND 99 */ | |
| 817 /* There are two options here. If the Critical Volume Agorithm setting is */ | |
| 818 /* "on" and the schedule is converged, or the Critical Volume Algorithm */ | |
| 819 /* setting was "off" in the first place, the program will re-assign variables */ | |
| 820 /* to their values at the start of ascent (end of bottom time) and process */ | |
| 821 /* a complete decompression schedule once again using all the same ascent */ | |
| 822 /* parameters and first stop depth. This decompression schedule will match */ | |
| 823 /* the last iteration of the Critical Volume Loop and the program will write */ | |
| 824 /* the final deco schedule to the output file. */ | |
| 825 | |
| 826 /* Note: if the Critical Volume Agorithm setting was "off", the final deco */ | |
| 827 /* schedule will be based on "Initial Allowable Supersaturation Gradients." */ | |
| 828 /* If it was "on", the final schedule will be based on "Adjusted Allowable */ | |
| 829 /* Supersaturation Gradients" (gradients that are "relaxed" as a result of */ | |
| 830 /* the Critical Volume Algorithm). */ | |
| 831 | |
| 832 /* If the Critical Volume Agorithm setting is "on" and the schedule is not */ | |
| 833 /* converged, the program will re-assign variables to their values at the */ | |
| 834 /* start of the deco zone and process another trial decompression schedule. */ | |
| 835 /* =============================================================================== */ | |
| 836 /* L70: */ | |
| 837 //Not more than 4 iteration allowed | |
| 838 count_critical_volume_iteration++; | |
| 839 if(count_critical_volume_iteration > 4) | |
| 840 { | |
| 841 //return CALC_FINAL_DECO; | |
| 842 if(calc_nulltime) | |
| 843 return CALC_FINAL_DECO; | |
| 844 else | |
| 845 return vpm_calc_final_deco(true); | |
| 846 } | |
| 847 if (schedule_converged || critical_volume_algorithm_off) | |
| 848 { | |
| 849 | |
| 850 //return CALC_FINAL_DECO; | |
| 851 if(calc_nulltime) | |
| 852 return CALC_FINAL_DECO; | |
| 853 else | |
| 854 return vpm_calc_final_deco(true); | |
| 855 /* final deco schedule */ | |
| 856 /* exit critical volume l */ | |
| 857 | |
| 858 /* =============================================================================== */ | |
| 859 /* IF SCHEDULE NOT CONVERGED, COMPUTE RELAXED ALLOWABLE SUPERSATURATION */ | |
| 860 /* GRADIENTS WITH VPM CRITICAL VOLUME ALGORITHM AND PROCESS ANOTHER */ | |
| 861 /* ITERATION OF THE CRITICAL VOLUME LOOP */ | |
| 862 /* =============================================================================== */ | |
| 863 | |
| 864 } else { | |
| 865 critical_volume(&deco_phase_volume_time); | |
| 866 deco_phase_volume_time = 0.; | |
| 867 run_time = run_time_start_of_deco_calc; | |
| 868 starting_depth = depth_start_of_deco_calc; | |
| 869 mix_number = mix_change[0]; | |
| 870 rate = rate_change[0]; | |
| 871 step_size = step_size_change[0]; | |
| 872 for (i = 1; i <= 16; ++i) | |
| 873 { | |
| 874 last_phase_volume_time[i - 1] = phase_volume_time[i - 1]; | |
| 875 helium_pressure[i - 1] = he_pressure_start_of_deco_zone[i - 1]; | |
| 876 nitrogen_pressure[i - 1] = n2_pressure_start_of_deco_zone[i - 1]; | |
| 877 } | |
| 878 if(calc_nulltime) | |
| 879 return CALC_CRITICAL; | |
| 880 else | |
| 881 return vpm_calc_critcal_volume(true, false); | |
| 882 } | |
| 883 /* end of critical volume decision */ | |
| 884 /* L100: */ | |
| 885 // }/* end of critical vol loop */ | |
| 886 } | |
| 887 | |
| 290 | 888 static void vpm_calc_deco_ceiling(void) |
| 38 | 889 { |
| 890 | |
| 891 short i; | |
| 892 // hw 1601209 float r1; | |
| 893 // hw 1601209 float stop_time; | |
| 894 // hw 1601209 int count = 0; | |
| 895 //static int dp_max; | |
| 896 //static float surfacetime; | |
| 897 // _Bool first_stop = false; | |
| 898 float tissue_He_saturation[16]; | |
| 899 float tissue_N2_saturation[16]; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
900 float vpm_buehlmann_safety_gradient = 1.0f - (((float)pDiveSettings->vpm_conservatism) / 40.0); |
| 38 | 901 //max_first_stop_depth = fmaxf(first_stop_depth,max_first_stop_depth); |
| 902 | |
| 903 /** CALC DECO Ceiling ******************************************************************/ | |
| 904 /** Not when Future stops */ | |
| 905 if(vpm_calc_what == DECOSTOPS) | |
| 906 { | |
| 907 | |
| 908 for (i = 1; i <= 16; ++i) | |
| 909 { | |
| 910 helium_pressure[i - 1] = he_pressure_start_of_deco_zone[i - 1]; | |
| 911 nitrogen_pressure[i - 1] = n2_pressure_start_of_deco_zone[i - 1]; | |
| 912 } | |
| 913 run_time = run_time_start_of_ascent;// run_time_start_of_ascent; | |
| 914 starting_depth = depth_change[0]; | |
| 915 mix_number = mix_change[0]; | |
| 916 rate = rate_change[0]; | |
| 917 //gas_loadings_ascent_descen(helium_pressure,nitrogen_pressure, starting_depth, depth_start_of_deco_calc, rate, true); | |
| 918 | |
| 919 float deco_ceiling_depth = 0.0f; | |
| 920 if(depth_start_of_deco_calc > max_deco_ceiling_depth) | |
| 921 { | |
| 922 calc_deco_ceiling(&deco_ceiling_depth, true); | |
| 923 } | |
| 924 if(buehlmannSafety) | |
| 925 { | |
| 926 for (i = 0; i < 16; i++) | |
| 927 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
928 tissue_He_saturation[i] = helium_pressure[i] / 10.0; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
929 tissue_N2_saturation[i] = nitrogen_pressure[i] / 10.0; |
| 38 | 930 } |
| 931 | |
| 932 if(!decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (deco_ceiling_depth / 10.0f) + pInput->pressure_surface_bar)) | |
| 933 { | |
| 934 vpm_violates_buehlmann = true; | |
| 935 do { | |
| 936 deco_ceiling_depth += 0.1f; | |
| 937 } while (!decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (deco_ceiling_depth / 10.0f) + pInput->pressure_surface_bar)); | |
| 938 } | |
| 939 } | |
| 940 | |
| 941 if (deco_ceiling_depth < depth_start_of_deco_calc) | |
| 942 { | |
| 943 projected_ascent(&depth_start_of_deco_calc, &rate, &deco_ceiling_depth, &step_size); | |
| 944 } | |
| 945 | |
| 946 max_deco_ceiling_depth = fmaxf(max_deco_ceiling_depth,deco_ceiling_depth); | |
| 947 | |
| 948 if(depth_start_of_deco_calc > deco_ceiling_depth) | |
| 949 { | |
| 950 gas_loadings_ascent_descen(helium_pressure,nitrogen_pressure, depth_start_of_deco_calc,deco_ceiling_depth, rate, true); | |
| 951 //surfacetime += segment_time; | |
| 952 } | |
| 953 | |
| 954 if(vpm_b) | |
| 955 { | |
| 956 BOYLES_LAW_COMPENSATION(&max_deco_ceiling_depth, &deco_ceiling_depth, &step_size); | |
| 957 } | |
| 958 calc_deco_ceiling(&deco_ceiling_depth, false); | |
| 959 | |
| 960 // buehlmann safety | |
| 961 if(vpm_violates_buehlmann) | |
| 962 { | |
| 963 for (i = 0; i < 16; i++) | |
| 964 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
965 tissue_He_saturation[i] = helium_pressure[i] / 10.0; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
966 tissue_N2_saturation[i] = nitrogen_pressure[i] / 10.0; |
| 38 | 967 } |
| 968 | |
| 969 if(!decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (deco_ceiling_depth / 10.0f) + pInput->pressure_surface_bar)) | |
| 970 { | |
| 971 vpm_violates_buehlmann = true; | |
| 972 do { | |
| 973 deco_ceiling_depth += 0.1f; | |
| 974 } while (!decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (deco_ceiling_depth / 10.0f) + pInput->pressure_surface_bar)); | |
| 975 } | |
| 976 } | |
| 977 // output_ceiling_meter | |
| 978 if(deco_ceiling_depth > first_stop_depth) | |
| 979 deco_ceiling_depth = first_stop_depth; | |
| 980 pDecoInfo->output_ceiling_meter = deco_ceiling_depth ; | |
| 981 } | |
| 982 else | |
| 983 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
984 pDecoInfo->output_ceiling_meter = 0.0; |
| 38 | 985 } |
| 986 | |
| 987 // fix hw 160627 | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
988 if(pDecoInfo->output_ceiling_meter < 0.0) |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
989 pDecoInfo->output_ceiling_meter = 0.0; |
| 38 | 990 |
| 991 /*** End CALC ceiling ***************************************************/ | |
| 992 } | |
| 993 | |
| 994 | |
| 995 /* =============================================================================== */ | |
| 996 /* DECO STOP LOOP BLOCK FOR FINAL DECOMPRESSION SCHEDULE */ | |
| 997 /* =============================================================================== */ | |
| 998 | |
| 290 | 999 static int vpm_calc_final_deco(_Bool begin) |
| 38 | 1000 { |
| 1001 short i; | |
| 1002 float r1; | |
| 1003 float stop_time; | |
| 1004 int count = 0; | |
| 1005 static int dp_max; | |
| 1006 static float surfacetime; | |
| 1007 _Bool first_stop = false; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1008 |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1009 short stop_time_seconds; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1010 |
| 38 | 1011 max_first_stop_depth = fmaxf(first_stop_depth,max_first_stop_depth); |
| 1012 if(begin) | |
| 1013 { | |
| 1014 gCNS_VPM = 0; | |
| 1015 dp_max = 0; | |
| 1016 for (i = 1; i <= 16; ++i) | |
| 1017 { | |
| 1018 helium_pressure[i - 1] = | |
| 1019 he_pressure_start_of_ascent[i - 1]; | |
| 1020 nitrogen_pressure[i - 1] = | |
| 1021 n2_pressure_start_of_ascent[i - 1]; | |
| 1022 } | |
| 1023 run_time = run_time_start_of_ascent;// run_time_start_of_ascent; | |
| 1024 starting_depth = depth_change[0]; | |
| 1025 mix_number = mix_change[0]; | |
| 1026 rate = rate_change[0]; | |
| 1027 step_size = step_size_change[0]; | |
| 1028 deco_stop_depth = first_stop_depth; | |
| 1029 max_first_stop_depth = fmaxf(first_stop_depth,max_first_stop_depth); | |
| 1030 last_run_time = 0.; | |
| 1031 | |
| 1032 | |
| 1033 | |
| 1034 /* =============================================================================== */ | |
| 1035 /* DECO STOP LOOP BLOCK FOR FINAL DECOMPRESSION SCHEDULE */ | |
| 1036 /* =============================================================================== */ | |
| 1037 surfacetime = 0; | |
| 1038 first_stop = true; | |
| 1039 } | |
| 1040 | |
| 1041 while(true) /* loop will run continuous until there is an break statement */ | |
| 1042 { | |
| 1043 if(starting_depth > deco_stop_depth) | |
| 1044 { | |
| 1045 gas_loadings_ascent_descen(helium_pressure,nitrogen_pressure, starting_depth,deco_stop_depth, rate, first_stop); | |
| 1046 surfacetime += segment_time; | |
| 1047 } | |
| 1048 | |
| 1049 /* =============================================================================== */ | |
| 1050 /* DURING FINAL DECOMPRESSION SCHEDULE PROCESS, COMPUTE MAXIMUM ACTUAL */ | |
| 1051 /* SUPERSATURATION GRADIENT RESULTING IN EACH COMPARTMENT */ | |
| 1052 /* If there is a repetitive dive, this will be used later in the VPM */ | |
| 1053 /* Repetitive Algorithm to adjust the values for critical radii. */ | |
| 1054 /* =============================================================================== */ | |
| 1055 if(vpm_calc_what == DECOSTOPS) | |
| 1056 calc_max_actual_gradient(&deco_stop_depth); | |
| 1057 | |
| 1058 if (deco_stop_depth <= 0.0f) { | |
| 1059 break; | |
| 1060 } | |
| 1061 if (number_of_changes > 1) | |
| 1062 { | |
| 1063 int i1 = number_of_changes; | |
| 1064 for (i = 2; i <= i1; ++i) | |
| 1065 { | |
| 1066 if (depth_change[i - 1] >= deco_stop_depth) | |
| 1067 { | |
| 1068 mix_number = mix_change[i - 1]; | |
| 1069 rate = rate_change[i - 1]; | |
| 1070 step_size = step_size_change[i - 1]; | |
| 1071 } | |
| 1072 } | |
| 1073 } | |
| 1074 | |
| 1075 if(first_stop) | |
| 1076 { | |
| 1077 run_time_first_stop = run_time; | |
| 1078 first_stop = false; | |
| 1079 } | |
| 1080 if(vpm_b) | |
| 1081 { | |
| 1082 BOYLES_LAW_COMPENSATION(&max_first_stop_depth, &deco_stop_depth, &step_size); | |
| 1083 } | |
| 1084 decompression_stop(&deco_stop_depth, &step_size, true); | |
| 1085 | |
| 1086 /* =============================================================================== */ | |
| 1087 /* This next bit justs rounds up the stop time at the first stop to be in */ | |
| 1088 /* whole increments of the minimum stop time (to make for a nice deco table). */ | |
| 1089 /* =============================================================================== */ | |
| 1090 | |
| 1091 if (last_run_time == 0.0f) | |
| 1092 { | |
| 1093 r1 = segment_time / minimum_deco_stop_time + 0.5f; | |
| 1094 stop_time = r_int(&r1) * minimum_deco_stop_time; | |
| 1095 } else { | |
| 1096 stop_time = run_time - last_run_time; | |
| 1097 } | |
| 1098 stop_time = segment_time; | |
| 1099 surfacetime += stop_time; | |
| 1100 if((vpm_calc_what == DECOSTOPS) || (vpm_calc_what == BAILOUTSTOPS)) | |
| 1101 { | |
| 1102 int dp = 0; | |
| 1103 if(deco_stop_depth == (float)pDiveSettings->last_stop_depth_bar * 10) | |
| 1104 { | |
| 1105 dp = 0; | |
| 1106 } | |
| 1107 else | |
| 1108 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1109 dp = 1 + (unsigned short)((deco_stop_depth - (pDiveSettings->input_second_to_last_stop_depth_bar * 10.0)) / step_size); |
| 38 | 1110 } |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1111 |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1112 //dp_max = (int)fmaxf(dp_max,dp); |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1113 if(dp > dp_max) |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1114 { |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1115 dp_max = dp; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1116 } |
| 38 | 1117 if(dp < DECOINFO_STRUCT_MAX_STOPS) |
| 1118 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1119 stop_time_seconds = (short) fminf((999.9 * 60.0), (stop_time *60.0)); |
| 38 | 1120 // |
| 1121 | |
| 1122 //if(vpm_calc_what == DECOSTOPS) | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1123 pDecoInfo->output_stop_length_seconds[dp] = stop_time_seconds; |
| 38 | 1124 //else |
| 1125 //decostop_bailout[dp] = (unsigned short)stop_time_seconds; | |
| 1126 } | |
| 1127 } | |
| 1128 | |
| 1129 | |
| 1130 /* =============================================================================== */ | |
| 1131 /* DURING FINAL DECOMPRESSION SCHEDULE, IF MINIMUM STOP TIME PARAMETER IS A */ | |
| 1132 /* WHOLE NUMBER (i.e. 1 minute) THEN WRITE DECO SCHEDULE USING short */ | |
| 1133 /* NUMBERS (looks nicer). OTHERWISE, USE DECIMAL NUMBERS. */ | |
| 1134 /* Note: per the request of a noted exploration diver(!), program now allows */ | |
| 1135 /* a minimum stop time of less than one minute so that total ascent time can */ | |
| 1136 /* be minimized on very long dives. In fact, with step size set at 1 fsw or */ | |
| 1137 /* 0.2 msw and minimum stop time set at 0.1 minute (6 seconds), a near */ | |
| 1138 /* continuous decompression schedule can be computed. */ | |
| 1139 /* =============================================================================== */ | |
| 1140 | |
| 1141 starting_depth = deco_stop_depth; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1142 if(deco_stop_depth == (float)pDiveSettings->last_stop_depth_bar * 10.0) |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1143 deco_stop_depth = 0.0; |
| 38 | 1144 else |
| 1145 { | |
| 1146 deco_stop_depth = deco_stop_depth - step_size; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1147 deco_stop_depth = fmaxf(deco_stop_depth,(float)pDiveSettings->last_stop_depth_bar * 10.0); |
| 38 | 1148 } |
| 1149 | |
| 1150 last_run_time = run_time; | |
| 1151 count++; | |
| 1152 //if(count > 14) | |
| 1153 //return CALC_FINAL_DECO2; | |
| 1154 /* L80: */ | |
| 1155 } /* for final deco sche */ | |
| 1156 | |
| 1157 if( (vpm_calc_what == DECOSTOPS) || (vpm_calc_what == BAILOUTSTOPS)) | |
| 1158 { | |
| 1159 for(int dp = dp_max +1;dp < DECOINFO_STRUCT_MAX_STOPS;dp++) | |
| 1160 { | |
| 1161 //if(vpm_calc_what == DECOSTOPS) | |
| 1162 pDecoInfo->output_stop_length_seconds[dp] = 0; | |
| 1163 //else | |
| 1164 //decostop_bailout[dp] = 0; | |
| 1165 } | |
| 1166 } | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1167 pDecoInfo->output_time_to_surface_seconds = (int)(surfacetime * 60.0); |
| 38 | 1168 pDecoInfo->output_ndl_seconds = 0; |
| 1169 | |
| 1170 vpm_calc_deco_ceiling(); | |
| 1171 /* end of deco stop lo */ | |
| 1172 return CALC_END; | |
| 1173 } | |
| 1174 | |
| 1175 /* =============================================================================== */ | |
| 1176 /* SUBROUTINE NUCLEAR_REGENERATION */ | |
| 1177 /* Purpose: This subprogram calculates the regeneration of VPM critical */ | |
| 1178 /* radii that takes place over the dive time. The regeneration time constant */ | |
| 1179 /* has a time scale of weeks so this will have very little impact on dives of */ | |
| 1180 /* normal length, but will have a major impact for saturation dives. */ | |
| 1181 /* =============================================================================== */ | |
| 1182 | |
| 290 | 1183 static int nuclear_regeneration(float *dive_time) |
| 38 | 1184 { |
| 1185 /* Local variables */ | |
| 1186 float crush_pressure_adjust_ratio_he, | |
| 1187 ending_radius_n2, | |
| 1188 ending_radius_he; | |
| 1189 short i; | |
| 1190 float crushing_pressure_pascals_n2, | |
| 1191 crushing_pressure_pascals_he, | |
| 1192 adj_crush_pressure_n2_pascals, | |
| 1193 adj_crush_pressure_he_pascals, | |
| 1194 crush_pressure_adjust_ratio_n2; | |
| 1195 | |
| 1196 /* loop */ | |
| 1197 /* =============================================================================== */ | |
| 1198 /* CALCULATIONS */ | |
| 1199 /* First convert the maximum crushing pressure obtained for each compartment */ | |
| 1200 /* to Pascals. Next, compute the ending radius for helium and nitrogen */ | |
| 1201 /* critical nuclei in each compartment. */ | |
| 1202 /* =============================================================================== */ | |
| 1203 | |
| 1204 for (i = 1; i <= 16; ++i) | |
| 1205 { | |
| 1206 crushing_pressure_pascals_he = | |
| 1207 pVpm->max_crushing_pressure_he[i - 1] / UNITS_FACTOR * 101325.0f; | |
| 1208 crushing_pressure_pascals_n2 = | |
| 1209 pVpm->max_crushing_pressure_n2[i - 1] / UNITS_FACTOR * 101325.0f; | |
| 1210 ending_radius_he = | |
| 1211 1.0f / (crushing_pressure_pascals_he / | |
| 1212 ((SKIN_COMPRESSION_GAMMAC - SURFACE_TENSION_GAMMA) * 2.0f) + | |
| 1213 1.0f / pVpm->adjusted_critical_radius_he[i - 1]); | |
| 1214 ending_radius_n2 = | |
| 1215 1.0f / (crushing_pressure_pascals_n2 / | |
| 1216 ((SKIN_COMPRESSION_GAMMAC - SURFACE_TENSION_GAMMA) * 2.0f) + | |
| 1217 1.0f / pVpm->adjusted_critical_radius_n2[i - 1]); | |
| 1218 | |
| 1219 /* =============================================================================== */ | |
| 1220 /* A "regenerated" radius for each nucleus is now calculated based on the */ | |
| 1221 /* regeneration time constant. This means that after application of */ | |
| 1222 /* crushing pressure and reduction in radius, a nucleus will slowly grow */ | |
| 1223 /* back to its original initial radius over a period of time. This */ | |
| 1224 /* phenomenon is probabilistic in nature and depends on absolute temperature. */ | |
| 1225 /* It is independent of crushing pressure. */ | |
| 1226 /* =============================================================================== */ | |
| 1227 | |
| 1228 regenerated_radius_he[i - 1] = | |
| 1229 pVpm->adjusted_critical_radius_he[i - 1] + | |
| 1230 (ending_radius_he - pVpm->adjusted_critical_radius_he[i - 1]) * | |
| 1231 expf(-(*dive_time) / REGENERATION_TIME_CONSTANT); | |
| 1232 regenerated_radius_n2[i - 1] = | |
| 1233 pVpm->adjusted_critical_radius_n2[i - 1] + | |
| 1234 (ending_radius_n2 - pVpm->adjusted_critical_radius_n2[i - 1]) * | |
| 1235 expf(-(*dive_time) / REGENERATION_TIME_CONSTANT); | |
| 1236 | |
| 1237 /* =============================================================================== */ | |
| 1238 /* In order to preserve reference back to the initial critical radii after */ | |
| 1239 /* regeneration, an "adjusted crushing pressure" for the nuclei in each */ | |
| 1240 /* compartment must be computed. In other words, this is the value of */ | |
| 1241 /* crushing pressure that would have reduced the original nucleus to the */ | |
| 1242 /* to the present radius had regeneration not taken place. The ratio */ | |
| 1243 /* for adjusting crushing pressure is obtained from algebraic manipulation */ | |
| 1244 /* of the standard VPM equations. The adjusted crushing pressure, in lieu */ | |
| 1245 /* of the original crushing pressure, is then applied in the VPM Critical */ | |
| 1246 /* Volume Algorithm and the VPM Repetitive Algorithm. */ | |
| 1247 /* =============================================================================== */ | |
| 1248 | |
| 1249 crush_pressure_adjust_ratio_he = | |
| 1250 ending_radius_he * (pVpm->adjusted_critical_radius_he[i - 1] - | |
| 1251 regenerated_radius_he[i - 1]) / | |
| 1252 (regenerated_radius_he[i - 1] * | |
| 1253 (pVpm->adjusted_critical_radius_he[i - 1] - | |
| 1254 ending_radius_he)); | |
| 1255 crush_pressure_adjust_ratio_n2 = | |
| 1256 ending_radius_n2 * (pVpm->adjusted_critical_radius_n2[i - 1] - | |
| 1257 regenerated_radius_n2[i - 1]) / | |
| 1258 (regenerated_radius_n2[i - 1] * | |
| 1259 (pVpm->adjusted_critical_radius_n2[i - 1] - | |
| 1260 ending_radius_n2)); | |
| 1261 adj_crush_pressure_he_pascals = | |
| 1262 crushing_pressure_pascals_he * crush_pressure_adjust_ratio_he; | |
| 1263 adj_crush_pressure_n2_pascals = | |
| 1264 crushing_pressure_pascals_n2 * crush_pressure_adjust_ratio_n2; | |
| 1265 pVpm->adjusted_crushing_pressure_he[i - 1] = | |
| 1266 adj_crush_pressure_he_pascals / 101325.0f * UNITS_FACTOR; | |
| 1267 pVpm->adjusted_crushing_pressure_n2[i - 1] = | |
| 1268 adj_crush_pressure_n2_pascals / 101325.0f * UNITS_FACTOR; | |
| 1269 } | |
| 1270 return 0; | |
| 1271 } /* nuclear_regeneration */ | |
| 1272 | |
| 1273 /* =============================================================================== */ | |
| 1274 /* SUBROUTINE CALC_INITIAL_ALLOWABLE_GRADIENT */ | |
| 1275 /* Purpose: This subprogram calculates the initial allowable gradients for */ | |
| 1276 /* helium and nitrogren in each compartment. These are the gradients that */ | |
| 1277 /* will be used to set the deco ceiling on the first pass through the deco */ | |
| 1278 /* loop. If the Critical Volume Algorithm is set to "off", then these */ | |
| 1279 /* gradients will determine the final deco schedule. Otherwise, if the */ | |
| 1280 /* Critical Volume Algorithm is set to "on", these gradients will be further */ | |
| 1281 /* "relaxed" by the Critical Volume Algorithm subroutine. The initial */ | |
| 1282 /* allowable gradients are referred to as "PssMin" in the papers by Yount */ | |
| 1283 /* and colleauges, i.e., the minimum supersaturation pressure gradients */ | |
| 1284 /* that will probe bubble formation in the VPM nuclei that started with the */ | |
| 1285 /* designated minimum initial radius (critical radius). */ | |
| 1286 | |
| 1287 /* The initial allowable gradients are computed directly from the */ | |
| 1288 /* "regenerated" radii after the Nuclear Regeneration subroutine. These */ | |
| 1289 /* gradients are tracked separately for helium and nitrogen. */ | |
| 1290 /* =============================================================================== */ | |
| 1291 | |
| 290 | 1292 static int calc_initial_allowable_gradient() |
| 38 | 1293 { |
| 1294 float initial_allowable_grad_n2_pa, | |
| 1295 initial_allowable_grad_he_pa; | |
| 1296 short i; | |
| 1297 | |
| 1298 /* loop */ | |
| 1299 /* =============================================================================== */ | |
| 1300 /* CALCULATIONS */ | |
| 1301 /* The initial allowable gradients are computed in Pascals and then converted */ | |
| 1302 /* to the diving pressure units. Two different sets of arrays are used to */ | |
| 1303 /* save the calculations - Initial Allowable Gradients and Allowable */ | |
| 1304 /* Gradients. The Allowable Gradients are assigned the values from Initial */ | |
| 1305 /* Allowable Gradients however the Allowable Gradients can be changed later */ | |
| 1306 /* by the Critical Volume subroutine. The values for the Initial Allowable */ | |
| 1307 /* Gradients are saved in a global array for later use by both the Critical */ | |
| 1308 /* Volume subroutine and the VPM Repetitive Algorithm subroutine. */ | |
| 1309 /* =============================================================================== */ | |
| 1310 | |
| 1311 for (i = 1; i <= 16; ++i) | |
| 1312 { | |
| 1313 initial_allowable_grad_n2_pa = | |
| 1314 SURFACE_TENSION_GAMMA * 2.0f * | |
| 1315 (SKIN_COMPRESSION_GAMMAC - SURFACE_TENSION_GAMMA) / | |
| 1316 (regenerated_radius_n2[i - 1] * SKIN_COMPRESSION_GAMMAC); | |
| 1317 | |
| 1318 initial_allowable_grad_he_pa = | |
| 1319 SURFACE_TENSION_GAMMA * 2.0f * | |
| 1320 (SKIN_COMPRESSION_GAMMAC - SURFACE_TENSION_GAMMA) / | |
| 1321 (regenerated_radius_he[i - 1] * SKIN_COMPRESSION_GAMMAC); | |
| 1322 | |
| 1323 pVpm->initial_allowable_gradient_n2[i - 1] = | |
| 1324 initial_allowable_grad_n2_pa / 101325.0f * UNITS_FACTOR; | |
| 1325 | |
| 1326 pVpm->initial_allowable_gradient_he[i - 1] = | |
| 1327 initial_allowable_grad_he_pa / 101325.0f * UNITS_FACTOR; | |
| 1328 | |
| 1329 allowable_gradient_he[i - 1] = | |
| 1330 pVpm->initial_allowable_gradient_he[i - 1]; | |
| 1331 | |
| 1332 allowable_gradient_n2[i - 1] = | |
| 1333 pVpm->initial_allowable_gradient_n2[i - 1]; | |
| 1334 } | |
| 1335 return 0; | |
| 1336 } /* calc_initial_allowable_gradient */ | |
| 1337 | |
| 1338 /* =============================================================================== */ | |
| 1339 /* SUBROUTINE CALC_DECO_CEILING */ | |
| 1340 /* Purpose: This subprogram calculates the deco ceiling (the safe ascent */ | |
| 1341 /* depth) in each compartment, based on the allowable gradients, and then */ | |
| 1342 /* finds the deepest deco ceiling across all compartments. This deepest */ | |
| 1343 /* value (Deco Ceiling Depth) is then used by the Decompression Stop */ | |
| 1344 /* subroutine to determine the actual deco schedule. */ | |
| 1345 /* =============================================================================== */ | |
| 1346 | |
| 290 | 1347 static int calc_deco_ceiling(float *deco_ceiling_depth,_Bool fallowable) |
| 38 | 1348 { |
| 1349 /* System generated locals */ | |
| 1350 float r1, r2; | |
| 1351 /* Local variables */ | |
| 1352 float weighted_allowable_gradient; | |
| 1353 short i; | |
| 1354 float compartment_deco_ceiling[16], | |
| 1355 gas_loading, | |
| 1356 tolerated_ambient_pressure; | |
| 1357 float gradient_he, gradient_n2; | |
| 1358 | |
| 1359 if(!vpm_b) | |
| 1360 fallowable = true; | |
| 1361 /* loop */ | |
| 1362 /* =============================================================================== */ | |
| 1363 /* CALCULATIONS */ | |
| 1364 /* Since there are two sets of allowable gradients being tracked, one for */ | |
| 1365 /* helium and one for nitrogen, a "weighted allowable gradient" must be */ | |
| 1366 /* computed each time based on the proportions of helium and nitrogen in */ | |
| 1367 /* each compartment. This proportioning follows the methodology of */ | |
| 1368 /* Buhlmann/Keller. If there is no helium and nitrogen in the compartment, */ | |
| 1369 /* such as after extended periods of oxygen breathing, then the minimum value */ | |
| 1370 /* across both gases will be used. It is important to note that if a */ | |
| 1371 /* compartment is empty of helium and nitrogen, then the weighted allowable */ | |
| 1372 /* gradient formula cannot be used since it will result in division by zero. */ | |
| 1373 /* =============================================================================== */ | |
| 1374 | |
| 1375 for (i = 1; i <= 16; ++i) | |
| 1376 { | |
| 1377 | |
| 1378 // abfrage raus und pointer stattdessen | |
| 1379 if(fallowable){ | |
| 1380 gradient_he = allowable_gradient_he[i-1]; | |
| 1381 gradient_n2 = allowable_gradient_n2[i-1]; | |
| 1382 } | |
| 1383 else{ | |
| 1384 gradient_he = deco_gradient_he[i-1]; | |
| 1385 gradient_n2 = deco_gradient_n2[i-1]; | |
| 1386 } | |
| 1387 | |
| 1388 gas_loading = helium_pressure[i - 1] + nitrogen_pressure[i - 1]; | |
| 1389 | |
| 1390 if (gas_loading > 0) | |
| 1391 { | |
| 1392 weighted_allowable_gradient = | |
| 1393 (gradient_he * helium_pressure[i - 1] + | |
| 1394 gradient_n2 * nitrogen_pressure[i - 1]) / | |
| 1395 (helium_pressure[i - 1] + nitrogen_pressure[i - 1]); | |
| 1396 | |
| 1397 tolerated_ambient_pressure = | |
| 1398 gas_loading + | |
| 1399 CONSTANT_PRESSURE_OTHER_GASES - | |
| 1400 weighted_allowable_gradient; | |
| 1401 } | |
| 1402 else | |
| 1403 { | |
| 1404 /* Computing MIN */ | |
| 1405 r1 = gradient_he; | |
| 1406 r2 = gradient_n2; | |
| 1407 weighted_allowable_gradient = fminf(r1,r2); | |
| 1408 | |
| 1409 tolerated_ambient_pressure = | |
| 1410 CONSTANT_PRESSURE_OTHER_GASES - weighted_allowable_gradient; | |
| 1411 } | |
| 1412 | |
| 1413 /* =============================================================================== */ | |
| 1414 /* The tolerated ambient pressure cannot be less than zero absolute, i.e., */ | |
| 1415 /* the vacuum of outer space! */ | |
| 1416 /* =============================================================================== */ | |
| 1417 | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1418 if (tolerated_ambient_pressure < 0.0) { |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1419 tolerated_ambient_pressure = 0.0; |
| 38 | 1420 } |
| 1421 compartment_deco_ceiling[i - 1] = | |
| 1422 tolerated_ambient_pressure - barometric_pressure; | |
| 1423 } | |
| 1424 | |
| 1425 /* =============================================================================== */ | |
| 1426 /* The Deco Ceiling Depth is computed in a loop after all of the individual */ | |
| 1427 /* compartment deco ceilings have been calculated. It is important that the */ | |
| 1428 /* Deco Ceiling Depth (max deco ceiling across all compartments) only be */ | |
| 1429 /* extracted from the compartment values and not be compared against some */ | |
| 1430 /* initialization value. For example, if MAX(Deco_Ceiling_Depth . .) was */ | |
| 1431 /* compared against zero, this could cause a program lockup because sometimes */ | |
| 1432 /* the Deco Ceiling Depth needs to be negative (but not less than zero */ | |
| 1433 /* absolute ambient pressure) in order to decompress to the last stop at zero */ | |
| 1434 /* depth. */ | |
| 1435 /* =============================================================================== */ | |
| 1436 | |
| 1437 *deco_ceiling_depth = compartment_deco_ceiling[0]; | |
| 1438 for (i = 2; i <= 16; ++i) | |
| 1439 { | |
| 1440 /* Computing MAX */ | |
| 1441 r1 = *deco_ceiling_depth; | |
| 1442 r2 = compartment_deco_ceiling[i - 1]; | |
| 1443 *deco_ceiling_depth = fmaxf(r1,r2); | |
| 1444 } | |
| 1445 return 0; | |
| 1446 } /* calc_deco_ceiling */ | |
| 1447 | |
| 1448 | |
| 1449 | |
| 1450 /* =============================================================================== */ | |
| 1451 /* SUBROUTINE CALC_MAX_ACTUAL_GRADIENT */ | |
| 1452 /* Purpose: This subprogram calculates the actual supersaturation gradient */ | |
| 1453 /* obtained in each compartment as a result of the ascent profile during */ | |
| 1454 /* decompression. Similar to the concept with crushing pressure, the */ | |
| 1455 /* supersaturation gradients are not cumulative over a multi-level, staged */ | |
| 1456 /* ascent. Rather, it will be the maximum value obtained in any one discrete */ | |
| 1457 /* step of the overall ascent. Thus, the program must compute and store the */ | |
| 1458 /* maximum actual gradient for each compartment that was obtained across all */ | |
| 1459 /* steps of the ascent profile. This subroutine is invoked on the last pass */ | |
| 1460 /* through the deco stop loop block when the final deco schedule is being */ | |
| 1461 /* generated. */ | |
| 1462 /* */ | |
| 1463 /* The max actual gradients are later used by the VPM Repetitive Algorithm to */ | |
| 1464 /* determine if adjustments to the critical radii are required. If the max */ | |
| 1465 /* actual gradient did not exceed the initial alllowable gradient, then no */ | |
| 1466 /* adjustment will be made. However, if the max actual gradient did exceed */ | |
| 1467 /* the intitial allowable gradient, such as permitted by the Critical Volume */ | |
| 1468 /* Algorithm, then the critical radius will be adjusted (made larger) on the */ | |
| 1469 /* repetitive dive to compensate for the bubbling that was allowed on the */ | |
| 1470 /* previous dive. The use of the max actual gradients is intended to prevent */ | |
| 1471 /* the repetitive algorithm from being overly conservative. */ | |
| 1472 /* =============================================================================== */ | |
| 1473 | |
| 290 | 1474 static int calc_max_actual_gradient(float *deco_stop_depth) |
| 38 | 1475 { |
| 1476 /* System generated locals */ | |
| 1477 float r1; | |
| 1478 | |
| 1479 /* Local variables */ | |
| 1480 short i; | |
| 1481 float compartment_gradient; | |
| 1482 | |
| 1483 /* loop */ | |
| 1484 /* =============================================================================== */ | |
| 1485 /* CALCULATIONS */ | |
| 1486 /* Note: negative supersaturation gradients are meaningless for this */ | |
| 1487 /* application, so the values must be equal to or greater than zero. */ | |
| 1488 /* =============================================================================== */ | |
| 1489 | |
| 1490 for (i = 1; i <= 16; ++i) | |
| 1491 { | |
| 1492 compartment_gradient = | |
| 1493 helium_pressure[i - 1] + | |
| 1494 nitrogen_pressure[i - 1] + | |
| 1495 CONSTANT_PRESSURE_OTHER_GASES - | |
| 1496 (*deco_stop_depth + barometric_pressure); | |
| 1497 if (compartment_gradient <= 0.0f) { | |
| 1498 compartment_gradient = 0.0f; | |
| 1499 } | |
| 1500 /* Computing MAX */ | |
| 1501 r1 = pVpm->max_actual_gradient[i - 1]; | |
| 1502 pVpm->max_actual_gradient[i - 1] = fmaxf(r1, compartment_gradient); | |
| 1503 } | |
| 1504 return 0; | |
| 1505 } /* calc_max_actual_gradient */ | |
| 1506 | |
| 1507 /* =============================================================================== */ | |
| 1508 /* SUBROUTINE CALC_SURFACE_PHASE_VOLUME_TIME */ | |
| 1509 /* Purpose: This subprogram computes the surface portion of the total phase */ | |
| 1510 /* volume time. This is the time factored out of the integration of */ | |
| 1511 /* supersaturation gradient x time over the surface interval. The VPM */ | |
| 1512 /* considers the gradients that allow bubbles to form or to drive bubble */ | |
| 1513 /* growth both in the water and on the surface after the dive. */ | |
| 1514 | |
| 1515 /* This subroutine is a new development to the VPM algorithm in that it */ | |
| 1516 /* computes the time course of supersaturation gradients on the surface */ | |
| 1517 /* when both helium and nitrogen are present. Refer to separate write-up */ | |
| 1518 /* for a more detailed explanation of this algorithm. */ | |
| 1519 /* =============================================================================== */ | |
| 1520 | |
| 290 | 1521 static int calc_surface_phase_volume_time() |
| 38 | 1522 { |
| 1523 /* Local variables */ | |
| 1524 float decay_time_to_zero_gradient; | |
| 1525 short i; | |
| 1526 float integral_gradient_x_time, | |
| 1527 surface_inspired_n2_pressure; | |
| 1528 | |
| 1529 /* loop */ | |
| 1530 /* =============================================================================== */ | |
| 1531 /* CALCULATIONS */ | |
| 1532 /* =============================================================================== */ | |
| 1533 | |
| 1534 surface_inspired_n2_pressure = | |
| 1535 (barometric_pressure - WATER_VAPOR_PRESSURE) * 0.79f; | |
| 1536 for (i = 1; i <= 16; ++i) | |
| 1537 { | |
| 1538 if (nitrogen_pressure[i - 1] > surface_inspired_n2_pressure) | |
| 1539 { | |
| 1540 surface_phase_volume_time[i - 1] = | |
| 1541 (helium_pressure[i - 1] / HELIUM_TIME_CONSTANT[i - 1] + | |
| 1542 (nitrogen_pressure[i - 1] - surface_inspired_n2_pressure) / | |
| 1543 NITROGEN_TIME_CONSTANT[i - 1]) / | |
| 1544 (helium_pressure[i - 1] + nitrogen_pressure[i - 1] - | |
| 1545 surface_inspired_n2_pressure); | |
| 1546 } else if (nitrogen_pressure[i - 1] <= surface_inspired_n2_pressure && | |
| 1547 helium_pressure[i - 1] + nitrogen_pressure[i - 1] >= surface_inspired_n2_pressure) | |
| 1548 { | |
| 1549 decay_time_to_zero_gradient = | |
| 1550 1.0f / (NITROGEN_TIME_CONSTANT[i - 1] - HELIUM_TIME_CONSTANT[i - 1]) * | |
| 1551 log((surface_inspired_n2_pressure - nitrogen_pressure[i - 1]) / | |
| 1552 helium_pressure[i - 1]); | |
| 1553 integral_gradient_x_time = | |
| 1554 helium_pressure[i - 1] / | |
| 1555 HELIUM_TIME_CONSTANT[i - 1] * | |
| 1556 (1.0f - expf(-HELIUM_TIME_CONSTANT[i - 1] * | |
| 1557 decay_time_to_zero_gradient)) + | |
| 1558 (nitrogen_pressure[i - 1] - surface_inspired_n2_pressure) / | |
| 1559 NITROGEN_TIME_CONSTANT[i - 1] * | |
| 1560 (1.0f - expf(-NITROGEN_TIME_CONSTANT[i - 1] * | |
| 1561 decay_time_to_zero_gradient)); | |
| 1562 surface_phase_volume_time[i - 1] = | |
| 1563 integral_gradient_x_time / | |
| 1564 (helium_pressure[i - 1] + | |
| 1565 nitrogen_pressure[i - 1] - | |
| 1566 surface_inspired_n2_pressure); | |
| 1567 } else { | |
| 1568 surface_phase_volume_time[i - 1] = 0.0f; | |
| 1569 } | |
| 1570 } | |
| 1571 return 0; | |
| 1572 } /* calc_surface_phase_volume_time */ | |
| 1573 | |
| 1574 /* =============================================================================== */ | |
| 1575 /* SUBROUTINE CRITICAL_VOLUME */ | |
| 1576 /* Purpose: This subprogram applies the VPM Critical Volume Algorithm. This */ | |
| 1577 /* algorithm will compute "relaxed" gradients for helium and nitrogen based */ | |
| 1578 /* on the setting of the Critical Volume Parameter Lambda. */ | |
| 1579 /* =============================================================================== */ | |
| 1580 | |
| 290 | 1581 static int critical_volume(float *deco_phase_volume_time) |
| 38 | 1582 { |
| 1583 /* System generated locals */ | |
| 1584 float r1; | |
| 1585 | |
| 1586 /* Local variables */ | |
| 1587 float initial_allowable_grad_n2_pa, | |
| 1588 initial_allowable_grad_he_pa, | |
| 1589 parameter_lambda_pascals, b, | |
| 1590 c; | |
| 1591 short i; | |
| 1592 float new_allowable_grad_n2_pascals, | |
| 1593 phase_volume_time[16], | |
| 1594 new_allowable_grad_he_pascals, | |
| 1595 adj_crush_pressure_n2_pascals, | |
| 1596 adj_crush_pressure_he_pascals; | |
| 1597 | |
| 1598 /* loop */ | |
| 1599 /* =============================================================================== */ | |
| 1600 /* CALCULATIONS */ | |
| 1601 /* Note: Since the Critical Volume Parameter Lambda was defined in units of */ | |
| 1602 /* fsw-min in the original papers by Yount and colleauges, the same */ | |
| 1603 /* convention is retained here. Although Lambda is adjustable only in units */ | |
| 1604 /* of fsw-min in the program settings (range from 6500 to 8300 with default */ | |
| 1605 /* 7500), it will convert to the proper value in Pascals-min in this */ | |
| 1606 /* subroutine regardless of which diving pressure units are being used in */ | |
| 1607 /* the main program - feet of seawater (fsw) or meters of seawater (msw). */ | |
| 1608 /* The allowable gradient is computed using the quadratic formula (refer to */ | |
| 1609 /* separate write-up posted on the Deco List web site). */ | |
| 1610 /* =============================================================================== */ | |
| 1611 | |
| 1612 /** | |
| 1613 ****************************************************************************** | |
| 1614 * @brief critical_volume comment by hw | |
| 1615 * @version V0.0.1 | |
| 1616 * @date 19-April-2014 | |
| 1617 * @retval global: allowable_gradient_he[i], allowable_gradient_n2[i] | |
| 1618 ****************************************************************************** | |
| 1619 */ | |
| 1620 | |
| 1621 parameter_lambda_pascals = | |
| 1622 CRIT_VOLUME_PARAMETER_LAMBDA / 33.0f * 101325.0f; | |
| 1623 for (i = 1; i <= 16; ++i) | |
| 1624 { | |
| 1625 phase_volume_time[i - 1] = | |
| 1626 *deco_phase_volume_time + surface_phase_volume_time[i - 1]; | |
| 1627 } | |
| 1628 for (i = 1; i <= 16; ++i) | |
| 1629 { | |
| 1630 | |
| 1631 adj_crush_pressure_he_pascals = | |
| 1632 pVpm->adjusted_crushing_pressure_he[i - 1] / UNITS_FACTOR * 101325.0f; | |
| 1633 | |
| 1634 initial_allowable_grad_he_pa = | |
| 1635 pVpm->initial_allowable_gradient_he[i - 1] / UNITS_FACTOR * 101325.0f; | |
| 1636 | |
| 1637 b = initial_allowable_grad_he_pa + parameter_lambda_pascals * | |
| 1638 SURFACE_TENSION_GAMMA / ( | |
| 1639 SKIN_COMPRESSION_GAMMAC * phase_volume_time[i - 1]); | |
| 1640 | |
| 1641 c = SURFACE_TENSION_GAMMA * ( | |
| 1642 SURFACE_TENSION_GAMMA * ( | |
| 1643 parameter_lambda_pascals * adj_crush_pressure_he_pascals)) / | |
| 1644 (SKIN_COMPRESSION_GAMMAC * | |
| 1645 (SKIN_COMPRESSION_GAMMAC * phase_volume_time[i - 1])); | |
| 1646 /* Computing 2nd power */ | |
| 1647 | |
| 1648 r1 = b; | |
| 1649 | |
| 1650 new_allowable_grad_he_pascals = | |
| 1651 (b + sqrtf(r1 * r1 - c * 4.0f)) / 2.0f; | |
| 1652 | |
| 1653 /* modify global variable */ | |
| 1654 allowable_gradient_he[i - 1] = | |
| 1655 new_allowable_grad_he_pascals / 101325.0f * UNITS_FACTOR; | |
| 1656 } | |
| 1657 | |
| 1658 for (i = 1; i <= 16; ++i) | |
| 1659 { | |
| 1660 adj_crush_pressure_n2_pascals = | |
| 1661 pVpm->adjusted_crushing_pressure_n2[i - 1] / UNITS_FACTOR * 101325.0f; | |
| 1662 | |
| 1663 initial_allowable_grad_n2_pa = | |
| 1664 pVpm->initial_allowable_gradient_n2[i - 1] / UNITS_FACTOR * 101325.0f; | |
| 1665 | |
| 1666 b = initial_allowable_grad_n2_pa + parameter_lambda_pascals * | |
| 1667 SURFACE_TENSION_GAMMA / ( | |
| 1668 SKIN_COMPRESSION_GAMMAC * phase_volume_time[i - 1]); | |
| 1669 | |
| 1670 c = SURFACE_TENSION_GAMMA * | |
| 1671 (SURFACE_TENSION_GAMMA * | |
| 1672 (parameter_lambda_pascals * adj_crush_pressure_n2_pascals)) / | |
| 1673 (SKIN_COMPRESSION_GAMMAC * | |
| 1674 (SKIN_COMPRESSION_GAMMAC * phase_volume_time[i - 1])); | |
| 1675 /* Computing 2nd power */ | |
| 1676 | |
| 1677 r1 = b; | |
| 1678 | |
| 1679 new_allowable_grad_n2_pascals = | |
| 1680 (b + sqrtf(r1 * r1 - c * 4.0f)) / 2.0f; | |
| 1681 | |
| 1682 /* modify global variable */ | |
| 1683 allowable_gradient_n2[i - 1] = | |
| 1684 new_allowable_grad_n2_pascals / 101325.0f * UNITS_FACTOR; | |
| 1685 } | |
| 1686 return 0; | |
| 1687 } /* critical_volume */ | |
| 1688 | |
| 1689 /* =============================================================================== */ | |
| 1690 /* SUBROUTINE CALC_START_OF_DECO_ZONE */ | |
| 1691 /* Purpose: This subroutine uses the Bisection Method to find the depth at */ | |
| 1692 /* which the leading compartment just enters the decompression zone. */ | |
| 1693 /* Source: "Numerical Recipes in Fortran 77", Cambridge University Press, */ | |
| 1694 /* 1992. */ | |
| 1695 /* =============================================================================== */ | |
| 1696 | |
| 290 | 1697 static int calc_start_of_deco_zone(float *starting_depth, |
| 38 | 1698 float *rate, |
| 1699 float *depth_start_of_deco_zone) | |
| 1700 { | |
| 1701 /* Local variables */ | |
| 1702 float last_diff_change, | |
| 1703 initial_helium_pressure, | |
| 1704 mid_range_nitrogen_pressure; | |
| 1705 short i, j; | |
| 1706 float initial_inspired_n2_pressure, | |
| 1707 cpt_depth_start_of_deco_zone, | |
| 1708 low_bound, | |
| 1709 initial_inspired_he_pressure, | |
| 1710 high_bound_nitrogen_pressure, | |
| 1711 nitrogen_rate, | |
| 1712 function_at_mid_range, | |
| 1713 function_at_low_bound, | |
| 1714 high_bound, | |
| 1715 mid_range_helium_pressure, | |
| 1716 mid_range_time, | |
| 1717 starting_ambient_pressure, | |
| 1718 initial_nitrogen_pressure, | |
| 1719 function_at_high_bound; | |
| 1720 | |
| 1721 float time_to_start_of_deco_zone, | |
| 1722 high_bound_helium_pressure, | |
| 1723 helium_rate, | |
| 1724 differential_change; | |
| 1725 float fraction_helium_begin; | |
| 1726 float fraction_helium_end; | |
| 1727 float fraction_nitrogen_begin; | |
| 1728 float fraction_nitrogen_end; | |
| 1729 float ending_ambient_pressure; | |
| 1730 float time_test; | |
| 1731 | |
| 1732 | |
| 1733 /* loop */ | |
| 1734 /* =============================================================================== */ | |
| 1735 /* CALCULATIONS */ | |
| 1736 /* First initialize some variables */ | |
| 1737 /* =============================================================================== */ | |
| 1738 | |
| 1739 *depth_start_of_deco_zone = 0.0f; | |
| 1740 starting_ambient_pressure = *starting_depth + barometric_pressure; | |
| 1741 | |
| 1742 //>>>>>>>>>>>>>>>>>>>> | |
| 1743 //Test depth to calculate helium_rate and nitrogen_rate | |
| 1744 ending_ambient_pressure = starting_ambient_pressure/2; | |
| 1745 | |
| 1746 time_test = (ending_ambient_pressure - starting_ambient_pressure) / *rate; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1747 decom_get_inert_gases(starting_ambient_pressure / 10.0, (&pDiveSettings->decogaslist[mix_number]), &fraction_nitrogen_begin, &fraction_helium_begin ); |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1748 decom_get_inert_gases(ending_ambient_pressure / 10.0, (&pDiveSettings->decogaslist[mix_number]), &fraction_nitrogen_end, &fraction_helium_end ); |
| 38 | 1749 initial_inspired_he_pressure = (starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_helium_begin; |
| 1750 initial_inspired_n2_pressure = (starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_nitrogen_begin; | |
| 1751 helium_rate = ((ending_ambient_pressure - WATER_VAPOR_PRESSURE)* fraction_helium_end - initial_inspired_he_pressure)/time_test; | |
| 1752 nitrogen_rate = ((ending_ambient_pressure - WATER_VAPOR_PRESSURE)* fraction_nitrogen_end - initial_inspired_n2_pressure)/time_test; | |
| 1753 //>>>>>>>>>>>>>>>>>>>>> | |
| 1754 /*initial_inspired_he_pressure = | |
| 1755 (starting_ambient_pressure - water_vapor_pressure) * | |
| 1756 fraction_helium[mix_number - 1]; | |
| 1757 initial_inspired_n2_pressure = | |
| 1758 (starting_ambient_pressure - water_vapor_pressure) * | |
| 1759 fraction_nitrogen[mix_number - 1]; | |
| 1760 helium_rate = *rate * fraction_helium[mix_number - 1]; | |
| 1761 nitrogen_rate = *rate * fraction_nitrogen[mix_number - 1];*/ | |
| 1762 | |
| 1763 /* =============================================================================== */ | |
| 1764 /* ESTABLISH THE BOUNDS FOR THE ROOT SEARCH USING THE BISECTION METHOD */ | |
| 1765 /* AND CHECK TO MAKE SURE THAT THE ROOT WILL BE WITHIN BOUNDS. PROCESS */ | |
| 1766 /* EACH COMPARTMENT INDIVIDUALLY AND FIND THE MAXIMUM DEPTH ACROSS ALL */ | |
| 1767 /* COMPARTMENTS (LEADING COMPARTMENT) */ | |
| 1768 /* In this case, we are solving for time - the time when the gas tension in */ | |
| 1769 /* the compartment will be equal to ambient pressure. The low bound for time */ | |
| 1770 /* is set at zero and the high bound is set at the time it would take to */ | |
| 1771 /* ascend to zero ambient pressure (absolute). Since the ascent rate is */ | |
| 1772 /* negative, a multiplier of -1.0 is used to make the time positive. The */ | |
| 1773 /* desired point when gas tension equals ambient pressure is found at a time */ | |
| 1774 /* somewhere between these endpoints. The algorithm checks to make sure that */ | |
| 1775 /* the solution lies in between these bounds by first computing the low bound */ | |
| 1776 /* and high bound function values. */ | |
| 1777 /* =============================================================================== */ | |
| 1778 | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1779 low_bound = 0.0; |
| 38 | 1780 high_bound = starting_ambient_pressure / *rate * -1.0f; |
| 1781 for (i = 1; i <= 16; ++i) | |
| 1782 { | |
| 1783 initial_helium_pressure = helium_pressure[i - 1]; | |
| 1784 initial_nitrogen_pressure = nitrogen_pressure[i - 1]; | |
| 1785 function_at_low_bound = | |
| 1786 initial_helium_pressure + | |
| 1787 initial_nitrogen_pressure + | |
| 1788 CONSTANT_PRESSURE_OTHER_GASES - | |
| 1789 starting_ambient_pressure; | |
| 1790 high_bound_helium_pressure = | |
| 1791 schreiner_equation__2(&initial_inspired_he_pressure, | |
| 1792 &helium_rate, | |
| 1793 &high_bound, | |
| 1794 &HELIUM_TIME_CONSTANT[i - 1], | |
| 1795 &initial_helium_pressure); | |
| 1796 high_bound_nitrogen_pressure = | |
| 1797 schreiner_equation__2(&initial_inspired_n2_pressure, | |
| 1798 &nitrogen_rate, | |
| 1799 &high_bound, | |
| 1800 &NITROGEN_TIME_CONSTANT[i - 1], | |
| 1801 &initial_nitrogen_pressure); | |
| 1802 function_at_high_bound = high_bound_helium_pressure + | |
| 1803 high_bound_nitrogen_pressure + | |
| 1804 CONSTANT_PRESSURE_OTHER_GASES; | |
| 1805 if (function_at_high_bound * function_at_low_bound >= 0.0f) | |
| 1806 { | |
| 1807 printf("\nERROR! ROOT IS NOT WITHIN BRACKETS"); | |
| 1808 } | |
| 1809 | |
| 1810 /* =============================================================================== */ | |
| 1811 /* APPLY THE BISECTION METHOD IN SEVERAL ITERATIONS UNTIL A SOLUTION WITH */ | |
| 1812 /* THE DESIRED ACCURACY IS FOUND */ | |
| 1813 /* Note: the program allows for up to 100 iterations. Normally an exit will */ | |
| 1814 /* be made from the loop well before that number. If, for some reason, the */ | |
| 1815 /* program exceeds 100 iterations, there will be a pause to alert the user. */ | |
| 1816 /* =============================================================================== */ | |
| 1817 | |
| 1818 if (function_at_low_bound < 0.0f) | |
| 1819 { | |
| 1820 time_to_start_of_deco_zone = low_bound; | |
| 1821 differential_change = high_bound - low_bound; | |
| 1822 } else { | |
| 1823 time_to_start_of_deco_zone = high_bound; | |
| 1824 differential_change = low_bound - high_bound; | |
| 1825 } | |
| 1826 for (j = 1; j <= 100; ++j) | |
| 1827 { | |
| 1828 last_diff_change = differential_change; | |
| 1829 differential_change = last_diff_change * 0.5f; | |
| 1830 mid_range_time = | |
| 1831 time_to_start_of_deco_zone + | |
| 1832 differential_change; | |
| 1833 mid_range_helium_pressure = | |
| 1834 schreiner_equation__2(&initial_inspired_he_pressure, | |
| 1835 &helium_rate, | |
| 1836 &mid_range_time, | |
| 1837 &HELIUM_TIME_CONSTANT[i - 1], | |
| 1838 &initial_helium_pressure); | |
| 1839 mid_range_nitrogen_pressure = | |
| 1840 schreiner_equation__2(&initial_inspired_n2_pressure, | |
| 1841 &nitrogen_rate, | |
| 1842 &mid_range_time, | |
| 1843 &NITROGEN_TIME_CONSTANT[i - 1], | |
| 1844 &initial_nitrogen_pressure); | |
| 1845 function_at_mid_range = | |
| 1846 mid_range_helium_pressure + | |
| 1847 mid_range_nitrogen_pressure + | |
| 1848 CONSTANT_PRESSURE_OTHER_GASES - | |
| 1849 (starting_ambient_pressure + *rate * mid_range_time); | |
| 1850 if (function_at_mid_range <= 0.0f) { | |
| 1851 time_to_start_of_deco_zone = mid_range_time; | |
| 1852 } | |
| 1853 if( fabs(differential_change) < 0.001f | |
| 1854 || function_at_mid_range == 0.0f) | |
| 1855 { | |
| 1856 goto L170; | |
| 1857 } | |
| 1858 /* L150: */ | |
| 1859 } | |
| 1860 printf("\nERROR! ROOT SEARCH EXCEEDED MAXIMUM ITERATIONS"); | |
| 1861 //pause(); | |
| 1862 | |
| 1863 /* =============================================================================== */ | |
| 1864 /* When a solution with the desired accuracy is found, the program jumps out */ | |
| 1865 /* of the loop to Line 170 and assigns the solution value for the individual */ | |
| 1866 /* compartment. */ | |
| 1867 /* =============================================================================== */ | |
| 1868 | |
| 1869 L170: | |
| 1870 cpt_depth_start_of_deco_zone = | |
| 1871 starting_ambient_pressure + | |
| 1872 *rate * time_to_start_of_deco_zone - | |
| 1873 barometric_pressure; | |
| 1874 | |
| 1875 /* =============================================================================== */ | |
| 1876 /* The overall solution will be the compartment with the maximum depth where */ | |
| 1877 /* gas tension equals ambient pressure (leading compartment). */ | |
| 1878 /* =============================================================================== */ | |
| 1879 | |
| 1880 *depth_start_of_deco_zone = | |
| 1881 fmaxf(*depth_start_of_deco_zone, cpt_depth_start_of_deco_zone); | |
| 1882 /* L200: */ | |
| 1883 } | |
| 1884 return 0; | |
| 1885 } /* calc_start_of_deco_zone */ | |
| 1886 | |
| 1887 /* =============================================================================== */ | |
| 1888 /* SUBROUTINE PROJECTED_ASCENT */ | |
| 1889 /* Purpose: This subprogram performs a simulated ascent outside of the main */ | |
| 1890 /* program to ensure that a deco ceiling will not be violated due to unusual */ | |
| 1891 /* gas loading during ascent (on-gassing). If the deco ceiling is violated, */ | |
| 1892 /* the stop depth will be adjusted deeper by the step size until a safe */ | |
| 1893 /* ascent can be made. */ | |
| 1894 /* =============================================================================== */ | |
| 1895 | |
| 290 | 1896 static int projected_ascent(float *starting_depth, |
| 38 | 1897 float *rate, |
| 1898 float *deco_stop_depth, | |
| 1899 float *step_size) | |
| 1900 { | |
| 1901 /* Local variables */ | |
| 1902 float weighted_allowable_gradient, | |
| 1903 ending_ambient_pressure, | |
| 1904 temp_gas_loading[16]; | |
| 1905 int i; | |
| 1906 float allowable_gas_loading[16]; | |
| 1907 float temp_nitrogen_pressure[16]; | |
| 1908 float temp_helium_pressure[16]; | |
| 1909 float run_time_save = 0; | |
| 1910 | |
| 1911 /* loop */ | |
| 1912 /* =============================================================================== */ | |
| 1913 /* CALCULATIONS */ | |
| 1914 /* =============================================================================== */ | |
| 1915 | |
| 1916 | |
| 1917 L665: | |
| 1918 ending_ambient_pressure = *deco_stop_depth + barometric_pressure; | |
| 1919 for (i = 1; i <= 16; ++i) { | |
| 1920 temp_helium_pressure[i - 1] = helium_pressure[i - 1]; | |
| 1921 temp_nitrogen_pressure[i - 1] = nitrogen_pressure[i - 1]; | |
| 1922 } | |
| 1923 run_time_save = run_time; | |
| 1924 gas_loadings_ascent_descen(temp_helium_pressure, temp_nitrogen_pressure, *starting_depth,*deco_stop_depth,*rate,true); | |
| 1925 run_time = run_time_save; | |
| 1926 | |
| 1927 for (i = 1; i <= 16; ++i) | |
| 1928 { | |
| 1929 temp_gas_loading[i - 1] = | |
| 1930 temp_helium_pressure[i - 1] + | |
| 1931 temp_nitrogen_pressure[i - 1]; | |
| 1932 if (temp_gas_loading[i - 1] > 0.0f) | |
| 1933 { | |
| 1934 weighted_allowable_gradient = | |
| 1935 (allowable_gradient_he[i - 1] * | |
| 1936 temp_helium_pressure[i - 1] + | |
| 1937 allowable_gradient_n2[i - 1] * | |
| 1938 temp_nitrogen_pressure[i - 1]) / temp_gas_loading[i - 1]; | |
| 1939 } else { | |
| 1940 /* Computing MIN */ | |
| 1941 weighted_allowable_gradient = fminf(allowable_gradient_he[i - 1],allowable_gradient_n2[i - 1]); | |
| 1942 } | |
| 1943 allowable_gas_loading[i - 1] = | |
| 1944 ending_ambient_pressure + | |
| 1945 weighted_allowable_gradient - | |
| 1946 CONSTANT_PRESSURE_OTHER_GASES; | |
| 1947 /* L670: */ | |
| 1948 } | |
| 1949 for (i = 1; i <= 16; ++i) { | |
| 1950 if (temp_gas_loading[i - 1] > allowable_gas_loading[i - 1]) { | |
| 1951 *deco_stop_depth += *step_size; | |
| 1952 goto L665; | |
| 1953 } | |
| 1954 /* L671: */ | |
| 1955 } | |
| 1956 return 0; | |
| 1957 } /* projected_ascent */ | |
| 1958 | |
| 1959 /* =============================================================================== */ | |
| 1960 /* SUBROUTINE DECOMPRESSION_STOP */ | |
| 1961 /* Purpose: This subprogram calculates the required time at each */ | |
| 1962 /* decompression stop. */ | |
| 1963 /* =============================================================================== */ | |
| 1964 | |
| 290 | 1965 static void decompression_stop(float *deco_stop_depth, |
| 38 | 1966 float *step_size, |
| 1967 _Bool final_deco_calculation) | |
| 1968 { | |
| 1969 /* Local variables */ | |
| 1970 float inspired_nitrogen_pressure; | |
| 1971 // short last_segment_number; | |
| 1972 // float weighted_allowable_gradient; | |
| 1973 float initial_helium_pressure[16]; | |
| 1974 /* by hw */ | |
|
51
8f8ea3a32e82
Resolved warnings pointing to possible invalid memory access
Ideenmodellierer
parents:
38
diff
changeset
|
1975 float initial_CNS = gCNS_VPM; |
| 38 | 1976 |
| 1977 //static float time_counter; | |
| 1978 short i; | |
| 1979 float ambient_pressure; | |
| 1980 float inspired_helium_pressure, | |
| 1981 next_stop; | |
| 1982 //last_run_time, | |
| 1983 //temp_segment_time; | |
| 1984 | |
| 1985 float deco_ceiling_depth, | |
| 1986 initial_nitrogen_pressure[16]; | |
| 1987 //round_up_operation; | |
| 1988 float fraction_helium_begin; | |
| 1989 float fraction_nitrogen_begin; | |
| 1990 int count = 0; | |
| 1991 _Bool buehlmann_wait = false; | |
| 1992 float tissue_He_saturation[16]; | |
| 1993 float tissue_N2_saturation[16]; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
1994 float vpm_buehlmann_safety_gradient = 1.0f - (((float)pDiveSettings->vpm_conservatism) / 40.0); |
| 38 | 1995 /* loop */ |
| 1996 /* =============================================================================== */ | |
| 1997 /* CALCULATIONS */ | |
| 1998 /* =============================================================================== */ | |
| 1999 | |
| 2000 segment_time = 0; | |
| 2001 // temp_segment_time = segment_time; | |
| 2002 ambient_pressure = *deco_stop_depth + barometric_pressure; | |
| 2003 //ending_ambient_pressure = ambient_pressure; | |
| 2004 decom_get_inert_gases(ambient_pressure / 10, (&pDiveSettings->decogaslist[mix_number]), &fraction_nitrogen_begin, &fraction_helium_begin ); | |
| 2005 | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2006 if(*deco_stop_depth == (pDiveSettings->last_stop_depth_bar * 10.0)) |
| 38 | 2007 next_stop = 0; |
| 2008 else | |
| 2009 { | |
| 2010 next_stop = *deco_stop_depth - *step_size; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2011 next_stop = fmaxf(next_stop, pDiveSettings->last_stop_depth_bar * 10.0); |
| 38 | 2012 } |
| 2013 | |
| 2014 inspired_helium_pressure = | |
| 2015 (ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_helium_begin; | |
| 2016 inspired_nitrogen_pressure = | |
| 2017 (ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_nitrogen_begin; | |
| 2018 | |
| 2019 /* =============================================================================== */ | |
| 2020 /* Check to make sure that program won't lock up if unable to decompress */ | |
| 2021 /* to the next stop. If so, write error message and terminate program. */ | |
| 2022 /* =============================================================================== */ | |
| 2023 | |
| 2024 //deco_ceiling_depth = next_stop +1; //deco_ceiling_depth = next_stop + 1; | |
| 2025 if(!vpm_violates_buehlmann) | |
| 149 | 2026 { |
| 38 | 2027 calc_deco_ceiling(&deco_ceiling_depth, false); //weg, weil auf jeden Fall schleife für safety und so konservativer |
| 149 | 2028 } |
| 38 | 2029 else |
| 149 | 2030 { |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2031 deco_ceiling_depth = next_stop + 1.0; |
| 149 | 2032 } |
| 38 | 2033 if(deco_ceiling_depth > next_stop) |
| 2034 { | |
| 2035 while (deco_ceiling_depth > next_stop) | |
| 2036 { | |
| 2037 | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2038 segment_time += 60.0; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2039 if(segment_time >= 999.0 ) |
| 38 | 2040 { |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2041 segment_time = 999.0 ; |
| 38 | 2042 run_time += segment_time; |
| 2043 return; | |
| 2044 } | |
| 2045 //goto L700; | |
| 2046 initial_CNS = gCNS_VPM; | |
| 2047 decom_oxygen_calculate_cns_exposure(60*60,&pDiveSettings->decogaslist[mix_number],ambient_pressure/10,&gCNS_VPM); | |
| 2048 for (i = 0; i < 16; i++) | |
| 2049 { | |
| 2050 initial_helium_pressure[i] = helium_pressure[i]; | |
| 2051 initial_nitrogen_pressure[i] = nitrogen_pressure[i]; | |
| 2052 helium_pressure[i] += (inspired_helium_pressure - helium_pressure[i]) * float_buehlmann_He_factor_expositon_one_hour[i]; | |
| 2053 nitrogen_pressure[i] += (inspired_nitrogen_pressure - nitrogen_pressure[i]) * float_buehlmann_N2_factor_expositon_one_hour[i]; | |
| 2054 } | |
| 2055 calc_deco_ceiling(&deco_ceiling_depth, false); | |
| 2056 } | |
| 2057 if(deco_ceiling_depth < next_stop) | |
| 2058 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2059 segment_time -= 60.0; |
| 38 | 2060 gCNS_VPM = initial_CNS; |
| 2061 for (i = 0; i < 16; i++) | |
| 2062 { | |
| 2063 helium_pressure[i] = initial_helium_pressure[i]; | |
| 2064 nitrogen_pressure[i] = initial_nitrogen_pressure[i]; | |
| 2065 } | |
| 2066 deco_ceiling_depth = next_stop +1; | |
| 2067 } | |
| 2068 count = 0; | |
| 2069 while (deco_ceiling_depth > next_stop && count < 13) | |
| 2070 { | |
| 2071 count++; | |
| 2072 segment_time += 5; | |
| 2073 //goto L700; | |
| 2074 initial_CNS = gCNS_VPM; | |
| 2075 decom_oxygen_calculate_cns_exposure(60*5,&pDiveSettings->decogaslist[mix_number],ambient_pressure/10,&gCNS_VPM); | |
| 2076 for (i = 0; i < 16; i++) | |
| 2077 { | |
| 2078 initial_helium_pressure[i] = helium_pressure[i]; | |
| 2079 initial_nitrogen_pressure[i] = nitrogen_pressure[i]; | |
| 2080 helium_pressure[i] += (inspired_helium_pressure - helium_pressure[i]) * float_buehlmann_He_factor_expositon_five_minutes[i]; | |
| 2081 nitrogen_pressure[i] += (inspired_nitrogen_pressure - nitrogen_pressure[i]) * float_buehlmann_N2_factor_expositon_five_minutes[i]; | |
| 2082 } | |
| 2083 calc_deco_ceiling(&deco_ceiling_depth, false); | |
| 2084 } | |
| 2085 if(deco_ceiling_depth < next_stop) | |
| 2086 { | |
| 2087 segment_time -= 5; | |
| 2088 gCNS_VPM = initial_CNS; | |
| 2089 for (i = 0; i < 16; i++) { | |
| 2090 helium_pressure[i] = initial_helium_pressure[i]; | |
| 2091 nitrogen_pressure[i] = initial_nitrogen_pressure[i]; | |
| 2092 } | |
| 2093 deco_ceiling_depth = next_stop +1; | |
| 2094 } | |
| 2095 buehlmann_wait = false; | |
| 2096 while (buehlmann_wait || (deco_ceiling_depth > next_stop)) | |
| 2097 { | |
| 2098 //time_counter = temp_segment_time; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2099 segment_time += 1.0; |
| 38 | 2100 |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2101 if(segment_time >= 999.0 ) |
| 38 | 2102 { |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2103 segment_time = 999.0 ; |
| 38 | 2104 run_time += segment_time; |
| 2105 return; | |
| 2106 } | |
| 2107 //goto L700; | |
| 2108 initial_CNS = gCNS_VPM; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2109 decom_oxygen_calculate_cns_exposure(60*1,&pDiveSettings->decogaslist[mix_number],ambient_pressure/10.0, &gCNS_VPM); |
| 38 | 2110 for (i = 0; i < 16; i++) |
| 2111 { | |
| 2112 initial_helium_pressure[i] = helium_pressure[i]; | |
| 2113 initial_nitrogen_pressure[i] = nitrogen_pressure[i]; | |
| 2114 helium_pressure[i] += (inspired_helium_pressure - helium_pressure[i]) * float_buehlmann_He_factor_expositon_one_minute[i]; | |
| 2115 nitrogen_pressure[i] += (inspired_nitrogen_pressure - nitrogen_pressure[i]) * float_buehlmann_N2_factor_expositon_one_minute[i]; | |
| 2116 } | |
| 2117 if(!buehlmann_wait) | |
| 2118 calc_deco_ceiling(&deco_ceiling_depth, false); | |
| 2119 | |
| 2120 if(buehlmannSafety && final_deco_calculation && !(deco_ceiling_depth > next_stop)) | |
| 2121 { | |
| 2122 for (i = 0; i < 16; i++) | |
| 2123 { | |
| 2124 tissue_He_saturation[i] = helium_pressure[i] / 10; | |
| 2125 tissue_N2_saturation[i] = nitrogen_pressure[i] / 10; | |
| 2126 } | |
| 2127 if( (fabsf(nitrogen_pressure[15] - inspired_nitrogen_pressure) < 0.00001f) && (fabsf(helium_pressure[15] - inspired_helium_pressure) < 0.00001f) | |
| 2128 && (fabsf(nitrogen_pressure[0] - inspired_nitrogen_pressure) < 0.00001f) && (fabsf(helium_pressure[0] - inspired_helium_pressure) < 0.00001f)) | |
| 2129 { | |
| 2130 buehlmann_wait_exceeded = true; | |
| 2131 break; | |
| 2132 } | |
| 2133 | |
| 2134 if(decom_tissue_test_tolerance(tissue_N2_saturation, tissue_He_saturation, vpm_buehlmann_safety_gradient, (next_stop / 10.0f) + pInput->pressure_surface_bar)) | |
| 2135 break; | |
| 2136 | |
| 2137 buehlmann_wait = true; | |
| 2138 } | |
| 2139 } | |
| 2140 if(buehlmann_wait) | |
| 149 | 2141 { |
| 38 | 2142 vpm_violates_buehlmann = true; |
| 149 | 2143 } |
| 2144 if(!buehlmann_wait) | |
| 38 | 2145 { |
| 2146 if(deco_ceiling_depth < next_stop) | |
| 2147 { | |
| 2148 segment_time -= 1; | |
| 2149 gCNS_VPM = initial_CNS; | |
| 2150 for (i = 0; i < 16; i++) { | |
| 2151 helium_pressure[i] = initial_helium_pressure[i]; | |
| 2152 nitrogen_pressure[i] = initial_nitrogen_pressure[i]; | |
| 2153 } | |
| 2154 deco_ceiling_depth = next_stop +1; | |
| 2155 } | |
| 2156 while (deco_ceiling_depth > next_stop) | |
| 2157 { | |
| 2158 //time_counter = temp_segment_time; | |
| 2159 segment_time += (float) 1.0f / 3.0f; | |
| 2160 //goto L700; | |
| 2161 initial_CNS = gCNS_VPM; | |
| 2162 decom_oxygen_calculate_cns_exposure(20,&pDiveSettings->decogaslist[mix_number],ambient_pressure/10,&gCNS_VPM); | |
| 2163 for (i = 0; i < 16; i++) | |
| 2164 { | |
| 2165 helium_pressure[i] += (inspired_helium_pressure - helium_pressure[i]) * float_buehlmann_He_factor_expositon_20_seconds[i]; | |
| 2166 nitrogen_pressure[i] += (inspired_nitrogen_pressure - nitrogen_pressure[i]) * float_buehlmann_N2_factor_expositon_20_seconds[i]; | |
| 2167 } | |
| 2168 calc_deco_ceiling(&deco_ceiling_depth, false); | |
| 2169 } | |
| 2170 } | |
| 2171 } | |
| 2172 | |
| 2173 /*float pressure_save =dive_data.pressure; | |
| 2174 dive_data.pressure = ambient_pressure/10; | |
| 2175 tissues_exposure_stage(st_deco_test,(int)(segment_time * 60), &dive_data, &gaslist); | |
| 2176 dive_data.pressure = pressure_save;*/ | |
| 2177 run_time += segment_time; | |
| 2178 return; | |
| 2179 } /* decompression_stop */ | |
| 2180 | |
| 2181 /* =============================================================================== */ | |
| 2182 // SUROUTINE BOYLES_LAW_COMPENSATION | |
| 2183 // Purpose: This subprogram calculates the reduction in allowable gradients | |
| 2184 // with decreasing ambient pressure during the decompression profile based | |
| 2185 // on Boyle's Law considerations. | |
| 2186 //=============================================================================== | |
| 290 | 2187 static void BOYLES_LAW_COMPENSATION (float* First_Stop_Depth, |
| 38 | 2188 float* Deco_Stop_Depth, |
| 2189 float* Step_Size) | |
| 2190 { | |
| 2191 short i; | |
| 2192 | |
| 2193 float Next_Stop; | |
| 2194 float Ambient_Pressure_First_Stop, Ambient_Pressure_Next_Stop; | |
| 2195 float Amb_Press_First_Stop_Pascals, Amb_Press_Next_Stop_Pascals; | |
| 2196 float A, B, C, Low_Bound, High_Bound, Ending_Radius; | |
| 2197 float Deco_Gradient_Pascals; | |
| 2198 float Allow_Grad_First_Stop_He_Pa, Radius_First_Stop_He; | |
| 2199 float Allow_Grad_First_Stop_N2_Pa, Radius_First_Stop_N2; | |
| 2200 | |
| 2201 //=============================================================================== | |
| 2202 // LO//AL ARRAYS | |
| 2203 //=============================================================================== | |
| 2204 // float Radius1_He[16], Radius2_He[16]; | |
| 2205 // float Radius1_N2[16], Radius2_N2[16]; | |
| 2206 float root_factor; | |
| 2207 | |
| 2208 //=============================================================================== | |
| 2209 // CALCULATIONS | |
| 2210 //=============================================================================== | |
| 2211 Next_Stop = *Deco_Stop_Depth - *Step_Size; | |
| 2212 | |
| 2213 Ambient_Pressure_First_Stop = *First_Stop_Depth + | |
| 2214 barometric_pressure; | |
| 2215 | |
| 2216 Ambient_Pressure_Next_Stop = Next_Stop + barometric_pressure; | |
| 2217 | |
| 2218 Amb_Press_First_Stop_Pascals = (Ambient_Pressure_First_Stop/UNITS_FACTOR) * 101325.0f; | |
| 2219 | |
| 2220 Amb_Press_Next_Stop_Pascals = | |
| 2221 (Ambient_Pressure_Next_Stop/UNITS_FACTOR) * 101325.0f; | |
| 2222 root_factor = powf(Amb_Press_First_Stop_Pascals/Amb_Press_Next_Stop_Pascals,1.0f / 3.0f); | |
| 2223 | |
| 2224 for( i = 0; i < 16;i++) | |
| 2225 { | |
| 2226 Allow_Grad_First_Stop_He_Pa = | |
| 2227 (allowable_gradient_he[i]/UNITS_FACTOR) * 101325.0f; | |
| 2228 | |
| 2229 Radius_First_Stop_He = (2.0f * SURFACE_TENSION_GAMMA) / | |
| 2230 Allow_Grad_First_Stop_He_Pa; | |
| 2231 | |
| 2232 // Radius1_He[i] = Radius_First_Stop_He; | |
| 2233 A = Amb_Press_Next_Stop_Pascals; | |
| 2234 B = -2.0f * SURFACE_TENSION_GAMMA; | |
| 2235 C = (Amb_Press_First_Stop_Pascals + (2.0f * SURFACE_TENSION_GAMMA)/ | |
| 2236 Radius_First_Stop_He)* Radius_First_Stop_He* | |
| 2237 (Radius_First_Stop_He*(Radius_First_Stop_He)); | |
| 2238 Low_Bound = Radius_First_Stop_He; | |
| 2239 High_Bound = Radius_First_Stop_He * root_factor; | |
| 2240 //*pow(Amb_Press_First_Stop_Pascals/Amb_Press_Next_Stop_Pascals,1.0/3.0); | |
| 2241 //*(Amb_Press_First_Stop_Pascals/Amb_Press_Next_Stop_Pascals)**(1.0/3.0); | |
| 2242 | |
| 2243 radius_root_finder(&A,&B,&C, &Low_Bound, &High_Bound, | |
| 2244 &Ending_Radius); | |
| 2245 | |
| 2246 // Radius2_He[i] = Ending_Radius; | |
| 2247 Deco_Gradient_Pascals = (2.0f * SURFACE_TENSION_GAMMA) / | |
| 2248 Ending_Radius; | |
| 2249 | |
| 2250 deco_gradient_he[i] = (Deco_Gradient_Pascals / 101325.0f)* | |
| 2251 UNITS_FACTOR; | |
| 2252 | |
| 2253 } | |
| 2254 | |
| 2255 for( i = 0; i < 16;i++) | |
| 2256 { | |
| 2257 Allow_Grad_First_Stop_N2_Pa = | |
| 2258 (allowable_gradient_n2[i]/UNITS_FACTOR) * 101325.0f; | |
| 2259 | |
| 2260 Radius_First_Stop_N2 = (2.0f * SURFACE_TENSION_GAMMA) / | |
| 2261 Allow_Grad_First_Stop_N2_Pa; | |
| 2262 | |
| 2263 // Radius1_N2[i] = Radius_First_Stop_N2; | |
| 2264 A = Amb_Press_Next_Stop_Pascals; | |
| 2265 B = -2.0f * SURFACE_TENSION_GAMMA; | |
| 2266 C = (Amb_Press_First_Stop_Pascals + (2.0f * SURFACE_TENSION_GAMMA)/ | |
| 2267 Radius_First_Stop_N2)* Radius_First_Stop_N2* | |
| 2268 (Radius_First_Stop_N2*(Radius_First_Stop_N2)); | |
| 2269 Low_Bound = Radius_First_Stop_N2; | |
| 2270 High_Bound = Radius_First_Stop_N2* root_factor;//pow(Amb_Press_First_Stop_Pascals/Amb_Press_Next_Stop_Pascals,1.0/3.0); | |
| 2271 | |
| 2272 //High_Bound = Radius_First_Stop_N2*exp(log(Amb_Press_First_Stop_Pascals/Amb_Press_Next_Stop_Pascals)/3); | |
| 2273 radius_root_finder(&A,&B,&C, &Low_Bound, &High_Bound, | |
| 2274 &Ending_Radius); | |
| 2275 | |
| 2276 // Radius2_N2[i] = Ending_Radius; | |
| 2277 Deco_Gradient_Pascals = (2.0f * SURFACE_TENSION_GAMMA) / | |
| 2278 Ending_Radius; | |
| 2279 | |
| 2280 deco_gradient_n2[i] = (Deco_Gradient_Pascals / 101325.0f)* | |
| 2281 UNITS_FACTOR; | |
| 2282 } | |
| 2283 } | |
| 2284 | |
| 2285 /* =============================================================================== */ | |
| 292 | 2286 // vpm_calc_ndl |
| 2287 // Purpose: This function computes NDL (time where no decostops are needed) | |
| 38 | 2288 //=============================================================================== |
|
291
24ff72e627f4
Deco Models: limit NDL to 240 minutes
Jan Mulder <jlmulder@xs4all.nl>
parents:
290
diff
changeset
|
2289 #define MAX_NDL 240 |
|
24ff72e627f4
Deco Models: limit NDL to 240 minutes
Jan Mulder <jlmulder@xs4all.nl>
parents:
290
diff
changeset
|
2290 |
| 292 | 2291 static int vpm_calc_ndl(void) |
| 38 | 2292 { |
| 2293 static float future_helium_pressure[16]; | |
| 2294 static float future_nitrogen_pressure[16]; | |
| 2295 static int temp_segment_time; | |
| 2296 static int mix_number; | |
| 2297 static float inspired_helium_pressure; | |
| 2298 static float inspired_nitrogen_pressure; | |
| 2299 | |
| 2300 float previous_helium_pressure[16]; | |
| 2301 float previous_nitrogen_pressure[16]; | |
| 2302 float ambient_pressure; | |
| 2303 float fraction_helium_begin; | |
| 2304 float fraction_nitrogen_begin; | |
| 2305 int i = 0; | |
| 2306 int count = 0; | |
| 2307 int status = CALC_END; | |
|
291
24ff72e627f4
Deco Models: limit NDL to 240 minutes
Jan Mulder <jlmulder@xs4all.nl>
parents:
290
diff
changeset
|
2308 |
| 38 | 2309 for(i = 0; i < 16;i++) |
| 2310 { | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2311 future_helium_pressure[i] = pInput->tissue_helium_bar[i] * 10.0;//tissue_He_saturation[st_dive][i] * 10; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2312 future_nitrogen_pressure[i] = pInput->tissue_nitrogen_bar[i] * 10.0; |
| 38 | 2313 } |
| 2314 temp_segment_time = 0; | |
| 2315 | |
| 2316 mix_number = 0; | |
|
863
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2317 ambient_pressure = pInput->pressure_ambient_bar * 10.0; |
|
0c89c6fa949c
Bugfix empty line in deco plan (VPM only):
Ideenmodellierer
parents:
830
diff
changeset
|
2318 decom_get_inert_gases( ambient_pressure / 10.0, (&pDiveSettings->decogaslist[mix_number]) , &fraction_nitrogen_begin, &fraction_helium_begin ); |
| 38 | 2319 inspired_helium_pressure =(ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_helium_begin; |
| 2320 inspired_nitrogen_pressure =(ambient_pressure - WATER_VAPOR_PRESSURE) *fraction_nitrogen_begin; | |
| 2321 | |
| 2322 status = CALC_END; | |
| 2323 while (status == CALC_END) | |
| 2324 { | |
| 2325 count++; | |
| 2326 temp_segment_time += 60; | |
|
291
24ff72e627f4
Deco Models: limit NDL to 240 minutes
Jan Mulder <jlmulder@xs4all.nl>
parents:
290
diff
changeset
|
2327 if(temp_segment_time >= MAX_NDL) |
| 38 | 2328 { |
| 2329 pDecoInfo->output_ndl_seconds = temp_segment_time * 60; | |
| 292 | 2330 return CALC_NDL; |
| 38 | 2331 } |
| 2332 run_time += 60; | |
| 2333 //goto L700; | |
| 2334 for (i = 1; i <= 16; ++i) { | |
| 2335 previous_helium_pressure[i-1] = future_helium_pressure[i - 1]; | |
| 2336 previous_nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1]; | |
| 2337 future_helium_pressure[i - 1] = future_helium_pressure[i - 1] + (inspired_helium_pressure - future_helium_pressure[i - 1]) * float_buehlmann_He_factor_expositon_one_hour[i-1]; | |
| 2338 future_nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1] + (inspired_nitrogen_pressure - future_nitrogen_pressure[i - 1]) * float_buehlmann_N2_factor_expositon_one_hour[i-1]; | |
| 2339 helium_pressure[i - 1] = future_helium_pressure[i - 1]; | |
| 2340 nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1]; | |
| 2341 } | |
| 2342 vpm_calc_deco(); | |
| 2343 while((status = vpm_calc_critcal_volume(true,true)) == CALC_CRITICAL); | |
| 2344 | |
| 2345 } | |
| 2346 | |
| 2347 temp_segment_time -= 60; | |
| 2348 run_time -= 60; | |
| 2349 for (i = 1; i <= 16; ++i) | |
| 2350 { | |
| 2351 future_helium_pressure[i - 1] = previous_helium_pressure[i-1]; | |
| 2352 future_nitrogen_pressure[i - 1] = previous_nitrogen_pressure[i - 1]; | |
| 2353 } | |
| 2354 | |
| 2355 status = CALC_END; | |
| 2356 if(temp_segment_time < 60) | |
| 2357 nullzeit_unter60 = true; | |
| 2358 | |
| 2359 while (status == CALC_END) | |
| 2360 { | |
| 2361 temp_segment_time += 5; | |
|
291
24ff72e627f4
Deco Models: limit NDL to 240 minutes
Jan Mulder <jlmulder@xs4all.nl>
parents:
290
diff
changeset
|
2362 if(temp_segment_time >= MAX_NDL) |
| 38 | 2363 { |
| 2364 pDecoInfo->output_ndl_seconds = temp_segment_time * 60; | |
| 292 | 2365 return CALC_NDL; |
| 38 | 2366 } |
| 2367 if(nullzeit_unter60 && temp_segment_time > 60) | |
| 2368 { | |
| 2369 nullzeit_unter60 = false; | |
| 292 | 2370 return CALC_NDL; |
| 38 | 2371 } |
| 2372 run_time += 5; | |
| 2373 //goto L700; | |
| 2374 for (i = 1; i <= 16; ++i) { | |
| 2375 previous_helium_pressure[i-1] = future_helium_pressure[i - 1]; | |
| 2376 previous_nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1]; | |
| 2377 future_helium_pressure[i - 1] = future_helium_pressure[i - 1] + (inspired_helium_pressure - future_helium_pressure[i - 1]) * float_buehlmann_He_factor_expositon_five_minutes[i-1]; | |
| 2378 future_nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1] + (inspired_nitrogen_pressure - future_nitrogen_pressure[i - 1]) * float_buehlmann_N2_factor_expositon_five_minutes[i-1]; | |
| 2379 helium_pressure[i - 1] = future_helium_pressure[i - 1]; | |
| 2380 nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1]; | |
| 2381 } | |
| 2382 vpm_calc_deco(); | |
| 2383 while((status =vpm_calc_critcal_volume(true,true)) == CALC_CRITICAL); | |
| 2384 } | |
| 2385 temp_segment_time -= 5; | |
| 2386 run_time -= 5; | |
| 2387 for (i = 1; i <= 16; ++i) { | |
| 2388 future_helium_pressure[i - 1] = previous_helium_pressure[i-1]; | |
| 2389 future_nitrogen_pressure[i - 1] = previous_nitrogen_pressure[i - 1]; | |
| 2390 } | |
| 2391 status = CALC_END; | |
|
291
24ff72e627f4
Deco Models: limit NDL to 240 minutes
Jan Mulder <jlmulder@xs4all.nl>
parents:
290
diff
changeset
|
2392 |
| 38 | 2393 if(temp_segment_time <= 20) |
| 2394 { | |
| 2395 while (status == CALC_END) | |
| 2396 { | |
| 2397 temp_segment_time += minimum_deco_stop_time; | |
| 2398 run_time += minimum_deco_stop_time; | |
| 2399 //goto L700; | |
| 2400 for (i = 1; i <= 16; ++i) { | |
| 2401 future_helium_pressure[i - 1] = future_helium_pressure[i - 1] + (inspired_helium_pressure - future_helium_pressure[i - 1]) * float_buehlmann_He_factor_expositon_one_minute[i-1]; | |
| 2402 future_nitrogen_pressure[i - 1] = future_nitrogen_pressure[i - 1] + (inspired_nitrogen_pressure - future_nitrogen_pressure[i - 1]) * float_buehlmann_N2_factor_expositon_one_minute[i-1]; | |
| 2403 helium_pressure[i - 1] = future_helium_pressure[i - 1]; | |
| 2404 nitrogen_pressure[i - 1] =future_nitrogen_pressure[i - 1]; | |
| 2405 | |
| 2406 } | |
| 2407 vpm_calc_deco(); | |
| 2408 while((status =vpm_calc_critcal_volume(true,true)) == CALC_CRITICAL); | |
| 2409 | |
| 2410 } | |
| 2411 } | |
| 2412 else | |
| 2413 temp_segment_time += 5; | |
| 2414 pDecoInfo->output_ndl_seconds = temp_segment_time * 60; | |
| 2415 if(temp_segment_time > 1) | |
| 292 | 2416 return CALC_NDL; |
| 38 | 2417 else |
| 2418 return CALC_BEGIN; | |
| 2419 } |