Mercurial > public > ostc4

comparison Discovery/Src/buehlmann.c @ 1052:3029dd4398a5 GasConsumption tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
1.7.5 beta release, fix Surface-GF
author heinrichsweikamp
date Mon, 15 Dec 2025 17:56:25 +0100
parents 88b6ab90c55a
children
comparison
equal deleted inserted replaced
1051:4b6afe5551e1 1052:3029dd4398a5
342 tissue_inertgas_saturation = gTissue_nitrogen_bar[ci] + gTissue_helium_bar[ci]; 342 tissue_inertgas_saturation = gTissue_nitrogen_bar[ci] + gTissue_helium_bar[ci];
343 inertgas_a = ( ( buehlmann_N2_a[ci] * gTissue_nitrogen_bar[ci]) + ( buehlmann_He_a[ci] * gTissue_helium_bar[ci]) ) / tissue_inertgas_saturation; 343 inertgas_a = ( ( buehlmann_N2_a[ci] * gTissue_nitrogen_bar[ci]) + ( buehlmann_He_a[ci] * gTissue_helium_bar[ci]) ) / tissue_inertgas_saturation;
344 inertgas_b = ( ( buehlmann_N2_b[ci] * gTissue_nitrogen_bar[ci]) + ( buehlmann_He_b[ci] * gTissue_helium_bar[ci]) ) / tissue_inertgas_saturation; 344 inertgas_b = ( ( buehlmann_N2_b[ci] * gTissue_nitrogen_bar[ci]) + ( buehlmann_He_b[ci] * gTissue_helium_bar[ci]) ) / tissue_inertgas_saturation;
345 } 345 }
346 346
347 M_surf = inertgas_a + inertgas_b * pres_surface; 347 M_surf = pres_surface / inertgas_b + inertgas_a;
348 if (M_surf > pres_surface) 348 if (M_surf > pres_surface)
349 { 349 {
350 gf_surf = (tissue_inertgas_saturation - pres_surface) / (M_surf - pres_surface); 350 gf_surf = (tissue_inertgas_saturation - pres_surface) / (M_surf - pres_surface);
351 if(M_surf > 9.99) 351 if(M_surf > 9.99)
352 { 352 {