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 (2020). 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