I want to set a non-erodible bottom for a part of the calculation domain (ex. Bed protection work etc.), however, the riverbed never evolves.
and what I did are below.
1. Create a private variable, NERO, in the geometry file. I used this variable to represent the thickness of the initial layers. In riverbeds where erosion does not occur at all, NERO=0, while in other cases, NERO=1 was assigned. Attached is a snapshot of the geometry file defining the variables.
2. Modify user_bed_init.f file as follow;
DO IPOIN=1,NPOIN
ESTRATUM(1,IPOIN) = PRIVE%ADR(1)%P%R(IPOIN)
ENDDO
3. Declare it by the statement of NAMES OF PRIVATE VARIABLES in the GAIA steering file.
I think the private variable is not correctly stored in ESTRATUM.
If anyone has had a similar experience, could you please share any solutions you found?