Changeset 2812


Ignore:
Timestamp:
Aug 14, 2019 3:51:08 PM (3 months ago)
Author:
gegorbet
Message:

density_match fix for floating percent valueand for use of smoothed s,D in other vector calcs

Location:
trunk/programs/us_density_match
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/programs/us_density_match/us_density_match.cpp

    r2811 r2812  
    13301330   {
    13311331      // *** Mi = si*R*T/(Di_avg*(1-vbari*rho))
    1332       double sedco     = alldis[ 0 ].bf_distro[ jj ].s * 1.0e-13;
    1333       double difco     = alldis[ 0 ].bf_distro[ jj ].d * 1.0e-7;
     1332//      double sedco     = alldis[ 0 ].bf_distro[ jj ].s * 1.0e-13;
     1333//      double difco     = alldis[ 0 ].bf_distro[ jj ].d * 1.0e-7;
     1334      double sedco     = v_sedcs[ 0 ][ jj ] * 1.0e-13;
     1335      double difco     = v_difcs[ 0 ][ jj ] * 1.0e-7;
    13341336      double vbari     = v_vbars[ jj ];
    13351337      double mmass     = sedco * R_GC * K20 / ( difco * ( 1.0 - vbari * DENS_20W ) );
     
    13511353      // *** fi/f_0i
    13521354      // *** ri = fi/(6 * pi * eta)   <-- hydrodynamic radius
    1353       double difco     = alldis[ 0 ].bf_distro[ jj ].d * 1.0e-7;
     1355//      double difco     = alldis[ 0 ].bf_distro[ jj ].d * 1.0e-7;
     1356      double difco     = v_difcs[ 0 ][ jj ] * 1.0e-7;
    13541357      double frico     = R_GC * K20 / ( difco * AVOGADRO );
    13551358      double hyrad     = frico / ( 6.0 * M_PI * VISC_20W );
     
    13721375      // *** fi/f_0i
    13731376      // *** ri = fi/(6 * pi * eta)   <-- hydrodynamic radius
    1374       double difco     = alldis[ 0 ].bf_distro[ jj ].d * 1.0e-7;
     1377//      double difco     = alldis[ 0 ].bf_distro[ jj ].d * 1.0e-7;
     1378      double difco     = v_difcs[ 0 ][ jj ] * 1.0e-7;
    13751379      double vbari     = v_vbars[ jj ];
    13761380//      double rzero     = pow( ( ( 0.75 / M_PI ) * vbari ), a_third );
  • trunk/programs/us_density_match/us_model_params.cpp

    r2811 r2812  
    269269         if ( !etext.isEmpty() )
    270270         {
    271             int pctval      = etext.toInt();
    272             if ( etext == "0" )
     271            double pctval   = etext.toDouble();
     272            if ( etext == "0"  ||  etext == "0.0" )
    273273            {  // Zero percent:  bump counts
    274274               nzd2pc++;
     
    276276DbgLv(1) << " vlf:    pctval" << pctval << "nzd2pc nad2pc" << nzd2pc << nad2pc;
    277277            }
    278             else if ( pctval != 0 )
     278            else if ( pctval > 0.0 )
    279279            {  // Non-zero numeric: bump percent count
    280280               nad2pc++;
     
    344344         if ( !etext.isEmpty() )
    345345         {
    346             int pctval      = etext.toInt();
    347             if ( etext == "0" )
     346            double pctval   = etext.toDouble();
     347            if ( etext == "0"  ||  etext == "0.0" )
    348348            {  // Count 0 percent and total percent
    349349               nzd2pc++;
    350350               nad2pc++;
    351351            }
    352             else if ( pctval != 0 )
     352            else if ( pctval > 0.0 )
    353353            {  // Count total non-empty percent values given
    354354               nad2pc++;
Note: See TracChangeset for help on using the changeset viewer.