Hello Charlotta
Finally i succeeded to obtain your results. You will find attached 2 results of the HM0 at time 550800 that shows that we get the value of 4.75 under some conditions :
- firstly : i didn't explain it well in my first message. In order to well discretized the set of frequencies around 0.118 for a number of frequencies of 29 you should take a minimal frequency of 0.06
- Then the main problem was the choice around the BOUNDARY ANGULAR DISTRIBUTION FUNCTION. For a value of 3, i found the good HM0 value of 4.75. If you take 2 like in your case then change the value of the BOUNDARY DIRECTIONAL SPREAD to 0.5. The default value is 2, and it is not good for this law as the integrale over angles is not 1 anymore for that value) If you do so you will get the results i attach.
- That was not the problem but the ";" is not the standard separator for fortran (read(*)), it might work with your compiler but in general it won't , the good separator for your txt file is ' '.
- Finally, just a detail but you don't have to do the loop on k each time you want a frequency for a new time step I join the fortran user file I modified a little bit to optimize it.
To be sure that there was no interaction with anything else I commented all the source terms.
I join
- the image of the 2 results at 550800 for s=0.5 and s=2
- The fortran user modified
- the steering file modified.
regards
T.