open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello users and moderators,
The issue is to increase the number of sources (more than 99 points)

In the debsce.f what about the following :
1/ FCT(3:3) , FCT(4:4) → I know that FCT take values of sources discharges Q(1) , Q(2),....but why 3:3 and 4:4 !!!
2/ In fortran FMT('I2') means character string defined by an integer of 2 digits → right?

If I increase the number of sources (more than 99 as this subroutine is programmed) HAVE I to modify debsce.f as follow :

ORIGINAL :
IF(OKDEBSCE(I).AND.T2D_FILES(T2DVEF)%NAME(1:1).NE.' ') THEN
FCT='Q( '
IF(I.LT.10) THEN
WRITE(FCT(3:3),FMT='(I1)') I
FCT(4:4)=')'
ELSEIF(I.LT.100) THEN
WRITE(FCT(3:4),FMT='(I2)') I
FCT(5:5)=')'
ELSE
WRITE(LU,*) 'DEBSCE NOT PROGRAMMED FOR MORE THAN 99 SOURCES'
CALL PLANTE(1)
STOP
ENDIF

MODIFIED :

IF(OKDEBSCE(I).AND.T2D_FILES(T2DVEF)%NAME(1:1).NE.' ') THEN
FCT='Q( '
IF(I.LT.10) THEN
WRITE(FCT(3:3),FMT='(I1)') I
FCT(4:4)=')'
ELSEIF(I.LT.100) THEN
WRITE(FCT(3:4),FMT='(I2)') I
FCT(5:5)=')'
!

---

modification starts here
ELSEIF(I.LT.1000) THEN
WRITE(FCT(3:5),FMT='(I3)') I
FCT(6:6)=')'
!

---

end of modification
ELSE
WRITE(LU,*) 'DEBSCE NOT PROGRAMMED FOR MORE THAN 99 SOURCES'
CALL PLANTE(1)
STOP
ENDIF
Any suggestion, remark or help will be appreciated
with regards

I just find the answer by myself. These modifications are needed to increase the number of sources up a thousand points.

可以詢問你後來是怎麼修正的嗎?
以下是我的debsce.f檔案的程式碼：
我只修改了這段
但依舊還是不能運行超過99個SOURCES

Hello Hsuan,

If you have a problem with the number of sources you may need to increase the number of columns read in the source file.
The procedure is described here: www.opentelemac.org/index.php/kunena/16-...the-number-set#42549


PS: Please edit your post to English so everyone can understand.

Best regards,
Jean-Paul

Sorry, I forgot to change the language. Can this website edit my post directly?

No, it might not be possible to edit it after some time. I suggest you try what was described in my last post. And if it is not working, you can try to give us more details on your problem with the error message. Or you can try to reproduce at first the code from others in the forum treating more than 100 sources.