Hi all,
An other question from my side.
I realized that once running in multidirectional spectral mode, ARTEMIS still uses the direction (and waveheight) defined in subroutine BORH. Not too suprising, but. Not defining it in this subroutine doesn't seem an option as the phase definition (which uses the defined angle) doesn't work anymore then (which is definitely necessary in the grid I'm using.
A solution might be to define THETAB as the angle calculated in dirale.f (Dale) which would then be used for the phase difference as well, but my problem is that I wouldn't know how to refer to this variable. Would anyone be able to help me with this?
Any other approach would be welcome as well.
Kind regards, Vera