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

use "../../data/processed/ocean_physics_1993_1998_reference.nc"
use "/Volumes/nascratch/OBS/HadISST/HadISST_sst.nc"

let tmask = if sst gt -2 then 1
let sst01=sst*tmask 
let sst02=sst01[l=@ave]

let msalOg=temp[d=1,gx=sst02,gy=sst02,k=1]
let dtemp=msalOg[l=1]-sst02[l=1]
set variable/title="MOMSO - HadISST Temperature" dtemp
set variable/units="degC" dtemp
save/clobber/file="../../data/processed/momsoMinusHadISSTemp.nc" dtemp