{ plot(ice, axes = FALSE) plot(grat, add = TRUE, lty = 3) text(labs, lab = parse(text = labs$lab), col = c("firebrick", "darkblue")[labs$islon + 1], cex = 0.85) title(sprintf("Sea ice concentration %s", gsub(".bin", "", basename(tfile))), cex.main = 0.8) title(sub = projection(ice), cex.sub = 0.6) }