In cvmix_kpp_coeffs_low, the indexing in variables like Gat1, DGat1, visc_at_OBL isn't very clear... I think I might put in a module-wide parameter index (CVMIX_KPP_MDIFF_IND, etc) to clarify.
Rather than add in module-wide index, I just renamed variables... so instead of visc_at_OBL(3) code has Mdiff_OBL, Tdiff_OBL, and Sdiff_OBL (similar for shape function variables)
In cvmix_kpp_coeffs_low, the indexing in variables like Gat1, DGat1, visc_at_OBL isn't very clear... I think I might put in a module-wide parameter index (CVMIX_KPP_MDIFF_IND, etc) to clarify.