V_s -(Q_s+em*sigma*(T_amb^4-T^4))/(rho_ice*L_ice) Sublimation interface velocity flc1hs(z-1e-5,1e-5) P_v 1[Pa]*exp(9.550426-5723.265[K]/T+3.53068*log(T[1/K])-0.00728332[1/K]*T) Vapor partial pressure Q_s down(fluxz_ht)*dnz+down(fluxy_ht)*dny+down(fluxx_ht)*dnx Normal heat flux c_init c_ice+(c_c-c_ice)*(Z-Z0)/(22[mm]-Z0) Initial vapor concentration T_init T_ice+(T0-T_ice)*(Z-Z0)/(22[mm]-Z0) Initial temperature distribution Q_in 5.5[kW/m^2]*((x/R0)^2+(y/R0)^2)*(x<0) Heat flux at vial bottom