Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
TELEMAC-2D
Using the IKLE connectivity table
Flat Threaded
  • Page:
  • 1
  • 2

TOPIC: Using the IKLE connectivity table

Using the IKLE connectivity table 9 months 3 weeks ago #42944

  • MartinPHYS
  • MartinPHYS's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
  • Thank you received: 1
Hi Christophe,

Yes, I have tried using

A = MESH%IKLE(IELEM,2)
B = MESH%IKLE(IELEM,3)
C = MESH%IKLE(IELEM,1)

And also I have tried using,

A = MESH%IKLE%I(IELEM,2)
B = MESH%IKLE%I(IELEM,3)
C = MESH%IKLE%I(IELEM,1)

as well as,

A = IKLE%I(IELEM,2)
B = IKLE%I(IELEM,3)
C = IKLE%I(IELEM,1)

The first gives me the error of,

126 | A = MESH%IKLE(IELEM,2)
| 1
Error: Unclassifiable statement at (1)

and the second and third both give me an error of,


126 | A = MESH%IKLE%I(IELEM,2)
| 1
Error: Rank mismatch in array reference at (1) (2/1)



126 | A = IKLE%I(IELEM,2)
| 1
Error: Rank mismatch in array reference at (1) (2/1)


I will attach my full fortran subroutine file in case this helps. Maybe there is another issue that I am not catching.
The administrator has disabled public write access.

Using the IKLE connectivity table 9 months 3 weeks ago #42945

  • MartinPHYS
  • MartinPHYS's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 24
  • Thank you received: 1
Here is my full fortran file
Attachments:
The administrator has disabled public write access.
The following user(s) said Thank You: jtravert

Using the IKLE connectivity table 9 months 2 weeks ago #42953

  • jtravert
  • jtravert's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 36
  • Thank you received: 23
Hello Martin,

Thank you for the file, it is easier to debug that way.
I do not remember everything about the 'theory' behind the IKLE connectivity table.
BUt here is a solution to your problem to get the nodes of all the elements following your implementation.

I tested it with your code on the 'break' example from TELEMAC-2D commenting all the things below, and it seems to work (with A, B and C the nodes number of the current element IELEM in the loop).


Here are the modifications (I just erased the following in the loop for clarity):

        DO IELEM = 1,NELEM
          A = IKLE%I(IELEM)
          B = IKLE%I(NELEM+IELEM)
          C = IKLE%I(2*NELEM+IELEM)
        ENDDO

Best regards,
Jean-Paul
The administrator has disabled public write access.
Flat Threaded
  • Page:
  • 1
  • 2
Moderators: pham
Forum
System Modules
TELEMAC-2D
Using the IKLE connectivity table

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