deco.ceiling {scuba} | R Documentation |
Decompression Ceiling for a Diver
Description
Given the state of a diver's tissues, find the diver's decompression ceiling depth (minimum permissible depth) or decompression ceiling pressure (minimum permissible ambient pressure).
Usage
deco.ceiling(state, model, what = c("depth", "pressure"))
Arguments
state |
Diver's tissue saturations. See Details. |
model |
Haldane type decompression model (object of class |
what |
What to calculate. See Details. |
Details
Given the current state of a diver's tissues, this command calculates the diver's decompression ceiling: the shallowest depth (or lowest ambient pressure) to which it is safe to ascend.
The argument state
gives the diver's current state
or the diver's progressive state during a dive.
It has typically been calculated using haldane
.
It should be either:
-
A numeric vector giving the diver's current nitrogen tensions for each tissue. The length of the vector matches the number of tissues in the
model
. -
A matrix giving the diver's current nitrogen and (if present) helium tensions for each tissue. The rows of the matrix correspond to tissues in the model. The first column (or the column labelled
"N2"
) contains the nitrogen tensions, and if present, the second column (or column labelled"He"
) contains the helium tensions. -
A three-dimensional array giving the progressive states of the diver over time. The first index corresponds to time. The second index corresponds to tissues in the model. The third coordinate corresponds to inert gases
"N2"
and"He"
.
Entries in the vector, matrix or array state
are pressures
in atmospheres absolute (ata).
The minimum tolerated ambient pressure Pamb.tol
is calculated
using the equivalent of Buehlmann's formula (Buehlmann et al, 1995,
section 6.5, page 117). If what="pressure"
, this minimum
tolerated ambient pressure is returned (negative pressures are reset
to zero). If what="depth"
this pressure is converted to
diving depth in metres assuming that the surface pressure is 1 ata
(negative depths are reset to zero).
Value
A numeric vector, matrix or array of the same type as state
and with the same dimensions as state
,
giving the decompression ceiling for each tissue,
each inert gas, and each time point as given.
If what="pressure"
the entries in the result are pressures
in ata. If what="depth"
the entries are depths in metres.
Author(s)
Adrian Baddeley Adrian.Baddeley@curtin.edu.au.
References
Buehlmann, A.B., Voellm, E.B. and Nussberger, P. (2002) Tauchmedizin. 5th edition. Springer-Verlag, Berlin. ISBN 3-540-42979-4.
See Also
hm
for decompression models.
haldane
to calculate tissue tensions
during or after a dive.
Examples
# diver state after saturation diving on EANx30 at 20 metres
s <- saturated.state("Z", 20, nitrox(0.30))
# decompression ceiling for each compartment
depthceiling <- deco.ceiling(s, "Z")
round(depthceiling, 2)
round(max(depthceiling), 2)
# decompress at 8 metres or deeper
# flying-after-diving with cabin pressure 0.8 ata
pressureceiling <- deco.ceiling(s, "Z", "pressure")
round(pressureceiling, 2)
any(pressureceiling > 0.8)
# immediate flying-after-diving is not safe