;
;jeff snider, univeristy of wyoming, march 2006
;
function satd_isentrope, in
;
common  constants, gascon, rho_h2o, mw_nh42so4, rho_nh42so4, sigma_0,      $
        sigma_t, dsigmadmolality_nh42so4, dsigmadmolality_nacl, mw_h2o,    $
        mw_nacl, rho_nacl, tkmelt, c_nh42so4, a_nh42so4, c_nacl, a_nacl,   $
        cp_air_o, cp_h2o_o, mw_air, gravity, p0, alpha, beta_local,        $
        gascon_h2o, cw_h2o_o, tk_o, lv_o, ew_o, tk_aerosol, salt_type,     $
        epsilon, c_twomey, k_twomey, aerosol_type, right_tail,             $
        gascon_air, tc_base, tk_base, p_base, h_start, hmax,               $
        r0, dt, rmin, rmax, coef, epsilon_aerosol, specific_volume_meas,   $
        nchan, mw_salt, rho_salt, a_salt, c_salt, specific_volume_base,    $
        dsigmadmolality_salt, rho_insoluble, volume_insoluble_to_soluble,  $
        mixrat_tot_1,r_1,geo_sigma_1,mixrat_tot_2,r_2,geo_sigma_2,         $
        tk_start,p_start,sratio_start,mixrat_tot_3,r_3,geo_sigma_3
;
common known_values_satd, rl_a, p_b
;
;print, p_base, tk_base, epsilon, gascon, cw_h2o_o, rl_a, p_b
;
return, $
[in(0) - epsilon * in(1) / (p_base - in(1)), $
 in(2) - epsilon * in(3) / (p_b - in(3)), $
 in(4) - (cp_air_o+(in(0)+rl_a)*cw_h2o_o)*alog(tk_base/tk_o)+gascon_air*alog(in(5))-in(0)*in(6)/tk_base, $
 in(7) - (cp_air_o+(in(2)+in(9))*cw_h2o_o)*alog(in(8)/tk_o)+gascon_air*alog(in(10))-in(2)*in(11)/in(8), $
 p_base   - in(5) - in(1), $
 p_b   - in(10) - in(3), $
 in(1) - ew_jeff(tk_base), $
 in(3) - ew_jeff(in(8)), $
 lv_o - in(6) - (cp_h2o_o-cw_h2o_o)*(tk_o-tk_base), $
 lv_o - in(11) - (cp_h2o_o-cw_h2o_o)*(tk_o-in(8)), $
 in(0) - rl_a - in(2) - in(9), $
 in(4) - in(7)]
;
end
;