BVF_MIXING upper and lower bounds swapped

Message

m.hadfield · Tue Feb 20, 2007 2:16 am

In mod_scalars.F the code that sets the parameters for the BVF_MIXING option is:

#ifdef BVF_MIXING
!
!-----------------------------------------------------------------------
!  Brunt-Vaisala frequency based vertical mixing variables.
!-----------------------------------------------------------------------
!
!    bvf_numax     Upper bound vertical diffusion (m2/s).
!    bvf_numin     Lower bound vertical diffusion (m2/s).
!    bvf_nu0       Proportionality constant (m2/s2).
!    bvf_nu0c      Convective diffusion (m2/s) in static unstable
!                    regime.
!
        real(r8) :: bvf_numax = 3.0E-5_r8      ! m2/s
        real(r8) :: bvf_numin = 4.0E-4_r8      ! m2/s
        real(r8) :: bvf_nu0 = 1.0E-7_r8        ! m2/s2
        real(r8) :: bvf_nu0c = 1.0_r8          ! m2/s
#endif

Surely that should be...

Code: Select all

        real(r8) :: bvf_numin = 3.0E-5_r8      ! m2/s
        real(r8) :: bvf_numax = 4.0E-4_r8      ! m2/s
[/code