mf <- par(mfrow=c(2,1))
ma <- par(mar=c(3,5,2,5))
dsal <- t(sw_sfac$del_sa[1,,])
dsal [dsal < -90]<-NA
image(sw_sfac$longs, sw_sfac$lats, dsal, col=femmecol(100),
asp=TRUE,xlab="dg",ylab="dg",
main="salinity conversion - p = 0 bar")
contour(sw_sfac$longs, sw_sfac$lats, dsal, asp=TRUE, add=TRUE)
dsal <- t(sw_sfac$del_sa[5,,]) # 5th depth level sw_sfac$p[5]
dsal [dsal < -90]<-NA
image(sw_sfac$longs, sw_sfac$lats, dsal, col=femmecol(100),
asp=TRUE, xlab="dg", ylab="dg",
main="salinity conversion - p = 4 bar")
contour(sw_sfac$longs, sw_sfac$lats, dsal, asp=TRUE, add=TRUE)
par("mfrow"=mf)
par("mar"=ma)
Run the code above in your browser using DataLab