Skip to content
This repository has been archived by the owner on Oct 21, 2024. It is now read-only.

Commit

Permalink
Source code of v0.1 release of HOS-ocean
Browse files Browse the repository at this point in the history 
Signed-off-by: Guillaume Ducrozet <[email protected]>
  • Loading branch information
@gducrozet
gducrozet committed Nov 19, 2014
0 parents commit 96267e6
Show file tree
Hide file tree
Showing 34 changed files with 16,774 additions and 0 deletions.
10 changes: 10 additions & 0 deletions .gitignore
Original file line number Diff line number Diff line change
@@ -0,0 +1,10 @@
# ignore obj directory
obj/
# ignore bin folder except input and waverf coef
bin/*
!bin/input_HOS.dat
!bin/waverf*
# ignore backup files + Mac specific files
*~
sources/*/*~
.DS_Store
674 changes: 674 additions & 0 deletions LICENSE.txt

Large diffs are not rendered by default.

59 changes: 59 additions & 0 deletions README.txt
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
HOS-ocean

=========

High-Order Spectral method for oceanic simulations

This readme file describes the different cases that might be computed with the periodic HOS model
and gives instructions to set the numerical parameters n1 to p2

Setting the value of integers n1, n2, M, p1 and p2 in variables_3D.f90

For a 2D simulation,
Compile with n2=1 AND p2=1 to adjust the memory allocation to minimum

