Adding Brooks-Corey SWRC Model – Issues with Defining Retention Curve Constants

[Hosseinan9675]

Hello Anura3D community,

I'm working on adding the Brooks-Corey Soil Water Retention Curve (SWRC) model to the Anura3D code, and I've encountered a few issues.

In the current codebase, I can see that the Van Genuchten and Linear SWRC models are implemented, but I couldn't find any explicit declarations for the retention curve identifiers, such as SWRC_VANGENUCHTEN or SWRC_LINEAR. These constants seem to be used within the conditional statements to differentiate between different SWRC models, but I don't see where they are defined in the source code.

To add the Brooks-Corey model, I have written a subroutine DynUpdateParticleDegreeOfSaturationBrooksCorey and integrated it into the main logic like so:

" else if (MatParams(ISet)%RetentionCurve == SWRC_BROOKSCOREY) then
call DynUpdateParticleDegreeOfSaturationBrooksCorey(ParticleIndex, ISet, Sr)
"
However, I'm getting the following errors:

Error #6303: The assignment operation or the binary expression operation is invalid for the data types of the two operands. [SWRC_BROOKSCOREY]
Error #6404: This name does not have a type and must have an explicit type. [SWRC_BROOKSCOREY]

Without any declarations for the existing SWRC models (like SWRC_VANGENUCHTEN), I'm unsure how to proceed with adding the Brooks-Corey model. Could someone clarify where these constants are typically defined in the Anura3D code, or if there's another approach to integrate a new SWRC model?

Any guidance on how to declare SWRC_BROOKSCOREY correctly or where to locate the declarations for the other SWRC models would be greatly appreciated!

Thank you in advance!

Best regards,
Hossein