p_water c*R_const*T Water partial pressure p P0+p_water Absolute pressure M_vap (P0*M_air+p_water*M_water)/p Vapor molar mass c_ice p_eq(T_ice)/(R_const*T_ice) Initial vapor concentration at ice/vapor interface c_c p_eq(T_c)/(R_const*T_c) Vapor concentration at condenser c_init c_ice+(c_c-c_ice)*(Zg-Z0)/(H0-Z0) Initial vapor concentration T_init T_ice+(T_c-T_ice)*(Zg-Z0)/(H0-Z0) Initial temperature distribution Q_in if(x>0,5.5[kW/m^2]*((x/R0)^2+(y/R0)^2),0) Heat flux at vial bottom V_s T_lm[W/m^2]/(rho_ice*L_ice) Sublimation interface velocity