open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi there,

I have implementing a vertical structure in a domain, and have recently parallelised my code. I think I need to better understand what the 'aire' variable is for, because when the domain is split into multiple regions, there are some in which there are no vertical structures, so the value of area remains zero and then causes the code to return an error and stop (because it would otherwise divide by zero).

Could someone please explain quite why we must divide by area? I want to know whether I could set it to default divide by one if the area is zero so I stop getting the exception.

Thanks,
Chris

Hello,

Yes, you are quite right, the treatment of parallelism is not done in subroutine dragfo in version 6.1. I'll do it in next version.

We apply a stress and first deal with forces, hence we must divide after by the area on which the force is exerted, so that when integrated over the area it gives back the force. You just need to :

Add :

DOUBLE PRECISION P_DSUM
EXTERNAL P_DSUM

at the beginning of the subroutine dragfo.

Before the test IF(AIRE.GT.1.D-6)... add the instruction :

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

I hope it works, with best regards,

Jean-Michel Hervouet