!calculating the difference between simulated and modeled sss
cancel da/all
cancel re/all
se mem/size=5000

use "/Volumes/nascratch/soeddy/data/processed/ocean_bling_1993_1998_reference.nc"
use "/Volumes/nascratch/OBS/woa_2009/dissolved_oxygen_annual_1deg.nc"


let msalOg=o2_b[d=1,gx=o_an[d=2],gy=o_an[d=2],gz=o_an[d=2]]
let do2=msalOg[l=1]*1e6-o_an[l=1]/22391.6*1000.0*1000.0
set variable/title="MOMSO - WOA09 O2" do2
set variable/units="[mmolO2/m^3]" do2


save/clobber/file="/Volumes/nascratch/soeddy/data/processed/momsoMinusWOAO2.nc" do2