#############################################################################
#                                                                           #
# ww3_tic1.1:Test script for WW III, one-dimensional propagation            #
#            through an ice layer                                           #
#  (this test was previously known as ice_test6)                            #
#                                                                           #
# For recommended switches, please refer to switch files in the input       #
# directories input_IC1, input_IC2, input_IC3, input_IC4, input_IC5         #
#                                                                           #
# Remarks :                                                                 #
#                                                                           #
# - Test case input (default):                                              #
#   * ww3_grid_2500m.inp                                                    #
#     + Spectral grid (applies to all)                                      #
#      Number of directions        :  24                                    #
#      Directional increment (deg.):  15.0                                  #
#      First direction       (deg.):   0.0                                  #
#      Number of frequencies       :   3                                    #
#      Frequency range        (Hz) :   0.0800-0.1250                        #
#      Increment factor            :   1.250                                #
#     + Spatial grid (applies to ww3_grid_2500m.inp only, many others       #
#                    are avaliable)                                         #
#      Dimensions                  :    156       3                         #
#      Increments             (km) :    2.50    2.50                        #
#      X range                (km) :   -2.50  385.00                        #
#      Y range                (km) :   -2.50    2.50                        #
#                                                                           #
#  Waves propagate from low-x to high-x. Boundary forcing is constant (not  #
#      updated) defined at x=0. SWH~1 m at boundary. Direction is 270 deg.  #
#      Waves propagate through ice, starting from rest (except at           #
#      boundary). Ice field is uniform. See tab50.ww3 for summary output.   #
#      This test case is provided with many different ww3_grid.inp files,   #
#      each with different resolution, to demonstrate that results are      #
#      sensitive to spatial resulution, but converge to a proper result     #
#      if adequate resolution is used. Also, there are tests for three      #
#      alternative routines to represent S_{ice}. These are identified by   #
#      their switches : IC1, IC2, IC3, IC4                                  #
#          IC1 = constant dissipation rate                                  #
#          IC2 = Liu and M-C (JPO 1988), Liu et al. (JGR 1991): dissipation #
#                is caused by interaction between waves and turbulence at   #
#                ice-water interface                                        #
#          IC3 = Wang and Shen (JGR 2010), dissipation is caused by         #
#                interaction between waves and ice, where ice layer is      #
#                represented as a continuous viscoelastic layer             #
#          IC4 = Choose from several empirical parameterizations            # 
#            IC4METHOD = 1 - Wadhams et al. (1988)                          #
#            IC4METHOD = 2 - Meylan et al. (2014)                           #
#            IC4METHOD = 3 - Kohout & Meylan (2008) in                      #
#                               Horvat & Tziperman (2015)                   #
#            IC4METHOD = 4 - Kohout et al. (2014)                           #
#            IC4METHOD = 5 - Simple ki step function                        #
#            IC4METHOD = 6 - Simple ki step function via namelist           #
#            IC4METHOD = 7 - Doble et al. (GRL 2015)                        #
#          IC5 = Mosig et al. (JGR 2015), dissipation is caused by          #
#                interaction between waves and ice, where ice layer is      #
#                represented as a continuous viscoelastic layer             #
#                                                                           #
#   Scattering of waves by ice floes will be treated separately from        #
#   IC1/IC2/IC3/IC4 and and is not included in this test case.              #
#                                                                           #
#   IC4 settings used in tests for IC4METHOD = 1,2,3,4 are given now :      #
#      ------------> &SIC4 IC4METHOD = 1 /                                  #
#            ...   ICECOEF1, ICECOEF2 are required:                         #
#                  T T     Ice parameter 1                                  #
#                  T T     Ice parameter 2                                  #
#            ...                                                            #
#                  'IC1' 19680606 000000   0.18                             #
#                  'IC2' 19680606 000000   7.3                              #
#                                                                           #
#      ------------> &SIC4 IC4METHOD = 2 /                                  #
#            ...ICECOEF1 to ICECOEF5 are required:                          #
#               T T     Ice parameter 1                                     #
#               T T     Ice parameter 2                                     #
#               T T     Ice parameter 3                                     #
#               T T     Ice parameter 4                                     #
#               T T     Ice parameter 5                                     #
#            ...                                                            #
#               'IC1' 19680606 000000   0.0                                 #
#               'IC2' 19680606 000000   0.0                                 #
#               'IC3' 19680606 000000   0.0212                              #
#               'IC4' 19680606 000000   0.0                                 #
#               'IC5' 19680606 000000   0.459                               #
#                                                                           #
#      ------------> &SIC4 IC4METHOD = 3 /                                  #
#                  and                                                      #
#      ------------> &SIC4 IC4METHOD = 7 /                                  #
#            ...ICECOEF1 is required, representing ice thickness.           #
#               T T     Ice parameter 1                                     #
#            ...                                                            #
#               'IC1' 19680606 000000   0.25                                #
#                                                                           #
#      ------------> &SIC4 IC4METHOD = 4 /                                  #
#            ...ICECOEF1, ICECOEF2 are required:                            #
#               T T     Ice parameter 1                                     #
#               T T     Ice parameter 2                                     #
#            ...                                                            #
#               'IC1' 19680606 000000   5.35E-6                             #
#               'IC2' 19680606 000000   16.05E-6                            #
#                                                                           #
#  Reference (w/plots): Rogers and Orzech, NRL Memorandum Report (2013)     #
#  available from http://www7320.nrlssc.navy.mil/pubs.php                   #
#  (This report only covers IC1 and IC2, not IC3, which is newer)           #
#                                                                           #
#  Example run_test commands:                                               #
#  (some details will vary by local system and configuration)               #
#                                                                           #
# ./bin/run_test -g 100m  -w w100m_IC1  -i input_IC1_156x3                  #
#                ../model ww3_tic1.1                                        #
# ./bin/run_test -g 350m  -w w350m_IC1  -i input_IC1_156x3                  #
#                ../model ww3_tic1.1                                        #
# ./bin/run_test -g 1000m -w w1000m_IC1 -i input_IC1_156x3                  #
#                ../model ww3_tic1.1                                        #
# ./bin/run_test -g 2500m -w w2500m_IC1 -i input_IC1_156x3                  #
#                ../model ww3_tic1.1                                        #
# ./bin/run_test -g 5km   -w w5km_IC1   -i input_IC1 ../model ww3_tic1.1    #
# ./bin/run_test -g 10km  -w w10km_IC1  -i input_IC1 ../model ww3_tic1.1    #
# ./bin/run_test -g 20km  -w w20km_IC1  -i input_IC1 ../model ww3_tic1.1    #
#                                                                           #
# ./bin/run_test -g 1000m_nondisp -w work_1000m_nondisp_IC2_ifr -i input_IC2_ifr -c gnu_debug  ../model ww3_tic1.1                  
# bin/run_test -g 1000m_nondisp -c gnu_debug                                #
#      -i input_IC2_ifr -o netcdf -s IC2IS2 -w workIS2  ../model ww3_tic1.1 #
# ./bin/run_test -g 1000m -i input_IC2_ifr -o netcdf ../model ww3_tic1.1    #
# ./bin/run_test -w work_IC2_nondisp_SMPL -i input_IC2_nondisp              #
#       ../model ww3_tic1.1                                                 #
# ./bin/run_test -g 1000m -w work_IC2_nrl_1000m -i input_IC2_nrl            #
#       ../model ww3_tic1.1                                                 #
# ./bin/run_test -g 1000m_nondisp -w work_IC2_nrl_1000m_nondisp             #
#       -i input_IC2_nrl ../model ww3_tic1.1                                #
#                                                                           #
# ./bin/run_test -g 350m  -w w350m_IC3  -i input_IC3 ../model ww3_tic1.1    #
# ./bin/run_test -g 1000m -w w1000m_IC3 -i input_IC3 ../model ww3_tic1.1    #
# ./bin/run_test -g 2500m -w w2500m_IC3 -i input_IC3 ../model ww3_tic1.1    #
#                                                                           #
# ./bin/run_test -w work_IC3_nondisp -i input_IC3_nondisp                   #
#                ../model ww3_tic1.1                                        #
#                                                                           #
# ./bin/run_test -w work_IC4_M3 -g 1000m -i input_IC4/M3                    #
#                ../model ww3_tic1.1                                        #
#                                                                           #
# ./bin/run_test -s MPI -p mpirun -n 3 -g CGICE_CHENG                       #
#        -w work_IC3_CGICE_CHENG -i input_IC3NL ../model ww3_tic1.1         #
# ./bin/run_test -s MPI -p mpirun -n 3 -g CGOW_CHENG                        #
#        -w work_IC3_CGOW_CHENG -i input_IC3NL ../model ww3_tic1.1          #
# ./bin/run_test -s MPI -p mpirun -n 3 -g CGOW_NOCHENG                      #
#        -w work_IC3_CGOW_NOCHENG -i input_IC3NL ../model ww3_tic1.1        #
# ./bin/run_test -g CGICE_NOCHENG -w work_IC3_CGICE_NOCHENG_TT              #
#        -i input_IC3 ../model ww3_tic1.1                                   #
# ./bin/run_test -g CGOW_NOCHENG -w work_IC3_CGOW_NOCHENG_TT                #
#        -i input_IC3 ../model ww3_tic1.1                                   #
#                                                                           #
# The following is removed, since the code no long supports this            #
#    combination (CGICE=TRUE + CHENG=FALSE + namelist input) :              #
# ./bin/run_test -s MPI -p mpirun -n 3 -g CGICE_NOCHENG                     #
#        -w work_IC3_CGICE_NOCHENG -i input_IC3NL ../model ww3_tic1.1       #
#                                                                           #
#  IC5:                                                                     #
# ./bin/run_test -g 350m  -w w350m_IC5  -i input_IC5 ../model ww3_tic1.1    #
# ./bin/run_test -g 1000m -w w1000m_IC5 -i input_IC5 ../model ww3_tic1.1    #
# ./bin/run_test -g 2500m -w w2500m_IC5 -i input_IC5 ../model ww3_tic1.1    #
# ./bin/run_test -g 1000m  -w w1000m_IC5  -i input_IC5 -o both              #
#       ../model ww3_tic1.1                                                 #
#                                                                           #
#                                 origination: Erick Rogers, Dec 2012       #
#                                     updated: Stefan Zieger, Oct 2013      #
#                                     updated: Erick Rogers, Nov 2013       #
#                                     updated: Fabrice Ardhuin, Jun 2014    #
#                                     updated: Erick Rogers, July 2014      #
#                                     updated: Fabrice Ardhuin, Sep 2014    #
#                                     updated: Clarence Collins, Dec 2015   #
#                                     updated: Erick Rogers, Apr 2016       #
#                                     updated: Jessica Meixner, May 2016    #
#                                last updated: Qingxiang Liu, Jul 2018      #
#    Copyright 2009-2014 National Weather Service (NWS),                    #
#       National Oceanic and Atmospheric Administration.  All rights        #
#       reserved.  WAVEWATCH III is a trademark of the NWS.                 #
#       No unauthorized use without permission.                             #
#                                                                           #
#############################################################################
