Westerschelde time series

## ============================================================================
## NIOZ westerschelde timeseries of selected monitoring variables
## data are in OceanView, called niozWS
## implemented by Karline Soetaert
## ============================================================================
 library(OceanView)
par(mfrow = c(5, 1), mar = c(3, 3, 3, 2))

varname <- c("SPMCHLA", "WNH4", "WNO3", "WPO4", "WCTEMP" )

Long <- WSnioz[as.character(WSnioz$VariableName) %in% varname & 
            WSnioz$Station == 10,
            c("SamplingDateTimeREAL", "VariableName", "DataValue")]

Cross <- db2cross(Long)

col <- ramp.col(c("lightblue", "black"))

scatter2D(x = Cross$x, y = Cross$z[,"SPMCHLA"], colvar = Cross$z[,"WCTEMP"],
          pch = 18, cex = 1.5, type = "b", main = "Chlorophyll", 
          colkey = FALSE, col = col, NAcol = "red")

scatter2D(x = Cross$x, y = Cross$z[,"WNH4"], colvar = Cross$z[,"WCTEMP"],
          pch = 18, cex = 1.5, type = "b", main = "Ammonium", 
          colkey = FALSE, col = col, NAcol = "red")

scatter2D(x = Cross$x, y = Cross$z[,"WNO3"], colvar = Cross$z[,"WCTEMP"],
          pch = 18, cex = 1.5, type = "b", main = "Nitrate", 
          colkey = FALSE, col = col, NAcol = "red")

scatter2D(x = Cross$x, y = Cross$z[,"WPO4"], colvar = Cross$z[,"WCTEMP"],
          pch = 18, cex = 1.5, type = "b", main = "Phosphate", 
          colkey = FALSE, col = col, NAcol = "red")

colkey (side = 3, clim = range(Cross$z[,"WCTEMP"], na.rm = TRUE), 
  clab = "Temperature", line.clab = 3, col = col, length = 0.5)

mtext(outer = TRUE, side = 1, "Station 10", line = -1)