Welcome, Guest
Username: Password: Remember me
  • Page:
  • 1
  • 2

TOPIC: Directional spectrum parameter

Directional spectrum parameter 11 years 9 months ago #8766

  • M.mire
  • M.mire's Avatar
Hello,

I've started to code the S function (of Goda). But I would like to have information regarding some variables... Where can I find the fortran file in which FREQ (discretized frequencies) and DIREC (discretized direction) tables are filled ?
Do you have a document in which all the fortran files are briefly explained ?

Regards,
Rémi
The administrator has disabled public write access.

Directional spectrum parameter 11 years 9 months ago #8768

  • Elodie
  • Elodie's Avatar
Hello,

FREQ and TETA are filled in the routine condiw.f.

I don't have a document in which all the fortran files are explained for the current version of Tomawac (or even for a recent version). I am sorry.

Regards,

Elodie
The administrator has disabled public write access.

Directional spectrum parameter 11 years 9 months ago #8779

  • M.mire
  • M.mire's Avatar
Thank you Elodie,

Ok, I cheked condiw.f file. And I have noticed that I was wrong in the last post, TETA is the discretized direction... So I'm wondering what is DIREC (the argument of the subroutine FSPRD3). Is it the same as TETA ? (in the comments of the fortran file DIREC is also the discretized direction).

Regards,

Rémi
The administrator has disabled public write access.

Directional spectrum parameter 11 years 9 months ago #8782

  • Elodie
  • Elodie's Avatar
Hello,

yes, DIREC in the subroutine fsprd3.f is the same as TETA in condiw.f or speini.f... You can see that easily if you look at the calls of the routine fsprd3.f in speini.f. The second argument is indeed TETA.

Regards,

Elodie
The administrator has disabled public write access.

Directional spectrum parameter 11 years 3 months ago #11345

  • M.mire
  • M.mire's Avatar
Hello,

It’s been a while that we’ve started this topic but I’ve finally had time to code the Goda method in order to have a spreading function changing with the frequencies. The results are quite satisfactory.
I’ve attached to this mail the fortran modifications and some results if you want to have look…

Regards,
Rémi

File Attachment:

File Name: h46a3999.f
File Size: 11 KB


h7222e8c.png


h956465b.png
The administrator has disabled public write access.

Directional spectrum parameter 11 years 2 months ago #11364

  • M.mire
  • M.mire's Avatar
Hi,

I made a small mistake in the limwac subroutine line 254:

ARGUM = ABS(COS(0.5D0*(DTETA)))and not ARGUM = ABS(COS(0.5D0*(DTETA-TETA1L)))

Sorry about that,

Regards,
Rémi

File Attachment:

File Name: hdadbcf7.f
File Size: 11 KB
The administrator has disabled public write access.
  • Page:
  • 1
  • 2
Moderators: tfouquet

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