Hello,
I'm still testing the boundary conditions for the sediment flux and this time I set the LIEBOR=2 in the SISYPHE boundary condition file. It seems give better results than other types of BC (KENT and KSORT). I also checked the code and I found the following lines in the subroutine DIFFIN:
IF(CLT(K1).EQ.KLOG.OR.CLT(K2).EQ.KLOG) THEN
! SEGMENTS OF TYPE NEUMANN
MASKTR%ADR(NEU)%P%R(K1)=1.D0
So apparently it masks the boundary nodes which have code KLOG(2) as the NEUMANN boundary. But I still don't understand what exactly it does after this masking procedure. Does it mean the DT/DN=0 on the boundary if I use KLOG as the boundary condition for sediment flux?
Thanks in advance!
Qilong