Vs -(Qs+em*sigma*(T_amb^4-T^4))/rho_ice/L_ice Sublimation interface velocity flc1hs(z-1e-5,1e-5) Pv 1[Pa]*exp(9.550426-5723.265[K]/T+3.53068*log(T*1[1/K])-0.00728332[1/K]*T) Vapor partial pressure Qs 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)/(0.022[m]-Z0) Initial vapor concentration T_init T_ice+(T0-T_ice)*(Z-Z0)/(0.022[m]-Z0) Initial temperature distribution Qin 5.5e3[W/m^2]*((x/R0)^2+(y/R0)^2)*(x>0) Heat flux at vial bottom