mifr/subsurface
1
0
Fork 0
mirror of https://github.com/subsurface/subsurface.git synced 2025-01-31 18:43:24 +00:00
Code Issues Projects Releases Packages Wiki Activity

Exclude surface intervals from sac rate calculation

Browse source 
We assume every sample with a depth of less than 10cm to be on the
surface.

This does not impact our interpolated pressures (one could assume that the
diver is not breathing from the regulator when on the surface - but
without air integration that's just an assumption).

It also doesn't change our tank pressure coloring by sac rate as that
always uses the momentary sac rate. Technically speaking this might impact
the actual colors printed (as those are relative to the total sac on the
dive which may go up due to this change).

Signed-off-by: Dirk Hohndel <dirk@hohndel.org>
This commit is contained in:
Dirk Hohndel 2011-11-21 13:23:13 -08:00
parent 0d1f8f9a5d
commit b577484998
1 changed files with 18 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

19
divelist.c
View file

@ -320,6 +320,7 @@ static double calculate_airuse(struct dive *dive)
static int calculate_sac(struct dive *dive)
{
double airuse, pressure, sac;
int duration, i;
airuse = calculate_airuse(dive);
if (!airuse)
@ -327,9 +328,25 @@ static int calculate_sac(struct dive *dive)
if (!dive->duration.seconds)
return 0;
/* find and eliminate long surface intervals */
duration = dive->duration.seconds;
for (i = 0; i < dive->samples; i++) {
if (dive->sample[i].depth.mm < 100) { /* less than 10cm */
int end = i + 1;
while (end < dive->samples && dive->sample[end].depth.mm < 100)
end++;
/* we only want the actual surface time during a dive */
if (end < dive->samples) {
end--;
duration -= dive->sample[end].time.seconds -
dive->sample[i].time.seconds;
i = end + 1;
}
}
}
/* Mean pressure in atm: 1 atm per 10m */
pressure = 1 + (dive->meandepth.mm / 10000.0);
sac = airuse / pressure * 60 / dive->duration.seconds;
sac = airuse / pressure * 60 / duration;
/* milliliters per minute.. */
return sac * 1000;