mass_domain at(0,intop1(hs.rho))-intop1(hs.rho) "Mass change in the domain"
mass_bulk intop2(Hk_Si_bulk*2329[kg/m^3]) "Total mass of bulk species"
mass_Si_surf intop2(Zk_Si_surf*hs.Gammat/hs.sigmak_Si_surf*0.028[kg/mol]) "Total mass of Si surface species"
mass_SiH_surf intop2(Zk_SiH_surf*hs.Gammat/hs.sigmak_SiH_surf*0.029[kg/mol]) "Total mass of SiH surface species"
dmass_Si_surf mass_Si_surf-at(0,mass_Si_surf) "Total mass minus initial mass of Si surface species"
dmass_SiH_surf mass_SiH_surf-at(0,mass_SiH_surf) "Total mass minus initial mass of SiH surface species"
mass_surf dmass_Si_surf+dmass_SiH_surf "Total mass minus initial mass of all surface species"