2018 Oct 16 11am MST

1. Carolyn now testing sloping surface

   a. Overview of code modifications

   i. Fixed dependency issue in initializing slope case: needed to gather number of grid points in x-y plane before initializing reference buoyancy profile.

   ii. Initializing slope case with full buoyancy: added reference salinity profile and salinity slope offset.

   iii. Applying offsets in advection scheme if cyclic lateral boundaries: added salinity offset.

   iv. Buoyancy term in momentum equation: calculate buoyancy using density field rather than temperature field.

   b. Test cases

   i. buoyancy forcing by temperature flux vs. buoyancy forcing by salinity flux

   ii. buoyancy forcing in atmosphere vs. equivalent buoyancy forcing in ocean

   iii. buoyancy forcing by salinity flux in low and high slope cases

   iv. Other recommendations?

2. Next step: new option for ice melting

   a. At each timestep, if constant flux layer or ice melting, call a module to calculate eddy diffusivities (exchange velocities) based on McPhee (1983) M-O Similarity Theory.

   b. Choose depth at which to define "mixed-layer/far-field" temperature and salinity. Choose at what length scale to average these quantities.

   c. Solve for melt rate, interface salinity and interface temperature using 3-equations. Only store heat flux and salt flux. Recalculate other values when outputting.