src/oceanbox.cpp
// 2 ways of solving for the Revelle factor
// keep track of last year pCO2 in the ocean and DIC
2 ways of solving for the Revelle factor
keep track of last year pCO2 in the ocean and DIC
2 ways of solving for the Revelle factor keep track of last year pCO2 in the ocean and DIC
unitval oceanbox::calc_revelle() {
    H_ASSERT( active_chemistry, "Active Chemistry required");
    unitval deltadic = mychemistry.convertToDIC(carbon) - dic_lastyear;
    H_ASSERT( deltadic.value( U_UMOL_KG) != 0, "DeltaDIC can not be zero");
     return unitval ( mychemistry.convertToDIC( carbon ) / mychemistry.CO3, U_UNITLESS ); 
}
unitval oceanbox :: calc_revelle ( ) { H_ASSERT ( active_chemistry , " " ) ; unitval deltadic = mychemistry . convertToDIC ( carbon ) - dic_lastyear ; H_ASSERT ( deltadic . value ( U_UMOL_KG ) != 0 , " " ) ; return unitval ( mychemistry . convertToDIC ( carbon ) / mychemistry . CO3 , U_UNITLESS ) ; }
2 ways of solving for the Revelle factor
keep track of last year pCO2 in the ocean and DIC
2 ways of solving for the Revelle factor keep track of last year pCO2 in the ocean and DIC
//    unitval deltapco2 = Ca - pco2_lastyear; // Revelle Factor can be calculated multiple ways.   // Based on changing atmospheric conditions as well approximated via DIC and CO3 // under high CO2, the HL box numbers are potentially unrealistic. 