open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi Fezai,

For cohesive sediments, please use init_compo_coh
Cath

Dear Cath,

Thank you for your prompt answer.

Actually I'm modeling non cohesive sediment transport. I told the steering file to show me D50 as graphic printout to check the distribution is correctly read but I'm getting the default values.

These are the code lines I implemented in init_compo subroutine :

CHARACTER(LEN=150):: SedPath
!
!

---

!
SedPath='...(dots replace the full path)\sim\granulo\'
!
OPEN(4444, FILE=TRIM(SedPath) // 'D50.txt', FORM='FORMATTED')
DO J=1, NPOIN
NCOUCHES(J) = 1
READ(4444, *) AVAIL(J,1,1)
ENDDO
CLOSE(4444)
OPEN(555, FILE=TRIM(SedPath)//'Log.txt', FORM='FORMATTED')
WRITE(555, 100) 'test d50'
CLOSE(555)
100 FORMAT(A)

Regards,
FEZAI

Well,

It needs some head scratching after following the "River-arc" example to find the solution. As it wasn't well explained in the documentation.

All the best,
FEZAI

Hello,

I'm working on a hydro-sedimentary model and I would like to spatialize the grain distribution because I work on a river with very different sediments/grain sizes (coarse sediments upstream and fine sediments downstream). I would like to modify only the minor river bed. If I do this have I to keep the keyword fraction for particular size class ?
I know that I have to modify the init_compo.f subroutine but I'm a beginner in fortran so can you give me some example please ?

Thanks,
Claire