set mode meta
go portrait1x2
set viewport top2
def vi/xli=0.154,0.846/yli=0.48,0.92 full1
set viewport full1

!go transparency2
go transparency
!go margins 2.8 0.8 2.2 2.7

use "./data/sat.cdf"
use "./data/tbar.cdf"


let satdegc=sat-273.15
plot/set/nolabel/line=14/vlimit=3:22/xlimit=1800:8000 satdegc[d=1,l=@sbx:51]
ppl axset,1,1,1,0
ppl ylab @P4T_o_c, @P2SAT (^oC)
ppl xlab TIME (year)
ppl plot
plot/ov/nolabel/line=16 tbar[d=2]-273.15
!!plot/ov/sym tbar[d=2,t=3380]-273.15
!!plot/ov/sym satdegc[d=1,l=@sbx:51,t=2790]
plot/ov/vs/line=1/dash/nolabel {2280,2280},{0,25}
plot/ov/vs/line=2/dash/nolabel {2790,2790},{0,25}
plot/ov/vs/line=4/dash/nolabel {3380,3380},{0,25}
list satdegc[d=1,l=1], satdegc[d=1,t=2790], satdegc[d=1,t=8000]
list tbar[d=2,l=1], tbar[d=2,t=3380], tbar[d=2,t=8000]


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

set viewport full2

ppl axset,0,0,0,1 !draw  only the right-side axis
ppl axlabp,,1 !labels on the right
plot/set/nolabel/col=13/hlimit=1800:8000/vlimit=260:2200 co2ccn
ppl ylab pCO_2 (@SGm@SRatm)
ppl plot
!!plot/ov/sym co2ccn[t=2280]

!

quit