def heat_stor(cp, spm, tdif):
    """
    Calculate the heat storage
    Args:
        cp: specific heat of layer (J/kg K)
        spm: layer specific mass (kg/m^2)
        tdif: temperature change (K)
    """

    return cp * spm * tdif