I have been trying to do some changes in BED1_SUSPENSION_ERODE and noticed that for pure sand case, erosion is calculated as erosion flux*settling velocity.
FLUER_PUR_SAND(ISAND,IPOIN) = CAE_ILAY(ISAND,ILAYER,IPOIN)*XWC(NUM_ISAND_ICLA(ISAND))
But this is the settling velocity (XWC) provided in the GAIA steering file by the user. I am wondering what is the purpose of the subroutine COMPUTE_SETTLING_VEL then?
Could someone guide me to use settling velocity (WCHU) from COMPUTE_SETTLING_VEL subroutine in the BED1_SUSPENSION_ERODE subroutine (so basically replace XWC with WCHU)?
Thanks a lot
Nithin