## JPM April 2015.
## 
## Blue whale voyage call locations.
## 
## Timestep is days since sampling began (DSSB).
## 
## Coloured points are sonobuoy locations.
## 
## Black points are all whale call locations.
## 
## Coloured vectors are whale calls from sonobuoy.
library(animation)
library(lattice)
par(mar = c(4, 4, 0.5, 0.5))
for (i in 1:(max.days + 1)) {
    j <- i - 1
    if (j %in% days.floor) {
        with(input, plot(whaleCallLongitude, whaleCallLatitude, 
            xlim = xl, ylim = yl, pch = 16, cex = 0.1, col = 1))
        with(input, points(sonobuoyLongitude, sonobuoyLatitude, 
            xlim = xl, ylim = yl, pch = 16, cex = 1, col = cols[input$dssf + 
              1]))
        legend("right", legend = min(input$dssf):max(input$dssf), 
            fill = cols, bty = "n", cex = 0.7, title = "DSSB")
        d <- input[input$dssf == j, ]
        with(d, points(sonobuoyLongitude, sonobuoyLatitude, 
            pch = 16, cex = 0.8, col = cols[d$dssf + 1]))
        with(d, segments(sonobuoyLongitude, sonobuoyLatitude, 
            whaleCallLongitude, whaleCallLatitude, col = cols[d$dssf + 
              1]))
    }
    ani.pause()
}
## R version 3.1.3 (2015-03-09)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Other packages: animation 2.3, lattice 0.20-31