Skip to content

Commit d82913b

Browse files
sbanihashsbanihash
authored
Addition of Regression Test (ww3_tic1.1/IC4_M10) (NOAA-EMC#1331)
1 parent d1fdda9 commit d82913b

16 files changed

+387
-0
lines changed

‎regtests/bin/matrix.base

+1
Original file line numberDiff line numberDiff line change
@@ -1957,6 +1957,7 @@
19571957
echo "$rtst -w work_IC4_M7 -i input_IC4_M7 $ww3 ww3_tic1.1" >> matrix.body
19581958
echo "$rtst -w work_IC4_M8 -i input_IC4_M8 $ww3 ww3_tic1.1" >> matrix.body
19591959
echo "$rtst -w work_IC4_M9 -i input_IC4_M9 $ww3 ww3_tic1.1" >> matrix.body
1960+
echo "$rtst -w work_IC4_M10 -i input_IC4_M10 $ww3 ww3_tic1.1" >> matrix.body
19601961
echo "$rtst -g 1000m -w work_IC5_M1 -i input_IC5_M1 $ww3 ww3_tic1.1" >> matrix.body
19611962
echo "$rtst -g 1000m -w work_IC5_M2 -i input_IC5_M2 $ww3 ww3_tic1.1" >> matrix.body
19621963
echo "$rtst -g 1000m -w work_IC5_M3 -i input_IC5_M3 $ww3 ww3_tic1.1" >> matrix.body

‎regtests/ww3_tic1.1/info

+9
Original file line numberDiff line numberDiff line change
@@ -54,6 +54,7 @@
5454
# IC4METHOD = 8 - Meylan et al. (2018) ; Liu et al. (2020) #
5555
# (NB: redundant with IC5+IC5VEMOD=3) #
5656
# IC4METHOD = 9 - RYW (2021) ; Yu et al. (2022) #
57+
# IC4METHOD = 10 - Meylan et al. (2021) #
5758
# IC5 = Choose from three different effective medium models #
5859
# IC5VEMOD = 1 - Extended Fox and Squire model (EFS) #
5960
# IC5VEMOD = 2 - Robinson and Palmer model (RP) #
@@ -101,6 +102,14 @@
101102
# 'IC1' 19680606 000000 5.35E-6 #
102103
# 'IC2' 19680606 000000 16.05E-6 #
103104
# #
105+
# ------------> &SIC4 IC4METHOD = 10 / #
106+
# ...ICECOEF1, ICECOEF5 are required: #
107+
# T T Ice parameter 1 #
108+
# T T Ice parameter 5 #
109+
# ... #
110+
# 'IC1' 19680606 000000 0.2 #
111+
# 'IC5' 19680606 000000 0.459 #
112+
# #
104113
# Reference (w/plots): Rogers and Orzech, NRL Memorandum Report (2013) #
105114
# available from http://www7320.nrlssc.navy.mil/pubs.php #
106115
# (This report only covers IC1 and IC2, not IC3, which is newer) #

‎regtests/ww3_tic1.1/input_IC4_M10/namelists_1-D.nml

+2
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,2 @@
1+
&SIC4 IC4METHOD = 10 /
2+
END OF NAMELISTS

‎regtests/ww3_tic1.1/input_IC4_M10/points.list

+16
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
0.00 0. 'Point 1 '
2+
1.00E3 0. 'Point 2 '
3+
2.00E3 0. 'Point 3 '
4+
3.00E3 0. 'Point 4 '
5+
4.00E3 0. 'Point 5 '
6+
5.00E3 0. 'Point 6 '
7+
6.00E3 0. 'Point 7 '
8+
7.00E3 0. 'Point 8 '
9+
8.00E3 0. 'Point 9 '
10+
9.00E3 0. 'Point 10 '
11+
10.00E3 0. 'Point 11 '
12+
11.00E3 0. 'Point 12 '
13+
12.00E3 0. 'Point 13 '
14+
13.00E3 0. 'Point 14 '
15+
14.00E3 0. 'Point 15 '
16+
15.00E3 0. 'Point 16 '

‎regtests/ww3_tic1.1/input_IC4_M10/switch

+1
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
NOGRB SHRD PR3 UQ FLX2 LN0 ST0 NL0 BT0 DB0 TR0 BS0 IC4 IS0 REF0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_grid.inp

+43
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,43 @@
1+
$ WAVEWATCH III Grid preprocessor input file
2+
$ ------------------------------------------
3+
'1-D parameterized ice test '
4+
$
5+
$ 1.1 0.04118 25 24 0.0
6+
1.1 0.0418 31 36 5.0
7+
$
8+
F T F F F T
9+
60. 60. 60. 60.
10+
$
11+
$ IC4METHOD determines calculation
12+
$ IC4METHOD = 1 - Wadhams et al. (1988)
13+
$ IC4METHOD = 2 - Meylan et al. (2014)
14+
$ IC4METHOD = 3 - Kohout & Meylan (2008) in Horvat & Tziperman (2015)
15+
$ IC4METHOD = 4 - Kohout et al. (2014)
16+
$ IC4METHOD = 5 - Simple ki step function
17+
$ IC4METHOD = 6 - Simple ki step function via namelist
18+
$ IC4METHOD = 7 - Doble et al. (GRL 2015)
19+
$ IC4METHOD = 8 - Meylan et al. (2018) ; Liu et al. (2020)
20+
$ IC4METHOD = 9 - RYW (2021) ; Yu et al. (2022)
21+
$ IC4M8 Fit to R21A L ChfM2=0.059
22+
$ IC4M10
23+
&SIC4 IC4METHOD = 10 , IC4CN = 0.059/
24+
END OF NAMELISTS
25+
$
26+
'RECT' F 'NONE'
27+
156 3
28+
1.0E3 1.0E3 1.
29+
-1.0E3 -1.0E3 1.
30+
$ dlim dmin file# scale layout# format# formatdescrip filetype# filenm
31+
-0.1 0.1 401 -1.0 1 1 '(....)' 'NAME' '../input_IC1/depth1d.flat'
32+
$
33+
10 1 1 '(....)' 'PART' 'input'
34+
$
35+
$ First grid
36+
$
37+
2 2 F
38+
$
39+
0 0 F
40+
0 0 F
41+
0 0
42+
$
43+
0. 0. 0. 0. 0

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_grid.nml

+81
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,81 @@
1+
! -------------------------------------------------------------------- !
2+
! WAVEWATCH III - ww3_grid.nml - Grid pre-processing !
3+
! -------------------------------------------------------------------- !
4+
5+
! -------------------------------------------------------------------- !
6+
! Define the spectrum parameterization via SPECTRUM_NML namelist
7+
! -------------------------------------------------------------------- !
8+
&SPECTRUM_NML
9+
SPECTRUM%XFR = 1.1
10+
SPECTRUM%FREQ1 = 0.0418
11+
SPECTRUM%NK = 31
12+
SPECTRUM%NTH = 36
13+
SPECTRUM%THOFF = 5.0
14+
/
15+
16+
! -------------------------------------------------------------------- !
17+
! Define the run parameterization via RUN_NML namelist
18+
! -------------------------------------------------------------------- !
19+
&RUN_NML
20+
RUN%FLCX = T
21+
RUN%FLSOU = T
22+
/
23+
24+
! -------------------------------------------------------------------- !
25+
! Define the timesteps parameterization via TIMESTEPS_NML namelist
26+
! -------------------------------------------------------------------- !
27+
&TIMESTEPS_NML
28+
TIMESTEPS%DTMAX = 60.
29+
TIMESTEPS%DTXY = 60.
30+
TIMESTEPS%DTKTH = 60.
31+
TIMESTEPS%DTMIN = 60.
32+
/
33+
34+
! -------------------------------------------------------------------- !
35+
! Define the grid to preprocess via GRID_NML namelist
36+
! -------------------------------------------------------------------- !
37+
&GRID_NML
38+
GRID%NAME = '1-D parameterized ice test'
39+
GRID%NML = '../input_IC4_M10/namelists_1-D.nml'
40+
GRID%TYPE = 'RECT'
41+
GRID%COORD = 'CART'
42+
GRID%CLOS = 'NONE'
43+
GRID%ZLIM = -0.1
44+
GRID%DMIN = 0.1
45+
/
46+
47+
! -------------------------------------------------------------------- !
48+
! Define the rectilinear grid type via RECT_NML namelist
49+
! -------------------------------------------------------------------- !
50+
&RECT_NML
51+
RECT%NX = 156
52+
RECT%NY = 3
53+
RECT%SX = 1.0E3
54+
RECT%SY = 1.0E3
55+
RECT%X0 = -1.0E3
56+
RECT%Y0 = -1.0E3
57+
/
58+
59+
! -------------------------------------------------------------------- !
60+
! Define the depth to preprocess via DEPTH_NML namelist
61+
! -------------------------------------------------------------------- !
62+
&DEPTH_NML
63+
DEPTH%SF = -1.0
64+
DEPTH%FILENAME = '../input_IC1/depth1d.flat'
65+
/
66+
67+
! -------------------------------------------------------------------- !
68+
! Define the input boundary points via INBND_COUNT_NML and
69+
! INBND_POINT_NML namelist
70+
! -------------------------------------------------------------------- !
71+
&INBND_COUNT_NML
72+
INBND_COUNT%N_POINT = 1
73+
/
74+
75+
&INBND_POINT_NML
76+
INBND_POINT(1) = 2 2 F
77+
/
78+
79+
! -------------------------------------------------------------------- !
80+
! WAVEWATCH III - end of namelist !
81+
! -------------------------------------------------------------------- !

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_ounf.inp

+20
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
$ WAVEWATCH III Grid output post-processing (netcdf)
2+
$---------------------------------------------------
3+
19680606 000000 3600. 99
4+
N
5+
$ Options: DPT CUR WND DT WLV ICE HS L T02 T01 TM1 FP DIR SPR DP EF
6+
$ TH1M STH1M PHS PTP PLP PDIR PSP WSF TWS PNR UST CHA CGE FAW
7+
$ TAW TWA WCC WCF WCH WCM SXY TWO BHD FOC TUS USS P2S WN USF
8+
$ P2L ABR UBR BED FBB TBB MSS MSC DTD FCT CFX CFT CFK US1 US2
9+
DPT WLV HS DIR
10+
$
11+
3 4
12+
0 1 2
13+
F
14+
ww3.
15+
4
16+
1 999 1 999 3 2
17+
$
18+
$ -------------------------------------------------------------------- $
19+
$ End of input file $
20+
$ -------------------------------------------------------------------- $

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_ounf.nml

+29
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,29 @@
1+
! -------------------------------------------------------------------- !
2+
! WAVEWATCH III ww3_ounf.nml - Grid output post-processing !
3+
! -------------------------------------------------------------------- !
4+
5+
! -------------------------------------------------------------------- !
6+
! Define the output fields to postprocess via FIELD_NML namelist
7+
! -------------------------------------------------------------------- !
8+
&FIELD_NML
9+
FIELD%TIMESTART = '19680606 000000'
10+
FIELD%TIMESTRIDE = '3600.'
11+
FIELD%TIMECOUNT = '99'
12+
FIELD%TIMESPLIT = 4
13+
FIELD%LIST = 'DPT WLV HS DIR'
14+
FIELD%PARTITION = '0 1 2'
15+
FIELD%SAMEFILE = F
16+
FIELD%TYPE = 4
17+
/
18+
19+
! -------------------------------------------------------------------- !
20+
! Define the content of the output file via FILE_NML namelist
21+
! -------------------------------------------------------------------- !
22+
&FILE_NML
23+
FILE%IXN = 999
24+
FILE%IYN = 999
25+
/
26+
27+
! -------------------------------------------------------------------- !
28+
! WAVEWATCH III - end of namelist !
29+
! -------------------------------------------------------------------- !

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_outf.inp

+13
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
$ WAVEWATCH III Grid output post-processing
2+
$ -----------------------------------------
3+
19680606 000000 3600. 99
4+
N
5+
$ Options: DPT CUR WND DT WLV ICE HS L T02 T01 TM1 FP DIR SPR DP EF
6+
$ TH1M STH1M PHS PTP PLP PDIR PSP WSF TWS PNR UST CHA CGE FAW
7+
$ TAW TWA WCC WCF WCH WCM SXY TWO BHD FOC TUS USS P2S WN USF
8+
$ P2L ABR UBR BED FBB TBB MSS MSC DTD FCT CFX CFT CFK US1 US2
9+
DPT WLV HS DIR
10+
$
11+
3 0
12+
$
13+
1 999 1 999 1 1

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_outp_spec.inp

+19
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,19 @@
1+
$ WAVEWATCH III Point output post-processing
2+
$ ------------------------------------------
3+
19680606 120000 3600. 1
4+
$
5+
1
6+
2
7+
3
8+
4
9+
5
10+
6
11+
7
12+
8
13+
9
14+
10
15+
11
16+
-1
17+
$
18+
1
19+
2 -1. 0. 33 F

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_outp_tab50.inp

+19
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,19 @@
1+
$ WAVEWATCH III Point output post-processing
2+
$ ------------------------------------------
3+
19680606 000000 600. 9999
4+
$
5+
1
6+
2
7+
3
8+
4
9+
5
10+
6
11+
7
12+
8
13+
9
14+
10
15+
11
16+
-1
17+
$
18+
2
19+
2 50

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_outp_tab51.inp

+10
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,10 @@
1+
$ WAVEWATCH III Point output post-processing
2+
$ ------------------------------------------
3+
19680606 000000 900. 49
4+
$
5+
$ 1
6+
11
7+
-1
8+
$
9+
2
10+
2 51

‎regtests/ww3_tic1.1/input_IC4_M10/ww3_prep_icecon.inp

+38
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,38 @@
1+
$ -------------------------------------------------------------------- $
2+
$ WAVEWATCH III Field preprocessor input file $
3+
$ -------------------------------------------------------------------- $
4+
$ Mayor types of field and time flag
5+
$ Field types : IC1, IC2, IC3, IC4, IC5 => Ice parameters (5)
6+
$ MDN => Mud densities
7+
$ MTH => Mud thicknesses
8+
$ MVS => Mud viscosities
9+
$ ICE => Ice concentrations.
10+
$ LEV => Water levels.
11+
$ WND => Winds.
12+
$ WNS => Winds (including air-sea temp. dif.)
13+
$ CUR => Currents.
14+
$ Format types : AI Transfer field 'as is'.
15+
$ LL Field defined on longitude-latitude grid.
16+
$ F1 Arbitrary grid, longitude and latitude of
17+
$ each grid point given in separate file.
18+
$ F2 Like F1, composite of 2 fields.
19+
$ Time flag : If true, time is included in file.
20+
$ Header flag : If true, write header on "*.ww3" data file
21+
$
22+
'ICE' 'AI' T T
23+
$
24+
$ Additional time input ---------------------------------------------- $
25+
$ If time flag is .FALSE., give time of field in yyyymmdd hhmmss format.
26+
$
27+
$ 19680606 000000
28+
$
29+
$ Define data files -------------------------------------------------- $
30+
$ The first input line identifies the file format with FROM, IDLA and
31+
$ IDFM, the second (third) lines give the file unit number and name.
32+
$
33+
'NAME' 1 2 '(I10,1x,I10)' '(1000(F6.2))'
34+
2345 '../input_IC2_nondisp/icecon.156x3.txt'
35+
$
36+
$ -------------------------------------------------------------------- $
37+
$ End of input file $
38+
$ -------------------------------------------------------------------- $

0 commit comments

Comments
 (0)