data(drag, package="ocedata")
par(mar=c(3,3,1,1), mgp=c(2, 0.7, 0))
plot(drag$U, drag$Cd, pch=ifelse(drag$method == "profile", 1, 19),
xlab="U [m/s]", ylab=expression(C[D]),
xlim=c(0, 22), ylim=c(0, 3e-3))
legend("topleft", pch=c(1,19),
legend=c("profile method", "eddy covariance method"))
Run the code above in your browser using DataLab