Now able to use namelist variables to control what gets run with CVMix and what POP handles natively; can toggle for all 4 wscale calls, computing boundary layer depth, computing shape functions, and computing the non-local term. Note that there is some dependence on the shape functions and the nonlocal term that I need to work out.