open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

No, I think it is OK now. You could check it by putting the time into power and see if it pops up in the results file. Just write for every point M power(M)=AT.

If it is OK then your problem is not the output and you should check your way of computing the power: in your subroutine dragfo I do not see how dragfo_called is initialised to 0.

You can declare it in this way:

INTEGER dragfo_called
DATA dragfo_called/0/

and it will be set to 0 at the beginning and saved from one call to the other, so it will be a real counter. You can add also a SAVE after all the declarations, so that all declared local variables are saved.

Another thing is that I do not know if you can combine two commands on the same line, separated by a ";", it is strange that your compiler does not complain.

Regards,

Jean-Michel Hervouet

Hi JMH,

Thanks very much for your help, the code does seem to be working to an extent. The issue now is that the power is being calculated in the wrong region of the channel, as shown by attachment 1, image 1. Instead of being calculated at the start of the channel at each node, it should be at the nodes in the square region in the center of the channel. This is defined in the dragfo file through the use of a do loop: do m=1,n_inpoly.

I find it strange that it is not doing this because the drag force calculation seems to be correct in that when I simulate flow velocity in Blue Kenue, flow is affected in the correct region (see attachment, image 2). However, by making the drag force equal to PRIVE%ADR(1)%P%R(m) and removing the calculation for power seems to give the same problem, in that when displaying variable 23 which is now the drag force in Blue Kenue the only nodes that have drag force that varies with time are the ones at the start of the channel, not in the central region. (So it looks the same as image 1).

Could this be because there is a different file that is manipulating which nodes are used? I have attached the dragfo file again.

Thanks

Danny

Apologies, its seems I cannot attach a docx fiel, here is attachment 1 in its correct format

Hello,

As I do not see what is the structure turbine_array, I am not sure to understand what you do, as a matter of fact I do not understand :

ta(1)%box_x = (/4750.00,4750.00,5250.00,5250.00/)

Which seems to code for :

ta(1)%box_x = 4750.D0
ta(2)%box_x = 4750.00
ta(3)%box_x = 5250.00
ta(4)%box_x = 5250.00 and same thing with box_y...

Otherwise you do not fill the values of ta for J=2 to 4 ?

Either this is wrong or this is a possibility of Fortran 90 that I am not aware of (or it works in C but maybe not in Fortran...).

With best regards,

Jean-Michel Hervouet