|
0
|
1 /*
|
|
|
2 * p2_deco_main.c
|
|
|
3 *
|
|
|
4 * Created on: 31.08.2009
|
|
|
5 * Author: christian.w @ heinrichsweikamp.com
|
|
|
6 *
|
|
|
7 */
|
|
|
8
|
|
|
9 //#include <p2_deco_header_c_v102d.h>
|
|
|
10
|
|
|
11 // OSTC - diving computer code
|
|
|
12 // Copyright (C) 2009 HeinrichsWeikamp GbR
|
|
|
13
|
|
|
14 // This program is free software: you can redistribute it and/or modify
|
|
|
15 // it under the terms of the GNU General Public License as published by
|
|
|
16 // the Free Software Foundation, either version 3 of the License, or
|
|
|
17 // (at your option) any later version.
|
|
|
18
|
|
|
19 // This program is distributed in the hope that it will be useful,
|
|
|
20 // but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
21 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
|
22 // GNU General Public License for more details.
|
|
|
23
|
|
|
24 // You should have received a copy of the GNU General Public License
|
|
|
25 // along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
|
26
|
|
|
27
|
|
|
28 // *****************************
|
|
|
29 // ** I N T R O D U C T I O N **
|
|
|
30 // *****************************
|
|
|
31 //
|
|
|
32 // OSTC
|
|
|
33 //
|
|
|
34 // code:
|
|
|
35 // p2_deco_main_c_v101.c
|
|
|
36 // part2 of the OSTC code
|
|
|
37 // code with constant O2 partial pressure routines
|
|
|
38 // under construction !!
|
|
|
39 //
|
|
|
40 // summary:
|
|
|
41 // decompression routines
|
|
|
42 // for the OSTC experimental project
|
|
|
43 // written by Christian Weikamp
|
|
|
44 // last revision __________
|
|
|
45 // comments added _________
|
|
|
46 //
|
|
|
47 // additional files:
|
|
|
48 // p2_tables_v100.romdata (other files)
|
|
|
49 // 18f4685_ostc_v100.lkr (linker script)
|
|
|
50 //
|
|
|
51 // history:
|
|
|
52 // 01/03/08 v100: first release candidate
|
|
|
53 // 03/13/08 v101: start of programming ppO2 code
|
|
|
54 // 03/13/25 v101a: backup of interrim version with ppO2 calculation
|
|
|
55 // 03/13/25 v101: open circuit gas change during deco
|
|
|
56 // 03/13/25 v101: CNS_fraction calculation
|
|
|
57 // 03/13/26 v101: optimization of tissue calc routines
|
|
|
58 // 07/xx/08 v102a: debug of bottom time routine
|
|
|
59 // 09/xx/08 v102d: Gradient Factor Model implemenation
|
|
|
60 // 10/10/08 v104: renamed to build v103 for v118 stable
|
|
|
61 // 10/14/08 v104: integration of temp_depth_last_deco for Gradient Model
|
|
|
62 // 03/31/090 v107: integration of FONT Incon24
|
|
|
63
|
|
|
64 //
|
|
|
65 // literature:
|
|
|
66 // B"uhlmann, Albert: Tauchmedizin; 4. Auflage;
|
|
|
67 // Schr"oder, Kai & Reith, Steffen; 2000; S"attigungsvorg"ange beim Tauchen, das Modell ZH-L16, Funktionsweise von Tauchcomputern; http://www.achim-und-kai.de/kai/tausim/saett_faq
|
|
|
68 // Morrison, Stuart; 2000; DIY DECOMPRESSION; http://www.lizardland.co.uk/DIYDeco.html
|
|
|
69 // Balthasar, Steffen; Dekompressionstheorie I: Neo Haldane Modelle; http://www.txfreak.de/dekompressionstheorie_1.pdf
|
|
|
70 // Baker, Erik C.; Clearing Up The Confusion About "Deep Stops"
|
|
|
71 // Baker, Erik C.; Understanding M-values; http://www.txfreak.de/understanding_m-values.pdf
|
|
|
72
|
|
|
73
|
|
|
74 // *********************
|
|
|
75 // ** I N C L U D E S **
|
|
|
76 // *********************
|
|
|
77 #include <p18f4685.h>
|
|
|
78 #include <math.h>
|
|
|
79
|
|
|
80 // ********************************
|
|
|
81 // ** C O N F I G U R A T I O N **
|
|
|
82 // ** for simulation without asm **
|
|
|
83 // ********************************
|
|
|
84 #pragma config OSC = IRCIO67
|
|
|
85 #pragma config FCMEN = OFF
|
|
|
86 #pragma config IESO = OFF
|
|
|
87 #pragma config PWRT = ON
|
|
|
88 #pragma config BOREN = OFF
|
|
|
89 #pragma config WDT = OFF
|
|
|
90 #pragma config WDTPS = 128
|
|
|
91 #pragma config MCLRE = ON
|
|
|
92 #pragma config LPT1OSC = OFF
|
|
|
93 #pragma config PBADEN = OFF
|
|
|
94 #pragma config DEBUG = OFF
|
|
|
95 #pragma config XINST = OFF
|
|
|
96 #pragma config LVP = OFF
|
|
|
97 #pragma config STVREN = OFF
|
|
|
98
|
|
|
99 // ****************************
|
|
|
100 // ** D E F I N E S **
|
|
|
101 // ** missing in p18f4685.h **
|
|
|
102 // ****************************
|
|
|
103 #define INT0IF 1
|
|
|
104 #define INT1IF 0
|
|
|
105 #define TMR1IF 0
|
|
|
106
|
|
|
107
|
|
|
108 # define DBG_c_gas 0b0000000000000001
|
|
|
109 # define DBG_c_ppO2 0b0000000000000010
|
|
|
110 # define DBG_RUN 0b0000000000000100
|
|
|
111 # define DBG_RESTART 0b0000000000001000
|
|
|
112
|
|
|
113 # define DBG_CdeSAT 0b0000000000010000
|
|
|
114 # define DBG_C_MODE 0b0000000000100000
|
|
|
115 # define DBG_C_SURF 0b0000000001000000
|
|
|
116 # define DBG_HEwoHE 0b0000000010000000
|
|
|
117
|
|
|
118 # define DBG_C_DPPO2 0b0000000100000000
|
|
|
119 # define DBG_C_DGAS 0b0000001000000000
|
|
|
120 # define DBG_C_DIST 0b0000010000000000
|
|
|
121 # define DBG_C_LAST 0b0000100000000000
|
|
|
122
|
|
|
123 # define DBG_C_GF 0b0001000000000000
|
|
|
124 # define DBG_ZH16ERR 0b0010000000000000
|
|
|
125 # define DBG_PHIGH 0b0100000000000000
|
|
|
126 # define DBG_PLOW 0b1000000000000000
|
|
|
127
|
|
|
128
|
|
|
129 # define DBS_mode 0b0000000000000001
|
|
|
130 # define DBS_ppO2 0b0000000000000010
|
|
|
131 # define DBS_HE_sat 0b0000000000000100
|
|
|
132 # define DBS_ppO2chg 0b0000000000001000
|
|
|
133
|
|
|
134 # define DBS_SAT2l 0b0000000000010000
|
|
|
135 # define DBS_SAT2h 0b0000000000100000
|
|
|
136 # define DBS_GFLOW2l 0b0000000001000000
|
|
|
137 # define DBS_GFLOW2h 0b0000000010000000
|
|
|
138
|
|
|
139 # define DBS_GFHGH2l 0b0000000100000000
|
|
|
140 # define DBS_GFHGH2h 0b0000001000000000
|
|
|
141 # define DBS_GASO22l 0b0000010000000000
|
|
|
142 # define DBS_GASO22h 0b0000100000000000
|
|
|
143
|
|
|
144 # define DBS_DIST2h 0b0001000000000000
|
|
|
145 # define DBS_LAST2h 0b0010000000000000
|
|
|
146 # define DBS_DECOO2l 0b0100000000000000
|
|
|
147 # define DBS_DECOO2h 0b1000000000000000
|
|
|
148
|
|
|
149
|
|
|
150 # define DBS2_PRES2h 0b0000000000000001
|
|
|
151 # define DBS2_PRES2l 0b0000000000000010
|
|
|
152 # define DBS2_SURF2l 0b0000000000000100
|
|
|
153 # define DBS2_SURF2h 0b0000000000001000
|
|
|
154
|
|
|
155 # define DBS2_DESAT2l 0b0000000000010000
|
|
|
156 # define DBS2_DESAT2h 0b0000000000100000
|
|
|
157 # define DBS2_GFDneg 0b0000000001000000
|
|
|
158 # define DBS2_ 0b000000000000000
|
|
|
159
|
|
|
160 # define DBS2_ 0b000000000000000
|
|
|
161 # define DBS2_ 0b000000000000000
|
|
|
162 # define DBS2_ 0b000000000000000
|
|
|
163 # define DBS2_ 0b000000000000000
|
|
|
164
|
|
|
165 // NDL_at_20mtr
|
|
|
166
|
|
|
167
|
|
|
168
|
|
|
169 // ***********************
|
|
|
170 // ** V A R I A B L E S **
|
|
|
171 // ***********************
|
|
|
172 // prefixes etc:
|
|
|
173 // _O_ = output for use in the assembler code
|
|
|
174 // _I_ = input from the assembler code for the c code
|
|
|
175 // char_ and int_ = used to identify output and input size
|
|
|
176 // var = variable (from b"uhlmann)
|
|
|
177 // pres = pressure
|
|
|
178 // gtissue = guiding tissue, the one limiting the ascent
|
|
|
179 // e2secs = exp of the b"uhlmann formula precalculated for a 2 second step
|
|
|
180 // e1min = same for 1 minute step
|
|
|
181 // sim = used in simulating the ascent to the surface
|
|
|
182 // nullzeit = remaining ground/bottom time for "no deco"
|
|
|
183 // hauptroutine = main
|
|
|
184
|
|
|
185 #define WP_FONT_SMALL_HEIGHT 24
|
|
|
186 #define WP_FONT_MEDIUM_HEIGHT 32
|
|
|
187 #define WP_FONT_LARGE_HEIGHT 58
|
|
|
188
|
|
|
189 #define oled_rw PORTA,2,0
|
|
|
190 #define oled_rs PORTE,0,0
|
|
|
191
|
|
|
192 #define U8 unsigned char
|
|
|
193 #define U16 unsigned int
|
|
|
194
|
|
|
195 // IO assembler
|
|
|
196 #pragma udata bank0a=0x060
|
|
|
197 volatile unsigned char wp_stringstore[26];
|
|
|
198 volatile U8 wp_color1;
|
|
|
199 volatile U8 wp_color2;
|
|
|
200 volatile U8 wp_top;
|
|
|
201 volatile U8 wp_leftx2;
|
|
|
202 volatile U8 wp_font;
|
|
|
203 volatile U8 wp_invert;
|
|
|
204 volatile U8 wp_temp_U8;
|
|
|
205 // internal C
|
|
|
206 #pragma udata bank0b=0x081
|
|
|
207 volatile U8 wp_txtptr;
|
|
|
208 volatile unsigned char wp_char;
|
|
|
209 volatile U8 wp_command;
|
|
|
210 volatile U16 wp_data_16bit;
|
|
|
211 volatile U8 wp_data_8bit_one;
|
|
|
212 volatile U8 wp_data_8bit_two;
|
|
|
213 volatile U16 wp_start;
|
|
|
214 volatile U16 wp_end;
|
|
|
215 volatile U16 wp_i;
|
|
|
216 volatile U8 wp_black;
|
|
|
217 // some spare
|
|
|
218 volatile U8 wp_debug_U8;
|
|
|
219
|
|
|
220 // asm only
|
|
|
221 #pragma udata bank0c=0x0D0
|
|
|
222 #define LENGTH_kf_bank0 48
|
|
|
223 volatile unsigned char keep_free_bank0[LENGTH_kf_bank0];
|
|
|
224
|
|
|
225
|
|
|
226 #pragma udata bank1=0x100
|
|
|
227 #define LENGTH_kf_bank1 256
|
|
|
228 volatile unsigned char keep_free_bank1[LENGTH_kf_bank1]; // used by the assembler code
|
|
|
229
|
|
|
230 #pragma udata bank2a=0x200
|
|
|
231 // output:
|
|
|
232 static unsigned int int_O_tissue_for_debug[32];
|
|
|
233 static unsigned int int_O_GF_spare____; // 0x240
|
|
|
234 static unsigned int int_O_GF_step; // 0x242
|
|
|
235 static unsigned int int_O_gtissue_limit; // 0x244
|
|
|
236 static unsigned int int_O_gtissue_press; // 0x246
|
|
|
237 static unsigned int int_O_limit_GF_low; // 0x248
|
|
|
238 static unsigned int int_O_gtissue_press_at_GF_low; // 0x24A
|
|
|
239 // 0x24C + 0x24D noch unbenutzt
|
|
|
240
|
|
|
241 #pragma udata bank2b=0x24E
|
|
|
242 static unsigned char char_O_GF_low_pointer; // 0x24E
|
|
|
243 static unsigned char char_O_actual_pointer; // 0x24F
|
|
|
244 #pragma udata bank2c=0x250
|
|
|
245 static unsigned char char_O_deco_table[32]; // 0x250
|
|
|
246 #pragma udata bank2d=0x270
|
|
|
247 static unsigned char char_I_table_deco_done[32];
|
|
|
248 #pragma udata bank2e=0x290
|
|
|
249 static unsigned int int_O_calc_tissue_call_counter; // 0x290
|
|
|
250 // internal:
|
|
|
251 unsigned char lock_GF_depth_list;
|
|
|
252 static float temp_limit;
|
|
|
253 static float GF_low;
|
|
|
254 static float GF_high;
|
|
|
255 static float GF_delta;
|
|
|
256 static float GF_temp;
|
|
|
257 static float GF_step;
|
|
|
258 static float GF_step2;
|
|
|
259 static float temp_pres_gtissue;
|
|
|
260 static float temp_pres_gtissue_diff;
|
|
|
261 static float temp_pres_gtissue_limit_GF_low;
|
|
|
262 static float temp_pres_gtissue_limit_GF_low_below_surface;
|
|
|
263 static unsigned int temp_depth_limit;
|
|
|
264 static unsigned char temp_decotime;
|
|
|
265 static unsigned char temp_gtissue_no;
|
|
|
266 static unsigned int temp_depth_last_deco; // new in v.101
|
|
|
267
|
|
|
268 static unsigned char temp_depth_GF_low_meter;
|
|
|
269 static unsigned char temp_depth_GF_low_number;
|
|
|
270 static unsigned char internal_deco_pointer;
|
|
|
271 static unsigned char internal_deco_table[32];
|
|
|
272 static float temp_pres_deco_GF_low;
|
|
|
273
|
|
|
274 static unsigned int debug_temp;
|
|
|
275
|
|
|
276
|
|
|
277 #pragma udata bank3a=0x300
|
|
|
278 static char output[32];
|
|
|
279 // used by the math routines
|
|
|
280 #pragma udata bank3b=0x380
|
|
|
281 volatile float pres_tissue_vault[32];
|
|
|
282 #pragma udata bank4a=0x400
|
|
|
283 // internal:
|
|
|
284 unsigned char ci ; // don't move - used in _asm routines - if moved then modify movlb commands
|
|
|
285 unsigned char x;
|
|
|
286 unsigned int main_i;
|
|
|
287 unsigned int int_temp;
|
|
|
288 unsigned int int_temp_decostatus;
|
|
|
289 static float pres_respiration;
|
|
|
290 static float pres_surface;
|
|
|
291 static float temp1;
|
|
|
292 static float temp2;
|
|
|
293 static float temp3;
|
|
|
294 static float temp4;
|
|
|
295 static float temp_deco;
|
|
|
296 static float temp_atem;
|
|
|
297 static float temp2_atem;
|
|
|
298 static float temp_tissue;
|
|
|
299 static float temp_surface;
|
|
|
300 static float N2_ratio;
|
|
|
301 static float He_ratio;
|
|
|
302 static float temp_ratio;
|
|
|
303 static float var_a;
|
|
|
304 static float var2_a;
|
|
|
305 static float var_b;
|
|
|
306 static float var2_b;
|
|
|
307 static float var_t05nc;
|
|
|
308 static float var2_t05nc;
|
|
|
309 static float var_e2secs;
|
|
|
310 static float var2_e2secs;
|
|
|
311 static float var_e1min;
|
|
|
312 static float var2_e1min;
|
|
|
313 static float var_halftimes;
|
|
|
314 static float var2_halftimes;
|
|
|
315 static float pres_gtissue_limit;
|
|
|
316 static float temp_pres_gtissue_limit;
|
|
|
317 static float actual_ppO2; // new in v.102
|
|
|
318 #pragma udata bank4b=0x480
|
|
|
319 static float pres_tissue[32];
|
|
|
320
|
|
|
321 #pragma udata bank5=0x500
|
|
|
322 // don't move positions in this bank, the registers are addressed directly from assembler code
|
|
|
323 // input:
|
|
|
324 static unsigned int int_I_pres_respiration; // 0x500
|
|
|
325 static unsigned int int_I_pres_surface; // 0x502
|
|
|
326 static unsigned int int_I_temp; // 0x504 new in v101
|
|
|
327 static unsigned char char_I_temp; // 0x506 new in v101
|
|
|
328 static unsigned char char_I_actual_ppO2; // 0x507
|
|
|
329 static unsigned int int_I_spare_3;
|
|
|
330 static unsigned int int_I_spare_4;
|
|
|
331 static unsigned int int_I_spare_5;
|
|
|
332 static unsigned int int_I_spare_6;
|
|
|
333 static unsigned char char_I_N2_ratio; // 0x510
|
|
|
334 static unsigned char char_I_He_ratio; // 0x511
|
|
|
335 static unsigned char char_I_saturation_multiplier; // for conservatism/safety values 1.0 (no conservatism) to 1.5 (50% faster saturation
|
|
|
336 static unsigned char char_I_desaturation_multiplier; // for conservatism/safety values 0.66 (50% slower desaturation) to 1.0 (no conservatism)// consveratism used in calc_tissue(), calc_tissue_step_1_min() and sim_tissue_1min()
|
|
|
337 static unsigned char char_I_GF_High_percentage; // 0x514 new in v.102
|
|
|
338 static unsigned char char_I_GF_Low_percentage; // 0x515 new in v.102
|
|
|
339 static unsigned char char_I_spare; // 0x516
|
|
|
340 static unsigned char char_I_deco_distance; // 0x517
|
|
|
341 static unsigned char char_I_const_ppO2; // 0x518 new in v.101
|
|
|
342 static unsigned char char_I_deco_ppO2_change; // 0x519 new in v.101
|
|
|
343 static unsigned char char_I_deco_ppO2; // 0x51A new in v.101
|
|
|
344 static unsigned char char_I_deco_gas_change; // 0x51B new in v.101
|
|
|
345 static unsigned char char_I_deco_N2_ratio; // 0x51C new in v.101
|
|
|
346 static unsigned char char_I_deco_He_ratio; // 0x51D new in v.101
|
|
|
347 static unsigned char char_I_depth_last_deco; // 0x51E new in v.101 unit: [m]
|
|
|
348 static unsigned char char_I_deco_model; // 0x51F new in v.102 ( 1 = MultiGraF, sonst Std. mit (de-)saturation_multiplier)
|
|
|
349 // output:
|
|
|
350 static unsigned int int_O_desaturation_time; // 0x520
|
|
|
351 static unsigned char char_O_nullzeit; // 0x522
|
|
|
352 static unsigned char char_O_deco_status; // 0x523
|
|
|
353 static unsigned char char_O_array_decotime[7]; // 0x524
|
|
|
354 static unsigned char char_O_array_decodepth[6]; // 0x52B
|
|
|
355 static unsigned char char_O_ascenttime; // 0x531
|
|
|
356 static unsigned char char_O_gradient_factor; // 0x532
|
|
|
357 static unsigned char char_O_tissue_saturation[32]; // 0x533
|
|
|
358 static unsigned char char_O_array_gradient_weighted[16]; // 0x553
|
|
|
359 static unsigned char char_O_gtissue_no; // 0x563
|
|
|
360 static unsigned char char_O_diluent; // 0x564 new in v.101
|
|
|
361 static unsigned char char_O_CNS_fraction; // 0x565 new in v.101
|
|
|
362 static unsigned char char_O_relative_gradient_GF; // 0x566 new in v.102
|
|
|
363
|
|
|
364 // internal:
|
|
|
365 static float pres_tissue_limit[16];
|
|
|
366 static float sim_pres_tissue_limit[16];
|
|
|
367 static float pres_diluent; // new in v.101
|
|
|
368 static float deco_diluent; // new in v.101
|
|
|
369 static float const_ppO2; // new in v.101
|
|
|
370 static float deco_ppO2_change; // new in v.101
|
|
|
371 static float deco_ppO2; // new in v.101
|
|
|
372
|
|
|
373
|
|
|
374
|
|
|
375 #pragma udata bank6=0x600
|
|
|
376 // internal:
|
|
|
377 static float sim_pres_tissue[32];
|
|
|
378 static float sim_pres_tissue_backup[32];
|
|
|
379
|
|
|
380 //#pragma udata bank7=0x700
|
|
|
381 //const unsigned char keep_free_bank7[256]; // used by the assembler code (DD font2display)
|
|
|
382
|
|
|
383 #pragma udata bank8=0x800
|
|
|
384 static char md_pi_subst[256];
|
|
|
385
|
|
|
386 #pragma udata bank9a=0x900
|
|
|
387 // output:
|
|
|
388 static char md_state[48]; // DONT MOVE !! // has to be at the beginning of bank 9 for the asm code!!!
|
|
|
389 #pragma udata bank9b=0x930
|
|
|
390 // output:
|
|
|
391 static unsigned int int_O_DBS_bitfield; // 0x930 new in v.108
|
|
|
392 static unsigned int int_O_DBS2_bitfield; // 0x932 new in v.108
|
|
|
393 static unsigned int int_O_DBG_pre_bitfield; // 0x934 new in v.108
|
|
|
394 static unsigned int int_O_DBG_post_bitfield; // 0x936 new in v.108
|
|
|
395 static char char_O_NDL_at_20mtr; // 0x938 new in v.108 // 0xFF == undefined, max. 254
|
|
|
396 // internal:
|
|
|
397 static char md_t;
|
|
|
398 static char md_buffer[16];
|
|
|
399 static char md_cksum[16];
|
|
|
400 static char md_i;
|
|
|
401 static char md_j;
|
|
|
402 static char md_temp;
|
|
|
403 static unsigned int md_pointer;
|
|
|
404 static float deco_N2_ratio; // new in v.101
|
|
|
405 static float deco_He_ratio; // new in v.101
|
|
|
406 static float calc_N2_ratio; // new in v.101
|
|
|
407 static float calc_He_ratio; // new in v.101
|
|
|
408 static float deco_gas_change; // new in v.101
|
|
|
409 static float CNS_fraction; // new in v.101
|
|
|
410 static float float_saturation_multiplier; // new in v.101
|
|
|
411 static float float_desaturation_multiplier; // new in v.101
|
|
|
412 static float float_deco_distance; // new in v.101
|
|
|
413 // internal, dbg:
|
|
|
414 static unsigned char DBG_char_I_deco_model; // new in v.108
|
|
|
415 static unsigned char DBG_char_I_depth_last_deco; // new in v.108
|
|
|
416 static float DBG_pres_surface; // new in v.108
|
|
|
417 static float DBG_GF_low; // new in v.108
|
|
|
418 static float DBG_GF_high; // new in v.108
|
|
|
419 static float DBG_const_ppO2; // new in v.108
|
|
|
420 static float DBG_deco_ppO2_change; // new in v.108
|
|
|
421 static float DBG_deco_ppO2; // new in v.108
|
|
|
422 static float DBG_deco_N2_ratio; // new in v.108
|
|
|
423 static float DBG_deco_He_ratio; // new in v.108
|
|
|
424 static float DBG_deco_gas_change; // new in v.108
|
|
|
425 static float DBG_float_saturation_multiplier; // new in v.108
|
|
|
426 static float DBG_float_desaturation_multiplier; // new in v.108
|
|
|
427 static float DBG_float_deco_distance; // new in v.108
|
|
|
428 static float DBG_deco_N2_ratio; // new in v.108
|
|
|
429 static float DBG_deco_He_ratio; // new in v.108
|
|
|
430 static float DBG_N2_ratio; // new in v.108
|
|
|
431 static float DBG_He_ratio; // new in v.108
|
|
|
432 static char flag_in_divemode; // new in v.108
|
|
|
433 static int int_dbg_i; // new in v.108
|
|
|
434 unsigned int temp_DBS;
|
|
|
435
|
|
|
436 // *************************
|
|
|
437 // ** P R O T O T Y P E S **
|
|
|
438 // *************************
|
|
|
439 void main_calc_hauptroutine(void);
|
|
|
440 void main_calc_without_deco(void);
|
|
|
441 void main_clear_tissue(void);
|
|
|
442 void main_calc_percentage(void);
|
|
|
443 void main_calc_wo_deco_step_1_min(void);
|
|
|
444 void main_debug(void);
|
|
|
445 void main_gradient_array(void);
|
|
|
446 void main_hash(void);
|
|
|
447
|
|
|
448 void calc_hauptroutine(void);
|
|
|
449 void calc_tissue(void);
|
|
|
450 void calc_nullzeit(void);
|
|
|
451 void backup_sim_pres_tissue(void);
|
|
|
452 void restore_sim_pres_tissue(void);
|
|
|
453
|
|
|
454 void calc_without_deco(void);
|
|
|
455 void clear_tissue(void);
|
|
|
456 void calc_ascenttime(void);
|
|
|
457 void update_startvalues(void);
|
|
|
458 void clear_decoarray(void);
|
|
|
459 void update_decoarray(void);
|
|
|
460 void sim_tissue_1min(void);
|
|
|
461 void sim_tissue_10min(void);
|
|
|
462 void calc_gradient_factor(void);
|
|
|
463 void calc_gradient_array_only(void);
|
|
|
464 void calc_desaturation_time(void);
|
|
|
465 void calc_wo_deco_step_1_min(void);
|
|
|
466 void calc_tissue_step_1_min(void);
|
|
|
467 void hash(void);
|
|
|
468 void clear_CNS_fraction(void);
|
|
|
469 void calc_CNS_fraction(void);
|
|
|
470 void calc_CNS_decrease_15min(void);
|
|
|
471 void calc_percentage(void);
|
|
|
472 void main(void);
|
|
|
473 void calc_hauptroutine_data_input(void);
|
|
|
474 void calc_hauptroutine_update_tissues(void);
|
|
|
475 void calc_hauptroutine_calc_deco(void);
|
|
|
476 void calc_hauptroutine_calc_ascend_to_deco(void);
|
|
|
477 void calc_nextdecodepth_GF(void);
|
|
|
478 void copy_deco_table_GF(void);
|
|
|
479 void clear_internal_deco_table_GF(void);
|
|
|
480 void update_internal_deco_table_GF(void);
|
|
|
481 void DD2_write(void);
|
|
|
482 void DD2_write_incon42(void);
|
|
|
483 void DD2_get_pointer_to_char(void);
|
|
|
484 void DD2_set_column(void);
|
|
|
485 void DD2_load_background(void);
|
|
|
486 void DD2_build_one_line_of_char(void);
|
|
|
487 void DD2_print_column(void);
|
|
|
488 void DD2_CmdWrite(void);
|
|
|
489 void DD2_DataWrite(void);
|
|
|
490 void push_tissues_to_vault(void);
|
|
|
491 void pull_tissues_from_vault(void);
|
|
|
492 void main_push_tissues_to_vault(void);
|
|
|
493 void main_pull_tissues_from_vault(void);
|
|
|
494 void wordprocessor(void);
|
|
|
495
|
|
|
496 // *******************************
|
|
|
497 // ** start **
|
|
|
498 // ** necessary for compilation **
|
|
|
499 // *******************************
|
|
|
500 #pragma romdata der_code = 0x0000
|
|
|
501 #pragma code der_start = 0x0000
|
|
|
502 void der_start(void)
|
|
|
503 {
|
|
|
504 _asm
|
|
|
505 goto main
|
|
|
506 _endasm
|
|
|
507 }
|
|
|
508
|
|
|
509 // ***********************************
|
|
|
510 // ** main code for simulation / **
|
|
|
511 // ** tests without assembler code **
|
|
|
512 // ** is NOT a part of the OSTC **
|
|
|
513 // ***********************************
|
|
|
514 #pragma code main = 0x9000
|
|
|
515 void main(void)
|
|
|
516 {
|
|
|
517 for(wp_temp_U8=0;wp_temp_U8<LENGTH_kf_bank0 - 1;wp_temp_U8++)
|
|
|
518 keep_free_bank0[wp_temp_U8] = 7;
|
|
|
519 keep_free_bank0[LENGTH_kf_bank0 - 1] = 7;
|
|
|
520
|
|
|
521 for(wp_temp_U8=0;wp_temp_U8<LENGTH_kf_bank1 - 1;wp_temp_U8++)
|
|
|
522 keep_free_bank1[wp_temp_U8] = 7;
|
|
|
523 keep_free_bank1[LENGTH_kf_bank1 - 1] = 7;
|
|
|
524
|
|
|
525 #if 1
|
|
|
526 // new main to test DR-5
|
|
|
527
|
|
|
528 wp_top = 10;
|
|
|
529 wp_leftx2 = 10;
|
|
|
530 wp_color1 = 255;
|
|
|
531 wp_color2 = 255;
|
|
|
532 wp_font = 0;
|
|
|
533 wp_invert = 0;
|
|
|
534 wp_stringstore[0] = ' ';
|
|
|
535 wp_stringstore[1] = '!';
|
|
|
536 wp_stringstore[2] = '"';
|
|
|
537 wp_stringstore[3] = ':';
|
|
|
538 wp_stringstore[4] = 0;
|
|
|
539 wordprocessor();
|
|
|
540
|
|
|
541 GF_low = 1.0;
|
|
|
542 GF_high = 1.0;
|
|
|
543
|
|
|
544 GF_temp = GF_low * GF_high;
|
|
|
545
|
|
|
546 clear_CNS_fraction();
|
|
|
547 //char_I_const_ppO2 = 100;
|
|
|
548 //for (main_i=0;main_i<255;main_i++)
|
|
|
549 //{
|
|
|
550 //calc_CNS_fraction();
|
|
|
551 //} //for
|
|
|
552
|
|
|
553
|
|
|
554
|
|
|
555
|
|
|
556 int_I_pres_respiration = 1000;//980;
|
|
|
557 int_I_pres_surface = 1000;//980;
|
|
|
558 char_I_N2_ratio = 39; //38;
|
|
|
559 char_I_He_ratio = 40; //50;
|
|
|
560 char_I_deco_distance = 0; // 10 = 1 meter
|
|
|
561 char_I_depth_last_deco = 3; // values below 3 (meter) are ignored
|
|
|
562
|
|
|
563 char_I_const_ppO2 = 0;
|
|
|
564 char_I_deco_ppO2_change = 0; // [dm] 10 = 1 meter
|
|
|
565 char_I_deco_ppO2 = 0;
|
|
|
566
|
|
|
567 char_I_deco_gas_change = 0; // [m] 1 = 1 meter
|
|
|
568 char_I_deco_N2_ratio = 0;
|
|
|
569 char_I_deco_He_ratio = 0;
|
|
|
570
|
|
|
571 //char_I_actual_ppO2; // 0x507
|
|
|
572 char_I_GF_High_percentage = 100; // 0x514 new in v.102
|
|
|
573 char_I_GF_Low_percentage = 100; // 0x515 new in v.102
|
|
|
574
|
|
|
575 char_I_saturation_multiplier = 110;
|
|
|
576 char_I_desaturation_multiplier = 90;
|
|
|
577
|
|
|
578 char_I_deco_model = 0;
|
|
|
579
|
|
|
580 main_clear_tissue();
|
|
|
581
|
|
|
582 int_I_pres_respiration = 1000 + int_I_pres_surface;
|
|
|
583 main_calc_wo_deco_step_1_min();
|
|
|
584 int_I_pres_respiration = 3000 + int_I_pres_surface;
|
|
|
585 main_calc_wo_deco_step_1_min();
|
|
|
586 int_I_pres_respiration = 5000 + int_I_pres_surface;
|
|
|
587 main_calc_wo_deco_step_1_min();
|
|
|
588
|
|
|
589 /*
|
|
|
590 int_I_pres_respiration = 6000 + int_I_pres_surface;
|
|
|
591 for (main_i=0;main_i<27;main_i++)
|
|
|
592 main_calc_wo_deco_step_1_min();
|
|
|
593 */
|
|
|
594
|
|
|
595 char_O_deco_status = 255;
|
|
|
596 while (char_O_deco_status)
|
|
|
597 main_calc_hauptroutine();
|
|
|
598 _asm
|
|
|
599 nop
|
|
|
600 _endasm
|
|
|
601
|
|
|
602 for (main_i=0;main_i<50;main_i++)
|
|
|
603 {
|
|
|
604 main_calc_hauptroutine();
|
|
|
605 }
|
|
|
606 int_I_pres_respiration = 10000;
|
|
|
607 for (main_i=0;main_i<1500;main_i++)
|
|
|
608 {
|
|
|
609 main_calc_hauptroutine();
|
|
|
610 }
|
|
|
611
|
|
|
612 _asm
|
|
|
613 nop
|
|
|
614 _endasm
|
|
|
615
|
|
|
616
|
|
|
617 int_I_pres_respiration = 3000;
|
|
|
618 for (main_i=0;main_i<150;main_i++)
|
|
|
619 {
|
|
|
620 calc_hauptroutine_data_input();
|
|
|
621 calc_hauptroutine_update_tissues();
|
|
|
622 } //for
|
|
|
623
|
|
|
624 update_startvalues();
|
|
|
625 clear_decoarray();
|
|
|
626 clear_internal_deco_table_GF();
|
|
|
627 calc_hauptroutine_calc_ascend_to_deco();
|
|
|
628 if (char_O_deco_status > 15) // can't go up to first deco, too deep to calculate in the given time slot
|
|
|
629 {
|
|
|
630 char_O_deco_status = 2;
|
|
|
631 // char_O_lock_depth_list = 255;
|
|
|
632 }
|
|
|
633 else
|
|
|
634 {
|
|
|
635 // char_O_lock_depth_list = lock_GF_depth_list;
|
|
|
636 calc_hauptroutine_calc_deco();
|
|
|
637 }
|
|
|
638 // build_debug_output();
|
|
|
639
|
|
|
640 _asm
|
|
|
641 nop
|
|
|
642 _endasm
|
|
|
643 while (char_O_deco_status == 1)
|
|
|
644 {
|
|
|
645 char_O_deco_status = 0;
|
|
|
646 // char_O_lock_depth_list = 255;
|
|
|
647 calc_hauptroutine_calc_deco();
|
|
|
648 // build_debug_output();
|
|
|
649 _asm
|
|
|
650 nop
|
|
|
651 _endasm
|
|
|
652 };
|
|
|
653 debug_temp = 60; // [mtr Aufstieg in 10 mtr/min (30steps'2sec/min]
|
|
|
654 int_I_pres_respiration = 9980;
|
|
|
655 for (main_i=0;main_i<debug_temp;main_i++)
|
|
|
656 {
|
|
|
657 int_I_pres_respiration = int_I_pres_respiration - 33;
|
|
|
658 calc_hauptroutine_data_input();
|
|
|
659 calc_hauptroutine_update_tissues();
|
|
|
660 int_I_pres_respiration = int_I_pres_respiration - 33;
|
|
|
661 calc_hauptroutine_data_input();
|
|
|
662 calc_hauptroutine_update_tissues();
|
|
|
663 int_I_pres_respiration = int_I_pres_respiration - 34;
|
|
|
664 calc_hauptroutine_data_input();
|
|
|
665 calc_hauptroutine_update_tissues();
|
|
|
666 } //for
|
|
|
667 _asm
|
|
|
668 nop
|
|
|
669 _endasm
|
|
|
670
|
|
|
671 update_startvalues();
|
|
|
672 clear_decoarray();
|
|
|
673 clear_internal_deco_table_GF();
|
|
|
674 calc_hauptroutine_calc_ascend_to_deco();
|
|
|
675 if (char_O_deco_status > 15) // can't go up to first deco, too deep to calculate in the given time slot
|
|
|
676 {
|
|
|
677 char_O_deco_status = 2;
|
|
|
678 // char_O_lock_depth_list = 255;
|
|
|
679 }
|
|
|
680 else
|
|
|
681 {
|
|
|
682 // char_O_lock_depth_list = lock_GF_depth_list;
|
|
|
683 calc_hauptroutine_calc_deco();
|
|
|
684 }
|
|
|
685 // build_debug_output();
|
|
|
686
|
|
|
687 _asm
|
|
|
688 nop
|
|
|
689 _endasm
|
|
|
690 while (char_O_deco_status == 1)
|
|
|
691 {
|
|
|
692 char_O_deco_status = 0;
|
|
|
693 // char_O_lock_depth_list = 255;
|
|
|
694 calc_hauptroutine_calc_deco();
|
|
|
695 // build_debug_output();
|
|
|
696 _asm
|
|
|
697 nop
|
|
|
698 _endasm
|
|
|
699 };
|
|
|
700 _asm
|
|
|
701 nop
|
|
|
702 _endasm
|
|
|
703 debug_temp = 60; // [mtr Aufstieg in 10 mtr/min (30steps'2sec/min]
|
|
|
704 int_I_pres_respiration = 9980;
|
|
|
705 debug_temp = debug_temp * 3;
|
|
|
706 for (main_i=0;main_i<debug_temp;main_i++)
|
|
|
707 {
|
|
|
708 calc_hauptroutine_data_input();
|
|
|
709 calc_hauptroutine_update_tissues();
|
|
|
710 } //for
|
|
|
711 _asm
|
|
|
712 nop
|
|
|
713 _endasm
|
|
|
714 #endif
|
|
|
715 // -----------------------
|
|
|
716
|
|
|
717 } // main
|
|
|
718
|
|
|
719 // ******************************************************
|
|
|
720 // ******************************************************
|
|
|
721 // ** THE FOLLOWING CODE HAS TO BE COPPIED TO THE OSTC **
|
|
|
722 // ******************************************************
|
|
|
723 // ******************************************************
|
|
|
724
|
|
|
725 // ***************
|
|
|
726 // ***************
|
|
|
727 // ** THE FONTS **
|
|
|
728 // ***************
|
|
|
729 // ***************
|
|
|
730 // all new for bigscreen
|
|
|
731
|
|
|
732 #pragma romdata font_data_large = 0x09A00
|
|
|
733 rom const rom U16 wp_large_data[] =
|
|
|
734 {
|
|
|
735 #include "ostc90.drx.txt"
|
|
|
736 };
|
|
|
737
|
|
|
738 #pragma romdata font_table_medium = 0x0A000
|
|
|
739 rom const rom U16 wp_medium_table[] =
|
|
|
740 {
|
|
|
741 #include "ostc48.tbl.txt" // length 0x22
|
|
|
742 };
|
|
|
743
|
|
|
744 #pragma romdata font_data_medium = 0x0A024
|
|
|
745 rom const rom U16 wp_medium_data[] =
|
|
|
746 {
|
|
|
747 #include "ostc48.drx.txt" // length 0x390
|
|
|
748 };
|
|
|
749
|
|
|
750 #pragma romdata font_table_small = 0x0A39A
|
|
|
751 rom const rom U16 wp_small_table[] =
|
|
|
752 {
|
|
|
753 #include "ostc28.tbl.txt"
|
|
|
754 };
|
|
|
755
|
|
|
756 #pragma romdata font_data_small = 0x0A484
|
|
|
757 rom const rom U16 wp_small_data[] =
|
|
|
758 {
|
|
|
759 #include "ostc28.drx.txt"
|
|
|
760 };
|
|
|
761
|
|
|
762 #pragma romdata font_table_large = 0x0BEE0
|
|
|
763 rom const rom U16 wp_large_table[] =
|
|
|
764 {
|
|
|
765 0x0000
|
|
|
766 //#include "ostc90.tbl.txt"
|
|
|
767 };
|
|
|
768
|
|
|
769 // ***********************
|
|
|
770 // ***********************
|
|
|
771 // ** THE SUBROUTINES 2 **
|
|
|
772 // ***********************
|
|
|
773 // ***********************
|
|
|
774 // all new in v.102
|
|
|
775 // moved from 0x0D000 to 0x0C000 in v.108
|
|
|
776
|
|
|
777 #pragma code subroutines2 = 0x0C000 // can be adapted to fit the romdata tables ahead
|
|
|
778
|
|
|
779 // -------------------------------
|
|
|
780 // DBS - debug on start of dive //
|
|
|
781 // -------------------------------
|
|
|
782 void create_dbs_set_dbg_and_ndl20mtr(void)
|
|
|
783 {
|
|
|
784 int_O_DBS_bitfield = 0;
|
|
|
785 int_O_DBS2_bitfield = 0;
|
|
|
786 if(int_O_DBG_pre_bitfield & DBG_RUN)
|
|
|
787 int_O_DBG_pre_bitfield = DBG_RESTART;
|
|
|
788 else
|
|
|
789 int_O_DBG_pre_bitfield = DBG_RUN;
|
|
|
790 int_O_DBG_post_bitfield = 0;
|
|
|
791 char_O_NDL_at_20mtr = 255;
|
|
|
792
|
|
|
793 DBG_N2_ratio = N2_ratio;
|
|
|
794 DBG_He_ratio = He_ratio;
|
|
|
795 DBG_char_I_deco_model = char_I_deco_model;
|
|
|
796 DBG_char_I_depth_last_deco = char_I_depth_last_deco;
|
|
|
797 DBG_pres_surface = pres_surface;
|
|
|
798 DBG_GF_low = GF_low;
|
|
|
799 DBG_GF_high = GF_high;
|
|
|
800 DBG_const_ppO2 = const_ppO2;
|
|
|
801 DBG_deco_ppO2_change = deco_ppO2_change;
|
|
|
802 DBG_deco_ppO2 = deco_ppO2;
|
|
|
803 DBG_deco_N2_ratio = deco_N2_ratio;
|
|
|
804 DBG_deco_He_ratio = deco_He_ratio;
|
|
|
805 DBG_deco_gas_change = deco_gas_change;
|
|
|
806 DBG_float_saturation_multiplier = float_saturation_multiplier;
|
|
|
807 DBG_float_desaturation_multiplier = float_desaturation_multiplier;
|
|
|
808 DBG_float_deco_distance = float_deco_distance;
|
|
|
809
|
|
|
810 if(char_I_deco_model)
|
|
|
811 int_O_DBS_bitfield |= DBS_mode;
|
|
|
812 if(const_ppO2)
|
|
|
813 int_O_DBS_bitfield |= DBS_ppO2;
|
|
|
814 for(int_dbg_i = 16; int_dbg_i < 32; int_dbg_i++)
|
|
|
815 if(pres_tissue[int_dbg_i])
|
|
|
816 int_O_DBS_bitfield |= DBS_HE_sat;
|
|
|
817 if(deco_ppO2_change)
|
|
|
818 int_O_DBS_bitfield |= DBS_ppO2chg;
|
|
|
819 if(float_saturation_multiplier < 0.99)
|
|
|
820 int_O_DBS_bitfield |= DBS_SAT2l;
|
|
|
821 if(float_saturation_multiplier > 1.3)
|
|
|
822 int_O_DBS_bitfield |= DBS_SAT2h;
|
|
|
823 if(GF_low < 0.19)
|
|
|
824 int_O_DBS_bitfield |= DBS_GFLOW2l;
|
|
|
825 if(GF_low > 1.01)
|
|
|
826 int_O_DBS_bitfield |= DBS_GFLOW2h;
|
|
|
827 if(GF_high < 0.6)
|
|
|
828 int_O_DBS_bitfield |= DBS_GFHGH2l;
|
|
|
829 if(GF_high > 1.01)
|
|
|
830 int_O_DBS_bitfield |= DBS_GFHGH2h;
|
|
|
831 if((N2_ratio + He_ratio) > 0.95)
|
|
|
832 int_O_DBS_bitfield |= DBS_GASO22l;
|
|
|
833 if((N2_ratio + He_ratio) < 0.05)
|
|
|
834 int_O_DBS_bitfield |= DBS_GASO22h;
|
|
|
835 if(float_deco_distance > 0.25)
|
|
|
836 int_O_DBS_bitfield |= DBS_DIST2h;
|
|
|
837 if(char_I_depth_last_deco > 8)
|
|
|
838 int_O_DBS_bitfield |= DBS_LAST2h;
|
|
|
839 if(DBG_deco_gas_change && ((deco_N2_ratio + deco_He_ratio) > 0.95))
|
|
|
840 int_O_DBS_bitfield |= DBS_DECOO2l;
|
|
|
841 if(DBG_deco_gas_change && ((deco_N2_ratio + deco_He_ratio) < 0.05))
|
|
|
842 int_O_DBS_bitfield |= DBS_DECOO2h;
|
|
|
843 if(pres_respiration > 3.0)
|
|
|
844 int_O_DBS2_bitfield |= DBS2_PRES2h;
|
|
|
845 if(pres_surface - pres_respiration > 0.2)
|
|
|
846 int_O_DBS2_bitfield |= DBS2_PRES2l;
|
|
|
847 if(pres_surface < 0.75)
|
|
|
848 int_O_DBS2_bitfield |= DBS2_SURF2l;
|
|
|
849 if(pres_surface > 1.11)
|
|
|
850 int_O_DBS2_bitfield |= DBS2_SURF2h;
|
|
|
851 if(float_desaturation_multiplier < 0.70)
|
|
|
852 int_O_DBS2_bitfield |= DBS2_DESAT2l;
|
|
|
853 if(float_desaturation_multiplier > 1.01)
|
|
|
854 int_O_DBS2_bitfield |= DBS2_DESAT2h;
|
|
|
855 if(GF_low > GF_high)
|
|
|
856 int_O_DBS2_bitfield |= DBS2_GFDneg;
|
|
|
857 }
|
|
|
858
|
|
|
859 // -------------------------------
|
|
|
860 // DBG - set DBG to end_of_dive //
|
|
|
861 // -------------------------------
|
|
|
862 void set_dbg_end_of_dive(void)
|
|
|
863 {
|
|
|
864 int_O_DBG_pre_bitfield &= (~DBG_RUN);
|
|
|
865 int_O_DBG_post_bitfield &= (~DBG_RUN);
|
|
|
866 }
|
|
|
867
|
|
|
868 // -------------------------------
|
|
|
869 // DBG - NDL at first 20 m. hit //
|
|
|
870 // -------------------------------
|
|
|
871 void check_ndl(void)
|
|
|
872 {
|
|
|
873 if((char_O_NDL_at_20mtr == -1) && (int_I_pres_respiration > 3000))
|
|
|
874 {
|
|
|
875 char_O_NDL_at_20mtr = char_O_nullzeit;
|
|
|
876 if(char_O_NDL_at_20mtr == 255)
|
|
|
877 char_O_NDL_at_20mtr == 254;
|
|
|
878 }
|
|
|
879 }
|
|
|
880
|
|
|
881 // -------------------------------
|
|
|
882 // DBG - multi main during dive //
|
|
|
883 // -------------------------------
|
|
|
884 void check_dbg(char is_post_check)
|
|
|
885 {
|
|
|
886 temp_DBS = 0;
|
|
|
887 if( (DBG_N2_ratio != N2_ratio) || (DBG_He_ratio != He_ratio) )
|
|
|
888 temp_DBS |= DBG_c_gas;
|
|
|
889 if(DBG_const_ppO2 != const_ppO2)
|
|
|
890 temp_DBS |= DBG_c_ppO2;
|
|
|
891 if((DBG_float_saturation_multiplier != float_saturation_multiplier) || (DBG_float_desaturation_multiplier != float_desaturation_multiplier))
|
|
|
892 temp_DBS |= DBG_CdeSAT;
|
|
|
893 if(DBG_char_I_deco_model != char_I_deco_model)
|
|
|
894 temp_DBS |= DBG_C_MODE;
|
|
|
895 if(DBG_pres_surface != pres_surface)
|
|
|
896 temp_DBS |= DBG_C_SURF;
|
|
|
897 if((!DBS_HE_sat) && (!He_ratio))
|
|
|
898 for(int_dbg_i = 16; int_dbg_i < 32; int_dbg_i++)
|
|
|
899 if(pres_tissue[int_dbg_i])
|
|
|
900 temp_DBS |= DBG_HEwoHE;
|
|
|
901 if(DBG_deco_ppO2 != deco_ppO2)
|
|
|
902 temp_DBS |= DBG_C_DPPO2;
|
|
|
903 if((DBG_deco_gas_change != deco_gas_change) || (DBG_deco_N2_ratio != deco_N2_ratio) || (DBG_deco_He_ratio != deco_He_ratio))
|
|
|
904 temp_DBS |= DBG_C_DGAS;
|
|
|
905 if(DBG_float_deco_distance != float_deco_distance)
|
|
|
906 temp_DBS |= DBG_C_DIST;
|
|
|
907 if(DBG_char_I_depth_last_deco != char_I_depth_last_deco)
|
|
|
908 temp_DBS |= DBG_C_LAST;
|
|
|
909 if((DBG_GF_low != GF_low) || (DBG_GF_high != GF_high))
|
|
|
910 temp_DBS |= DBG_C_GF;
|
|
|
911 if(pres_respiration > 13.0)
|
|
|
912 temp_DBS |= DBG_PHIGH;
|
|
|
913 if(pres_surface - pres_respiration > 0.2)
|
|
|
914 temp_DBS |= DBG_PLOW;
|
|
|
915 /*
|
|
|
916 if()
|
|
|
917 temp_DBS |= ;
|
|
|
918 if()
|
|
|
919 temp_DBS |= ;
|
|
|
920 */
|
|
|
921 if(is_post_check)
|
|
|
922 int_O_DBG_post_bitfield |= temp_DBS;
|
|
|
923 else
|
|
|
924 int_O_DBG_pre_bitfield |= temp_DBS;
|
|
|
925 }
|
|
|
926
|
|
|
927 // -------------------------------
|
|
|
928 // DBG - prior to calc. of dive //
|
|
|
929 // -------------------------------
|
|
|
930 void check_pre_dbg(void)
|
|
|
931 {
|
|
|
932 check_dbg(0);
|
|
|
933 }
|
|
|
934
|
|
|
935 // -------------------------------
|
|
|
936 // DBG - after decocalc of dive //
|
|
|
937 // -------------------------------
|
|
|
938 void check_post_dbg(void)
|
|
|
939 {
|
|
|
940 check_dbg(1);
|
|
|
941 }
|
|
|
942
|
|
|
943
|
|
|
944
|
|
|
945 // -------------------------
|
|
|
946 // calc_next_decodepth_GF //
|
|
|
947 // -------------------------
|
|
|
948 // new in v.102
|
|
|
949 void calc_nextdecodepth_GF(void)
|
|
|
950 {
|
|
|
951 // INPUT, changing during dive:
|
|
|
952 // temp_pres_gtissue_limit_GF_low
|
|
|
953 // temp_pres_gtissue_limit_GF_low_below_surface
|
|
|
954 // temp_pres_gtissue
|
|
|
955 // temp_pres_gtissue_diff
|
|
|
956 // lock_GF_depth_list
|
|
|
957
|
|
|
958 // INPUT, fixed during dive:
|
|
|
959 // pres_surface
|
|
|
960 // GF_delta
|
|
|
961 // GF_high
|
|
|
962 // GF_low
|
|
|
963 // temp_depth_last_deco
|
|
|
964 // float_deco_distance
|
|
|
965
|
|
|
966 // OUTPUT
|
|
|
967 // GF_step
|
|
|
968 // temp_deco
|
|
|
969 // temp_depth_limt
|
|
|
970 // lock_GF_depth_list
|
|
|
971
|
|
|
972 // USES
|
|
|
973 // temp1
|
|
|
974 // temp2
|
|
|
975 // int_temp
|
|
|
976
|
|
|
977 char_I_table_deco_done[0] = 0; // safety if changed somewhere else. Needed for exit
|
|
|
978 if (char_I_deco_model == 1)
|
|
|
979 {
|
|
|
980 if (lock_GF_depth_list == 0)
|
|
|
981 {
|
|
|
982 temp2 = temp_pres_gtissue_limit_GF_low_below_surface / 0.29985; // = ... / 99.95 / 0.003;
|
|
|
983 int_temp = (int) (temp2 + 0.99);
|
|
|
984 if (int_temp > 31)
|
|
|
985 int_temp = 31; // deepest deco at 93 meter (31 deco stops)
|
|
|
986 if (int_temp < 0)
|
|
|
987 int_temp = 0;
|
|
|
988 temp_depth_GF_low_number = int_temp;
|
|
|
989 temp_depth_GF_low_meter = 3 * temp_depth_GF_low_number;
|
|
|
990 temp2 = (float)temp_depth_GF_low_meter * 0.09995;
|
|
|
991 temp_pres_deco_GF_low = temp2 + float_deco_distance + pres_surface;
|
|
|
992 if (temp_depth_GF_low_number == 0)
|
|
|
993 GF_step = 0;
|
|
|
994 else
|
|
|
995 GF_step = GF_delta / (float)temp_depth_GF_low_number;
|
|
|
996 if (GF_step < 0)
|
|
|
997 GF_step = 0;
|
|
|
998 if (GF_step > GF_delta)
|
|
|
999 GF_step = GF_delta;
|
|
|
1000 int_O_GF_step = (int)(GF_step * 10000);
|
|
|
1001 int_O_limit_GF_low = (int)(temp_pres_deco_GF_low * 1000);
|
|
|
1002 int_O_gtissue_press_at_GF_low = (int)(temp_pres_gtissue * 1000);
|
|
|
1003 char_O_GF_low_pointer = temp_depth_GF_low_number;
|
|
|
1004 lock_GF_depth_list = 1;
|
|
|
1005 internal_deco_pointer = 0;
|
|
|
1006 }
|
|
|
1007 if (internal_deco_pointer == 0) // new run
|
|
|
1008 {
|
|
|
1009 internal_deco_pointer = temp_depth_GF_low_number;
|
|
|
1010 GF_temp = GF_high - ((float)internal_deco_pointer * GF_step);
|
|
|
1011 int_temp = char_I_table_deco_done[internal_deco_pointer];
|
|
|
1012 output[8] = int_temp;
|
|
|
1013 output[9] = 33;
|
|
|
1014 }
|
|
|
1015 else
|
|
|
1016 {
|
|
|
1017 int_temp = 1;
|
|
|
1018 }
|
|
|
1019 while (int_temp == 1)
|
|
|
1020 {
|
|
|
1021 int_temp = internal_deco_pointer - 1;
|
|
|
1022 if (int_temp == 1) // new in v104
|
|
|
1023 {
|
|
|
1024 temp2 = (float)(temp_depth_last_deco * int_temp) * 0.09995;
|
|
|
1025 GF_step2 = GF_step/3.0 * ((float)(6 - temp_depth_last_deco));
|
|
|
1026 }
|
|
|
1027 else
|
|
|
1028 if (int_temp == 0)
|
|
|
1029 {
|
|
|
1030 temp2 = 0.0;
|
|
|
1031 GF_step2 = GF_high - GF_temp;
|
|
|
1032 }
|
|
|
1033 else
|
|
|
1034 {
|
|
|
1035 temp2 = (float)(3 *int_temp) * 0.09995;
|
|
|
1036 GF_step2 = GF_step;
|
|
|
1037 }
|
|
|
1038 temp2 = temp2 + pres_surface; // next deco stop to be tested
|
|
|
1039 temp1 = ((GF_temp + GF_step2)* temp_pres_gtissue_diff) + temp_pres_gtissue; // upper limit (lowest pressure allowed) // changes GF_step2 in v104
|
|
|
1040 if (temp1 > temp2) // check if ascent to next deco stop is ok
|
|
|
1041 {
|
|
|
1042 int_temp = 0; // no
|
|
|
1043 }
|
|
|
1044 else
|
|
|
1045 {
|
|
|
1046 internal_deco_pointer = int_temp;
|
|
|
1047 GF_temp = GF_temp + GF_step2; // changed in v104
|
|
|
1048 int_temp = char_I_table_deco_done[internal_deco_pointer]; // yes and check for ascent to even next stop if deco_done is set
|
|
|
1049 }
|
|
|
1050 } // while
|
|
|
1051 if (internal_deco_pointer > 0)
|
|
|
1052 {
|
|
|
1053 temp2 = (float)(0.29985 * internal_deco_pointer);
|
|
|
1054 temp_deco = temp2 + float_deco_distance + pres_surface;
|
|
|
1055 if (internal_deco_pointer == 1) // new in v104
|
|
|
1056 temp_depth_limit = temp_depth_last_deco;
|
|
|
1057 else
|
|
|
1058 temp_depth_limit = 3 * internal_deco_pointer;
|
|
|
1059 if (output[9] == 33)
|
|
|
1060 {
|
|
|
1061 output[9] = internal_deco_pointer;
|
|
|
1062 output[10] = char_I_table_deco_done[internal_deco_pointer];
|
|
|
1063 output[12] = output[12] + 1;
|
|
|
1064 if (output[12] == 100)
|
|
|
1065 output[12] = 0;
|
|
|
1066 }
|
|
|
1067 }
|
|
|
1068 else // if (char_I_deco_model == 1)
|
|
|
1069 {
|
|
|
1070 temp_deco = pres_surface;
|
|
|
1071 temp_depth_limit = 0;
|
|
|
1072 }
|
|
|
1073 }
|
|
|
1074 else
|
|
|
1075 {
|
|
|
1076 // calc_nextdecodepth - original
|
|
|
1077 // optimized in v.101
|
|
|
1078 // depth_last_deco included in v.101
|
|
|
1079
|
|
|
1080 temp1 = temp_pres_gtissue_limit - pres_surface;
|
|
|
1081 if (temp1 >= 0)
|
|
|
1082 {
|
|
|
1083 temp1 = temp1 / 0.29985; // = temp1 / 99.95 / 0.003;
|
|
|
1084 temp_depth_limit = (int) (temp1 + 0.99);
|
|
|
1085 temp_depth_limit = 3 * temp_depth_limit; // depth for deco [m]
|
|
|
1086 if (temp_depth_limit == 0)
|
|
|
1087 temp_deco = pres_surface;
|
|
|
1088 else
|
|
|
1089 {
|
|
|
1090 if (temp_depth_limit < temp_depth_last_deco)
|
|
|
1091 temp_depth_limit = temp_depth_last_deco;
|
|
|
1092 temp1 = (float)temp_depth_limit * 0.09995;
|
|
|
1093 temp_deco = temp1 + float_deco_distance + pres_surface; // depth for deco [bar]
|
|
|
1094 } // if (temp_depth_limit == 0)
|
|
|
1095 } // if (temp1 >= 0)
|
|
|
1096 else
|
|
|
1097 {
|
|
|
1098 temp_deco = pres_surface;
|
|
|
1099 temp_depth_limit = 0;
|
|
|
1100 } // if (temp1 >= 0)
|
|
|
1101 } // calc_nextdecodepth original
|
|
|
1102 } // calc_nextdecodepth_GF
|
|
|
1103
|
|
|
1104
|
|
|
1105 #if 0
|
|
|
1106 void build_debug_output(void)
|
|
|
1107 {
|
|
|
1108 output[0] = 0; // not used in asm PLED output
|
|
|
1109 output[1] = (int) (GF_low * 100);
|
|
|
1110 output[2] = (int) (GF_high * 100);
|
|
|
1111 output[3] = (int) (GF_step * 100);
|
|
|
1112 output[4] = (int) temp_depth_GF_low_number;
|
|
|
1113 output[5] = (int) temp_depth_GF_low_meter;
|
|
|
1114 //output[6]
|
|
|
1115 output[7] = (int) internal_deco_pointer;
|
|
|
1116 //output[8] = char_I_table_deco_done[temp_depth_GF_low_number]
|
|
|
1117 //output[9] = internal_deco_pointer @ new run
|
|
|
1118 //output[10] = char_I_table_deco_done[internal_deco_pointer] @ new run
|
|
|
1119 output [11] = (int) (temp_pres_deco_GF_low * 10);
|
|
|
1120 } // build_debug_output
|
|
|
1121 #endif
|
|
|
1122
|
|
|
1123 // ---------------------
|
|
|
1124 // copy_deco_table_GF //
|
|
|
1125 // ---------------------
|
|
|
1126 // new in v.102
|
|
|
1127 void copy_deco_table_GF(void)
|
|
|
1128 {
|
|
|
1129 if (char_I_deco_model == 1)
|
|
|
1130 {
|
|
|
1131 int_temp = 32;
|
|
|
1132 for (ci=0;ci<int_temp;ci++)
|
|
|
1133 char_O_deco_table[ci] = internal_deco_table[ci];
|
|
|
1134 }
|
|
|
1135 } // copy_deco_table_GF
|
|
|
1136
|
|
|
1137
|
|
|
1138 // ------------------------------
|
|
|
1139 // clear_internal_deco_table_GF//
|
|
|
1140 // ------------------------------
|
|
|
1141 // new in v.102
|
|
|
1142 void clear_internal_deco_table_GF(void)
|
|
|
1143 {
|
|
|
1144 if (char_I_deco_model == 1)
|
|
|
1145 {
|
|
|
1146 for (ci=0;ci<32;ci++) // cycle through the 16 b"uhlmann tissues for Helium
|
|
|
1147 {
|
|
|
1148 internal_deco_table[ci] = 0;
|
|
|
1149 }
|
|
|
1150 }
|
|
|
1151 } // clear_internal_deco_table_GF
|
|
|
1152
|
|
|
1153
|
|
|
1154 // --------------------------------
|
|
|
1155 // update_internal_deco_table_GF //
|
|
|
1156 // --------------------------------
|
|
|
1157 // new in v.102
|
|
|
1158 void update_internal_deco_table_GF(void)
|
|
|
1159 {
|
|
|
1160 if ((char_I_deco_model == 1) && (internal_deco_table[internal_deco_pointer] < 255))
|
|
|
1161 internal_deco_table[internal_deco_pointer] = internal_deco_table[internal_deco_pointer] + 1;
|
|
|
1162 } // update_internal_deco_table_GF
|
|
|
1163
|
|
|
1164
|
|
|
1165 // ---------------------
|
|
|
1166 // temp_tissue_safety //
|
|
|
1167 // ---------------------
|
|
|
1168 // outsourced in v.102
|
|
|
1169 void temp_tissue_safety(void)
|
|
|
1170 {
|
|
|
1171 if (char_I_deco_model == 1)
|
|
|
1172 {
|
|
|
1173 }
|
|
|
1174 else
|
|
|
1175 {
|
|
|
1176 if (temp_tissue < 0.0)
|
|
|
1177 temp_tissue = temp_tissue * float_desaturation_multiplier;
|
|
|
1178 else
|
|
|
1179 temp_tissue = temp_tissue * float_saturation_multiplier;
|
|
|
1180 }
|
|
|
1181 } // temp_tissue_safety
|
|
|
1182
|
|
|
1183 // -----------
|
|
|
1184 // dd2 OLD //
|
|
|
1185 // -----------
|
|
|
1186 void DD2_write(void)
|
|
|
1187 {
|
|
|
1188 _asm
|
|
|
1189 nop
|
|
|
1190 _endasm
|
|
|
1191 }
|
|
|
1192 void DD2_write_incon42(void)
|
|
|
1193 {
|
|
|
1194 DD2_write();
|
|
|
1195 }
|
|
|
1196
|
|
|
1197 void DD2_write_incon24(void)
|
|
|
1198 {
|
|
|
1199 DD2_write();
|
|
|
1200 }
|
|
|
1201 void DD2_get_pointer_to_char(void)
|
|
|
1202 {
|
|
|
1203 DD2_write();
|
|
|
1204 }
|
|
|
1205 void DD2_set_column(void)
|
|
|
1206 {
|
|
|
1207 DD2_write();
|
|
|
1208 }
|
|
|
1209 void DD2_load_background(void)
|
|
|
1210 {
|
|
|
1211 DD2_write();
|
|
|
1212 }
|
|
|
1213 void DD2_build_one_line_of_char(void)
|
|
|
1214 {
|
|
|
1215 DD2_write();
|
|
|
1216 }
|
|
|
1217 void DD2_print_column(void)
|
|
|
1218 {
|
|
|
1219 DD2_write();
|
|
|
1220 }
|
|
|
1221 void DD2_CmdWrite(void)
|
|
|
1222 {
|
|
|
1223 DD2_write();
|
|
|
1224 }
|
|
|
1225 void DD2_DataWrite(void)
|
|
|
1226 {
|
|
|
1227 DD2_write();
|
|
|
1228 }
|
|
|
1229
|
|
|
1230 // **********************
|
|
|
1231 // **********************
|
|
|
1232 // ** THE JUMP-IN CODE **
|
|
|
1233 // ** for the asm code **
|
|
|
1234 // **********************
|
|
|
1235 // **********************
|
|
|
1236 #pragma code main_calc_hauptroutine = 0x10000
|
|
|
1237 void main_calc_hauptroutine(void)
|
|
|
1238 {
|
|
|
1239 calc_hauptroutine();
|
|
|
1240 int_O_desaturation_time = 65535;
|
|
|
1241 } // divemode
|
|
|
1242 #pragma code main_without_deco = 0x10020
|
|
|
1243 void main_calc_without_deco(void)
|
|
|
1244 {
|
|
|
1245 calc_without_deco();
|
|
|
1246 calc_desaturation_time();
|
|
|
1247 }
|
|
|
1248
|
|
|
1249 #pragma code main_clear_CNS_fraction = 0x10030
|
|
|
1250 void main_clear_CNS_fraction(void)
|
|
|
1251 {
|
|
|
1252 clear_CNS_fraction();
|
|
|
1253 }
|
|
|
1254
|
|
|
1255 #pragma code main_calc_CNS_decrease_15min = 0x10034
|
|
|
1256 void main_calc_CNS_decrease_15min(void)
|
|
|
1257 {
|
|
|
1258 calc_CNS_decrease_15min();
|
|
|
1259 }
|
|
|
1260
|
|
|
1261 #pragma code main_calc_percentage = 0x10038
|
|
|
1262 void main_calc_percentage (void)
|
|
|
1263 {
|
|
|
1264 calc_percentage();
|
|
|
1265 }
|
|
|
1266
|
|
|
1267 #pragma code main_clear_tissue = 0x10040
|
|
|
1268 void main_clear_tissue(void)
|
|
|
1269 {
|
|
|
1270 clear_tissue();
|
|
|
1271 char_I_depth_last_deco = 0; // for compatibility with v.101pre_no_last_deco
|
|
|
1272 }
|
|
|
1273
|
|
|
1274 #pragma code main_calc_CNS_fraction = 0x10050
|
|
|
1275 void main_calc_CNS_fraction(void)
|
|
|
1276 {
|
|
|
1277 calc_CNS_fraction();
|
|
|
1278 }
|
|
|
1279
|
|
|
1280 #pragma code main_calc_desaturation_time = 0x10060
|
|
|
1281 void main_calc_desaturation_time(void)
|
|
|
1282 {
|
|
|
1283 calc_desaturation_time();
|
|
|
1284 }
|
|
|
1285
|
|
|
1286 #pragma code main_calc_wo_deco_step_1_min = 0x10080
|
|
|
1287 void main_calc_wo_deco_step_1_min(void)
|
|
|
1288 {
|
|
|
1289 calc_wo_deco_step_1_min();
|
|
|
1290 char_O_deco_status = 3; // surface new in v.102 overwrites value of calc_wo_deco_step_1_min
|
|
|
1291 calc_desaturation_time();
|
|
|
1292 } // surface mode
|
|
|
1293
|
|
|
1294 #pragma code main_debug = 0x100A0
|
|
|
1295 void main_debug(void)
|
|
|
1296 {
|
|
|
1297 //debug();
|
|
|
1298 }
|
|
|
1299
|
|
|
1300 #pragma code main_DD2_write_incon42 = 0x100B0
|
|
|
1301 void main_DD2_write_incon42(void)
|
|
|
1302 {
|
|
|
1303 DD2_write_incon42();
|
|
|
1304 }
|
|
|
1305
|
|
|
1306 #pragma code main_DD2_write_incon24 = 0x100B4
|
|
|
1307 void main_DD2_write_incon24(void)
|
|
|
1308 {
|
|
|
1309 DD2_write_incon24();
|
|
|
1310 }
|
|
|
1311
|
|
|
1312 #pragma code main_wordprocessor = 0x100B8
|
|
|
1313 void main_wordprocessor(void)
|
|
|
1314 {
|
|
|
1315 wordprocessor();
|
|
|
1316 }
|
|
|
1317
|
|
|
1318 #pragma code main_gradient_array = 0x100C0
|
|
|
1319 void main_gradient_array(void)
|
|
|
1320 {
|
|
|
1321 calc_gradient_array_only();
|
|
|
1322 }
|
|
|
1323 #pragma code main_push_tissues = 0x100C4
|
|
|
1324 void main_push_tissues_to_vault(void)
|
|
|
1325 {
|
|
|
1326 push_tissues_to_vault();
|
|
|
1327 }
|
|
|
1328 #pragma code main_pull_tissues = 0x100C8
|
|
|
1329 void main_pull_tissues_from_vault(void)
|
|
|
1330 {
|
|
|
1331 pull_tissues_from_vault();
|
|
|
1332 }
|
|
|
1333
|
|
|
1334 #pragma code main_hash = 0x100E0
|
|
|
1335 void main_hash(void)
|
|
|
1336 {
|
|
|
1337 hash();
|
|
|
1338 }
|
|
|
1339
|
|
|
1340 // ***********************
|
|
|
1341 // ***********************
|
|
|
1342 // ** THE LOOKUP TABLES **
|
|
|
1343 // ***********************
|
|
|
1344 // ***********************
|
|
|
1345
|
|
|
1346 #pragma romdata tables = 0x10200
|
|
|
1347 #include <p2_tables.romdata> // new table for deco_main_v.101 (var_a modified)
|
|
|
1348
|
|
|
1349 #pragma romdata tables2 = 0x10600
|
|
|
1350 rom const rom unsigned int md_pi[] =
|
|
|
1351 {
|
|
|
1352 0x292E, 0x43C9, 0xA2D8, 0x7C01, 0x3D36, 0x54A1, 0xECF0, 0x0613
|
|
|
1353 , 0x62A7, 0x05F3, 0xC0C7, 0x738C, 0x9893, 0x2BD9, 0xBC4C, 0x82CA
|
|
|
1354 , 0x1E9B, 0x573C, 0xFDD4, 0xE016, 0x6742, 0x6F18, 0x8A17, 0xE512
|
|
|
1355 , 0xBE4E, 0xC4D6, 0xDA9E, 0xDE49, 0xA0FB, 0xF58E, 0xBB2F, 0xEE7A
|
|
|
1356 , 0xA968, 0x7991, 0x15B2, 0x073F, 0x94C2, 0x1089, 0x0B22, 0x5F21
|
|
|
1357 , 0x807F, 0x5D9A, 0x5A90, 0x3227, 0x353E, 0xCCE7, 0xBFF7, 0x9703
|
|
|
1358 , 0xFF19, 0x30B3, 0x48A5, 0xB5D1, 0xD75E, 0x922A, 0xAC56, 0xAAC6
|
|
|
1359 , 0x4FB8, 0x38D2, 0x96A4, 0x7DB6, 0x76FC, 0x6BE2, 0x9C74, 0x04F1
|
|
|
1360 , 0x459D, 0x7059, 0x6471, 0x8720, 0x865B, 0xCF65, 0xE62D, 0xA802
|
|
|
1361 , 0x1B60, 0x25AD, 0xAEB0, 0xB9F6, 0x1C46, 0x6169, 0x3440, 0x7E0F
|
|
|
1362 , 0x5547, 0xA323, 0xDD51, 0xAF3A, 0xC35C, 0xF9CE, 0xBAC5, 0xEA26
|
|
|
1363 , 0x2C53, 0x0D6E, 0x8528, 0x8409, 0xD3DF, 0xCDF4, 0x4181, 0x4D52
|
|
|
1364 , 0x6ADC, 0x37C8, 0x6CC1, 0xABFA, 0x24E1, 0x7B08, 0x0CBD, 0xB14A
|
|
|
1365 , 0x7888, 0x958B, 0xE363, 0xE86D, 0xE9CB, 0xD5FE, 0x3B00, 0x1D39
|
|
|
1366 , 0xF2EF, 0xB70E, 0x6658, 0xD0E4, 0xA677, 0x72F8, 0xEB75, 0x4B0A
|
|
|
1367 , 0x3144, 0x50B4, 0x8FED, 0x1F1A, 0xDB99, 0x8D33, 0x9F11, 0x8314
|
|
|
1368 };
|
|
|
1369
|
|
|
1370 // *********************
|
|
|
1371 // *********************
|
|
|
1372 // ** THE SUBROUTINES **
|
|
|
1373 // *********************
|
|
|
1374 // *********************
|
|
|
1375
|
|
|
1376 #pragma code subroutines = 0x10700 // can be adapted to fit the romdata tables ahead
|
|
|
1377
|
|
|
1378
|
|
|
1379 // ---------------
|
|
|
1380 // CLEAR tissue //
|
|
|
1381 // ---------------
|
|
|
1382 // optimized in v.101 (var_a)
|
|
|
1383
|
|
|
1384 void clear_tissue(void) // preload tissues with standard pressure for the given ambient pressure
|
|
|
1385 {
|
|
|
1386
|
|
|
1387 flag_in_divemode = 0;
|
|
|
1388 int_O_DBS_bitfield = 0;
|
|
|
1389 int_O_DBS2_bitfield = 0;
|
|
|
1390 int_O_DBG_pre_bitfield = 0;
|
|
|
1391 int_O_DBG_post_bitfield = 0;
|
|
|
1392 char_O_NDL_at_20mtr = 255;
|
|
|
1393
|
|
|
1394 _asm
|
|
|
1395 lfsr 1, 0x300 // C math routines shall use this variable bank
|
|
|
1396 movlw 0x01
|
|
|
1397 movwf TBLPTRU,0
|
|
|
1398 _endasm
|
|
|
1399
|
|
|
1400 // N2_ratio = (float)char_I_N2_ratio; // the 0.0002 of 0.7902 are missing with standard air
|
|
|
1401 N2_ratio = 0.7902; // N2_ratio / 100.0;
|
|
|
1402 pres_respiration = (float)int_I_pres_respiration / 1000.0;
|
|
|
1403 for (ci=0;ci<16;ci++) // cycle through the 16 b"uhlmann tissues
|
|
|
1404 {
|
|
|
1405 pres_tissue[ci] = N2_ratio * (pres_respiration - 0.0627) ;
|
|
|
1406 _asm
|
|
|
1407 movlw 0x02
|
|
|
1408 movwf TBLPTRH,0
|
|
|
1409 movlb 4 // fuer ci
|
|
|
1410 movf ci,0,1
|
|
|
1411 addwf ci,0,1
|
|
|
1412 addwf ci,0,1
|
|
|
1413 addwf ci,0,1
|
|
|
1414 addlw 0x80
|
|
|
1415 movwf TBLPTRL,0
|
|
|
1416 TBLRDPOSTINC
|
|
|
1417 movff TABLAT,var_a+1
|
|
|
1418 TBLRDPOSTINC
|
|
|
1419 movff TABLAT,var_a
|
|
|
1420 TBLRDPOSTINC
|
|
|
1421 movff TABLAT,var_a+3
|
|
|
1422 TBLRD
|
|
|
1423 movff TABLAT,var_a+2
|
|
|
1424 addlw 0x80
|
|
|
1425 movwf TBLPTRL,0
|
|
|
1426 incf TBLPTRH,1,0
|
|
|
1427 TBLRDPOSTINC
|
|
|
1428 movff TABLAT,var_b+1
|
|
|
1429 TBLRDPOSTINC
|
|
|
1430 movff TABLAT,var_b
|
|
|
1431 TBLRDPOSTINC
|
|
|
1432 movff TABLAT,var_b+3
|
|
|
1433 TBLRD
|
|
|
1434 movff TABLAT,var_b+2
|
|
|
1435 _endasm
|
|
|
1436
|
|
|
1437 pres_tissue_limit[ci] = (pres_tissue[ci] - var_a) * var_b ;
|
|
|
1438 // now update the guiding tissue
|
|
|
1439 if (pres_tissue_limit[ci] < 0)
|
|
|
1440 pres_tissue_limit[ci] = 0;
|
|
|
1441 } // for 0 to 16
|
|
|
1442
|
|
|
1443 for (ci=16;ci<32;ci++) // cycle through the 16 b"uhlmann tissues for Helium
|
|
|
1444 {
|
|
|
1445 pres_tissue[ci] = 0.0;
|
|
|
1446 } // for
|
|
|
1447
|
|
|
1448 clear_decoarray();
|
|
|
1449 char_O_deco_status = 0;
|
|
|
1450 char_O_nullzeit = 0;
|
|
|
1451 char_O_ascenttime = 0;
|
|
|
1452 char_O_gradient_factor = 0;
|
|
|
1453 char_O_relative_gradient_GF = 0;
|
|
|
1454 } // clear_tissue(void)
|
|
|
1455
|
|
|
1456
|
|
|
1457 // --------------------
|
|
|
1458 // calc_without_deco //
|
|
|
1459 // fixed N2_ratio ! //
|
|
|
1460 // --------------------
|
|
|
1461 // optimized in v.101 (float_..saturation_multiplier)
|
|
|
1462
|
|
|
1463 void calc_without_deco(void)
|
|
|
1464 {
|
|
|
1465 _asm
|
|
|
1466 lfsr 1, 0x300
|
|
|
1467 _endasm
|
|
|
1468 N2_ratio = 0.7902; // FIXED RATIO !! sum as stated in b"uhlmann
|
|
|
1469 pres_respiration = (float)int_I_pres_respiration / 1000.0; // assembler code uses different digit system
|
|
|
1470 pres_surface = (float)int_I_pres_surface / 1000.0;
|
|
|
1471 temp_atem = N2_ratio * (pres_respiration - 0.0627); // 0.0627 is the extra pressure in the body
|
|
|
1472 temp2_atem = 0.0;
|
|
|
1473 temp_surface = pres_surface; // the b"uhlmann formula using temp_surface does apply to the pressure without any inert ratio
|
|
|
1474 float_desaturation_multiplier = char_I_desaturation_multiplier / 100.0;
|
|
|
1475 float_saturation_multiplier = char_I_saturation_multiplier / 100.0;
|
|
|
1476
|
|
|
1477 calc_tissue(); // update the pressure in the 16 tissues in accordance with the new ambient pressure
|
|
|
1478
|
|
|
1479 clear_decoarray();
|
|
|
1480 char_O_deco_status = 0;
|
|
|
1481 char_O_nullzeit = 0;
|
|
|
1482 char_O_ascenttime = 0;
|
|
|
1483 calc_gradient_factor();
|
|
|
1484
|
|
|
1485 } // calc_without_deco
|
|
|
1486
|
|
|
1487
|
|
|
1488 // --------------------
|
|
|
1489 // calc_hauptroutine //
|
|
|
1490 // --------------------
|
|
|
1491 // this is the major code in dive mode
|
|
|
1492 // calculates:
|
|
|
1493 // the tissues,
|
|
|
1494 // the bottom time
|
|
|
1495 // and simulates the ascend with all deco stops
|
|
|
1496
|
|
|
1497 void calc_hauptroutine(void)
|
|
|
1498 {
|
|
|
1499 calc_hauptroutine_data_input();
|
|
|
1500
|
|
|
1501 if(!flag_in_divemode)
|
|
|
1502 {
|
|
|
1503 flag_in_divemode = 1;
|
|
|
1504 create_dbs_set_dbg_and_ndl20mtr();
|
|
|
1505 }
|
|
|
1506 else
|
|
|
1507 check_pre_dbg();
|
|
|
1508
|
|
|
1509 calc_hauptroutine_update_tissues();
|
|
|
1510 calc_gradient_factor();
|
|
|
1511
|
|
|
1512
|
|
|
1513 switch (char_O_deco_status) // toggle between calculation for nullzeit (bottom time), deco stops and more deco stops (continue)
|
|
|
1514 {
|
|
|
1515 case 0:
|
|
|
1516 update_startvalues();
|
|
|
1517 calc_nullzeit();
|
|
|
1518 check_ndl();
|
|
|
1519 char_O_deco_status = 255; // calc deco next time
|
|
|
1520 break;
|
|
|
1521 case 1:
|
|
|
1522 if (char_O_deco_status == 3)
|
|
|
1523 break;
|
|
|
1524 char_O_deco_status = 0;
|
|
|
1525 // char_O_lock_depth_list = 255;
|
|
|
1526 calc_hauptroutine_calc_deco();
|
|
|
1527 // build_debug_output();
|
|
|
1528 break;
|
|
|
1529 case 3: // new dive
|
|
|
1530 clear_decoarray();
|
|
|
1531 clear_internal_deco_table_GF();
|
|
|
1532 copy_deco_table_GF();
|
|
|
1533 internal_deco_pointer = 0;
|
|
|
1534 lock_GF_depth_list = 0;
|
|
|
1535 update_startvalues();
|
|
|
1536 calc_nextdecodepth_GF();
|
|
|
1537 char_O_deco_status = 0;
|
|
|
1538 break;
|
|
|
1539 default:
|
|
|
1540 update_startvalues();
|
|
|
1541 clear_decoarray();
|
|
|
1542 clear_internal_deco_table_GF();
|
|
|
1543 output[6] = 1;
|
|
|
1544 calc_hauptroutine_calc_ascend_to_deco();
|
|
|
1545 if (char_O_deco_status > 15) // can't go up to first deco, too deep to calculate in the given time slot
|
|
|
1546 {
|
|
|
1547 char_O_deco_status = 2;
|
|
|
1548 // char_O_lock_depth_list = 255;
|
|
|
1549 }
|
|
|
1550 else
|
|
|
1551 {
|
|
|
1552 // char_O_lock_depth_list = lock_GF_depth_list;
|
|
|
1553 calc_hauptroutine_calc_deco();
|
|
|
1554 }
|
|
|
1555 // build_debug_output();
|
|
|
1556 break;
|
|
|
1557 }
|
|
|
1558 calc_ascenttime();
|
|
|
1559 check_post_dbg();
|
|
|
1560 }
|
|
|
1561
|
|
|
1562 void calc_hauptroutine_data_input(void)
|
|
|
1563 {
|
|
|
1564 pres_respiration = (float)int_I_pres_respiration / 1000.0;
|
|
|
1565 pres_surface = (float)int_I_pres_surface / 1000.0;
|
|
|
1566
|
|
|
1567 N2_ratio = (float)char_I_N2_ratio / 100.0;; // the 0.0002 of 0.7902 are missing with standard air
|
|
|
1568 He_ratio = (float)char_I_He_ratio / 100.0;;
|
|
|
1569 deco_N2_ratio = (float)char_I_deco_N2_ratio / 100.0;
|
|
|
1570 deco_He_ratio = (float)char_I_deco_He_ratio / 100.0;
|
|
|
1571 float_deco_distance = (float)char_I_deco_distance / 100.0;
|
|
|
1572 if(char_I_deco_gas_change)
|
|
|
1573 {
|
|
|
1574 deco_gas_change = (float)char_I_deco_gas_change / 9.995 + pres_surface;
|
|
|
1575 deco_gas_change = deco_gas_change + float_deco_distance;
|
|
|
1576 }
|
|
|
1577 else
|
|
|
1578 deco_gas_change = 0;
|
|
|
1579 const_ppO2 = (float)char_I_const_ppO2 / 100.0;
|
|
|
1580 deco_ppO2_change = (float)char_I_deco_ppO2_change / 99.95 + pres_surface;
|
|
|
1581 deco_ppO2_change = deco_ppO2_change + float_deco_distance;
|
|
|
1582 deco_ppO2 = (float)char_I_deco_ppO2 / 100.0;
|
|
|
1583 float_desaturation_multiplier = char_I_desaturation_multiplier / 100.0;
|
|
|
1584 float_saturation_multiplier = char_I_saturation_multiplier / 100.0;
|
|
|
1585 GF_low = (float)char_I_GF_Low_percentage / 100.0;
|
|
|
1586 GF_high = (float)char_I_GF_High_percentage / 100.0;
|
|
|
1587 GF_delta = GF_high - GF_low;
|
|
|
1588
|
|
|
1589 temp2 = (pres_respiration - pres_surface) / 0.29985;
|
|
|
1590 int_temp = (int)(temp2);
|
|
|
1591 if (int_temp < 0)
|
|
|
1592 int_temp = 0;
|
|
|
1593 if (int_temp > 255)
|
|
|
1594 int_temp = 255;
|
|
|
1595 char_O_actual_pointer = int_temp;
|
|
|
1596
|
|
|
1597 temp_depth_last_deco = (int)char_I_depth_last_deco;
|
|
|
1598 }
|
|
|
1599
|
|
|
1600 void calc_hauptroutine_update_tissues(void)
|
|
|
1601 {
|
|
|
1602 int_O_calc_tissue_call_counter = int_O_calc_tissue_call_counter + 1;
|
|
|
1603 if (char_I_const_ppO2 == 0) // new in v.101
|
|
|
1604 pres_diluent = pres_respiration; // new in v.101
|
|
|
1605 else // new in v.101
|
|
|
1606 pres_diluent = ((pres_respiration - const_ppO2)/(N2_ratio + He_ratio)); // new in v.101
|
|
|
1607 if (pres_diluent > pres_respiration) // new in v.101
|
|
|
1608 pres_diluent = pres_respiration; // new in v.101
|
|
|
1609 if (pres_diluent > 0.0627) // new in v.101
|
|
|
1610 {
|
|
|
1611 temp_atem = N2_ratio * (pres_diluent - 0.0627); // changed in v.101
|
|
|
1612 temp2_atem = He_ratio * (pres_diluent - 0.0627); // changed in v.101
|
|
|
1613 char_O_diluent = (char)(pres_diluent/pres_respiration*100.0);
|
|
|
1614 }
|
|
|
1615 else // new in v.101
|
|
|
1616 {
|
|
|
1617 temp_atem = 0.0; // new in v.101
|
|
|
1618 temp2_atem = 0.0; // new in v.101
|
|
|
1619 char_O_diluent = 0;
|
|
|
1620 }
|
|
|
1621 temp_surface = pres_surface;
|
|
|
1622 calc_tissue();
|
|
|
1623 int_O_gtissue_limit = (int)(pres_tissue_limit[char_O_gtissue_no] * 1000);
|
|
|
1624 int_O_gtissue_press = (int)((pres_tissue[char_O_gtissue_no] + pres_tissue[char_O_gtissue_no+16]) * 1000);
|
|
|
1625 if (char_I_deco_model == 1)
|
|
|
1626 {
|
|
|
1627 temp1 = temp1 * GF_high;
|
|
|
1628 }
|
|
|
1629 else
|
|
|
1630 {
|
|
|
1631 temp1 = temp_surface;
|
|
|
1632 }
|
|
|
1633 if (pres_gtissue_limit > temp1 && char_O_deco_status == 0) // if guiding tissue can not be exposed to surface pressure immediately
|
|
|
1634 {
|
|
|
1635 char_O_nullzeit = 0; // deco necessary
|
|
|
1636 char_O_deco_status = 255; // calculate deco skip nullzeit calculation
|
|
|
1637 }
|
|
|
1638 } // calc_hauptroutine_update_tissues
|
|
|
1639 void calc_hauptroutine_calc_deco(void)
|
|
|
1640 {
|
|
|
1641 do
|
|
|
1642 {
|
|
|
1643 int_temp_decostatus = 0;
|
|
|
1644 calc_nextdecodepth_GF();
|
|
|
1645 if (temp_depth_limit > 0)
|
|
|
1646 {
|
|
|
1647 if (char_I_const_ppO2 == 0) // new in v.101
|
|
|
1648 {
|
|
|
1649 deco_diluent = temp_deco; // new in v.101
|
|
|
1650 if (temp_deco > deco_gas_change)
|
|
|
1651 {
|
|
|
1652 calc_N2_ratio = N2_ratio;
|
|
|
1653 calc_He_ratio = He_ratio;
|
|
|
1654 }
|
|
|
1655 else
|
|
|
1656 {
|
|
|
1657 calc_N2_ratio = deco_N2_ratio;
|
|
|
1658 calc_He_ratio = deco_He_ratio;
|
|
|
1659 }
|
|
|
1660 }
|
|
|
1661 else // new in v.101
|
|
|
1662 {
|
|
|
1663 calc_N2_ratio = N2_ratio;
|
|
|
1664 calc_He_ratio = He_ratio;
|
|
|
1665 if (temp_deco > deco_ppO2_change)
|
|
|
1666 {
|
|
|
1667 deco_diluent = ((temp_deco - const_ppO2)/(N2_ratio + He_ratio)); // new in v.101
|
|
|
1668 }
|
|
|
1669 else
|
|
|
1670 {
|
|
|
1671 deco_diluent = ((temp_deco - deco_ppO2)/(N2_ratio + He_ratio)); // new in v.101
|
|
|
1672 }
|
|
|
1673 }
|
|
|
1674 if (deco_diluent > temp_deco) // new in v.101
|
|
|
1675 deco_diluent = temp_deco; // new in v.101
|
|
|
1676 if (deco_diluent > 0.0627) // new in v.101
|
|
|
1677 {
|
|
|
1678 temp_atem = calc_N2_ratio * (deco_diluent - 0.0627); // changed in v.101
|
|
|
1679 temp2_atem = calc_He_ratio * (deco_diluent - 0.0627); // changed in v.101
|
|
|
1680 }
|
|
|
1681 else // new in v.101
|
|
|
1682 {
|
|
|
1683 temp_atem = 0.0; // new in v.101
|
|
|
1684 temp2_atem = 0.0; // new in v.101
|
|
|
1685 }
|
|
|
1686 sim_tissue_1min();
|
|
|
1687 update_internal_deco_table_GF();
|
|
|
1688 temp_decotime = 1;
|
|
|
1689 update_decoarray();
|
|
|
1690 char_O_deco_status = char_O_deco_status + 1;
|
|
|
1691 if (char_O_deco_status < 16)
|
|
|
1692 int_temp_decostatus = 1;
|
|
|
1693 }
|
|
|
1694 else // if (temp_depth_limit > 0)
|
|
|
1695 {
|
|
|
1696 char_O_deco_status = 0;
|
|
|
1697 }
|
|
|
1698 } while (int_temp_decostatus == 1);
|
|
|
1699 if (char_O_deco_status > 15)
|
|
|
1700 {
|
|
|
1701 char_O_deco_status = 1;
|
|
|
1702 }
|
|
|
1703 else
|
|
|
1704 {
|
|
|
1705 copy_deco_table_GF();
|
|
|
1706 char_O_deco_status = 0;
|
|
|
1707 }
|
|
|
1708 }
|
|
|
1709
|
|
|
1710 void calc_hauptroutine_calc_ascend_to_deco(void)
|
|
|
1711 {
|
|
|
1712 update_startvalues();
|
|
|
1713 char_O_deco_status = 0;
|
|
|
1714 temp_deco = pres_respiration;
|
|
|
1715 lock_GF_depth_list = 1; // new in v.102
|
|
|
1716 do // go up to first deco
|
|
|
1717 {
|
|
|
1718 int_temp_decostatus = 0;
|
|
|
1719 temp_deco = temp_deco - 1.0;
|
|
|
1720 if ( char_I_deco_model == 1) // new in v.102 , 4 = deep stops
|
|
|
1721 temp_limit = temp_pres_gtissue_limit_GF_low;
|
|
|
1722 else
|
|
|
1723 temp_limit = temp_pres_gtissue_limit;
|
|
|
1724 if ((temp_deco > temp_limit) && (temp_deco > pres_surface)) // changes in v.102
|
|
|
1725 {
|
|
|
1726 lock_GF_depth_list = 0; // new in v.102, distance to first stop > 10 mtr.
|
|
|
1727 output[6] = 0;
|
|
|
1728 if (char_I_const_ppO2 == 0) // new in v.101 // calculate at half of the ascent
|
|
|
1729 {
|
|
|
1730 deco_diluent = temp_deco + 0.5; // new in v.101
|
|
|
1731 if (temp_deco + 0.5 > deco_gas_change)
|
|
|
1732 {
|
|
|
1733 calc_N2_ratio = N2_ratio;
|
|
|
1734 calc_He_ratio = He_ratio;
|
|
|
1735 }
|
|
|
1736 else
|
|
|
1737 {
|
|
|
1738 calc_N2_ratio = deco_N2_ratio;
|
|
|
1739 calc_He_ratio = deco_He_ratio;
|
|
|
1740 }
|
|
|
1741 }
|
|
|
1742 else // new in v.101
|
|
|
1743 {
|
|
|
1744 calc_N2_ratio = N2_ratio;
|
|
|
1745 calc_He_ratio = He_ratio;
|
|
|
1746 if (temp_deco + 0.5 > deco_ppO2_change)
|
|
|
1747 deco_diluent = ((temp_deco + 0.5 - const_ppO2)/(N2_ratio + He_ratio)); // new in v.101 // calculate at half of the ascent
|
|
|
1748 else
|
|
|
1749 deco_diluent = ((temp_deco + 0.5 - deco_ppO2)/(N2_ratio + He_ratio)); // new in v.101 // calculate at half of the ascent
|
|
|
1750 if (deco_diluent > (temp_deco +0.5)) // new in v.101
|
|
|
1751 deco_diluent = temp_deco + 0.5; // new in v.101 // calculate at half of the ascent
|
|
|
1752 }
|
|
|
1753 if (deco_diluent > 0.0627) // new in v.101
|
|
|
1754 {
|
|
|
1755 temp_atem = calc_N2_ratio * (deco_diluent - 0.0627); // changed in v.101
|
|
|
1756 temp2_atem = calc_He_ratio * (deco_diluent - 0.0627); // changed in v.101
|
|
|
1757 }
|
|
|
1758 else // new in v.101
|
|
|
1759 {
|
|
|
1760 temp_atem = 0.0; // new in v.101
|
|
|
1761 temp2_atem = 0.0; // new in v.101
|
|
|
1762 }
|
|
|
1763 sim_tissue_1min();
|
|
|
1764 char_O_deco_status = char_O_deco_status + 1;
|
|
|
1765 if (char_O_deco_status < 16) // 16 is the limit of calculations for one time slot
|
|
|
1766 int_temp_decostatus = 1;
|
|
|
1767 }
|
|
|
1768 } while (int_temp_decostatus == 1);
|
|
|
1769 } // calc_hauptroutine_calc_ascend_to_deco
|
|
|
1770
|
|
|
1771 // --------------
|
|
|
1772 // calc_tissue //
|
|
|
1773 // --------------
|
|
|
1774 // optimized in v.101
|
|
|
1775
|
|
|
1776 void calc_tissue(void)
|
|
|
1777 {
|
|
|
1778 _asm
|
|
|
1779 lfsr 1, 0x300
|
|
|
1780 movlw 0x01
|
|
|
1781 movwf TBLPTRU,0
|
|
|
1782 _endasm
|
|
|
1783
|
|
|
1784 char_O_gtissue_no = 255;
|
|
|
1785 pres_gtissue_limit = 0.0;
|
|
|
1786
|
|
|
1787 for (ci=0;ci<16;ci++)
|
|
|
1788 {
|
|
|
1789 _asm
|
|
|
1790 movlw 0x02
|
|
|
1791 movwf TBLPTRH,0
|
|
|
1792 movlb 4 // fuer ci
|
|
|
1793 movf ci,0,1
|
|
|
1794 addwf ci,0,1
|
|
|
1795 addwf ci,0,1
|
|
|
1796 addwf ci,0,1
|
|
|
1797 movwf TBLPTRL,0
|
|
|
1798 TBLRDPOSTINC
|
|
|
1799 movff TABLAT,var_e2secs+1 // the order is confussing
|
|
|
1800 TBLRDPOSTINC
|
|
|
1801 movff TABLAT,var_e2secs // low byte first, high afterwards
|
|
|
1802 TBLRDPOSTINC
|
|
|
1803 movff TABLAT,var_e2secs+3
|
|
|
1804 TBLRD
|
|
|
1805 movff TABLAT,var_e2secs+2
|
|
|
1806 addlw 0x40
|
|
|
1807 movwf TBLPTRL,0
|
|
|
1808 TBLRDPOSTINC
|
|
|
1809 movff TABLAT,var2_e2secs+1
|
|
|
1810 TBLRDPOSTINC
|
|
|
1811 movff TABLAT,var2_e2secs
|
|
|
1812 TBLRDPOSTINC
|
|
|
1813 movff TABLAT,var2_e2secs+3
|
|
|
1814 TBLRD
|
|
|
1815 movff TABLAT,var2_e2secs+2
|
|
|
1816 addlw 0x40
|
|
|
1817 movwf TBLPTRL,0
|
|
|
1818 TBLRDPOSTINC
|
|
|
1819 movff TABLAT,var_a+1
|
|
|
1820 TBLRDPOSTINC
|
|
|
1821 movff TABLAT,var_a
|
|
|
1822 TBLRDPOSTINC
|
|
|
1823 movff TABLAT,var_a+3
|
|
|
1824 TBLRD
|
|
|
1825 movff TABLAT,var_a+2
|
|
|
1826 addlw 0x40
|
|
|
1827 movwf TBLPTRL,0
|
|
|
1828 TBLRDPOSTINC
|
|
|
1829 movff TABLAT,var2_a+1
|
|
|
1830 TBLRDPOSTINC
|
|
|
1831 movff TABLAT,var2_a
|
|
|
1832 TBLRDPOSTINC
|
|
|
1833 movff TABLAT,var2_a+3
|
|
|
1834 TBLRD
|
|
|
1835 movff TABLAT,var2_a+2
|
|
|
1836 addlw 0x40
|
|
|
1837 movwf TBLPTRL,0
|
|
|
1838 incf TBLPTRH,1,0
|
|
|
1839 TBLRDPOSTINC
|
|
|
1840 movff TABLAT,var_b+1
|
|
|
1841 TBLRDPOSTINC
|
|
|
1842 movff TABLAT,var_b
|
|
|
1843 TBLRDPOSTINC
|
|
|
1844 movff TABLAT,var_b+3
|
|
|
1845 TBLRD
|
|
|
1846 movff TABLAT,var_b+2
|
|
|
1847 addlw 0x40
|
|
|
1848 movwf TBLPTRL,0
|
|
|
1849 TBLRDPOSTINC
|
|
|
1850 movff TABLAT,var2_b+1
|
|
|
1851 TBLRDPOSTINC
|
|
|
1852 movff TABLAT,var2_b
|
|
|
1853 TBLRDPOSTINC
|
|
|
1854 movff TABLAT,var2_b+3
|
|
|
1855 TBLRD
|
|
|
1856 movff TABLAT,var2_b+2
|
|
|
1857 _endasm
|
|
|
1858 // the start values are the previous end values // write new values in temp
|
|
|
1859
|
|
|
1860 if( (var_e2secs < 0.0000363)
|
|
|
1861 || (var_e2secs > 0.00577)
|
|
|
1862 || (var2_e2secs < 0.0000961)
|
|
|
1863 || (var2_e2secs > 0.150)
|
|
|
1864 || (var_a < 0.231)
|
|
|
1865 || (var_a > 1.27)
|
|
|
1866 || (var_b < 0.504)
|
|
|
1867 || (var_b > 0.966)
|
|
|
1868 || (var2_a < 0.510)
|
|
|
1869 || (var2_a > 1.75)
|
|
|
1870 || (var2_b < 0.423)
|
|
|
1871 || (var2_b > 0.927)
|
|
|
1872 )
|
|
|
1873 int_O_DBG_pre_bitfield |= DBG_ZH16ERR;
|
|
|
1874
|
|
|
1875 // N2
|
|
|
1876 temp_tissue = (temp_atem - pres_tissue[ci]) * var_e2secs;
|
|
|
1877 temp_tissue_safety();
|
|
|
1878 pres_tissue[ci] = pres_tissue[ci] + temp_tissue;
|
|
|
1879
|
|
|
1880 // He
|
|
|
1881 temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var2_e2secs;
|
|
|
1882 temp_tissue_safety();
|
|
|
1883 pres_tissue[ci+16] = pres_tissue[ci+16] + temp_tissue;
|
|
|
1884
|
|
|
1885 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16];
|
|
|
1886
|
|
|
1887 var_a = (var_a * pres_tissue[ci] + var2_a * pres_tissue[ci+16]) / temp_tissue;
|
|
|
1888 var_b = (var_b * pres_tissue[ci] + var2_b * pres_tissue[ci+16]) / temp_tissue;
|
|
|
1889 pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
|
|
|
1890 if (pres_tissue_limit[ci] < 0)
|
|
|
1891 pres_tissue_limit[ci] = 0;
|
|
|
1892 if (pres_tissue_limit[ci] > pres_gtissue_limit)
|
|
|
1893 {
|
|
|
1894 pres_gtissue_limit = pres_tissue_limit[ci];
|
|
|
1895 char_O_gtissue_no = ci;
|
|
|
1896 }//if
|
|
|
1897 } // for
|
|
|
1898 }//calc_tissue(void)
|
|
|
1899
|
|
|
1900 // ----------------
|
|
|
1901 // calc_nullzeit //
|
|
|
1902 // ----------------
|
|
|
1903 // calculates the remaining bottom time
|
|
|
1904
|
|
|
1905 // unchanged in v.101
|
|
|
1906
|
|
|
1907 void calc_nullzeit(void)
|
|
|
1908 {
|
|
|
1909 char_O_nullzeit = 0;
|
|
|
1910 int_temp = 1;
|
|
|
1911 do
|
|
|
1912 {
|
|
|
1913 backup_sim_pres_tissue();
|
|
|
1914 sim_tissue_10min();
|
|
|
1915 char_O_nullzeit = char_O_nullzeit + 10;
|
|
|
1916 int_temp = int_temp + 1;
|
|
|
1917 if (char_I_deco_model == 1)
|
|
|
1918 temp1 = GF_high * temp_pres_gtissue_diff + temp_pres_gtissue;
|
|
|
1919 else
|
|
|
1920 temp1 = temp_pres_gtissue_limit;
|
|
|
1921 if (temp1 > temp_surface) // changed in v.102 , if guiding tissue can not be exposed to surface pressure immediately
|
|
|
1922 int_temp = 255;
|
|
|
1923 } while (int_temp < 17);
|
|
|
1924 if (int_temp == 255)
|
|
|
1925 {
|
|
|
1926 restore_sim_pres_tissue();
|
|
|
1927 char_O_nullzeit = char_O_nullzeit - 10;
|
|
|
1928 } //if int_temp == 255]
|
|
|
1929 int_temp = 1;
|
|
|
1930 if (char_O_nullzeit < 60)
|
|
|
1931 {
|
|
|
1932 do
|
|
|
1933 {
|
|
|
1934 sim_tissue_1min();
|
|
|
1935 char_O_nullzeit = char_O_nullzeit + 1;
|
|
|
1936 int_temp = int_temp + 1; // new in v.102a
|
|
|
1937 if (char_I_deco_model == 1)
|
|
|
1938 temp1 = GF_high * temp_pres_gtissue_diff + temp_pres_gtissue;
|
|
|
1939 else
|
|
|
1940 temp1 = temp_pres_gtissue_limit;
|
|
|
1941 if (temp1 > temp_surface) // changed in v.102 , if guiding tissue can not be exposed to surface pressure immediately
|
|
|
1942 int_temp = 255;
|
|
|
1943 } while (int_temp < 10);
|
|
|
1944 if (int_temp == 255)
|
|
|
1945 char_O_nullzeit = char_O_nullzeit - 1;
|
|
|
1946 } // if char_O_nullzeit < 60
|
|
|
1947 } //calc_nullzeit
|
|
|
1948
|
|
|
1949 // -------------------------
|
|
|
1950 // backup_sim_pres_tissue //
|
|
|
1951 // -------------------------
|
|
|
1952 void backup_sim_pres_tissue(void)
|
|
|
1953 {
|
|
|
1954 for (x = 0;x<16;x++)
|
|
|
1955 {
|
|
|
1956 sim_pres_tissue_backup[x] = sim_pres_tissue[x];
|
|
|
1957 sim_pres_tissue_backup[x+16] = sim_pres_tissue[x+16];
|
|
|
1958 }
|
|
|
1959 } // backup_sim
|
|
|
1960
|
|
|
1961 // --------------------------
|
|
|
1962 // restore_sim_pres_tissue //
|
|
|
1963 // --------------------------
|
|
|
1964 void restore_sim_pres_tissue(void)
|
|
|
1965 {
|
|
|
1966 for (x = 0;x<16;x++)
|
|
|
1967 {
|
|
|
1968 sim_pres_tissue[x] = sim_pres_tissue_backup[x];
|
|
|
1969 sim_pres_tissue[x+16] = sim_pres_tissue_backup[x+16];
|
|
|
1970 }
|
|
|
1971 } // restore_sim
|
|
|
1972
|
|
|
1973 // ------------------
|
|
|
1974 // calc_ascenttime //
|
|
|
1975 // ------------------
|
|
|
1976
|
|
|
1977 void calc_ascenttime(void)
|
|
|
1978 {
|
|
|
1979 if (pres_respiration > pres_surface)
|
|
|
1980 {
|
|
|
1981 switch (char_O_deco_status)
|
|
|
1982 {
|
|
|
1983 case 2:
|
|
|
1984 char_O_ascenttime = 255;
|
|
|
1985 break;
|
|
|
1986 case 1:
|
|
|
1987 break;
|
|
|
1988 default:
|
|
|
1989 temp1 = pres_respiration - pres_surface + 0.6; // + 0.6 hence 1 minute ascent time from a depth of 4 meter on
|
|
|
1990 if (temp1 < 0)
|
|
|
1991 temp1 = 0;
|
|
|
1992 if (temp1 > 255)
|
|
|
1993 temp1 = 255;
|
|
|
1994 char_O_ascenttime = (char)temp1;
|
|
|
1995
|
|
|
1996 for(ci=0;ci<7;ci++)
|
|
|
1997 {
|
|
|
1998 x = char_O_ascenttime + char_O_array_decotime[ci];
|
|
|
1999 if (x < char_O_ascenttime)
|
|
|
2000 char_O_ascenttime = 255;
|
|
|
2001 else
|
|
|
2002 char_O_ascenttime = x;
|
|
|
2003 }
|
|
|
2004 }
|
|
|
2005 }
|
|
|
2006 else
|
|
|
2007 char_O_ascenttime = 0;
|
|
|
2008 } // calc_ascenttime()
|
|
|
2009
|
|
|
2010
|
|
|
2011 // ---------------------
|
|
|
2012 // update_startvalues //
|
|
|
2013 // ---------------------
|
|
|
2014 // updated in v.102
|
|
|
2015
|
|
|
2016 void update_startvalues(void)
|
|
|
2017 {
|
|
|
2018 temp_pres_gtissue_limit = pres_gtissue_limit;
|
|
|
2019 temp_pres_gtissue = pres_tissue[char_O_gtissue_no] + pres_tissue[char_O_gtissue_no+16];
|
|
|
2020 temp_pres_gtissue_diff = temp_pres_gtissue_limit - temp_pres_gtissue; // negative number
|
|
|
2021 temp_pres_gtissue_limit_GF_low = GF_low * temp_pres_gtissue_diff + temp_pres_gtissue;
|
|
|
2022 temp_pres_gtissue_limit_GF_low_below_surface = temp_pres_gtissue_limit_GF_low - pres_surface;
|
|
|
2023 if (temp_pres_gtissue_limit_GF_low_below_surface < 0)
|
|
|
2024 temp_pres_gtissue_limit_GF_low_below_surface = 0;
|
|
|
2025
|
|
|
2026 temp_gtissue_no = char_O_gtissue_no;
|
|
|
2027 for (x = 0;x<16;x++)
|
|
|
2028 {
|
|
|
2029 sim_pres_tissue[x] = pres_tissue[x];
|
|
|
2030 sim_pres_tissue[x+16] = pres_tissue[x+16];
|
|
|
2031 sim_pres_tissue_limit[x] = pres_tissue_limit[x];
|
|
|
2032 }
|
|
|
2033 } // update_startvalues
|
|
|
2034
|
|
|
2035
|
|
|
2036 // ------------------
|
|
|
2037 // sim_tissue_1min //
|
|
|
2038 // ------------------
|
|
|
2039 // optimized in v.101
|
|
|
2040
|
|
|
2041 void sim_tissue_1min(void)
|
|
|
2042 {
|
|
|
2043 temp_pres_gtissue_limit = 0.0;
|
|
|
2044 temp_gtissue_no = 255;
|
|
|
2045
|
|
|
2046 _asm
|
|
|
2047 lfsr 1, 0x300
|
|
|
2048 movlw 0x01
|
|
|
2049 movwf TBLPTRU,0
|
|
|
2050 _endasm
|
|
|
2051
|
|
|
2052
|
|
|
2053 for (ci=0;ci<16;ci++)
|
|
|
2054 {
|
|
|
2055 _asm
|
|
|
2056 movlw 0x02
|
|
|
2057 movwf TBLPTRH,0
|
|
|
2058 movlb 4 // fuer ci
|
|
|
2059 movf ci,0,1
|
|
|
2060 addwf ci,0,1
|
|
|
2061 addwf ci,0,1
|
|
|
2062 addwf ci,0,1
|
|
|
2063 addlw 0x80
|
|
|
2064 movwf TBLPTRL,0
|
|
|
2065 TBLRDPOSTINC
|
|
|
2066 movff TABLAT,var_a+1
|
|
|
2067 TBLRDPOSTINC
|
|
|
2068 movff TABLAT,var_a
|
|
|
2069 TBLRDPOSTINC
|
|
|
2070 movff TABLAT,var_a+3
|
|
|
2071 TBLRD
|
|
|
2072 movff TABLAT,var_a+2
|
|
|
2073 addlw 0x40
|
|
|
2074 movwf TBLPTRL,0
|
|
|
2075 TBLRDPOSTINC
|
|
|
2076 movff TABLAT,var2_a+1
|
|
|
2077 TBLRDPOSTINC
|
|
|
2078 movff TABLAT,var2_a
|
|
|
2079 TBLRDPOSTINC
|
|
|
2080 movff TABLAT,var2_a+3
|
|
|
2081 TBLRD
|
|
|
2082 movff TABLAT,var2_a+2
|
|
|
2083 addlw 0x40
|
|
|
2084 movwf TBLPTRL,0
|
|
|
2085 incf TBLPTRH,1,0
|
|
|
2086 TBLRDPOSTINC
|
|
|
2087 movff TABLAT,var_b+1
|
|
|
2088 TBLRDPOSTINC
|
|
|
2089 movff TABLAT,var_b
|
|
|
2090 TBLRDPOSTINC
|
|
|
2091 movff TABLAT,var_b+3
|
|
|
2092 TBLRD
|
|
|
2093 movff TABLAT,var_b+2
|
|
|
2094 addlw 0x40
|
|
|
2095 movwf TBLPTRL,0
|
|
|
2096 TBLRDPOSTINC
|
|
|
2097 movff TABLAT,var2_b+1
|
|
|
2098 TBLRDPOSTINC
|
|
|
2099 movff TABLAT,var2_b
|
|
|
2100 TBLRDPOSTINC
|
|
|
2101 movff TABLAT,var2_b+3
|
|
|
2102 TBLRD
|
|
|
2103 movff TABLAT,var2_b+2
|
|
|
2104 addlw 0xC0
|
|
|
2105 movwf TBLPTRL,0
|
|
|
2106 incf TBLPTRH,1,0
|
|
|
2107 TBLRDPOSTINC
|
|
|
2108 movff TABLAT,var_e1min+1
|
|
|
2109 TBLRDPOSTINC
|
|
|
2110 movff TABLAT,var_e1min
|
|
|
2111 TBLRDPOSTINC
|
|
|
2112 movff TABLAT,var_e1min+3
|
|
|
2113 TBLRD
|
|
|
2114 movff TABLAT,var_e1min+2
|
|
|
2115 addlw 0x40
|
|
|
2116 movwf TBLPTRL,0
|
|
|
2117 TBLRDPOSTINC
|
|
|
2118 movff TABLAT,var2_e1min+1
|
|
|
2119 TBLRDPOSTINC
|
|
|
2120 movff TABLAT,var2_e1min
|
|
|
2121 TBLRDPOSTINC
|
|
|
2122 movff TABLAT,var2_e1min+3
|
|
|
2123 TBLRD
|
|
|
2124 movff TABLAT,var2_e1min+2
|
|
|
2125 _endasm
|
|
|
2126 // N2
|
|
|
2127 temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_e1min;
|
|
|
2128 temp_tissue_safety();
|
|
|
2129 sim_pres_tissue[ci] = sim_pres_tissue[ci] + temp_tissue;
|
|
|
2130 // He
|
|
|
2131 temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var2_e1min;
|
|
|
2132 temp_tissue_safety();
|
|
|
2133 sim_pres_tissue[ci+16] = sim_pres_tissue[ci+16] + temp_tissue;
|
|
|
2134 // pressure limit
|
|
|
2135 temp_tissue = sim_pres_tissue[ci] + sim_pres_tissue[ci+16];
|
|
|
2136 var_a = (var_a * sim_pres_tissue[ci] + var2_a * sim_pres_tissue[ci+16]) / temp_tissue;
|
|
|
2137 var_b = (var_b * sim_pres_tissue[ci] + var2_b * sim_pres_tissue[ci+16]) / temp_tissue;
|
|
|
2138 sim_pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
|
|
|
2139
|
|
|
2140 if (sim_pres_tissue_limit[ci] < 0)
|
|
|
2141 sim_pres_tissue_limit[ci] = 0;
|
|
|
2142 if (sim_pres_tissue_limit[ci] > temp_pres_gtissue_limit)
|
|
|
2143 {
|
|
|
2144 temp_pres_gtissue = temp_tissue;
|
|
|
2145 temp_pres_gtissue_limit = sim_pres_tissue_limit[ci];
|
|
|
2146 temp_gtissue_no = ci;
|
|
|
2147 }
|
|
|
2148 } // for
|
|
|
2149 temp_pres_gtissue_diff = temp_pres_gtissue_limit - temp_pres_gtissue;
|
|
|
2150 temp_pres_gtissue_limit_GF_low = GF_low * temp_pres_gtissue_diff + temp_pres_gtissue;
|
|
|
2151 temp_pres_gtissue_limit_GF_low_below_surface = temp_pres_gtissue_limit_GF_low - pres_surface;
|
|
|
2152 if (temp_pres_gtissue_limit_GF_low_below_surface < 0)
|
|
|
2153 temp_pres_gtissue_limit_GF_low_below_surface = 0;
|
|
|
2154 } //sim_tissue_1min()
|
|
|
2155
|
|
|
2156 //--------------------
|
|
|
2157 // sim_tissue_10min //
|
|
|
2158 //--------------------
|
|
|
2159
|
|
|
2160 // Attention!! uses var_e1min und var2_e1min to load 10min data !!!
|
|
|
2161 // is identical to sim_tissue_1min routine except for the different load of those variables
|
|
|
2162
|
|
|
2163 // optimized in v.101
|
|
|
2164
|
|
|
2165 void sim_tissue_10min(void)
|
|
|
2166 {
|
|
|
2167 temp_pres_gtissue_limit = 0.0;
|
|
|
2168 temp_gtissue_no = 255;
|
|
|
2169
|
|
|
2170 _asm
|
|
|
2171 lfsr 1, 0x300
|
|
|
2172 movlw 0x01
|
|
|
2173 movwf TBLPTRU,0
|
|
|
2174 _endasm
|
|
|
2175
|
|
|
2176 for (ci=0;ci<16;ci++)
|
|
|
2177 {
|
|
|
2178 _asm
|
|
|
2179 movlw 0x02
|
|
|
2180 movwf TBLPTRH,0
|
|
|
2181 movlb 4 // fuer ci
|
|
|
2182 movf ci,0,1
|
|
|
2183 addwf ci,0,1
|
|
|
2184 addwf ci,0,1
|
|
|
2185 addwf ci,0,1
|
|
|
2186 addlw 0x80
|
|
|
2187 movwf TBLPTRL,0
|
|
|
2188 TBLRDPOSTINC
|
|
|
2189 movff TABLAT,var_a+1
|
|
|
2190 TBLRDPOSTINC
|
|
|
2191 movff TABLAT,var_a
|
|
|
2192 TBLRDPOSTINC
|
|
|
2193 movff TABLAT,var_a+3
|
|
|
2194 TBLRD
|
|
|
2195 movff TABLAT,var_a+2
|
|
|
2196 addlw 0x40
|
|
|
2197 movwf TBLPTRL,0
|
|
|
2198 TBLRDPOSTINC
|
|
|
2199 movff TABLAT,var2_a+1
|
|
|
2200 TBLRDPOSTINC
|
|
|
2201 movff TABLAT,var2_a
|
|
|
2202 TBLRDPOSTINC
|
|
|
2203 movff TABLAT,var2_a+3
|
|
|
2204 TBLRD
|
|
|
2205 movff TABLAT,var2_a+2
|
|
|
2206 addlw 0x40
|
|
|
2207 movwf TBLPTRL,0
|
|
|
2208 incf TBLPTRH,1,0
|
|
|
2209 TBLRDPOSTINC
|
|
|
2210 movff TABLAT,var_b+1
|
|
|
2211 TBLRDPOSTINC
|
|
|
2212 movff TABLAT,var_b
|
|
|
2213 TBLRDPOSTINC
|
|
|
2214 movff TABLAT,var_b+3
|
|
|
2215 TBLRD
|
|
|
2216 movff TABLAT,var_b+2
|
|
|
2217 addlw 0x40
|
|
|
2218 movwf TBLPTRL,0
|
|
|
2219 TBLRDPOSTINC
|
|
|
2220 movff TABLAT,var2_b+1
|
|
|
2221 TBLRDPOSTINC
|
|
|
2222 movff TABLAT,var2_b
|
|
|
2223 TBLRDPOSTINC
|
|
|
2224 movff TABLAT,var2_b+3
|
|
|
2225 TBLRD
|
|
|
2226 movff TABLAT,var2_b+2
|
|
|
2227 addlw 0xC0 // different to 1 min
|
|
|
2228 movwf TBLPTRL,0
|
|
|
2229 incf TBLPTRH,1,0
|
|
|
2230 incf TBLPTRH,1,0 // different to 1 min
|
|
|
2231 TBLRDPOSTINC
|
|
|
2232 movff TABLAT,var_e1min+1
|
|
|
2233 TBLRDPOSTINC
|
|
|
2234 movff TABLAT,var_e1min
|
|
|
2235 TBLRDPOSTINC
|
|
|
2236 movff TABLAT,var_e1min+3
|
|
|
2237 TBLRD
|
|
|
2238 movff TABLAT,var_e1min+2
|
|
|
2239 addlw 0x40
|
|
|
2240 movwf TBLPTRL,0
|
|
|
2241 //incf TBLPTRH,1,0 // different to 1 min
|
|
|
2242 TBLRDPOSTINC
|
|
|
2243 movff TABLAT,var2_e1min+1
|
|
|
2244 TBLRDPOSTINC
|
|
|
2245 movff TABLAT,var2_e1min
|
|
|
2246 TBLRDPOSTINC
|
|
|
2247 movff TABLAT,var2_e1min+3
|
|
|
2248 TBLRD
|
|
|
2249 movff TABLAT,var2_e1min+2
|
|
|
2250 _endasm
|
|
|
2251 // N2
|
|
|
2252 temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_e1min;
|
|
|
2253 temp_tissue_safety();
|
|
|
2254 sim_pres_tissue[ci] = sim_pres_tissue[ci] + temp_tissue;
|
|
|
2255 // He
|
|
|
2256 temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var2_e1min;
|
|
|
2257 temp_tissue_safety();
|
|
|
2258 sim_pres_tissue[ci+16] = sim_pres_tissue[ci+16] + temp_tissue;
|
|
|
2259 // pressure limit
|
|
|
2260 temp_tissue = sim_pres_tissue[ci] + sim_pres_tissue[ci+16];
|
|
|
2261 var_a = (var_a * sim_pres_tissue[ci] + var2_a * sim_pres_tissue[ci+16]) / temp_tissue;
|
|
|
2262 var_b = (var_b * sim_pres_tissue[ci] + var2_b * sim_pres_tissue[ci+16]) / temp_tissue;
|
|
|
2263
|
|
|
2264 sim_pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
|
|
|
2265 if (sim_pres_tissue_limit[ci] < 0)
|
|
|
2266 sim_pres_tissue_limit[ci] = 0;
|
|
|
2267 if (sim_pres_tissue_limit[ci] > temp_pres_gtissue_limit)
|
|
|
2268 {
|
|
|
2269 temp_pres_gtissue = temp_tissue;
|
|
|
2270 temp_pres_gtissue_limit = sim_pres_tissue_limit[ci];
|
|
|
2271 temp_gtissue_no = ci;
|
|
|
2272 }
|
|
|
2273 } // for
|
|
|
2274 temp_pres_gtissue_diff = temp_pres_gtissue_limit - temp_pres_gtissue; // negative number
|
|
|
2275 temp_pres_gtissue_limit_GF_low = GF_low * temp_pres_gtissue_diff + temp_pres_gtissue;
|
|
|
2276 temp_pres_gtissue_limit_GF_low_below_surface = temp_pres_gtissue_limit_GF_low - pres_surface;
|
|
|
2277 if (temp_pres_gtissue_limit_GF_low_below_surface < 0)
|
|
|
2278 temp_pres_gtissue_limit_GF_low_below_surface = 0;
|
|
|
2279 } //sim_tissue_10min()
|
|
|
2280
|
|
|
2281
|
|
|
2282 // ------------------
|
|
|
2283 // clear_decoarray //
|
|
|
2284 // ------------------
|
|
|
2285 // unchanged in v.101
|
|
|
2286
|
|
|
2287 void clear_decoarray(void)
|
|
|
2288 {
|
|
|
2289 char_O_array_decodepth[0] = 0;
|
|
|
2290 char_O_array_decodepth[1] = 0;
|
|
|
2291 char_O_array_decodepth[2] = 0;
|
|
|
2292 char_O_array_decodepth[3] = 0;
|
|
|
2293 char_O_array_decodepth[4] = 0;
|
|
|
2294 char_O_array_decodepth[5] = 0;
|
|
|
2295 char_O_array_decotime[0] = 0;
|
|
|
2296 char_O_array_decotime[1] = 0;
|
|
|
2297 char_O_array_decotime[2] = 0;
|
|
|
2298 char_O_array_decotime[3] = 0;
|
|
|
2299 char_O_array_decotime[4] = 0;
|
|
|
2300 char_O_array_decotime[5] = 0;
|
|
|
2301 char_O_array_decotime[6] = 0;
|
|
|
2302 } // clear_decoarray
|
|
|
2303
|
|
|
2304
|
|
|
2305 // -------------------
|
|
|
2306 // update_decoarray //
|
|
|
2307 // -------------------
|
|
|
2308 // unchanged in v.101
|
|
|
2309
|
|
|
2310 void update_decoarray()
|
|
|
2311 {
|
|
|
2312 x = 0;
|
|
|
2313 do
|
|
|
2314 {
|
|
|
2315 if (char_O_array_decodepth[x] == temp_depth_limit)
|
|
|
2316 {
|
|
|
2317 int_temp = char_O_array_decotime[x] + temp_decotime;
|
|
|
2318 if (int_temp < 0)
|
|
|
2319 int_temp = 0;
|
|
|
2320 if (int_temp > 240)
|
|
|
2321 int_temp = 240;
|
|
|
2322 char_O_array_decotime[x] = int_temp;
|
|
|
2323 x = 10; // exit
|
|
|
2324 } // if
|
|
|
2325 else
|
|
|
2326 {
|
|
|
2327 if (char_O_array_decodepth[x] == 0)
|
|
|
2328 {
|
|
|
2329 if (temp_depth_limit > 255)
|
|
|
2330 char_O_array_decodepth[x] = 255;
|
|
|
2331 else
|
|
|
2332 char_O_array_decodepth[x] = (char)temp_depth_limit;
|
|
|
2333 int_temp = char_O_array_decotime[x] + temp_decotime;
|
|
|
2334 if (int_temp > 240)
|
|
|
2335 char_O_array_decotime[x] = 240;
|
|
|
2336 else
|
|
|
2337 char_O_array_decotime[x] = (char)int_temp;
|
|
|
2338 x = 10; // exit
|
|
|
2339 } // if
|
|
|
2340 else
|
|
|
2341 x++;
|
|
|
2342 } // else
|
|
|
2343 } while (x<6);
|
|
|
2344 if (x == 6)
|
|
|
2345 {
|
|
|
2346 int_temp = char_O_array_decotime[6] + temp_decotime;
|
|
|
2347 if (int_temp > 220)
|
|
|
2348 char_O_array_decotime[6] = 220;
|
|
|
2349 else
|
|
|
2350 char_O_array_decotime[6] = (char)int_temp;
|
|
|
2351 } // if x == 6
|
|
|
2352 } // update_decoarray
|
|
|
2353
|
|
|
2354
|
|
|
2355 // -----------------------
|
|
|
2356 // calc_gradient_factor //
|
|
|
2357 // -----------------------
|
|
|
2358 // optimized in v.101 (var_a)
|
|
|
2359 // new code in v.102
|
|
|
2360
|
|
|
2361 void calc_gradient_factor(void)
|
|
|
2362 {
|
|
|
2363 // tissue > respiration (entsaettigungsvorgang)
|
|
|
2364 // gradient ist wieviel prozent an limit mit basis tissue
|
|
|
2365 // dh. 0% = respiration == tissue
|
|
|
2366 // dh. 100% = respiration == limit
|
|
|
2367 temp_tissue = pres_tissue[char_O_gtissue_no] + pres_tissue[char_O_gtissue_no+16];
|
|
|
2368 temp1 = temp_tissue - pres_respiration;
|
|
|
2369 temp2 = temp_tissue - pres_tissue_limit[char_O_gtissue_no]; // changed in v.102
|
|
|
2370 temp2 = temp1/temp2;
|
|
|
2371 temp2 = temp2 * 100; // displayed in percent
|
|
|
2372 if (temp2 < 0)
|
|
|
2373 temp2 = 0;
|
|
|
2374 if (temp2 > 255)
|
|
|
2375 temp2 = 255;
|
|
|
2376 if (temp1 < 0)
|
|
|
2377 char_O_gradient_factor = 0;
|
|
|
2378 else
|
|
|
2379 char_O_gradient_factor = (char)temp2;
|
|
|
2380
|
|
|
2381 temp3 = temp2;
|
|
|
2382
|
|
|
2383 if (char_I_deco_model == 1) // calculate relative gradient factor
|
|
|
2384 {
|
|
|
2385 temp1 = (float)temp_depth_GF_low_meter * 0.09995;
|
|
|
2386 temp2 = pres_respiration - pres_surface;
|
|
|
2387 if (temp2 <= 0)
|
|
|
2388 temp1 = GF_high;
|
|
|
2389 else
|
|
|
2390 if (temp2 >= temp1)
|
|
|
2391 temp1 = GF_low;
|
|
|
2392 else
|
|
|
2393 temp1 = GF_low + (temp1 - temp2)/temp1*GF_delta;
|
|
|
2394 if (temp_depth_GF_low_meter == 0)
|
|
|
2395 temp1 = GF_high;
|
|
|
2396 temp2 = temp3 / temp1; // temp3 is already in percent
|
|
|
2397 if (temp2 < 0)
|
|
|
2398 temp2 = 0;
|
|
|
2399 if (temp2 > 255)
|
|
|
2400 temp2 = 255;
|
|
|
2401 char_O_relative_gradient_GF = (char)temp2;
|
|
|
2402 } // calc relative gradient factor
|
|
|
2403 else
|
|
|
2404 {
|
|
|
2405 char_O_relative_gradient_GF = char_O_gradient_factor;
|
|
|
2406 }
|
|
|
2407 } // calc_gradient
|
|
|
2408
|
|
|
2409 // ---------------------------
|
|
|
2410 // calc_gradient_array_only //
|
|
|
2411 // ---------------------------
|
|
|
2412 // optimized in v.101 (var_a)
|
|
|
2413 // new code in v.102
|
|
|
2414
|
|
|
2415 void calc_gradient_array_only()
|
|
|
2416 {
|
|
|
2417 pres_respiration = (float)int_I_pres_respiration / 1000.0; // assembler code uses different digit system
|
|
|
2418 for (ci=0;ci<16;ci++)
|
|
|
2419 {
|
|
|
2420 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16];
|
|
|
2421 temp1 = temp_tissue - pres_respiration;
|
|
|
2422 temp2 = temp_tissue - pres_tissue_limit[ci];
|
|
|
2423 temp2 = temp1/temp2;
|
|
|
2424 temp2 = temp2 * 200; // because of output in (Double-)percentage
|
|
|
2425 if (temp2 < 0)
|
|
|
2426 temp2 = 0;
|
|
|
2427 if (temp2 > 255)
|
|
|
2428 temp2 = 255;
|
|
|
2429 if (temp1 < 0)
|
|
|
2430 char_O_array_gradient_weighted[ci] = 0;
|
|
|
2431 else
|
|
|
2432 char_O_array_gradient_weighted[ci] = (char)temp2;
|
|
|
2433 } // for
|
|
|
2434 } // calc_gradient_array_only
|
|
|
2435
|
|
|
2436
|
|
|
2437 // -------------------------
|
|
|
2438 // calc_desaturation_time //
|
|
|
2439 // -------------------------
|
|
|
2440 // FIXED N2_ratio
|
|
|
2441 // unchanged in v.101
|
|
|
2442
|
|
|
2443 void calc_desaturation_time(void)
|
|
|
2444 {
|
|
|
2445 _asm
|
|
|
2446 lfsr 1, 0x300
|
|
|
2447 movlw 0x01
|
|
|
2448 movwf TBLPTRU,0
|
|
|
2449 _endasm
|
|
|
2450 N2_ratio = 0.7902; // FIXED sum as stated in b"uhlmann
|
|
|
2451 pres_surface = (float)int_I_pres_surface / 1000.0;
|
|
|
2452 temp_atem = N2_ratio * (pres_surface - 0.0627);
|
|
|
2453 int_O_desaturation_time = 0;
|
|
|
2454 float_desaturation_multiplier = char_I_desaturation_multiplier / 142.0; // new in v.101 (70,42%/100.=142)
|
|
|
2455
|
|
|
2456 for (ci=0;ci<16;ci++)
|
|
|
2457 {
|
|
|
2458 _asm
|
|
|
2459 movlw 0x04
|
|
|
2460 movwf TBLPTRH,0
|
|
|
2461 movlb 4 // fuer ci
|
|
|
2462 movf ci,0,1
|
|
|
2463 addwf ci,0,1
|
|
|
2464 addwf ci,0,1
|
|
|
2465 addwf ci,0,1
|
|
|
2466 addlw 0x80
|
|
|
2467 movwf TBLPTRL,0
|
|
|
2468 TBLRDPOSTINC
|
|
|
2469 movff TABLAT,var_halftimes+1
|
|
|
2470 TBLRDPOSTINC
|
|
|
2471 movff TABLAT,var_halftimes
|
|
|
2472 TBLRDPOSTINC
|
|
|
2473 movff TABLAT,var_halftimes+3
|
|
|
2474 TBLRD
|
|
|
2475 movff TABLAT,var_halftimes+2
|
|
|
2476 addlw 0x40
|
|
|
2477 movwf TBLPTRL,0
|
|
|
2478 TBLRDPOSTINC
|
|
|
2479 movff TABLAT,var2_halftimes+1
|
|
|
2480 TBLRDPOSTINC
|
|
|
2481 movff TABLAT,var2_halftimes
|
|
|
2482 TBLRDPOSTINC
|
|
|
2483 movff TABLAT,var2_halftimes+3
|
|
|
2484 TBLRD
|
|
|
2485 movff TABLAT,var2_halftimes+2
|
|
|
2486 _endasm
|
|
|
2487
|
|
|
2488 // saturation_time (for flight) and N2_saturation in multiples of halftime
|
|
|
2489 // version v.100: 1.1 = 10 percent distance to totally clean (totally clean is not possible, would take infinite time )
|
|
|
2490 // new in version v.101: 1.07 = 7 percent distance to totally clean (totally clean is not possible, would take infinite time )
|
|
|
2491 // changes in v.101: 1.05 = 5 percent dist to totally clean is new desaturation point for display and noFly calculations
|
|
|
2492 // N2
|
|
|
2493 temp1 = 1.05 * temp_atem;
|
|
|
2494 temp1 = temp1 - pres_tissue[ci];
|
|
|
2495 temp2 = temp_atem - pres_tissue[ci];
|
|
|
2496 if (temp2 >= 0.0)
|
|
|
2497 {
|
|
|
2498 temp1 = 0;
|
|
|
2499 temp2 = 0;
|
|
|
2500 }
|
|
|
2501 else
|
|
|
2502 temp1 = temp1 / temp2;
|
|
|
2503 if (temp1 > 0.0)
|
|
|
2504 {
|
|
|
2505 temp1 = log(1.0 - temp1);
|
|
|
2506 temp1 = temp1 / -0.6931; // temp1 is the multiples of half times necessary.
|
|
|
2507 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested.
|
|
|
2508 // minus because log is negative
|
|
|
2509 temp2 = var_halftimes * temp1 / float_desaturation_multiplier; // time necessary (in minutes ) for complete desaturation (see comment about 10 percent) , new in v.101: float_desaturation_multiplier
|
|
|
2510 }
|
|
|
2511 else
|
|
|
2512 {
|
|
|
2513 temp1 = 0;
|
|
|
2514 temp2 = 0;
|
|
|
2515 }
|
|
|
2516
|
|
|
2517 // He
|
|
|
2518 temp3 = 0.1 - pres_tissue[ci+16];
|
|
|
2519 if (temp3 >= 0.0)
|
|
|
2520 {
|
|
|
2521 temp3 = 0;
|
|
|
2522 temp4 = 0;
|
|
|
2523 }
|
|
|
2524 else
|
|
|
2525 temp3 = -1.0 * temp3 / pres_tissue[ci+16];
|
|
|
2526 if (temp3 > 0.0)
|
|
|
2527 {
|
|
|
2528 temp3 = log(1.0 - temp3);
|
|
|
2529 temp3 = temp3 / -0.6931; // temp1 is the multiples of half times necessary.
|
|
|
2530 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested.
|
|
|
2531 // minus because log is negative
|
|
|
2532 temp4 = var2_halftimes * temp3 / float_desaturation_multiplier; // time necessary (in minutes ) for "complete" desaturation, new in v.101 float_desaturation_multiplier
|
|
|
2533 }
|
|
|
2534 else
|
|
|
2535 {
|
|
|
2536 temp3 = 0;
|
|
|
2537 temp4 = 0;
|
|
|
2538 }
|
|
|
2539
|
|
|
2540 // saturation_time (for flight)
|
|
|
2541 if (temp4 > temp2)
|
|
|
2542 int_temp = (int)temp4;
|
|
|
2543 else
|
|
|
2544 int_temp = (int)temp2;
|
|
|
2545 if(int_temp > int_O_desaturation_time)
|
|
|
2546 int_O_desaturation_time = int_temp;
|
|
|
2547
|
|
|
2548 // N2 saturation in multiples of halftime for display purposes
|
|
|
2549 temp2 = temp1 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8
|
|
|
2550 temp2 = temp2 + 80.0; // set center
|
|
|
2551 if (temp2 < 0.0)
|
|
|
2552 temp2 = 0.0;
|
|
|
2553 if (temp2 > 255.0)
|
|
|
2554 temp2 = 255.0;
|
|
|
2555 char_O_tissue_saturation[ci] = (char)temp2;
|
|
|
2556 // He saturation in multiples of halftime for display purposes
|
|
|
2557 temp4 = temp3 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8
|
|
|
2558 temp4 = temp4 + 80.0; // set center
|
|
|
2559 if (temp4 < 0.0)
|
|
|
2560 temp4 = 0.0;
|
|
|
2561 if (temp4 > 255.0)
|
|
|
2562 temp4 = 255.0;
|
|
|
2563 char_O_tissue_saturation[ci+16] = (char)temp4;
|
|
|
2564 } // for
|
|
|
2565 } // calc_desaturation_time
|
|
|
2566
|
|
|
2567
|
|
|
2568 // --------------------------
|
|
|
2569 // calc_wo_deco_step_1_min //
|
|
|
2570 // --------------------------
|
|
|
2571 // FIXED N2 Ratio
|
|
|
2572 // optimized in v.101 (...saturation_multiplier)
|
|
|
2573 // desaturation slowed down to 70,42%
|
|
|
2574
|
|
|
2575 void calc_wo_deco_step_1_min(void)
|
|
|
2576 {
|
|
|
2577 if(flag_in_divemode)
|
|
|
2578 {
|
|
|
2579 flag_in_divemode = 0;
|
|
|
2580 set_dbg_end_of_dive();
|
|
|
2581 }
|
|
|
2582 _asm
|
|
|
2583 lfsr 1, 0x300
|
|
|
2584 _endasm
|
|
|
2585 N2_ratio = 0.7902; // FIXED, sum lt. buehlmann
|
|
|
2586 pres_respiration = (float)int_I_pres_respiration / 1000.0; // assembler code uses different digit system
|
|
|
2587 pres_surface = (float)int_I_pres_surface / 1000.0;
|
|
|
2588 temp_atem = N2_ratio * (pres_respiration - 0.0627); // 0.0627 is the extra pressure in the body
|
|
|
2589 temp2_atem = 0.0;
|
|
|
2590 temp_surface = pres_surface; // the b"uhlmann formula using temp_surface does not use the N2_ratio
|
|
|
2591 float_desaturation_multiplier = char_I_desaturation_multiplier / 142.0; // new in v.101 (70,42%/100.=142)
|
|
|
2592 float_saturation_multiplier = char_I_saturation_multiplier / 100.0;
|
|
|
2593
|
|
|
2594 calc_tissue_step_1_min(); // update the pressure in the 16 tissues in accordance with the new ambient pressure
|
|
|
2595 clear_decoarray();
|
|
|
2596 char_O_deco_status = 0;
|
|
|
2597 char_O_nullzeit = 0;
|
|
|
2598 char_O_ascenttime = 0;
|
|
|
2599 calc_gradient_factor();
|
|
|
2600
|
|
|
2601 } // calc_wo_deco_step_1_min(void)
|
|
|
2602
|
|
|
2603
|
|
|
2604 // -------------------------
|
|
|
2605 // calc_tissue_step_1_min //
|
|
|
2606 // -------------------------
|
|
|
2607 // optimized in v.101
|
|
|
2608
|
|
|
2609 void calc_tissue_step_1_min(void)
|
|
|
2610 {
|
|
|
2611 _asm
|
|
|
2612 lfsr 1, 0x300
|
|
|
2613 movlw 0x01
|
|
|
2614 movwf TBLPTRU,0
|
|
|
2615 _endasm
|
|
|
2616
|
|
|
2617 char_O_gtissue_no = 255;
|
|
|
2618 pres_gtissue_limit = 0.0;
|
|
|
2619
|
|
|
2620 for (ci=0;ci<16;ci++)
|
|
|
2621 {
|
|
|
2622 _asm
|
|
|
2623 movlw 0x02
|
|
|
2624 movwf TBLPTRH,0
|
|
|
2625 movlb 4 // fuer ci
|
|
|
2626 movf ci,0,1
|
|
|
2627 addwf ci,0,1
|
|
|
2628 addwf ci,0,1
|
|
|
2629 addwf ci,0,1
|
|
|
2630 addlw 0x80
|
|
|
2631 movwf TBLPTRL,0
|
|
|
2632 TBLRDPOSTINC
|
|
|
2633 movff TABLAT,var_a+1
|
|
|
2634 TBLRDPOSTINC
|
|
|
2635 movff TABLAT,var_a
|
|
|
2636 TBLRDPOSTINC
|
|
|
2637 movff TABLAT,var_a+3
|
|
|
2638 TBLRD
|
|
|
2639 movff TABLAT,var_a+2
|
|
|
2640 addlw 0x40
|
|
|
2641 movwf TBLPTRL,0
|
|
|
2642 TBLRDPOSTINC
|
|
|
2643 movff TABLAT,var2_a+1
|
|
|
2644 TBLRDPOSTINC
|
|
|
2645 movff TABLAT,var2_a
|
|
|
2646 TBLRDPOSTINC
|
|
|
2647 movff TABLAT,var2_a+3
|
|
|
2648 TBLRD
|
|
|
2649 movff TABLAT,var2_a+2
|
|
|
2650 addlw 0x40
|
|
|
2651 movwf TBLPTRL,0
|
|
|
2652 incf TBLPTRH,1,0
|
|
|
2653 TBLRDPOSTINC
|
|
|
2654 movff TABLAT,var_b+1
|
|
|
2655 TBLRDPOSTINC
|
|
|
2656 movff TABLAT,var_b
|
|
|
2657 TBLRDPOSTINC
|
|
|
2658 movff TABLAT,var_b+3
|
|
|
2659 TBLRD
|
|
|
2660 movff TABLAT,var_b+2
|
|
|
2661 addlw 0x40
|
|
|
2662 movwf TBLPTRL,0
|
|
|
2663 TBLRDPOSTINC
|
|
|
2664 movff TABLAT,var2_b+1
|
|
|
2665 TBLRDPOSTINC
|
|
|
2666 movff TABLAT,var2_b
|
|
|
2667 TBLRDPOSTINC
|
|
|
2668 movff TABLAT,var2_b+3
|
|
|
2669 TBLRD
|
|
|
2670 movff TABLAT,var2_b+2
|
|
|
2671 addlw 0xC0
|
|
|
2672 movwf TBLPTRL,0
|
|
|
2673 incf TBLPTRH,1,0
|
|
|
2674 TBLRDPOSTINC
|
|
|
2675 movff TABLAT,var_e1min+1
|
|
|
2676 TBLRDPOSTINC
|
|
|
2677 movff TABLAT,var_e1min
|
|
|
2678 TBLRDPOSTINC
|
|
|
2679 movff TABLAT,var_e1min+3
|
|
|
2680 TBLRD
|
|
|
2681 movff TABLAT,var_e1min+2
|
|
|
2682 addlw 0x40
|
|
|
2683 movwf TBLPTRL,0
|
|
|
2684 TBLRDPOSTINC
|
|
|
2685 movff TABLAT,var2_e1min+1
|
|
|
2686 TBLRDPOSTINC
|
|
|
2687 movff TABLAT,var2_e1min
|
|
|
2688 TBLRDPOSTINC
|
|
|
2689 movff TABLAT,var2_e1min+3
|
|
|
2690 TBLRD
|
|
|
2691 movff TABLAT,var2_e1min+2
|
|
|
2692 _endasm
|
|
|
2693
|
|
|
2694 // N2 1 min
|
|
|
2695 temp_tissue = (temp_atem - pres_tissue[ci]) * var_e1min;
|
|
|
2696 temp_tissue_safety();
|
|
|
2697 pres_tissue[ci] = pres_tissue[ci] + temp_tissue;
|
|
|
2698
|
|
|
2699 // He 1 min
|
|
|
2700 temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var2_e1min;
|
|
|
2701 temp_tissue_safety();
|
|
|
2702 pres_tissue[ci+16] = pres_tissue[ci+16] + temp_tissue;
|
|
|
2703
|
|
|
2704 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16];
|
|
|
2705 var_a = (var_a * pres_tissue[ci] + var2_a * pres_tissue[ci+16]) / temp_tissue;
|
|
|
2706 var_b = (var_b * pres_tissue[ci] + var2_b * pres_tissue[ci+16]) / temp_tissue;
|
|
|
2707 pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
|
|
|
2708 if (pres_tissue_limit[ci] < 0)
|
|
|
2709 pres_tissue_limit[ci] = 0;
|
|
|
2710 if (pres_tissue_limit[ci] > pres_gtissue_limit)
|
|
|
2711 {
|
|
|
2712 pres_gtissue_limit = pres_tissue_limit[ci];
|
|
|
2713 char_O_gtissue_no = ci;
|
|
|
2714 }//if
|
|
|
2715
|
|
|
2716 // gradient factor array for graphical display
|
|
|
2717 // display range is 0 to 250! in steps of 5 for 1 pixel
|
|
|
2718 // the display is divided in 6 blocks
|
|
|
2719 // -> double the gradient 100% = 200
|
|
|
2720 // tissue > respiration (entsaettigungsvorgang)
|
|
|
2721 // gradient ist wieviel prozent an limit von tissue aus
|
|
|
2722 // dh. 0% = respiration == tissue
|
|
|
2723 // dh. 100% = respiration == limit
|
|
|
2724 temp1 = temp_tissue - pres_respiration;
|
|
|
2725 temp2 = temp_tissue - pres_tissue_limit[ci]; // changed in v.102
|
|
|
2726 temp2 = temp1/temp2;
|
|
|
2727 temp2 = temp2 * 200; // because of output in (Double-)percentage
|
|
|
2728 if (temp2 < 0)
|
|
|
2729 temp2 = 0;
|
|
|
2730 if (temp2 > 255)
|
|
|
2731 temp2 = 255;
|
|
|
2732 if (temp1 < 0)
|
|
|
2733 char_O_array_gradient_weighted[ci] = 0;
|
|
|
2734 else
|
|
|
2735 char_O_array_gradient_weighted[ci] = (char)temp2;
|
|
|
2736
|
|
|
2737 } // for
|
|
|
2738 } // calc wo deco 1min
|
|
|
2739
|
|
|
2740 #if 0
|
|
|
2741 // --------
|
|
|
2742 // debug //
|
|
|
2743 // --------
|
|
|
2744 void debug(void)
|
|
|
2745 {
|
|
|
2746 for (ci=0;ci<32;ci++)
|
|
|
2747 {
|
|
|
2748 int_O_tissue_for_debug[ci] = (unsigned int)(pres_tissue[ci] *1000);
|
|
|
2749 }
|
|
|
2750 } // void debug(void)
|
|
|
2751 #endif
|
|
|
2752
|
|
|
2753 // ----------
|
|
|
2754 // md hash //
|
|
|
2755 // ----------
|
|
|
2756 void hash(void)
|
|
|
2757 {
|
|
|
2758 // init
|
|
|
2759 for (md_i=0;md_i<16;md_i++)
|
|
|
2760 {
|
|
|
2761 md_state[md_i] = 0;
|
|
|
2762 md_cksum[md_i] = 0;
|
|
|
2763 } // for md_i 16
|
|
|
2764
|
|
|
2765 _asm
|
|
|
2766 movlw 0x01
|
|
|
2767 movwf TBLPTRU,0
|
|
|
2768 movlw 0x06
|
|
|
2769 movwf TBLPTRH,0
|
|
|
2770 movlw 0x00
|
|
|
2771 movwf TBLPTRL,0
|
|
|
2772 _endasm
|
|
|
2773 for (md_i=0;md_i<127;md_i++)
|
|
|
2774 {
|
|
|
2775 _asm
|
|
|
2776 TBLRDPOSTINC
|
|
|
2777 movff TABLAT,md_temp
|
|
|
2778 _endasm
|
|
|
2779 md_pi_subst[md_i] = md_temp;
|
|
|
2780 } // for md_i 256
|
|
|
2781 _asm
|
|
|
2782 TBLRDPOSTINC
|
|
|
2783 movff TABLAT,md_temp
|
|
|
2784 _endasm
|
|
|
2785 md_pi_subst[127] = md_temp;
|
|
|
2786 for (md_i=0;md_i<127;md_i++)
|
|
|
2787 {
|
|
|
2788 _asm
|
|
|
2789 TBLRDPOSTINC
|
|
|
2790 movff TABLAT,md_temp
|
|
|
2791 _endasm
|
|
|
2792 md_pi_subst[md_i+128] = md_temp;
|
|
|
2793 } // for md_i 256
|
|
|
2794 _asm
|
|
|
2795 TBLRD
|
|
|
2796 movff TABLAT,md_temp
|
|
|
2797 _endasm
|
|
|
2798 md_pi_subst[255] = md_temp;
|
|
|
2799
|
|
|
2800 _asm
|
|
|
2801 movlw 0x00
|
|
|
2802 movwf TBLPTRU,0
|
|
|
2803 movlw 0x00
|
|
|
2804 movwf TBLPTRH,0
|
|
|
2805 movlw 0x00
|
|
|
2806 movwf TBLPTRL,0
|
|
|
2807 _endasm
|
|
|
2808 // cycle buffers
|
|
|
2809 for (md_pointer=0x0000;md_pointer<0x17f3;md_pointer++)
|
|
|
2810 {
|
|
|
2811 md_t = 0;
|
|
|
2812 for (md_i=0;md_i<16;md_i++)
|
|
|
2813 {
|
|
|
2814 if(md_pointer == 9)
|
|
|
2815 md_temp = md_cksum[md_i];
|
|
|
2816 else
|
|
|
2817 {
|
|
|
2818 _asm
|
|
|
2819 TBLRDPOSTINC
|
|
|
2820 movff TABLAT,md_temp
|
|
|
2821 _endasm
|
|
|
2822 } // else
|
|
|
2823 md_buffer[md_i] = md_temp;
|
|
|
2824 md_state[md_i+16] = md_buffer[md_i];
|
|
|
2825 md_state[md_i+32] = (unsigned char)(md_buffer[md_i] ^ md_state[md_i]);
|
|
|
2826 } // for md_i 16
|
|
|
2827
|
|
|
2828 for (md_i=0;md_i<18;md_i++)
|
|
|
2829 {
|
|
|
2830 for (md_j=0;md_j<48;md_j++)
|
|
|
2831 {
|
|
|
2832 md_state[md_j] = (unsigned char)(md_state[md_j] ^ md_pi_subst[md_t]);
|
|
|
2833 md_t = md_state[md_j];
|
|
|
2834 } // for md_j 48
|
|
|
2835 md_t = (unsigned char)(md_t+1);
|
|
|
2836 } // for md_i 18
|
|
|
2837 md_t = md_cksum[15];
|
|
|
2838
|
|
|
2839 for (md_i=0;md_i<16;md_i++)
|
|
|
2840 {
|
|
|
2841 md_cksum[md_i] = (unsigned char)(md_cksum[md_i] ^ md_pi_subst[(md_buffer[md_i] ^ md_t)]);
|
|
|
2842 md_t = md_cksum[md_i];
|
|
|
2843 } // for md_i 16
|
|
|
2844 } // for md_pointer
|
|
|
2845 } // void hash(void)
|
|
|
2846
|
|
|
2847 // ---------------------
|
|
|
2848 // clear_CNS_fraction //
|
|
|
2849 // ---------------------
|
|
|
2850 // new in v.101
|
|
|
2851
|
|
|
2852 void clear_CNS_fraction(void)
|
|
|
2853 {
|
|
|
2854 CNS_fraction = 0.0;
|
|
|
2855 char_O_CNS_fraction = 0;
|
|
|
2856 } // void clear_CNS_fraction(void)
|
|
|
2857
|
|
|
2858
|
|
|
2859 // --------------------
|
|
|
2860 // calc_CNS_fraction //
|
|
|
2861 // --------------------
|
|
|
2862 // new in v.101
|
|
|
2863 // optimized in v.102 : with new variables char_I_actual_ppO2 and actual_ppO2
|
|
|
2864
|
|
|
2865 // Input: char_I_actual_ppO2
|
|
|
2866 // Output: char_O_CNS_fraction
|
|
|
2867 // Uses and Updates: CNS_fraction
|
|
|
2868 // Uses: acutal_ppO2
|
|
|
2869
|
|
|
2870 void calc_CNS_fraction(void)
|
|
|
2871 {
|
|
|
2872 actual_ppO2 = (float)char_I_actual_ppO2 / 100.0;
|
|
|
2873
|
|
|
2874 if (char_I_actual_ppO2 < 50)
|
|
|
2875 CNS_fraction = CNS_fraction;// no changes
|
|
|
2876 else if (char_I_actual_ppO2 < 60)
|
|
|
2877 CNS_fraction = 1/(-54000.0 * actual_ppO2 + 54000.0) + CNS_fraction;
|
|
|
2878 else if (char_I_actual_ppO2 < 70)
|
|
|
2879 CNS_fraction = 1/(-45000.0 * actual_ppO2 + 48600.0) + CNS_fraction;
|
|
|
2880 else if (char_I_actual_ppO2 < 80)
|
|
|
2881 CNS_fraction = 1/(-36000.0 * actual_ppO2 + 42300.0) + CNS_fraction;
|
|
|
2882 else if (char_I_actual_ppO2 < 90)
|
|
|
2883 CNS_fraction = 1/(-27000.0 * actual_ppO2 + 35100.0) + CNS_fraction;
|
|
|
2884 else if (char_I_actual_ppO2 < 110)
|
|
|
2885 CNS_fraction = 1/(-18000.0 * actual_ppO2 + 27000.0) + CNS_fraction;
|
|
|
2886 else if (char_I_actual_ppO2 < 150)
|
|
|
2887 CNS_fraction = 1/(-9000.0 * actual_ppO2 + 17100.0) + CNS_fraction;
|
|
|
2888 else if (char_I_actual_ppO2 < 160)
|
|
|
2889 CNS_fraction = 1/(-22500.0 * actual_ppO2 + 37350.0) + CNS_fraction;
|
|
|
2890 else if (char_I_actual_ppO2 < 165)
|
|
|
2891 CNS_fraction = 0.000755 + CNS_fraction; // Arieli et all.(2002): Modeling pulmonary and CNS O2 toxicity... Formula (A1) based on value for 1.55 and c=20
|
|
|
2892 else if (char_I_actual_ppO2 < 170)
|
|
|
2893 CNS_fraction = 0.00102 + CNS_fraction; // example calculation: Sqrt((1.7/1.55)^20)*0.000404
|
|
|
2894 else if (char_I_actual_ppO2 < 175)
|
|
|
2895 CNS_fraction = 0.00136 + CNS_fraction;
|
|
|
2896 else if (char_I_actual_ppO2 < 180)
|
|
|
2897 CNS_fraction = 0.00180 + CNS_fraction;
|
|
|
2898 else if (char_I_actual_ppO2 < 185)
|
|
|
2899 CNS_fraction = 0.00237 + CNS_fraction;
|
|
|
2900 else if (char_I_actual_ppO2 < 190)
|
|
|
2901 CNS_fraction = 0.00310 + CNS_fraction;
|
|
|
2902 else if (char_I_actual_ppO2 < 195)
|
|
|
2903 CNS_fraction = 0.00401 + CNS_fraction;
|
|
|
2904 else if (char_I_actual_ppO2 < 200)
|
|
|
2905 CNS_fraction = 0.00517 + CNS_fraction;
|
|
|
2906 else if (char_I_actual_ppO2 < 230)
|
|
|
2907 CNS_fraction = 0.0209 + CNS_fraction;
|
|
|
2908 else
|
|
|
2909 CNS_fraction = 0.0482 + CNS_fraction; // value for 2.5
|
|
|
2910
|
|
|
2911 if (CNS_fraction > 2.5)
|
|
|
2912 CNS_fraction = 2.5;
|
|
|
2913 if (CNS_fraction < 0.0)
|
|
|
2914 CNS_fraction = 0.0;
|
|
|
2915 char_O_CNS_fraction = (char)((CNS_fraction + 0.005)* 100.0);
|
|
|
2916 } // void calc_CNS_fraction(void)
|
|
|
2917
|
|
|
2918 // --------------------------
|
|
|
2919 // calc_CNS_decrease_15min //
|
|
|
2920 // --------------------------
|
|
|
2921 // new in v.101
|
|
|
2922
|
|
|
2923 // calculates the half time of 90 minutes in 6 steps of 15 min
|
|
|
2924
|
|
|
2925 // Output: char_O_CNS_fraction
|
|
|
2926 // Uses and Updates: CNS_fraction
|
|
|
2927
|
|
|
2928 void calc_CNS_decrease_15min(void)
|
|
|
2929 {
|
|
|
2930 CNS_fraction = 0.890899 * CNS_fraction;
|
|
|
2931 char_O_CNS_fraction = (char)((CNS_fraction + 0.005)* 100.0);
|
|
|
2932 }// calc_CNS_decrease_15min(void)
|
|
|
2933
|
|
|
2934
|
|
|
2935 // ------------------
|
|
|
2936 // calc_percentage //
|
|
|
2937 // ------------------
|
|
|
2938 // new in v.101
|
|
|
2939
|
|
|
2940 // calculates int_I_temp * char_I_temp / 100
|
|
|
2941 // output is int_I_temp
|
|
|
2942
|
|
|
2943 void calc_percentage(void)
|
|
|
2944 {
|
|
|
2945 temp1 = (float)int_I_temp;
|
|
|
2946 temp2 = (float)char_I_temp / 100.0;
|
|
|
2947 temp3 = temp1 * temp2;
|
|
|
2948 int_I_temp = (int)temp3;
|
|
|
2949 }
|
|
|
2950 void push_tissues_to_vault(void)
|
|
|
2951 {
|
|
|
2952 for (ci=0;ci<32;ci++)
|
|
|
2953 pres_tissue_vault[ci] = pres_tissue[ci];
|
|
|
2954 }
|
|
|
2955 void pull_tissues_from_vault(void)
|
|
|
2956 {
|
|
|
2957 for (ci=0;ci<32;ci++)
|
|
|
2958 pres_tissue[ci] = pres_tissue_vault[ci];
|
|
|
2959 }
|
|
|
2960
|
|
|
2961 void wp_write_command(void)
|
|
|
2962 {
|
|
|
2963 _asm
|
|
|
2964 bcf oled_rs
|
|
|
2965 movff wp_command,PORTD
|
|
|
2966 bcf oled_rw
|
|
|
2967 bsf oled_rw
|
|
|
2968 _endasm
|
|
|
2969 }
|
|
|
2970
|
|
|
2971 void wp_write_data(void)
|
|
|
2972 {
|
|
|
2973 wp_data_8bit_one = wp_data_16bit >> 8;
|
|
|
2974 wp_data_8bit_two = wp_data_16bit;
|
|
|
2975 _asm
|
|
|
2976 bsf oled_rs
|
|
|
2977 movff wp_data_8bit_one,PORTD
|
|
|
2978 bcf oled_rw
|
|
|
2979 bsf oled_rw
|
|
|
2980 movff wp_data_8bit_two,PORTD
|
|
|
2981 bcf oled_rw
|
|
|
2982 bsf oled_rw
|
|
|
2983 _endasm
|
|
|
2984 }
|
|
|
2985
|
|
|
2986 void wp_write_black(void)
|
|
|
2987 {
|
|
|
2988 _asm
|
|
|
2989 movff wp_black,PORTD
|
|
|
2990 bcf oled_rw
|
|
|
2991 bsf oled_rw
|
|
|
2992 bcf oled_rw
|
|
|
2993 bsf oled_rw
|
|
|
2994 _endasm
|
|
|
2995 }
|
|
|
2996
|
|
|
2997 void wp_write_color(void)
|
|
|
2998 {
|
|
|
2999 _asm
|
|
|
3000 movff wp_color1,PORTD
|
|
|
3001 bcf oled_rw
|
|
|
3002 bsf oled_rw
|
|
|
3003 movff wp_color2,PORTD
|
|
|
3004 bcf oled_rw
|
|
|
3005 bsf oled_rw
|
|
|
3006 _endasm
|
|
|
3007 }
|
|
|
3008
|
|
|
3009 void wp_set_window(void)
|
|
|
3010 {
|
|
|
3011 // x axis start ( 0 - 319)
|
|
|
3012 wp_command = 0x35;
|
|
|
3013 wp_write_command();
|
|
|
3014 wp_data_16bit = ((U16)wp_leftx2) << 1;
|
|
|
3015 wp_write_data();
|
|
|
3016 // x axis end ( 0 - 319)
|
|
|
3017 wp_command = 0x36;
|
|
|
3018 wp_write_command();
|
|
|
3019 wp_data_16bit = 319;
|
|
|
3020 wp_write_data();
|
|
|
3021 // y axis start + end ( 0 - 239 )
|
|
|
3022 wp_command = 0x37;
|
|
|
3023 wp_write_command();
|
|
|
3024 // the bottom part
|
|
|
3025 wp_data_16bit = wp_top;
|
|
|
3026 if(wp_font == 2)
|
|
|
3027 wp_data_16bit += WP_FONT_LARGE_HEIGHT;
|
|
|
3028 else if(wp_font == 1)
|
|
|
3029 wp_data_16bit += WP_FONT_MEDIUM_HEIGHT;
|
|
|
3030 else
|
|
|
3031 wp_data_16bit += WP_FONT_SMALL_HEIGHT;
|
|
|
3032 wp_data_16bit--;
|
|
|
3033 if(wp_data_16bit > 239)
|
|
|
3034 wp_data_16bit = 239;
|
|
|
3035 // the top part
|
|
|
3036 wp_data_16bit |= ((U16)wp_top) << 8;
|
|
|
3037 // all together in one 16bit transfer
|
|
|
3038 wp_write_data();
|
|
|
3039
|
|
|
3040 // start
|
|
|
3041 wp_command = 0x20;
|
|
|
3042 wp_write_command();
|
|
|
3043 wp_data_16bit = wp_top;
|
|
|
3044 wp_write_data();
|
|
|
3045
|
|
|
3046 wp_command = 0x21;
|
|
|
3047 wp_write_command();
|
|
|
3048 wp_data_16bit = ((U16)wp_leftx2) << 1;
|
|
|
3049 wp_write_data();
|
|
|
3050 }
|
|
|
3051
|
|
|
3052 void wp_set_char_font_small(void)
|
|
|
3053 {
|
|
|
3054 // space is A1
|
|
|
3055 if (wp_char > 0x7E) // skip space between ~ and ¡
|
|
|
3056 wp_char -= 34;
|
|
|
3057
|
|
|
3058 if (wp_char == ' ')
|
|
|
3059 wp_char = 0xA1;
|
|
|
3060
|
|
|
3061 if((wp_char < '!') || (wp_char > 0xA1)) // font has 34 chars after ~ // ¾ + 4 chars limit to end of battery at the moment
|
|
|
3062 wp_char = 0x82; // ¤
|
|
|
3063
|
|
|
3064 wp_start = wp_small_table[wp_char - '!'];
|
|
|
3065 wp_end = wp_small_table[1 + wp_char - '!'];
|
|
|
3066 }
|
|
|
3067
|
|
|
3068 void wp_set_char_font_medium(void)
|
|
|
3069 {
|
|
|
3070 // space is 3E
|
|
|
3071 if (wp_char == 0x27) // 0x27 == '
|
|
|
3072 wp_char = 0x3B;
|
|
|
3073 if (wp_char == '"')
|
|
|
3074 wp_char = 0x3C;
|
|
|
3075 if (wp_char == 'm')
|
|
|
3076 wp_char = 0x3D;
|
|
|
3077 if (wp_char == ' ')
|
|
|
3078 wp_char = 0x3E;
|
|
|
3079
|
|
|
3080 if((wp_char < '.') || (wp_char > 0x3E))
|
|
|
3081 wp_char = 0x3E;
|
|
|
3082 wp_start = wp_medium_table[wp_char - '.'];
|
|
|
3083 wp_end = wp_medium_table[1 + wp_char - '.'];
|
|
|
3084 }
|
|
|
3085
|
|
|
3086 void wp_set_char_font_large(void)
|
|
|
3087 {
|
|
|
3088 // space is / = 0x2F
|
|
|
3089 if (wp_char == ' ')
|
|
|
3090 wp_char = 0x2F;
|
|
|
3091
|
|
|
3092 if((wp_char < '.') || (wp_char > '9'))
|
|
|
3093 wp_char = 0x2F;
|
|
|
3094 wp_start = wp_large_table[wp_char - '.'];
|
|
|
3095 wp_end = wp_large_table[1 + wp_char - '.'];
|
|
|
3096 }
|
|
|
3097
|
|
|
3098 void wordprocessor(void)
|
|
|
3099 {
|
|
|
3100 wp_set_window();
|
|
|
3101
|
|
|
3102 // access to GRAM
|
|
|
3103 wp_command = 0x22;
|
|
|
3104 wp_write_command();
|
|
|
3105 _asm
|
|
|
3106 bsf oled_rs
|
|
|
3107 _endasm
|
|
|
3108
|
|
|
3109 wp_txtptr = 0;
|
|
|
3110 wp_char = wp_stringstore[wp_txtptr];
|
|
|
3111
|
|
|
3112 while(wp_char)
|
|
|
3113 {
|
|
|
3114 if(wp_font == 2)
|
|
|
3115 wp_set_char_font_large();
|
|
|
3116 else if(wp_font == 1)
|
|
|
3117 wp_set_char_font_medium();
|
|
|
3118 else
|
|
|
3119 wp_set_char_font_small();
|
|
|
3120
|
|
|
3121 wp_black = 0;
|
|
|
3122
|
|
|
3123 for(wp_i = wp_start; wp_i<wp_end;wp_i++)
|
|
|
3124 {
|
|
|
3125 if(wp_font == 2)
|
|
|
3126 wp_data_16bit = wp_large_data[wp_i / 2];
|
|
|
3127 else if(wp_font == 1)
|
|
|
3128 wp_data_16bit = wp_medium_data[wp_i / 2];
|
|
|
3129 else
|
|
|
3130 wp_data_16bit = wp_small_data[wp_i / 2];
|
|
|
3131 if(wp_i & 1)
|
|
|
3132 wp_temp_U8 = wp_data_16bit & 0xFF;
|
|
|
3133 else
|
|
|
3134 wp_temp_U8 = wp_data_16bit >> 8;
|
|
|
3135 if((wp_temp_U8 & 128))
|
|
|
3136 {
|
|
|
3137 wp_temp_U8 -= 127;
|
|
|
3138 if(wp_invert)
|
|
|
3139 {
|
|
|
3140 while(wp_temp_U8 > 0)
|
|
|
3141 {
|
|
|
3142 wp_temp_U8--;
|
|
|
3143 wp_write_color();
|
|
|
3144 }
|
|
|
3145 }
|
|
|
3146 else
|
|
|
3147 {
|
|
|
3148 _asm
|
|
|
3149 movff wp_black,PORTD
|
|
|
3150 _endasm
|
|
|
3151 while(wp_temp_U8 > 0)
|
|
|
3152 {
|
|
|
3153 wp_temp_U8--;
|
|
|
3154 _asm
|
|
|
3155 bcf oled_rw
|
|
|
3156 bsf oled_rw
|
|
|
3157 bcf oled_rw
|
|
|
3158 bsf oled_rw
|
|
|
3159 _endasm
|
|
|
3160 }
|
|
|
3161 }
|
|
|
3162 }
|
|
|
3163 else
|
|
|
3164 {
|
|
|
3165 wp_temp_U8++;
|
|
|
3166 if(wp_invert)
|
|
|
3167 {
|
|
|
3168 _asm
|
|
|
3169 movff wp_black,PORTD
|
|
|
3170 _endasm
|
|
|
3171 while(wp_temp_U8 > 0)
|
|
|
3172 {
|
|
|
3173 wp_temp_U8--;
|
|
|
3174 _asm
|
|
|
3175 bcf oled_rw
|
|
|
3176 bsf oled_rw
|
|
|
3177 bcf oled_rw
|
|
|
3178 bsf oled_rw
|
|
|
3179 _endasm
|
|
|
3180 }
|
|
|
3181 }
|
|
|
3182 else
|
|
|
3183 {
|
|
|
3184 while(wp_temp_U8 > 0)
|
|
|
3185 {
|
|
|
3186 wp_temp_U8--;
|
|
|
3187 wp_write_color();
|
|
|
3188 }
|
|
|
3189 }
|
|
|
3190 }
|
|
|
3191 }
|
|
|
3192 wp_txtptr++;
|
|
|
3193 wp_char = wp_stringstore[wp_txtptr];
|
|
|
3194 }
|
|
|
3195 wp_command = 0x00;
|
|
|
3196 wp_write_command();
|
|
|
3197 wp_top = 0;
|
|
|
3198 wp_leftx2 = 0;
|
|
|
3199 wp_font = 0;
|
|
|
3200 wp_invert = 0;
|
|
|
3201 }
|
|
|
3202
|
