open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi
The velocities are small on the X-axis direction but stronger on Y-axis...
It's also the beginning of the simulation so maybe currents needs time to be well established...

Hope this helps

Thanks Christophe!

However, I have double-checked other files, the magnitude is also small..

There is another point I am confusing is that I got 8 files separately. So I want wonder why? were these files generated by the number of processors or layers?

Thanks you so much.

Best Regards
Huy

Hi
As you run in parallel, there is 1 formatted file per processor...
And if you put the node number directly, this mean you have the value of the node in the subdomain... If the node doesn't exist, you will crash!
You should modify your program to only print the node with the right global number...

Regards

Thank you a lot

So do you think it is possibly linked with the error I got below this post while running a coupled model?

www.opentelemac.org/index.php/kunena/12-...allel?start=30#25416

Sorry for bothering!

Best Regards
HUY

Hi
Yes it's possible, it depends if you hardcoded a node number somewhere in a fortran file...

Regards

Thanks Christophe!

By the way, could you guide me how to print the node with the right global number? I have followed Pilou's previous post:

IF(NCSIZE.GT.1) THEN
N = 1000
NLOC = GLOBAL_TO_LOCAL_POINT(N,MESH3D)
IF(NLOC.NE.0) THEN
......
ENDIF
ENDIF

This works for one point, but I actually want 3 points to be extracted at the same time (e.g. 8118, 8496 and 7935). So how could I add/ modify the code? Could I slip these 3 points in three different files?

Thanks a lot in advance.

Best Regards
Huy

hi, you can use these:
IF(GLOBAL_TO_LOCAL_POINT(your global ID,MESH3D).NE.0) THEN
open(unit=123,file='..\TEMP1_89506.data',Access = 'append')
WRITE(123,*) AT,your global ID,
& TA%ADR(1)%P%R(GLOBAL_TO_LOCAL_POINT(your global ID,MESH3D)),
END IF

Then you can have 1 file for each node. And the processors which do not have your ID will not generate files. So you will get 3 files instead of 3*# of processors

Good luck

Yuan

Thanks Yuan,

To tell you the truth, I am not good at coding

What does "open(unit=123,file='..\TEMP1_89506.data',Access = 'append')" mean?

As I mentioned I want to extract three points (8118, 8496 and 7935) in the mesh.. what could I write exactly?

Best Regards
Huy

Hi, Huy,

I am not that good at neither;)

I stored the output of the nodes separately. To get output from 3nodes, I generated 3 files. Each file, you need to give it an ID. In the above example case, the file ID is 123, name is TEMP1_89506.data, and at each time step, you want to open it and write the results into it, so we need to open it, and write the results. If you only need 3 point, it is easier to only have 3 files, so you can loop the codes. In this way, you can only ask the precessor who has your point t write out the result, so your rest processors will no generate files.This can save you some space.

Best regards
Yuan

Hi all,

Thanks a lot for your discussion. Now I want to export time series of 2D u, v, and tracer at each node from a Telemac 3D project using a scope file, would you please tell me how to acheive this aim? Thanks!