θ_conv_variable_Pr = (
                      Cᵂu★  = 5.18,
                      CᵂwΔ  = 7.37,
                      Cᴷu⁻  = 0.0939,
                      Cᴷc⁻  = 0.835,
                      Cᴷe⁻  = 1.04,
                      Cᴷuʳ  = -0.508,
                      Cᴷcʳ  = +0.256,
                      Cᴷeʳ  = -0.0696,
                      Cᴬc   = 0.544,
                      Cᴬe   = 9.77e-3,
                      Cᴰ    = 0.283,
                      Cᴸᵇ   = 0.0561,
                      Cᴬu   = 0.0,
                      CᴷRiʷ = 0.0,
                      CᴷRiᶜ = 0.0,
                      Cᵟu   = 0.5,
                      Cᵟc   = 0.5,
                      Cᵟe   = 0.5)

θ_variable_Pr = (
                 CᵂwΔ  = 14.1,
                 Cᵂu★  = 2.55,
                 Cᴰ    = 0.354,
                 Cᴸᵇ   = 0.159,
                 Cᴷu⁻  = 0.0963,
                 Cᴷc⁻  = 0.279,
                 Cᴷe⁻  = 1.03,
                 Cᴷuʳ  = -0.525,
                 Cᴷcʳ  = +0.670,
                 Cᴷeʳ  = +1.14,
                 CᴷRiʷ = 0.111,
                 CᴷRiᶜ = 1.98,
                 Cᴬu   = 0.0,
                 Cᴬc   = 0.0,
                 Cᴬe   = 0.0,
                 Cᵟu   = 0.5,
                 Cᵟc   = 0.5,
                 Cᵟe   = 0.5,
                )

θ_constant_Pr = (
                 CᵂwΔ  = 8.66,
                 Cᵂu★  = 2.89,
                 Cᴰ    = 0.489,
                 Cᴸᵇ   = 0.0286,
                 Cᴷu⁻  = 0.0776,
                 Cᴷc⁻  = 0.567,
                 Cᴷe⁻  = 9.00,
                 Cᴷuʳ  = 0.0,
                 Cᴷcʳ  = 0.0,
                 Cᴷeʳ  = 0.0,
                 CᴷRiʷ = 0.0,
                 CᴷRiᶜ = 0.0,
                 Cᴬu   = 0.0,
                 Cᴬc   = 0.0,
                 Cᴬe   = 0.0,
                 Cᵟu   = 0.5,
                 Cᵟc   = 0.5,
                 Cᵟe   = 0.5,
                )

