0001 function displaymaterials(md)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 disp(sprintf(' Materials:'));
0012 disp(sprintf(' rho_ice: %g [kg/m^3]',md.rho_ice));
0013 disp(sprintf(' rho_water: %g [kg/m^3]',md.rho_water));
0014 if isempty(md.B), disp(sprintf(' B: N/A (flow law parameter)')); else disp(sprintf(' B: (%i) (flow law parameter [Pa/s^(1/n)])',length(md.B))); end
0015 if isempty(md.n), disp(sprintf(' n: N/A (flow law coefficient)')); else disp(sprintf(' n: (%i) (flow law coefficient)',length(md.n))); end
0016 disp(sprintf(' heatcapacity: %g [J/kg/K]',md.heatcapacity));
0017 disp(sprintf(' thermalconductivity: %g [W/m/K]',md.thermalconductivity));
0018 disp(sprintf(' meltingpoint: %g [K] (melting point of ice at 1atm in K)',md.meltingpoint));
0019 disp(sprintf(' latentheat: %g [J/kg] (latent heat of fusion)',md.latentheat));
0020 disp(sprintf(' beta: %g [K/Pa]',md.beta));
0021 disp(sprintf(' mixed_layer_capacity: %g [W/kg/K]',md.mixed_layer_capacity));
0022 disp(sprintf(' thermal_exchange_velocity: %g [m/s]',md.thermal_exchange_velocity));