-
Notifications
You must be signed in to change notification settings - Fork 105
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Implementation of HDS in SUMMA #568
base: develop
Are you sure you want to change the base?
Conversation
@martynpclark, please review this. I can't add you as a reviewer to this PR because I don't have permission. |
I can include the test step I have been experimenting with (if that's helpful). |
real(rkind) :: depressionVol ! depression volume [m3] | ||
!***************************************** | ||
! initialize HDS pothole storage variables | ||
! GRU level |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Consider updating to new refactored code standards that are currently being developed
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@@ -391,7 +397,51 @@ subroutine read_icond(iconFile, & ! intent(in): name of | |||
! -------------------------------------------------------------------------------------------------------- | |||
! (2) now get the basin variable(s) | |||
! -------------------------------------------------------------------------------------------------------- | |||
|
|||
!******************************** |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Can this code be generalized?
This is the initial implementation of the Hysteretic Depressional Storage (HDS) algorithm for prairie potholes in the SUMMA model.
The model has been tested under different I/O configurations and works as expected on a hypothetical test case.
To activate the HDS model, the following has to be done:
a model decision must be defined to activate HDS as:
The following HDS parameters must be defined at the GRU (basin) level in the
zParamTrial.nc
file (if the parameters are spatially variable) or in theBasinParamInfo.txt
if they are spatially constant. The following is an example of dummy parameters value.pondVolFrac
at the GRU level in thezInitialCond.nc