Welcome, Guest
Username: Password: Remember me

TOPIC: To record iteration result at specific node

To record iteration result at specific node 4 years 11 months ago #35112

  • amyrhulbakar
  • amyrhulbakar's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 46
  • Thank you received: 4
Dear users,

Is there any way to record the iteration result only at a specific node? I used a small time step, e.g. 0.1 sec. and need to observe changes at every step. Yet printing out the result for the whole domain produced a very huge result file. Please advise or could you direct me to any existing topic that had been discussed.

Many thanks,
Amyrhul
The administrator has disabled public write access.

To record iteration result at specific node 4 years 11 months ago #35121

  • amyrhulbakar
  • amyrhulbakar's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 46
  • Thank you received: 4
Dear users,

Further to my question, I found that by assigning with the VARIABLES TO BE PRINTED is one of the way to record the iteration result at every time step - based on the previous discussion here:

www.opentelemac.org/index.php/kunena/16-...-to-be-printed#15769

Based on the discussion, user has defined the location of section in SECTIONS INPUT FILE with the result in SECTIONS OUTPUT FILE. For my case, where should I define the location of a specific node that I want to extract the values from and where this result should be recorded?

Many thanks,
Amyrhul
The administrator has disabled public write access.

To record iteration result at specific node 4 years 10 months ago #35122

  • amyrhulbakar
  • amyrhulbakar's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 46
  • Thank you received: 4
Dear users,

The variable (e.g. water depth) was printed for every time step in the listing-printout log files as in the attachments. As the model was run in parallel, every processor printed values at the local level, i.e. I have 250,963 nodes at the global domain. Which subroutine(s) should I approach to have the variable printed at the global level? Further, how to only print the variable from a specific node?

Many thanks,
Amyrhul

File Attachment:

File Name: PE00039-00001.log
File Size: 1,532 KB


File Attachment:

File Name: PE00039-00002.log
File Size: 1,545 KB


File Attachment:

File Name: PE00039-00003.log
File Size: 1,524 KB
The administrator has disabled public write access.

To record iteration result at specific node 4 years 10 months ago #35123

  • amyrhulbakar
  • amyrhulbakar's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 46
  • Thank you received: 4

File Attachment:

File Name: PE00039-00001.txt
File Size: 1,532 KB


File Attachment:

File Name: PE00039-00002.txt
File Size: 1,545 KB


File Attachment:

File Name: PE00039-00003.txt
File Size: 1,524 KB
The administrator has disabled public write access.

To record iteration result at specific node 4 years 10 months ago #35192

  • amyrhulbakar
  • amyrhulbakar's Avatar
  • OFFLINE
  • Junior Boarder
  • Posts: 46
  • Thank you received: 4
Dear all,

Could anybody who has experience on this do comments based on my previous questions please. FYI, I am printing this information at the listing printout with specifying the VARIABLES TO BE PRINTED.

Questions:
1. Which subroutine(s) should I approach to have the variable printed at the global numbering?
2. How to only print the variable for a specific node?

Kind Regards,
Amyrhul
The administrator has disabled public write access.

To record iteration result at specific node 4 years 10 months ago #35193

  • pilou1253
  • pilou1253's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 584
  • Thank you received: 106
Hello!

You can have a look at what is done in the test case Rouse from telemac3d and especially in the utimp.f subroutine.

The part IF(NCSIZE.GT.1) allows you to select a local node number corresponding to global ID 46. Note that there is a subroutine called global_to_local_point.f in sources/utils/bief that also allows to switch between global and local node numbers.

To get your results for water depth, the best would be to print them in a formatted results data file, see corresponding keyword. To write in this file, simply add a line like:

WRITE(LU_file,*) AT,H%R(your_local_node_ID)


LU_file is the logical unit corresponding to the formatted result file you use. See more info in the .dico file and I guess in the user manual.
AT is the current time step in seconds

Good luck!

Best regards
PL
The administrator has disabled public write access.
The following user(s) said Thank You: amyrhulbakar
Moderators: pham

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