melt.py

melt.LWC_correct(self)

If there is LWC in a layer after temperature diffusion and the temperature is less than zero, one option is to just balance the energy to increase the temperature and lower the LWC. It isn’t the best way to solve the problem but it is one way.

This should be vectorized but that is not a priority.

melt.bucketVV(self, iii)

Bucket scheme for meltwater percolation, retention, refreezing, and runoff

@author: verjans Used in Verjans et al. (2019)

melt.effectiveT(self, iii)

trying to see what happens if we raise the temperature to an ‘effective temperature’ that is the temperature plus the latent heat from the liquid. Puts the firn above T_melt. Potential issue with this method: there might be effective diffusion of mass because we are now diffusing with volumes warmer than T_melt –> adjacent volumes that were dry might end up as having liquid.

melt.percolation_bucket(self, iii)

This is the bucket scheme that allows liquid water to persist in the firn. It includes consideration of irreducible liquid water content (LWC) an maximum LWC. Water that encounters a slab of a certain density (impermeable_rho) will not percolate through.

LWC is in volume (m^3), and since we are working in one dimension we assume that it is m^3/m^2.

@author: maxstev