mdot kvap*rhoL*(p_star-pG)/pG*step2(thetaL)*step1(p_star/pG) "Evaporation rate"
p_star 10^(A-B/(C+T))*1e5[Pa] "Vapor pressure from Antoine equation"
thetaG 1-thetaL-thetaS "Volume fraction of gas"
DL alpha*(thetaL-thetaL_star)*step1(thetaL/thetaL_star) "Apparent moisture diffusivity"
thetaL0 thetaS*(rhoS/rhoL)*(wL0/(1-wL0)) "Initial volume fraction of liquid"