Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
TELEMAC-2D
DEBSCE.f
Flat Threaded

TOPIC: DEBSCE.f

DEBSCE.f 7 years 4 months ago #24592

  • mourad
  • mourad's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 150
  • Thank you received: 6
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
The administrator has disabled public write access.
The following user(s) said Thank You: erwanC

DEBSCE.f 7 years 4 months ago #24593

  • mourad
  • mourad's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 150
  • Thank you received: 6
I just find the answer by myself. These modifications are needed to increase the number of sources up a thousand points.
The administrator has disabled public write access.
The following user(s) said Thank You: erwanC

DEBSCE.f 7 months 2 weeks ago #43264

  • Hsuan
  • Hsuan's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 16
可以詢問你後來是怎麼修正的嗎?
以下是我的debsce.f檔案的程式碼：
      IF(OKDEBSCE(I).AND.T3D_FILES(T3DVEF)%NAME(1:1).NE.' ') THEN
!
!       FCT WILL BE Q(1), Q(2), ETC, Q(99), DEPENDING ON I
        FCT='Q(       '
        IF(I.LT.10) THEN
          WRITE(FCT(3:3),FMT='(I1)') I
          FCT(4:4)=')'
		ELSEIF(I.LT.1000) THEN
		  WRITE(FCT(3:5),FMT='(I3)') I
		  FCT(6:6)=')'
        ELSE
          WRITE(LU,*) 'T3D_DEBSCE NOT PROGRAMMED'
          WRITE(LU,*) 'FOR MORE THAN 99 SOURCES'
          CALL PLANTE(1)
          STOP
        ENDIF
        CALL READ_FIC_SOURCES(T3D_DEBSCE,FCT,AT,T3D_FILES(T3DVEF)%LU,
     &                        INFOGR,OKDEBSCE(I))
!
      ENDIF

我只修改了這段
		ELSEIF(I.LT.1000) THEN
		  WRITE(FCT(3:5),FMT='(I3)') I
		  FCT(6:6)=')'

但依舊還是不能運行超過99個SOURCES
The administrator has disabled public write access.

DEBSCE.f 7 months 2 weeks ago #43265

  • jtravert
  • jtravert's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 36
  • Thank you received: 23
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
The administrator has disabled public write access.

DEBSCE.f 7 months 2 weeks ago #43266

  • Hsuan
  • Hsuan's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 16
Sorry, I forgot to change the language. Can this website edit my post directly?
The administrator has disabled public write access.

DEBSCE.f 7 months 1 week ago #43319

  • jtravert
  • jtravert's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 36
  • Thank you received: 23
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.
The administrator has disabled public write access.
Flat Threaded
Moderators: pham
Forum
System Modules
TELEMAC-2D
DEBSCE.f

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.