B (T-T_m+dT)/(2*dT)*((T<=(T_m+dT))*(T>=(T_m-dT)))+(T>(T_m+dT)) Fraction of liquid phase fB (1-B)^2/(B^3+epsilon) Volume force prefactor Fr -fB*A_mush*u Volume force density, r-component Fz -fB*A_mush*(w-v_cast) Volume force density, z-component D exp(-(T-T_m)^2/(dT^2))/(sqrt(pi)*dT) Temperature dependence, latent heat Cp1 380[J/(kg*K)]+(dH/T_m)*flc2hs((T-T_m)[1/K],dT[1/K]) Heat capacity + latent heat contribution