Skip to content

Cbal v1 - new canopy snow routine#469

Closed
acebulsk wants to merge 120 commits intosrlabUsask:masterfrom
acebulsk:cbal_v1
Closed

Cbal v1 - new canopy snow routine#469
acebulsk wants to merge 120 commits intosrlabUsask:masterfrom
acebulsk:cbal_v1

Conversation

@acebulsk
Copy link
Collaborator

First go at implementing new canopy snow energy and mass budgets as new CRHM module. Caught some bugs in other modules / files along the way which should be assessed. The main changes are the addition of two new modules to represent the initial loading of snow in the canopy:

CanopyVectorBased which follows developments from Cebulski & Pomeroy (2025b)

the second module is:

CanSnobal which follows developments from Cebulski & Pomeroy (in prep.)

together these new modules replace the canopy snow routines previously incorporated in CanopyClearingGap. The original radiation routines in CanopyClearingGap were copied over to the CanopyVectorBased module which feeds into the surface snowpack routines so CanopyClearingGap is no longer needed.

Note: there are several commits that are non-relevant as initial attempts to this were to revive the PSP module from Parv et al., (2000) however I could not get this working so ended up creating something new.

…pears to work using new files w old algorithm now to modify.
…x canopy snow load for sublimation to constant
…ding switch. now we are getting exactly the same as before adding the switches
…. Originally was immediately segfaulting due to some coding bugs.
…realistic canopy and ice sphere temps but still need to validate and plot the temps and energy terms
…d from canopy module. also fixed infinite loop error.
… is different in PSPnew. changed wind reference height to work with fortress instruments where wind is at mid canopy
…ome variables to global, and some internal params to global
added Switch to change between cionco or log wind profiles

changed Sicart 2006 cloudy longwave to use last computed value when QdflatE == 0

Adjusted energy bal check
… profile currently only used for trajectory calculation, should also use this for unloading
also changed to pull ground temperature form Canopy module
acebulsk added 28 commits March 18, 2025 16:40
…le to use cionco wind profile for IP and ablation
…is still entered if no snow but still liquid
…r non canopy hrus. also casted switches to ints to avoid issues with parameter values which are longs when compared against the specified integer values
@acebulsk
Copy link
Collaborator Author

Closing after more recent PR #471

@acebulsk acebulsk closed this Oct 24, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant