open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I got an error after 9 minutes when calculating by coupling Telemac and Sisyphe.
my model has a trapezoid profile. The embankments are made of cement, only the sole in the main channel is mobable. ( image Cross section_Channel)





***************************
SUBROUTINE NOEROD
***************************
!

---

! RIGID BEDS POSITION
!

---

DO J=1,NPOIN

ZR(J) = ZF(J)

IF( Y(J) .GT. 0.3D0
& .AND. Y(J) .LT. 0.5D0 )THEN

ZR(J) =ZF(J) + (0.3D0-Y(J))

ELSEIF( Y(J) .GE. 0.5D0
& .AND. Y(J) .LE. 1.5D0 )THEN

ZR(J) =ZF(J) - 0.2D0

ELSEIF( Y(J) .GT. 1.5D0
& .AND. Y(J) .LT. 1.7D0 )THEN

ZR(J) =ZF(J) + (Y(J)-1.7D0)

ENDIF
ENDDO

////////////////////////////

In Sis_Steering_File:(the sohl material is defined as multigrain)
/

---

/ MULTI GRAIN
/

---

NUMBER OF SIZE-CLASSES OF BED MATERIAL = 3
SEDIMENT DIAMETERS = 0.00544D0; 0.00317D0; 0.0017D0
INITIAL FRACTION FOR PARTICULAR SIZE CLASS = 0.5D0; 0.3D0; 0.2D0

/

---

/ MULTI LAYER
/

---

NUMBER OF BED LOAD MODEL LAYERS = 3
CONSTANT ACTIVE LAYER THICKNESS = YES
ACTIVE LAYER THICKNESS = 0.20
SECONDARY CURRENTS = YES

////////////////////////////

The error in Listing.File:
================================================================
ITERATION 113800 TIME: 9 MN 29.0000 S ( 569.0000 S)
.
.
MAXIMAL EVOLUTION : 0.3159248E-04 NODE : 84191
MINIMAL EVOLUTION : -0.4070518E-04 NODE : 86027
TOTAL MAXIMAL EVOLUTION : 0.8143133E-01 NODE : 1070
TOTAL MINIMAL EVOLUTION : -0.1011259 NODE : 83933
ERROR IN LAYER CASE 1

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2
_____________
runcode::main:

////////////////////////////

I ask for help.

Many Thanks
Rima

Hello Rima,

Subroutine layer is a frequent cause of trouble... Some general remarks:

1. First, uncomment all checks the say "ACTIVATE BEFORE INVESTIGATING PROBLEMS IN LAYER...". They may give you more insight in the nature of your problem.
2. Then, try setting ZERO in layer.f to something smaller than 1.D-7. I have even used 1.D-4 on difficult situations.
3. Do you really need stratification? If not keep only one layer and set ACTIVE LAYER THICKNESS to something big. If again you definitely need it, start with one layer only and if you manage to complete the simulation without errors, then introduce the feature.

Best Regards,
Costas

Hi Costas,
Thank you very much for the support.This has finally worked
but I still have please a question about the results (Mean Diameter).



I have defined the Embankments as Rigid Bed in Subroutine Noerd.




Nevertheless, a classification for the bottom material on the embankments was shown. They do not have to be movable.
can you please explain this.

Thanks alot
Rima

Hello Rima,

This is what I get too in areas where non-erodible bed has been exposed. The logical nature of the problem is that there is no sediment there to calculate D50 for. The mathematical aspect is that D50 is calculated as SUM(D50*AVAIL for each class)/SUM(AVAIL for each class). If AVAIL for every class as well the SUM(AVAIL) is almost zero, you can see where this is going...
You can mask out (or set to zero) the areas where D50 is invalid.

Regards,
Costas