mk2: code_part1/OSTC_code_c_part2/p2

comparison code_part1/OSTC_code_c_part2/p2_deco.c @ 237:55178aa1f972

Hunting for NoFly bug... + Asserts for NoFly range. + Bugfix surface pressure in simulator mode.

author	JeanDo
date	Sun, 20 Mar 2011 10:19:19 +0100
parents	3dbeacf42e9e
children	d995e220ddac

comparison

equal deleted inserted replaced

-:3dbeacf42e9e
+:55178aa1f972
 //---- ZH-L16 + GRADIENT FACTOR model ------------------------------------
 	if (char_I_deco_model == 1)
 	{
 // Recompute leading gas limit, at current depth:
 overlay float depth = (temp_deco - pres_surface) / 0.09995;
-assert( depth >= -0.01 );       // -epsilon tolerance.
+assert( depth >= -0.2 );        // Allow for 200mbar of weather change.
 assert( low_depth < 255 );
 if( depth > low_depth )
 sim_limit( GF_low );
 else
 void calc_hauptroutine_update_tissues(void)
 {
 assert( 0.00 <= N2_ratio && N2_ratio <= 1.00 );
 assert( 0.00 <= He_ratio && He_ratio <= 1.00 );
 assert( (N2_ratio + He_ratio) <= 0.95 );
+assert( 0.800 < pres_respiration && pres_respiration < 14.0 );
 if (char_I_const_ppO2 == 0)													// new in v.101
 		pres_diluent = pres_respiration;										// new in v.101
 	else
 {
 // Apply the Eric Baker's varying gradient factor correction.
 // Note: the correction factor depends both on GF and b,
 //       Actual values are in the 1.5 .. 1.0 range (for a GF=30%),
 //       so that can change who is the leading gas...
 // Note: Also depends of the GF_current...
-//       *BUT* calc_tissue() is used to compute bottom time,
-//       hence what would happend at surface,
-//       hence at GF_high.
 if( char_I_deco_model == 1 )
 p = ( p - var_N2_a * GF_current) * var_N2_b
 / (GF_current + var_N2_b * (1.0 - GF_current));
 else
 p = (p - var_N2_a) * var_N2_b;
-if( p < 0 ) p = 0;
+if( p < 0.0 ) p = 0.0;
 sim_pres_tissue_limit[ci] = p;
 if( p > sim_lead_tissue_limit )
 {
 sim_lead_tissue_no = ci;
 sim_lead_tissue_limit = p;
 }
 } // for ci
 assert( sim_lead_tissue_no < 16 );
-assert( 0.0 <= sim_lead_tissue_limit && sim_lead_tissue_limit <= 14.0);
+assert( 0.0 <= sim_lead_tissue_limit && sim_lead_tissue_limit <= 14.0 );
 }
 //////////////////////////////////////////////////////////////////////////////
 // clear_deco_table
 //
 // new code in v.102
 //
 static void calc_gradient_factor(void)
 {
 overlay float gf;
 assert( char_O_gtissue_no < 16 );
+assert( 0.800 <= pres_respiration && pres_respiration < 14.0 );
 	// tissue > respiration (entsaettigungsvorgang)
 	// gradient ist wieviel prozent an limit mit basis tissue
 	// dh. 0% = respiration == tissue
 	// dh. 100% = respiration == limit
 RESET_C_STACK
 assert( 800 < int_I_pres_surface && int_I_pres_surface < 1100 );
 assert( 0 < char_I_desaturation_multiplier && char_I_desaturation_multiplier <= 100 );
-N2_ratio = 0.7902; // FIXED sum as stated in b"uhlmann
+N2_ratio = 0.7902; // FIXED sum as stated in b�hlmann
 pres_surface = int_I_pres_surface * 0.001;
 ppN2 = N2_ratio * (pres_surface - ppWVapour);
 int_O_desaturation_time = 0;
 float_desaturation_multiplier = char_I_desaturation_multiplier / 142.0; // new in v.101	(70,42%/100.=142)
 // new in v.101
 //
 void deco_clear_CNS_fraction(void)
 {
 RESET_C_STACK
 CNS_fraction = 0.0;
 char_O_CNS_fraction = 0;
 }
 //////////////////////////////////////////////////////////////////////////////
 //
 void deco_calc_CNS_fraction(void)
 {
 overlay float actual_ppO2;
 RESET_C_STACK
+assert( 0.0 <= CNS_fraction && CNS_fraction <= 2.5 );
+assert( char_I_actual_ppO2 > 15 );
 actual_ppO2 = (float)char_I_actual_ppO2 / 100.0;
 if (char_I_actual_ppO2 < 50)
 CNS_fraction = CNS_fraction;// no changes
 // Output: char_O_CNS_fraction
 // Uses and Updates: CNS_fraction
 //
 void deco_calc_CNS_decrease_15min(void)
 {
 RESET_C_STACK
+assert( 0.0 <= CNS_fraction && CNS_fraction <= 2.5 );
 CNS_fraction =  0.890899 * CNS_fraction;
 char_O_CNS_fraction = (char)(CNS_fraction * 100.0 + 0.5);
 }
 //////////////////////////////////////////////////////////////////////////////
 //
 void deco_calc_percentage(void)
 {
 RESET_C_STACK
+assert( 60 <= char_I_temp && char_I_temp <= 100 );
+assert(  0 <= int_I_temp  && int_I_temp  < 2880 );      // Less than 48h...
 int_I_temp = (unsigned short)(((float)int_I_temp * (float)char_I_temp) * 0.01 );
+assert( int_I_temp < 1440 );                            // Less than 24h...
 }
 //////////////////////////////////////////////////////////////////////////////
 // deco_gas_volumes

Mercurial > public > mk2

comparison code_part1/OSTC_code_c_part2/p2_deco.c @ 237:55178aa1f972