Maintenance scheduled for Thursday, February 27th at 15:00 MDT. Expected downtime <1 hour.

Commit e024ec6d authored by Asquith, William H.'s avatar Asquith, William H.

added graphical support for aquifer top drawing

parent 5f9bf402
......@@ -68,6 +68,7 @@ for(i in 1:num.sites) { # MASTER LOOP START: Mloop
refalt <- TMP$SPECIAL_REFALT_VA[1]
alt <- TMP$ALT_VA[1] # the altitude of the land surface in the alt_datum_cd
aqfrBttm <- TMP$AQFR_BOTTOM_ALT_VA[1] # the optional bottom of the aquifer
aqfrTop <- TMP$AQFR_TOP_ALT_VA[1] # the optional top of the aquifer
lev <- TMP$LEV_ALT_VA # SPEED FIX
casingTop <- alt - TMP$CSNG_TOP_VA[1] # this is the maximum of all casing tops (can be multiple in NWIS)
casingBttm <- alt - TMP$CSNG_BOTTOM_VA[1] # this is the minimum of all casing bottoms (can be multiple in NWIS)
......@@ -109,11 +110,13 @@ for(i in 1:num.sites) { # MASTER LOOP START: Mloop
ylim <- c(-ylim.dn.buffer*yf,ylim.up.buffer*yf) + range(lev)
if(TOGGLE.force.ylim.to.alt) ylim <- range(c(ylim,alt,refalt), na.rm=TRUE)
if(TOGGLE.show.aquifer.bottom) ylim <- range(c(ylim,aqfrBttm), na.rm=TRUE)
if(TOGGLE.show.aquifer.top ) ylim <- range(c(ylim,aqfrTop), na.rm=TRUE)
if(length(FLO$ALT_VA) != 0) ylim <- range(c(ylim,alt), na.rm=TRUE)
} else {
ylim <- range(c(alt, lev, casingTop, casingBttm, screenTop,
screenBttm, welldepth, refalt), na.rm=TRUE)
if(TOGGLE.show.aquifer.bottom) ylim <- range(c(ylim,aqfrBttm), na.rm=TRUE)
if(TOGGLE.show.aquifer.top ) ylim <- range(c(ylim,aqfrTop), na.rm=TRUE)
}
if(length(ylim.lock.to.this.interval) == 2) ylim <- ylim.lock.to.this.interval
options(opts) # resetting the warnings to on again
......@@ -213,7 +216,10 @@ for(i in 1:num.sites) { # MASTER LOOP START: Mloop
lines(c(usr[1]+0.2*(usr[2]-usr[1]), usr[2]), rep(aqfrBttm,2), lty=aqfr.bottom.lty, lwd=aqfr.bottom.lwd, col=aqfr.bottom.col)
text(usr[1], aqfrBttm, "aquifer bottom", cex=aqfr.bottom.cex, pos=4, col=aqfr.bottom.col)
}
if(TOGGLE.show.aquifer.top ) {
lines(c(usr[1]+0.2*(usr[2]-usr[1]), usr[2]), rep(aqfrTop,2), lty=aqfr.top.lty, lwd=aqfr.top.lwd, col=aqfr.top.col)
text(usr[1], aqfrTop, "aquifer top", cex=aqfr.top.cex, pos=4, col=aqfr.top.col)
}
if(showGRIDLINES) {
draw.decadeGridLines(TMP, ylim=usr[3:4], xlim=xlim)
draw.annualGridLines(TMP, ylim=usr[3:4], xlim=xlim)
......@@ -509,6 +515,10 @@ for(i in 1:num.sites) { # MASTER LOOP START: Mloop
lines(c(usr[1]+0.2*(usr[2]-usr[1]), usr[2]), rep(aqfrBttm,2), lty=aqfr.bottom.lty, lwd=aqfr.bottom.lwd, col=aqfr.bottom.col)
text(usr[1], aqfrBttm, "aquifer bottom", cex=aqfr.bottom.cex, pos=4, col=aqfr.bottom.col)
}
if(TOGGLE.show.aquifer.top ) {
lines(c(usr[1]+0.2*(usr[2]-usr[1]), usr[2]), rep(aqfrTop,2), lty=aqfr.top.lty, lwd=aqfr.top.lwd, col=aqfr.top.col)
text(usr[1], aqfrTop, "aquifer top", cex=aqfr.top.cex, pos=4, col=aqfr.top.col)
}
if(showGRIDLINES) {
draw.decadeGridLines(TMP, ylim=usr[3:4], xlim=xlim)
draw.annualGridLines(TMP, ylim=usr[3:4], xlim=xlim)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment