File Exchange

image thumbnail

Steady state snow and firn density model

version 1.0.0.0 (3.46 KB) by Aslak Grinsted
Calculates snow/firn depth density profiles using Herron-Langway type models.

1 Download

Updated 29 Jul 2014

View License

Calculates steady state snow/firn depth density profiles using Herron-Langway type models.

[rho,zieq]=densitymodel(Tavg,bdot,rhos,z,model)

Herron-Langway type models. (Arthern et al. 2010 formulation).

INPUT:
Tavg: 10m temperature in celcius
bdot: accumulation rate in mwe/yr or (kg/m2/yr)
rhos: surface density in kg/m3
z: depth in true_metres

model can be: {'HLJohnsen' 'HerronLangway' 'LiZwally' 'Helsen' 'NabarroHerring'}
default is herronlangway. (The other models are tuned for non-stationary modelling (Read Arthern et al.2010 before applying in steady state).


OUTPUT:
rho: density (kg/m3) for all z-values.
zieq: ice equivalent depth for all z-values.
t: age for all z-values (only taking densification into account.)

Example usage:
z=0:300;
[rho,zieq,t]=densitymodel(-31.5,177,340,z,'HerronLangway');
plot(z,rho)

Aslak Grinsted 2010

Cite As

Aslak Grinsted (2021). Steady state snow and firn density model (https://www.mathworks.com/matlabcentral/fileexchange/47386-steady-state-snow-and-firn-density-model), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (0)

MATLAB Release Compatibility
Created with R2010a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!