mifr/subsurface
1
0
Fork 0
mirror of https://github.com/subsurface/subsurface.git synced 2024-11-28 05:00:20 +00:00
Code Issues Projects Releases Packages Wiki Activity

Fix up SAC calculations for ATM/bar confusion

Browse source 
We even documented that we did SAC in bar*l/min, but the "S" in SAC
stands for "Surface".  So we should normalize SAC rate to surface
pressure, not one bar.

It's a tiny 1% difference, and doesn't actually matter in practice, but
it's noticeable when you want to explicitly test for SAC-rate by
creating a test-dive that averages exactly 10m.  Suddenly you don't get
the round numbers you expect.

[ Side note: 10m is not _exactly_ one extra atmosphere according to our
  calculations, but it's darn close in sea water: the standard salinity
  of 1.03 kg/l together with the standard acceleration of 9.81m/s^2
  gives an additional pressure of 1.01 bar, which is within a fraction
  of a percent of one ATM.

  Of course, divers have likely chosen that value exactly for the math
  to come out that way, since the true average salinity of seawater is
  actually slightly lower ]

So here's a few test-dives, along with the SAC rate fixup to make them
look right.

(There's also a one-liner to dive.c that makes the duration come out
right if the last sample has a non-zero depth, and the previous sample
did not: one of my original test-dives did the "average 10m depth" by
starting at 0 and ending at 20m, and dive.c got a tiny bit confused
about that ;)

[ The rationale for me testing our SAC rate calculations in the first
  place was that on snorkkeli.net user "Poltsi" reported that our SAC rate
  calculations differ from the ones that Suunto DM4 reports. So I wanted
  to verify that we did things right.

  Note that Poltsi reported differences larger than the difference of
  BAR/ATM, so this is not the cause. I'll continue to look at this. ]

Signed-off-by: Linus Torvalds <torvalds@linux-foundation.org>
Signed-off-by: Dirk Hohndel <dirk@hohndel.org>
This commit is contained in:
Linus Torvalds 2013-02-24 10:01:18 -08:00 committed by Dirk Hohndel
parent 1f20d5feed
commit f4bf16d5db
3 changed files with 53 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

1
dive.c
View file

@ -607,6 +607,7 @@ static void fixup_dive_dc(struct dive *dive, struct divecomputer *dc)
end = time;
if (depth > SURFACE_THRESHOLD) {
end = time;
if (start < 0)
start = lasttime;
if (depth > maxdepth)

4
divelist.c
View file

@ -749,8 +749,8 @@ static int calculate_sac(struct dive *dive)
}
}
}
/* Mean pressure in bar (SAC calculations are in bar*l/min) */
pressure = depth_to_mbar(dc->meandepth.mm, dive) / 1000.0;
/* Mean pressure in ATM (SAC calculations are in atm*l/min) */
pressure = (double) depth_to_mbar(dc->meandepth.mm, dive) / SURFACE_PRESSURE;
sac = airuse / pressure * 60 / duration;
/* milliliters per minute.. */

50
dives/sac-test.xml Normal file
View file

@ -0,0 +1,50 @@
<divelog program='subsurface' version='2'>
<settings>
</settings>
<dives>
<dive number='1' date='2013-02-24' time='09:26:04' duration='25:00 min'>
<location>SAC test dive 1</location>
<notes>SAC should be 20l/min (0.7063 cuft/min)
We have exactly 100 ATM of gas in a 10l cylinder, average 10m for 25min</notes>
<cylinder size='10 l' description='10l' start='202.6 bar' end='101.3 bar' />
<divecomputer>
<depth max='10.0 m' mean='10.0 m' />
<sample time='0:00 min' depth='10.0 m' />
<sample time='5:00 min' depth='10.0 m' />
<sample time='10:00 min' depth='10.0 m' />
<sample time='15:00 min' depth='10.0 m' />
<sample time='20:00 min' depth='10.0 m' />
<sample time='25:00 min' depth='10.0 m' />
</divecomputer>
</dive>
<dive number='2' date='2013-02-24' time='10:26:04' duration='25:00 min'>
<location>SAC test dive 2</location>
<notes>SAC should be 20l/min (0.7063 cuft/min)
We have exactly 100 ATM of gas in a 10l cylinder, average 10m for 25min</notes>
<cylinder size='10 l' description='10l' start='202.6 bar' end='101.3 bar' />
<divecomputer>
<depth max='20.0 m' mean='10.0 m' />
<sample time='0:00 min' depth='0.0 m' />
<sample time='12:30 min' depth='20.0 m' />
<sample time='25:00 min' depth='0.0 m' />
</divecomputer>
</dive>
<dive number='3' date='2013-02-24' time='11:26:04' duration='25:00 min'>
<location>SAC test dive 3</location>
<notes>SAC should be 20l/min (0.7063 cuft/min)
We have exactly 100 ATM of gas in 15l of cylinder, average 20m for 25min</notes>
<cylinder size='10 l' description='10l' start='202.6 bar' end='101.3 bar' />
<cylinder size='5 l' description='10l' start='202.6 bar' end='101.3 bar' />
<divecomputer>
<event time='17:00 min' name='gaschange' />
<depth max='20.0 m' mean='10.0 m' />
<sample time='0:00 min' depth='0.0 m' />
<sample time='5:00 min' depth='20.0 m' />
<sample time='10:00 min' depth='40.0 m' />
<sample time='15:00 min' depth='0.0 m' />
<sample time='20:00 min' depth='40.0 m' />
<sample time='25:00 min' depth='0.0 m' />
</divecomputer>
</dive>
</dives>
</divelog>