Approximation of Melet 2013 Fig 2, Using CESM Case P1

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