Welcome, Guest
Username: Password: Remember me

TOPIC: Specification of input-current

Specification of input-current 6 years 7 months ago #29508

  • Papadimitriou
  • Papadimitriou's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 80
  • Thank you received: 3
Dear all, I would like to examine the wave-current intaractions specified in TOMAWAC.
I was wondering if there exists a modified version of the COUUTI.f subroutine,or some general directions to modify it, because I might not be able to use the .slf file standard for the input current velocities.

Kind regards,
Andreas
The administrator has disabled public write access.

Specification of input-current 6 years 7 months ago #29511

  • Leballeur
  • Leballeur's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 434
  • Thank you received: 163
Hi Andreas,

What is exactly your issue with the use of the "BINARY CURRENTS FILE" ?
The subroutine couuti.f is used to read the current data from a user-defined format, so there are a lot possibilities to program into.

Regards,
Laurent
The administrator has disabled public write access.

Specification of input-current 6 years 7 months ago #29513

  • Papadimitriou
  • Papadimitriou's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 80
  • Thank you received: 3
I was hoping to insert the current velocities at cetrain points of the mesh like this for example: (X,Y,U,V).I was considering the use of the Coouti.f subroutine kinda like the way meteo.f interpolates wind field in the whole mesh.
I was wondering if there is an example of a modified routine that performs a similar role for TOMAWAC or alternatively, if there is a possibility to set this input at a Telemac-2d simulation and obtain the required binary .slf currents file through the hydrodynamic model (like meteo.f does for wind)

Thank you very much for the assistance
Kind regards,
Andreas
The administrator has disabled public write access.

Specification of input-current 6 years 7 months ago #29514

  • Leballeur
  • Leballeur's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 434
  • Thank you received: 163
The best way is to first run a Telemac-2d simulation on the same grid and then you can give the results file to Tomawac, without programming anything.

Otherwise you can code what you want directly in couuti.f by filling the arrays UC1,UC2,VC1,VC2 for the NPOIN of your mesh at each time step AT. You can take inspiration from meteo.f, but I think spatial interpolations are not done... If you plan to do it by yourself, you can use the subroutine fasp.f for instance, but be aware that your resulting currents may no longer have physical coherence.

Regards,
Laurent
The administrator has disabled public write access.
The following user(s) said Thank You: Papadimitriou

Specification of input-current 6 years 6 months ago #30058

  • Papadimitriou
  • Papadimitriou's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 80
  • Thank you received: 3
Dear all,

I am trying to reproduce the example in the Tomawac test case "opposing current" using the couuti.f subroutine for user specified current formats.

I have exported the coordinates and velocities to a text file "LAI.txt" and using the Meteo.f subroutine as blueprint I have tried to read it using the couuti.f

however i have two problems

1) How can i declare the KP1BOR variable in my couuti.f subroutine. Even if i use USE BIEF, it still says it is a dummy variable
2) If I set the variable kP1BOR in the fasp.f subroutine as NBOR in order for the case to run, I get velocities UV equal to zero

I am attaching my files, in hope that you can assist me

Kind regards,
Andreas
The administrator has disabled public write access.

Specification of input-current 6 years 6 months ago #30059

  • Papadimitriou
  • Papadimitriou's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 80
  • Thank you received: 3
Actually after some cleanup, I am receiving a segmentation fault.
I have checked and I think i save all the temporal arays from my currents file.
Here are the modified versions of the files.

Kind regards,
Andreas
The administrator has disabled public write access.

Specification of input-current 6 years 6 months ago #30066

  • tfouquet
  • tfouquet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 294
  • Thank you received: 112
Hello

I modified your user fortran and this one seems to work

regards
Attachments:
The administrator has disabled public write access.
The following user(s) said Thank You: Papadimitriou

Specification of input-current 6 years 6 months ago #30070

  • Papadimitriou
  • Papadimitriou's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 80
  • Thank you received: 3
Thank you very much Mr. Fouquet you have been as helpful as always.

For future reference I took inspiration from meteo.f and modified couuti.f to read spatially and temporally varying currents from a formatted file and perform the interpolation on the mesh using fasp.f

I am also attaching an example file, which corresponds to the Lai experiment (1989) test case

Best regards,
Andreas
Attachments:
The administrator has disabled public write access.
Moderators: tfouquet

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