I have a very basic Sisyphe question
I have been running Sisyphe and Telemac coupled mode with a specified SSC boundary concentration that works perfectly using the keyword
CONCENTRATION PER CLASS AT BOUNDARIES = 0.001
I am now trying to run sisyphe in uncoupled mode. However it now no longer applies the specified boundary value - it looks like it is applying some sort of free condition e.g. zero gradient
I have also tried specifying the condition in the boundary file with the following at all boundary points
5 6 6 0.000 0.000 0.000 0.000 5 0.001 0.000 0.000 231937 1
I've also tried setting the value in conlit.f see below but that causes a crash
!
!
IF(SUSP) THEN
!
DO K=1,NPTFR
!
! SO FAR LICBOR=LIEBOR (WITH KADH CHANGED INTO KLOG, SEE ABOVE,
! BUT CAN BE CHANGED)
!
LICBOR%I(K) = 0.001 !LIEBOR%I(K)
!
I'd be very grateful if anyone could help with this
Thanks
John