set mode meta
go portrait1x2
set viewport top2
go transparency
!go margins 2.8 0.8 2.2 2.7
def vi/xli=0.154,0.846/yli=0.48,0.92 full1
set viewport full1

use "./data/suboxic.cdf"

plot/set/nolabel/line=13/vlimit=0:15/xlimit=1800:8000 suboxic[d=1,l=@sbx:11]/1e15
ppl axset,1,1,1,0
ppl ylab V_s_u_b_o_x (10^1^5 m^3)
ppl xlab TIME (year)
ppl plot

plot/ov/vs/line=1/dash/nolabel {3410,3410},{0,25}
plot/ov/vs/line=3/dash/nolabel {5687,5687},{0,25}

can data/all

use "./data/no3bar.cdf"
def vi/xli=0.154,0.846/yli=0.48,0.92 full2
set viewport full2

let ovol=1.358*1E+18
let ovol1=1.35769E+18

ppl axset,0,0,0,1 !draw  only the right-side axis
ppl axlabp,,1 !labels on the right
plot/set/nolabel/col=15/hlimit=1800:8000 no3bar*ovol/1e15
ppl ylab @P9NO_3 (Pmol)
ppl plot

list no3bar[t=5600:5900]*1000-23.80

quit