use mean evaporation rate for q_flux (Kg/m^2/s)