VIC icon indicating copy to clipboard operation
VIC copied to clipboard

Is snow_melt.c handeling 'RefreezeEnergy' correctly?

Open jhamman opened this issue 11 years ago • 1 comments

In these lines in snow_melt.c:

    /* Refreeze liquid water in the pack.
       variable 'RefreezeEnergy' is the heat released to the snow pack
       if all liquid water were refrozen.
       if RefreezeEnergy < PackCC then all water IS refrozen
       PackCC always <=0.0

       WORK IN PROGRESS: This energy is NOT added to MeltEnergy, since this does
       not involve energy transported to the pixel.  Instead heat from the snow
       pack is used to refreeze water */

    snow->pack_water += melt[0]; /* add surface layer outflow to pack
                                        liquid water*/
    PackRefreezeEnergy = snow->pack_water * Lf * RHO_W;

Does anyone think this is true or is this just an old comment?

jhamman avatar Nov 23 '14 16:11 jhamman

This is also echoed in snow_intercept.c:

    /*** FIX THIS ***/
    *MeltEnergy = RefreezeEnergy / Dt;

jhamman avatar Nov 23 '14 16:11 jhamman