TELEMAC-MASCARET Forum: Using the DRAGFO subroutine to model TEC (2/2)

Welcome, Guest

Page:
1
2

TOPIC: Using the DRAGFO subroutine to model TEC

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31659

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	What is the error message?
	Christophe The administrator has disabled public write access.

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31663

r.clayton	Hi, When I try to compute the code, the following command window opens (see attached) and it seems like it is computing then it just all closes and nothing happens. Thanks, Attachments: Errormessage.png (65KB)
	The administrator has disabled public write access.

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31666

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	You should open a cmd windows manually, go into your working directory and run telemac manually in order to have the full listing This will help us to understand where the problem is regards
	Christophe The administrator has disabled public write access.

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31676

r.clayton	Hi, Please find full error message attached...any suggestions of what I need to change? Thanks, Attachments: Errormessage_v2.png (107KB)
	The administrator has disabled public write access.

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31682

Leballeur OFFLINE openTELEMAC Guru Posts: 434 Thank you received: 163	Hello, It seems that you're using a version of drago.f subroutine more recent that your telemac installation. As v6p3 version is quite old, I advise you to upgrade your telemac installation. Regards, Laurent
	The administrator has disabled public write access.

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31727

s_gegenleithner
OFFLINE
Fresh Boarder
Posts: 18
Thank you received: 1

Hello,

I have a question regarding the use of dragfo.f in parallel, please correct me if I am wrong.

To make the dragfo.f subroutine work in TELEMAC, the following code was added:

IF(NCSIZE.GT.0) AIRE=P_DSUM(AIRE)

However, taking a look at:

      AIRE=0.D0
      DO I=1,BIEF_NBPTS(11,MESH)
!
        IF(INPOLY(X(I),Y(I),XSOM,YSOM,NSOM)) THEN
          UNORM = SQRT(UN%R(I)**2+VN%R(I)**2)
          FUDRAG%R(I) =  - 0.5D0 * N * DIAM * CD * UNORM
          FVDRAG%R(I) =  - 0.5D0 * N * DIAM * CD * UNORM
          AIRE = AIRE + T1%R(I)
        ENDIF
!
      ENDDO

rises some questions for me. The Mesh object, at least in my case represents one part of the decomposed mesh, not the full one. Now checking if the coordinates are inside of the polygon allows to check only for this specific part. In my case the polygon is outside of this specific decomposed domain, and thus TELEMAC never jumps inside of the IF statement.

Any hint is appreciated.

Cheers,
Sebastian

The administrator has disabled public write access.

Using the DRAGFO subroutine to model TEC 6 years 1 month ago #31736

c.coulet
OFFLINE
Moderator
Posts: 3722
Thank you received: 1031

Hi Sebastian
if the polygon doesn't match with the subdomain, then AIRE (area) stay to 0.
But in your model, there is at least one or more subdomain which intersect the polygon, then those subdomain provide some positive values to AIRE which are globally summed with the call to P_DSUM...
This is required if a node inside the polygon belongs to more than one subdomain? Otherwise the area will be only the area of element inside the subdomain, not the real area around the node...

Hope this helps

Christophe

The administrator has disabled public write access.

The following user(s) said Thank You: s_gegenleithner

Using the DRAGFO subroutine to model TEC 6 years 3 weeks ago #31835

s_gegenleithner OFFLINE Fresh Boarder Posts: 18 Thank you received: 1	Hello Christophe, Thank you for your reply. I implemented a dragforce for several subdomains - it works fine. I got confused because I didn't know that only processor 0 contributes to the listing printout. Sorry for the confusion, Sebastian
	The administrator has disabled public write access.