If partial dealiasing is used, compile with p1 set to maximal required value (total dealiasing is obtained
with p1=M but it can be reduced if p1 is further set to a value below M

For a 3D simulation,
Compile with n2\=1 AND p2 set to required value

If partial dealiasing is used in x-direction, compile with p1 set to maximal required value (total dealiasing is obtained
with p1=M but it can be reduced if p1 is further set to a value below M

If partial dealiasing is used in y-direction, compile with p2 set to maximal required value (total dealiasing is obtained
with p2=M but it can be reduced if p2 is further set to a value below M

Setting the value of integer i_case.

i_case = 1
starts from rest
i_case = 2 and 21
starts with a natural mode, either progressive (2) or stationary (21).
The number of the mode, its amplitude and phase have to be chosen in the module initial_condition.f90.
The x-component may be positive or negative for the progressive case.
Should be positive for the stationary case.
The y-component may be positive or negative for both cases.
In the progressive case, the potential on the free surface is set accordingly to linear theory.
In the stationary case, it is set to zero (no velocity at t=0).
i_case = 3 and 31 and 32
Irregular directional sea-state (linear) from spectrum specified with
input file parameters (3)
WAVEWATCH III® spectrum file (31)
previous HOS-ocean simulation (32)
i_case = 8...
Rienecker and Fenton test
this means that xlen and ylen are interpreted as integer and represent the number of wavelength
in x- and y- direction respectively
i_case = 81
steepness 0.1
i_case = 82
steepness 0.2
i_case = 83
steepness 0.3
i_case = 84
steepness 0.4
i_case = 89
steepness 0.09

Further details about input file, output of the code... may be find at the Wiki page of HOS-ocean project: https://github.com/LHEEA/HOS-ocean/wiki
26 changes: 26 additions & 0 deletions bin/input_HOS.dat
Original file line number Diff line number Diff line change
@@ -0,0 +1,26 @@
Choice of computed case :: i_case :: 3
--- Geometry of the horizontal domain
Length in x-direction :: xlen :: 80.0
Length in y-direction :: ylen :: 20.0
--- Time stuff
Duration of the simulation :: T_stop :: 5.0
Sampling frequency (output) :: f_out :: 5.0
Tolerance of RK scheme :: toler :: 1.0e-7
Dommermuth initialisation :: n :: 4
Dommermuth initialisation :: Ta :: 0.0
--- Physical dimensional parameters
Gravity :: grav :: 9.81
Water depth :: depth :: 35.
--- Irregular waves (i_case=3)
Peak period in s :: Tp_real :: 10.0
Significant wave height in m :: Hs_real :: 4.5
JONSWAP Spectrum :: gamma :: 3.3
Directionality (Dysthe) :: beta :: 0.785398
--- Output files
Tecplot version :: tecplot :: 11
Output: 1-dim. ; 0-nondim. :: i_out_dim :: 1
3d free surface quantities :: i_3d :: 0
3d modes :: i_a_3d :: 1
2d free surface, center line :: i_2d :: 0
Wave probes in domain :: i_prob :: 0
Swense output 1='yes',0='no' :: i_sw :: 0
155 changes: 155 additions & 0 deletions bin/waverf_L628_inf_ka009_N50_100.cof
Original file line number Diff line number Diff line change
@@ -0,0 +1,155 @@
6.28318530717959 0.180000000000000 1.00000000000000
6.25778973731953 1.00405823316635 0.000000000000000E+000
0.000000000000000E+000 50 100
0 -1.00405823316635
1 8.962625141103669E-002
2 3.289225957487150E-005
3 5.165858754824442E-007
4 8.387440315440270E-009
5 1.293778976869871E-010
6 1.919583610661396E-012
7 2.822834311237738E-014
8 4.221893575099491E-016
9 5.747892737504613E-018
10 5.522926235216216E-019
11 1.877928647612090E-018
12 -3.240722812129962E-018
13 -1.417506518532076E-019
14 1.860508384829791E-018
15 -2.675875654203579E-019
16 -1.331555314774992E-018
17 1.371985747084411E-018
18 1.102662446845281E-019
19 -1.086784557417955E-018
20 1.371251422913937E-018
21 -1.886110862454749E-018
22 2.238561281591652E-018
23 -2.095895144973650E-018
24 1.598556002462387E-018
25 -6.012183194674783E-019
26 4.373987104435328E-019
27 -1.096344893602576E-018
28 2.108673138650750E-018
29 -1.917443644281366E-018
30 -5.574775601975155E-019
31 3.911069159997590E-018
32 -5.569027665180324E-018
33 4.374661466073393E-018
34 -1.500740536545008E-018
35 -1.259659760801246E-018
36 1.987490316384722E-018
37 -2.566062976033432E-019
38 -2.036131789814112E-018
39 3.062056883928202E-018
40 -3.088637442355000E-018
41 3.160733147274003E-018
42 -3.203120111095265E-018
43 2.017731412427283E-018
44 6.652505744134810E-019
45 -3.112895402546005E-018
46 3.663599529981498E-018
47 -2.437000263303890E-018
48 8.345827121050447E-019
49 -2.928818411446611E-020
50 -5.194142394169160E-020
0 0.000000000000000E+000
1 8.971999761003183E-002
2 4.071719650066735E-003
3 2.779775834786763E-004
4 2.252077243749956E-005
5 2.005889481097034E-006
6 1.897597704408893E-007
7 1.871654957880109E-008
8 1.903462209518203E-009
9 1.981687858254588E-010
10 2.101650492525309E-011
11 2.262502407646139E-012
12 2.466077342333506E-013
13 2.714661828662202E-014
14 3.029798634202052E-015
15 3.425038030968608E-016
16 7.188694084447888E-017
17 2.581268532253489E-017
18 -1.915134717478395E-017
19 -2.831068712794149E-017
20 -4.801714581503802E-017
21 1.582067810090848E-017
22 -8.049116928532384E-018
23 -3.330669073875469E-017
24 4.718447854656916E-018
25 -1.387778780781446E-017
26 -2.886579864025407E-017
27 1.110223024625157E-017
28 2.234323837058128E-017
29 -1.526556658859590E-017
30 -1.179611963664229E-017
31 -1.401656568589260E-017
32 -4.982125823005390E-017
33 -4.024558464266192E-018
34 2.872702076217592E-017
35 -2.248201624865942E-017
36 4.718447854656916E-018
37 1.554312234475219E-017
38 3.066991105526995E-017
39 3.774758283725533E-017
40 -4.843347944927246E-017
41 3.663735981263017E-017
42 -8.992806499463769E-017
43 3.302913498259841E-017
44 6.016021014687567E-017
45 5.134781488891349E-017
46 -4.302114220422482E-018
47 2.994132719535969E-017
48 3.150257832373882E-017
49 5.466981034540908E-017
50 -1.042029802807706E-016
51 2.556115041851825E-017
52 -1.994931997373328E-018
53 -1.734723475976807E-017
54 -3.420874694626263E-017
55 3.851086116668512E-018
56 -5.370703881624195E-017
57 1.387778780781446E-017
58 3.469446951953614E-018
59 -2.213507155346406E-017
60 4.475586568020162E-017
61 -6.938893903907229E-019
62 -1.233735336114705E-016
63 -1.486311074216928E-016
64 -3.053113317719181E-017
65 -5.315192730392937E-017
66 1.151856388048600E-017
67 8.035239140724571E-017
68 -1.561251128379126E-016
69 9.436895709313831E-018
70 7.091549569793188E-017
71 -3.316791286067655E-017
72 -1.817990202823694E-017
73 5.023759186428833E-017
74 -9.520162436160718E-017
75 -3.358424649491098E-017
76 5.107025913275720E-017
77 -1.948441408217150E-016
78 -3.858025010572419E-017
79 -2.525757381022231E-017
80 -4.329869796038111E-017
81 2.498001805406602E-018
82 7.494005416219807E-017
83 -1.673661209622424E-016
84 2.775557561562891E-017
85 -2.256528297550631E-016
86 -2.331468351712829E-017
87 -3.247402347028583E-017
88 -4.996003610813205E-017
89 -8.326672684688674E-017
90 2.220446049250313E-017
91 -5.079270337660091E-017
92 -7.771561172376097E-018
93 1.199040866595169E-016
94 2.220446049250313E-018
95 2.942091015256665E-017
96 -9.853229343548265E-017
97 -4.884981308350689E-017
98 -1.609823385706477E-017
99 -1.110223024625157E-018
100 -4.163336342344337E-018
50 changes: 50 additions & 0 deletions bin/waverf_L628_inf_ka01_N15_30.cof
Original file line number Diff line number Diff line change
@@ -0,0 +1,50 @@
6.28318530717959 0.200000000000000 1.00000000000000
6.25184754973936 1.00501255943798 0.000000000000000E+000
0.000000000000000E+000 15 30
0 -1.00501255943798
1 9.948427798901736E-002
2 5.016359814956224E-005
3 8.847110715062713E-007
4 1.624492305254340E-008
5 2.862936655948687E-010
6 4.911462881176951E-012
7 8.440483376220313E-014
8 1.486426726474920E-015
9 2.485328365286742E-017
10 1.701945934163557E-018
11 -1.771726142400146E-018
12 2.341093350754148E-018
13 7.727598558724145E-019
14 -3.220070811180185E-018
15 2.210949699793718E-018
0 0.000000000000000E+000
1 9.961371978412105E-002
2 5.033047475278240E-003
3 3.828087046639302E-004
4 3.456241508610625E-005
5 3.431196168595293E-006
6 3.618274090168757E-007
7 3.978398054327314E-008
8 4.510552742420515E-009
9 5.235239274499944E-010
10 6.189946379538588E-011
11 7.429291441300261E-012
12 9.028211049126138E-013
13 1.108097410125926E-013
14 1.378758218706366E-014
15 1.701421053388107E-015
16 1.864249495516409E-016
17 5.088522196198634E-018
18 6.962023550253586E-017
19 -1.295260195396016E-017
20 4.903485025427775E-017
21 1.017704439239727E-017
22 -6.476300976980080E-018
23 1.017704439239727E-017
24 2.775557561562892E-018
25 4.718447854656915E-017
26 7.494005416219807E-017
27 -1.174986034394957E-016
28 4.810966440042345E-017
29 5.458596537740353E-017
30 -4.625929269271485E-018
65 changes: 65 additions & 0 deletions bin/waverf_L628_inf_ka021_N20_40.cof
Original file line number Diff line number Diff line change
@@ -0,0 +1,65 @@
6.28318530717959 0.420000000000000 1.00000000000000
6.14614517779847 1.02229692358653 0.000000000000000E+000
0.000000000000000E+000 20 40
0 -1.02229692358653
1 0.204678319297992
2 9.854084431457394E-004
3 4.374188203232261E-005
4 2.218977600400424E-006
5 1.196605000730274E-007
6 6.812427795733324E-009
7 4.065109517719706E-010
8 2.517341718150921E-011
9 1.603740523375740E-012
10 1.044569612410480E-013
11 6.920956059730814E-015
12 4.685732865113373E-016
13 3.542355437589990E-017
14 -7.090251224258199E-018
15 1.221416650682497E-017
16 -1.259843755664276E-017
17 1.128753077166092E-017
18 -7.346960099930986E-018
19 3.015384905576172E-018
20 -6.296544979803001E-019
0 0.000000000000000E+000
1 0.206019267587262
2 2.266843591635342E-002
3 3.802362504788413E-003
4 7.611875607910365E-004
5 1.680165953315119E-004
6 3.945698465245146E-005
7 9.671282309953178E-006
8 2.445968383976705E-006
9 6.335877376764665E-007
10 1.672453130860396E-007
11 4.482545639117319E-008
12 1.216639128470698E-008
13 3.337233379918647E-009
14 9.236693375436200E-010
15 2.576384910968344E-010
16 7.234912495490686E-011
17 2.043746595570273E-011
18 5.803607941445854E-012
19 1.655884630802351E-012
20 4.743315753547074E-013
21 1.366118156098661E-013
22 3.942089710218255E-014
23 1.129304982860901E-014
24 3.122502256758253E-015
25 1.087324674742263E-015
26 1.915134717478395E-016
27 9.159339953157541E-017
28 1.873501354054952E-017
29 4.718447854656915E-017
30 1.679212324745549E-016
31 -6.522560269672795E-017
32 3.191891195797325E-017
33 9.298117831235686E-017
34 4.996003610813205E-017
35 -1.804112415015879E-017
36 -1.068589661201713E-016
37 1.221245327087672E-016
38 -1.068589661201713E-016
39 -4.440892098500626E-017
40 1.804112415015879E-017
Loading

0 comments on commit 96267e6

Please sign in to comment.