Approximation of Melet 2013 Fig 1, Using CESM Case P1

Note: Case P1 = "2DPM" with time-varying E(x,y,t) = f(Nb)(giaf.c112.tidal.2DPM.test1.Nb.bugfix)
          

Note: Data for interior ocean mean potential temperature are unavailable from case J3