GLS MIXING

Wilkin: Here we would write a more verbose description of the GLS scheme (but not a full algorithm/science description). Basically just enough to help a user with a meaningful choice among the following suboptions:

The following CPP options modify the configuration of the Generic Length-Scale vertical turbulence closure:

CANUTO_A CANUTO_B CHARNOK CRAIG_BANNER KANTHA_CLAYSON K_C2ADVECTION K_C4ADVECTION N2S2_HORAVG

There are parameters in ocean.in that determine whether GLS_MIXING implements the turbulence closure options of k-kl, k-$\epsilon$, or k-$\omega$