2 simple tests from Section 7.5.1 of the CVMix manual:

* computing w_m when Bf = 0 m^2/s^3 and u* ~= 0.014 m/s ("neutral forcing" from 7.5.1.3)
* computing w_m and w_s when u* = 0 m/s, Bf ~= -5.9e-8 m^2 / s^3, and h = 6000 m (convective limit test from 7.5.1.4)

None of the infrastructure around these tests has been built (no namelist, no output beyond what is printed to screen), but both tests produce turbulent scale values that match the manual. Here's the additional output (this is "Test 6" in the kpp driver):

 6a: Bf = 0 m^2/s^3 and u* = sqrt(1/5175) m/s
     => w_m = w_s ~= 0.0056 m/s
 w_m =   5.56038437485532715E-003
 
 6b: u* = 0 m/s and Bf ~= -5.9e-8 m^2/s^3
     => w_m = 0.4 * (-Bf * OBL_depth)^1/3 * (c_m*0.04)^1/3 m/s
        w_s = 0.4 * (-Bf * OBL_depth)^1/3 * (c_s*0.04)^1/3 m/s
     (so w_m ~= 0.020 m/s and w_s ~= 0.045 m/s)
 w_m =   1.96908794643690342E-002
 w_s =   4.48353646932846961E-002
 
I think I'll have all the physical parameters from the manual (tao, rho0, g, alpha, etc) as namelist variables and then print out both the analytic and numeric values. First, though, I want to code up test 6c, which is finding the minimum value of the unresolved velocity scale needed to compute the bulk Richardson number (Section 7.5.6.5 in the manual)



git-svn-id: https://cvmix.googlecode.com/svn/trunk@264 ec264166-b127-864d-6cb2-e100125a5c76
Realized that the unresolved velocity scale test was unrelated to the turbulent scales tests, so I went ahead and added namelist support for the two tests from the previous commit. Also tweaked the output a little bit:

 Test 6: 2 simple tests for velocity scale
 ----------
 6a: Bf = 0 m^2/s^3 and u* = sqrt(tao/rho0)
                           =  1.39009609371383187E-002
     => w_m = w_s ~= vonkarman*u*
                   =  5.56038437485532801E-003
 w_m =   5.56038437485532801E-003
 w_s =   5.56038437485532801E-003
 
 6b: u* = 0 m/s and Bf = (grav*alpha/(rho0*Cp0))*Qnonpen
                       = -5.92973386386299205E-008
     => w_m = vonkarman * (-Bf * OBL_depth)^1/3 * (c_m*0.1*vonkarman)^1/3 m/s
            =   1.96908794643690342E-002
     => w_s = vonkarman * (-Bf * OBL_depth)^1/3 * (c_s*0.1*vonkarman)^1/3 m/s
            =   4.48353646932846961E-002
 w_m =   1.96908794643690342E-002
 w_s =   4.48353646932846961E-002



git-svn-id: https://cvmix.googlecode.com/svn/trunk@265 ec264166-b127-864d-6cb2-e100125a5c76