us_model: delete vbar updating. correct the equation of calculation o…

…f D from f_f0.
ehb54 · Jan 25, 2025 · 4b74572 · 4b74572
1 parent cd1079a
commit 4b74572
Showing 1 changed file with 0 additions and 7 deletions.
diff --git a/utils/us_model.cpp b/utils/us_model.cpp
@@ -197,12 +197,6 @@ bool US_Model::calc_coefficients( SimulationComponent& component )
    // Unless custom grid is used, we set the vbar to the reciprocal of the positive s-value vbar.
 
    buoyancyb      = 1.0 - ( vbar * DENS_20W );
-   // if ( ( s < 0.0  &&  buoyancyb > 0.0 )  ||
-   //      ( s > 0.0  &&  buoyancyb < 0.0 ) )
-   // {
-   //    buoyancyb      = 1.0 - ( DENS_20W / vbar );
-   //    vbar           = 1.0 / vbar;
-   // }
 
    // Start with already calculated s if possible
    if ( s != 0.0 )
@@ -215,7 +209,6 @@ bool US_Model::calc_coefficients( SimulationComponent& component )
       {                                          ///////////////
          // Please see calc_coefficient_documentation.pdf in this directory,
          // equation 14, with adjustments for units in poise.
-         // double numer   = 0.02 * vbar * s * VISC_20W * f_f0;
          double numer   = 0.02 * vbar * s * VISC_20W;
          f0             = 0.09 * VISC_20W * M_PI * sqrt( numer / buoyancyb );
          fv             = f_f0 * f0;