oce (version 0.8-3)

ctdTrim: Trim start/end portions of a CTD cast

Description

Trim start/end portions of a CTD cast.

Usage

ctdTrim(x, method=c("downcast","index","range"), 
inferWaterDepth=TRUE, removeDepthInversions=FALSE,
parameters, debug=getOption("oceDebug"))

Arguments

x

A ctd object, e.g. as read by read.ctd.

method

Various methods exist, some of which use parameters: [object Object]

"index"{Select values only in the list of indices specified in parameters. The indices may be integers, e.g. parameters=10:30<

Value

• An object of class "ctd", with data having been trimmed in some way.

item

• inferWaterDepth
• removeDepthInversions
• parameters
• debug

code

ctdDecimate

Details

For normal CTD profiling, the goal is to isolate only the downcast, discarding measurements made in the air, in an equilibration phase in which the device is held below the water surface, and then the upcast phase that follows the downcast. This is handled reasonably well by ctdTrim with method="downcast". (The datasets provided with oce were produced this way.)

However, for detailed work it makes sense to do things semi-manually. The eye is simply better at handling exceptional cases. The process is simple: use plotScan() to get an idea of the scan indices of the downacast, and then use ctdTrim with method="index". A few trials will normally identify the downcast very well.

References

The Seabird CTD instrument is described at http://www.seabird.com/products/spec_sheets/19plusdata.htm.

See Also

The documentation for ctd-class explains the structure of CTD objects, and also outlines the other functions dealing with them.

Examples

library(oce)
data(ctdRaw) 
plot(ctdRaw) # barely recognizable, due to pre- and post-cast junk
plot(ctdTrim(ctdRaw)) # looks like a real profile ...
plot(ctdDecimate(ctdTrim(ctdRaw),method="boxcar")) # ... smoothed