Introduce IC4M8 and IC4M9 to WW3 (#1176)

manual/eqs/ICE4.tex

@@ -52,6 +52,20 @@ \subsubsection{~$S_{ice}$: Empirical/parametric damping by sea ice} \label{sec:I
 
 {\code IC4M7}: This is a formula for dissipation from \cite{art:Dob15}, developed for a mixture of pancake and frazil ice, using data collected in the Weddell Sea (Antarctica). The formula depends on wave frequency and ice thickness:
 \begin{equation}\label{eq:ice7}
-  {\alpha=0.2T^{-2.13}h} \:\:\: .
+  {\alpha=2k_i=0.2h^1f^{2.13}} \:\:\: .
 \end{equation}
 This method is described in \cite{rep:RPLA18}.
+
+{\code IC4M8}: Like {\code IC4M7}, this method is in the general form of 
+\begin{equation}\label{eq:ice8}
+  {k_i=C_{hf}h^mf^n} \:\:\: .
+\end{equation}
+The formula is taken from \cite{Meylan2018}, where it is described as a ``Model with Order 3 Power Law''. It is applied by \cite{Liu2020}, where it is referred to as the ``M2'' model. The model specifies $m=1$ and $n=3$, and $C_{hf}$ is a user-specified calibration coefficient. \cite{Liu2020} provide calibration to two field cases and \cite{rep:RYW2021} provides a calibration to a third field case, \cite{art:RMK2021}. The third calibration is set as the default for {\code IC4M8}, $C_{hf}=0.059$, but can be changed in using the namelist parameter (constant and uniform) {\code IC4CN}, or using the spatially and/or temporally variable parameter  ${C_{ice,2}}$ . Further details on the calibrations are available in the inline documentation in {\file w3sic4md.F90}.  This method is functionally the same as the ``{\code M2}'' model in {\code IC5} (i.e., {\code IC5} with {\code IC5VEMOD=3}) and is redundantly included here as {\code IC4M8} because it is in the same ``family'' as {\code IC4M7} and {\code IC4M9}, being in the form of Eq. (\ref{eq:ice8}). 
+
+For an example of setting the namelist parameter, see {\file /regtests/ww3\_tic1.1/input\_IC4\_M8}.
+
+{\code IC4M9}: This formula is taken from the ``monomial power fit'' given in section 2.2.3 of \cite{rep:RYW2021}. Like {\code IC4M7} and {\code IC4M8}, it is a specific case of the general form of Eq. (\ref{eq:ice8}). The specificity is the constraint that $m=n/2-1$. This constraint is derived by \cite{rep:RYW2021} by invoking the scaling from \cite{art:YRW2019}, which is based on Reynolds number with ice thickness as the relevant length scale. This is also given as equation 2 in \cite{art:YRW2022}. The default namelist settings are $C_{hf}=2.9$ and $n=4.5$, from calibration by \cite{rep:RYW2021} to  \cite{art:RMK2021}. Further details, including alternative calibrations such as \cite{art:Yu2022}, are available in the inline documentation in {\file w3sic4md.F90}. Constant values can be set using namelist parameters, where $C_{hf}$ and $n$ are {\code IC4CN(1)} and {\code IC4CN(2)}, respectively. Spatially and/or temporally versions of the same can be specified as ${C_{ice,2}}$ and ${C_{ice,3}}$, respectively.
+
+The namelist default $C_{hf}$ values in {\code IC4M8} and {\code IC4M9} are consistent with those of identical formulae implemented in \cite{man:SWAN4145A}.
+
+

manual/eqs/ICE5.tex

@@ -25,7 +25,7 @@ \subsubsection{~$S_{ice}$: Damping by sea ice (effective medium models)} \label{
 \begin{align}
 k_i^{EFS} &\propto \eta h_i^3 \sigma^{11},\label{eq:fspw}\\ k_i^{RP} &\propto \frac{\eta}{\rho_w g^2} \sigma^3,\label{eq:rppw}
 \end{align}
-whereas previous field measurements \citep[e.g.,][]{Meylan2018, Rogers2021} support a power law $k_i \propto \sigma^n$, with $n$ between 2 and 4. Eqs.~(\ref{eq:fspw}) and (\ref{eq:rppw}) indicate at certain regimes (i.e., $k_r \approx k_0$ and low $k_i$), $k_i$ of the EFS model is too sensitive to wave frequency and $k_i$ of the RP model shows no dependence on ice thickness.
+whereas previous field measurements \citep[e.g.,][]{Meylan2018, RMK21} support a power law $k_i \propto \sigma^n$, with $n$ between 2 and 4. Eqs.~(\ref{eq:fspw}) and (\ref{eq:rppw}) indicate at certain regimes (i.e., $k_r \approx k_0$ and low $k_i$), $k_i$ of the EFS model is too sensitive to wave frequency and $k_i$ of the RP model shows no dependence on ice thickness.
 
 The third model included in the {\code IC5} module is based on the ``Model with Order 3 Power Law'' proposed by \citet[][their section 6.2; hereafter the M2 model]{Meylan2018}, which assumes the loss of wave energy is proportional to the horizontal ice velocity squared times the ice thickness. The attenuation rate is given by
 \begin{equation}
@@ -52,4 +52,4 @@ \subsubsection{~$S_{ice}$: Damping by sea ice (effective medium models)} \label{
 %
 \cit{IC5VEMOD} {the sea ice model to be selected: 1 - {\code EFS}, 2 - {\code RP}, 3 - {\code M2}; Default=3 (i.e., \textbf{the {\code M2} model is chosen}).}
 \end{clist}
-The first 6 parameters were introduced to improve the stability of the numerical solver for the EFS model \citep[the solver may fail for small wave periods in some rare cases, particularly for shallow water depth $d$ and low $G$; see][]{Liu2020}. Nonetheless, since version 7.12, the M2 model becomes the default option and these limiters are therefore not used by default.
+The first 6 parameters were introduced to improve the stability of the numerical solver for the EFS model \citep[the solver may fail for small wave periods in some rare cases, particularly for shallow water depth $d$ and low $G$; see][]{Liu2020}. Nonetheless, since version 7.12, the M2 model becomes the default option and these limiters are therefore not used by default.

manual/manual.bib

@@ -524,7 +524,7 @@ @TECHREPORT{rep:CR17
   INSTITUTION = "{N}aval {R}esearch {L}aboratory, {S}tennis {S}pace {C}enter, {MS}",
   TYPE        = "NRL Memorandum Report",
   NUMBER      = "NRL/MR/7320--17-9726",
-  NOTE       = "25 pp., www7320.nrlssc.navy.mil/pubs.php" }
+  NOTE       = "25 pp., www7320.nrlssc.navy.mil/pubs" }
 
 % item art:CRT17
 
@@ -1764,7 +1764,7 @@ @INPROCEEDINGS{pro:RZ14
   TITLE        = "New wave-ice interaction physics in {WAVEWATCH III}",
   BOOKTITLE    = Ice14,
   PUBLISHER    = "IAHR",
-  NOTE        = "8 pp., www7320.nrlssc.navy.mil/pubs.php" }
+  NOTE        = "8 pp., www7320.nrlssc.navy.mil/pubs" }
 
 % item rep:RPLA18
 
@@ -1775,7 +1775,18 @@ @TECHREPORT{rep:RPLA18
   INSTITUTION = "{N}aval {R}esearch {L}aboratory, {S}tennis {S}pace {C}enter, {MS}",
   TYPE        = "NRL Memorandum Report",
   NUMBER      = "NRL/MR/7320--18-9786",
-  NOTE       = "179 pp., www7320.nrlssc.navy.mil/pubs.php" }
+  NOTE       = "179 pp., www7320.nrlssc.navy.mil/pubs" }
+
+% item rep:RYW2021
+
+@TECHREPORT{rep:RYW2021,
+  AUTHOR      = "W. E. Rogers and J. Yu and D. W. Wang",
+  YEAR        = "2021",
+  TITLE       = "Incorporating dependencies on ice thickness in empirical parameterizations of wave dissipation by sea ice",
+  INSTITUTION = "{N}aval {R}esearch {L}aboratory, {S}tennis {S}pace {C}enter, {MS}",
+  TYPE        = "NRL Technical Report",
+  NUMBER      = "NRL/OT/7320-21-5145",
+  NOTE       = "35 pp., https://arxiv.org/abs/2104.01246" }
 
 % item rep:RMK18
 
@@ -1786,7 +1797,7 @@ @TECHREPORT{rep:RMK18
   INSTITUTION = "{N}aval {R}esearch {L}aboratory, {S}tennis {S}pace {C}enter, {MS}",
   TYPE        = "NRL Memorandum Report",
   NUMBER      = "NRL/MR/7320--18-9801",
-  NOTE       = "25 pp., www7320.nrlssc.navy.mil/pubs.php" }
+  NOTE       = "25 pp., www7320.nrlssc.navy.mil/pubs" }
 
 % item art:RH09
 
@@ -1811,6 +1822,33 @@ @ARTICLE{art:RTS16
   doi = "doi:10.1002/2016JC012251"
   }
 
+% item art:YRW2019
+
+@ARTICLE{art:YRW2019,
+  AUTHOR = "J. Yu and W. E. Rogers and D. W. Wang",
+  YEAR   = 2019,
+  TITLE  = "A Scaling for Wave Dispersion Relationships in Ice-Covered Waters",
+  JOURNAL = JGR,
+  VOLUME = "124",
+  PAGES  = "8429--8438" ,
+  doi = "doi:10.1029/2018JC014870"
+  }
+
+% item art:Yu2022
+
+@Article{art:Yu2022,
+AUTHOR = {Yu, Jie},
+TITLE = {Wave Boundary Layer at the Ice-Water Interface},
+JOURNAL = {Journal of Marine Science and Engineering},
+VOLUME = {10},
+YEAR = {2022},
+NUMBER = {10},
+ARTICLE-NUMBER = {1472},
+URL = {https://www.mdpi.com/2077-1312/10/10/1472},
+ISSN = {2077-1312},
+DOI = {10.3390/jmse10101472}
+}
+
 % item art:CFSRR10
 
 @ARTICLE{art:CFSRR10,
@@ -2346,7 +2384,7 @@ @TECHREPORT{rep:RC09
   INSTITUTION = "{N}aval {R}esearch {L}aboratory, {S}tennis {S}pace {C}enter, {MS}",
   TYPE        = "NRL Memorandum Report",
   NUMBER      = "NRL/MR/7320--09-9193",
-  NOTE       = "42 pp., www7320.nrlssc.navy.mil/pubs.php" }
+  NOTE       = "42 pp., www7320.nrlssc.navy.mil/pubs" }
 
 % item rep:RO13
 
@@ -2357,7 +2395,7 @@ @TECHREPORT{rep:RO13
   INSTITUTION = "{N}aval {R}esearch {L}aboratory, {S}tennis {S}pace {C}enter, {MS}",
   TYPE        = "NRL Memorandum Report",
   NUMBER      = "NRL/MR/7320--13-9462",
-  NOTE       = "31 pp., www7320.nrlssc.navy.mil/pubs.php" }
+  NOTE       = "31 pp., www7320.nrlssc.navy.mil/pubs" }
 
 % item rep:Roland2008
 
@@ -2606,6 +2644,17 @@ @MANUAL{man:SWAN3
   ADDRESS     = "P.O. Box 5048, 2600 GA Delft, The Netherlands",
   NOTE        = "see http://swan.ct.tudelft.nl" }
 
+% item man:SWAN4145A
+
+@MANUAL{man:SWAN4145A,
+  AUTHOR      = "{SWAN team}",
+  YEAR        = "2023",
+  TITLE       = "{SWAN Cycle III} version 41.45A User Manual",
+  ORGANIZATION = "Delft University of Technology,
+                  Faculty of Civil Engineering and Geosciences",
+  ADDRESS     = "P.O. Box 5048, 2600 GA Delft, The Netherlands",
+  NOTE        = "see https://swanmodel.sourceforge.io/" }
+
 % item man:Jones98
 
 @MANUAL{man:Jones98,
@@ -3485,10 +3534,14 @@ @article{Liu2021
 	title={{Global Wave Hindcasts Using the Observation-based Source Terms: Description and Validation}},
 	author={Liu, Qingxiang and Babanin, Alexander and Rogers, W Erick and Zieger, Stefan and Young, Ian and Bidlot, Jean-Raymond and Durrant, Tom and Ewans, Kevin and Guan, Changlong and Kirezci, Cagil and Lemos, Gil and MacHutchon, Keith and Moon, Il-Ju and Rapizo, Henrique and Ribal, Agustinus and Semedo, Alvaro and Wang, Juanjuan},
 	journal={Journal of Advances in Modeling Earth Systems (JAMES)},
-	year={submitted}
+	year = {2021},
+	volume = {13},
+	number = {8},
+	pages = {e2021MS002493},
+	doi = {https://doi.org/10.1029/2021MS002493},
 }
 
-@article{Rogers2021,
+@article{art:RMK2021,
 	title = {Estimates of spectral wave attenuation in Antarctic sea ice, using model/data inversion},
 	journal = {Cold Regions Science and Technology},
 	volume = {182},
@@ -3499,6 +3552,18 @@ @article{Rogers2021
 	author = {W. Erick Rogers and Michael H. Meylan and Alison L. Kohout}
 }
 
+@article{art:YRW2022,
+	title = {A new method for parameterization of wave dissipation by sea ice},
+	journal = {Cold Regions Science and Technology},
+	volume = {199},
+	pages = {103582},
+	year = {2022},
+	issn = {0165-232X},
+	doi = {https://doi.org/10.1016/j.coldregions.2022.103582},
+	url = {https://www.sciencedirect.com/science/article/pii/S0165232X2200101X},
+	author = {Jie Yu and W. Erick Rogers and David W. Wang},
+}
+
 @article{Forristall1981,
 	author = {Forristall, George Z.},
 	doi = {10.1029/JC086iC09p08075},

model/nml/ww3_shel.nml

@@ -319,9 +319,9 @@
 ! * the number of each homogeneous input is defined by HOMOG_COUNT
 ! * the total number of homogeneous input is automatically calculated
 ! * the homogeneous input must start from index 1 to N
-! * if VALUE1 is equal 0, then the homogeneous input is desactivated
+! * if VALUE1 is equal 0, then the homogeneous input is deactivated
 ! * NAME can be IC1, IC2, IC3, IC4, IC5, MDN, MTH, MVS, LEV, CUR, WND, ICE, MOV
-! * each homogeneous input is defined over a maximum of 3 values detailled below :
+! * each homogeneous input is defined over a maximum of 3 values detailed below :
 !     - IC1 is defined by thickness
 !     - IC2 is defined by viscosity
 !     - IC3 is defined by density

model/src/w3gdatmd.F90

@@ -199,10 +199,14 @@ MODULE W3GDATMD
   !                               Default is 1.0, meaning that 100% ice
   !                               concentration result in zero source term
   !                               If set to 0.0, then ice has no direct impact on Sln / Sin / Snl / Sds
-  !      IC3PARS   R.A.  Public   various parameters for use in IC4, handled as
+  !      IC3PARS   R.A.  Public   various parameters for use in IC3, handled as
   !                               an array for simplicity
-  !      IC4_KI    R.A.  Public   KI (dissipation rate) values for use in IC4
-  !      IC4_FC    R.A.  Public   FC (frequency bin separators) for use in IC4
+  !      IC4_KI    R.A.  Public   KI (dissipation rate) values for use in IC4M6
+  !      IC4_FC    R.A.  Public   FC (frequency bin separators) for use in IC4M6
+  !      IC4_CN    R.A.  Public   Coefficients for use in IC4M2
+  !      IC4_FMIN  Real  Public   Minimum frequency below which ki is set to 
+  !                               some background level dissipation (for S_ice)
+  !      IC4_KIBK  Real  Public   Low, background level dissipation (for S_ice)
   !      PFMOVE    Real  Public   Tunable parameter in GSE correction
   !                               for moving grids.
   !      GRIDSHIFT Real  Public   Grid offset for multi-grid w/SCRIP
@@ -615,7 +619,7 @@ MODULE W3GDATMD
        IPARS = -1, NAUXGR
   !
 #ifdef W3_IC4
-  INTEGER, PARAMETER      :: NIC4=10
+  INTEGER, PARAMETER      :: NIC4=16 ,  NIC42=5
 #endif
   INTEGER, PARAMETER      :: RLGTYPE = 1
   INTEGER, PARAMETER      :: CLGTYPE = 2
@@ -732,6 +736,8 @@ MODULE W3GDATMD
     INTEGER, POINTER      :: IC4PARS(:)
     REAL, POINTER         :: IC4_KI(:)
     REAL, POINTER         :: IC4_FC(:)
+    REAL, POINTER         :: IC4_CN(:)
+    REAL                  :: IC4_FMIN, IC4_KIBK
 #endif
 #ifdef W3_IC5
     REAL,    POINTER      :: IC5PARS(:)
@@ -1144,6 +1150,8 @@ MODULE W3GDATMD
   INTEGER, POINTER        :: IC4PARS(:)
   REAL, POINTER           :: IC4_KI(:)
   REAL, POINTER           :: IC4_FC(:)
+  REAL, POINTER           :: IC4_CN(:)
+  REAL, POINTER           :: IC4_FMIN, IC4_KIBK
 #endif
 #ifdef W3_IC5
   REAL,    POINTER        :: IC5PARS(:)
@@ -1840,6 +1848,8 @@ SUBROUTINE W3DIMX  ( IMOD, MX, MY, MSEA, NDSE, NDST   &
     CHECK_ALLOC_STATUS ( ISTAT )
     ALLOCATE ( GRIDS(IMOD)%IC4_FC(NIC4), STAT=ISTAT )
     CHECK_ALLOC_STATUS ( ISTAT )
+    ALLOCATE ( GRIDS(IMOD)%IC4_CN(NIC42), STAT=ISTAT )
+    CHECK_ALLOC_STATUS ( ISTAT )
 #endif
 #ifdef W3_IC5
     ALLOCATE ( GRIDS(IMOD)%IC5PARS(9), STAT=ISTAT )
@@ -2314,6 +2324,9 @@ SUBROUTINE W3SETG ( IMOD, NDSE, NDST )
     IC4PARS => GRIDS(IMOD)%IC4PARS
     IC4_KI => GRIDS(IMOD)%IC4_KI
     IC4_FC => GRIDS(IMOD)%IC4_FC
+    IC4_CN => GRIDS(IMOD)%IC4_CN
+    IC4_FMIN => GRIDS(IMOD)%IC4_FMIN
+    IC4_KIBK => GRIDS(IMOD)%IC4_KIBK
 #endif
 #ifdef W3_IC5
     IC5PARS => GRIDS(IMOD)%IC5PARS

model/src/w3gridmd.F90

@@ -114,6 +114,7 @@ MODULE W3GRIDMD
   !/    07-Jun-2021 : S_{nl} GKE NL5 (Q. Liu)             ( version 7.13 )
   !/    19-Jul-2021 : Momentum and air density support    ( version 7.14 )
   !/    28-Feb-2023 : GQM as an alternative for NL1       ( version 7.15 )
+  !/    11-Jan-2024 : New namelist parameters for IC4     ( version 7.15 )
   !/
   !/    Copyright 2009-2013 National Weather Service (NWS),
   !/       National Oceanic and Atmospheric Administration.  All rights
@@ -764,9 +765,10 @@ MODULE W3GRIDMD
 
 #ifdef W3_IC4
   INTEGER                 :: IC4METHOD
-  REAL                    :: IC4KI(NIC4), IC4FC(NIC4)
+  REAL                    :: IC4KI(NIC4), IC4FC(NIC4),            &
+                             IC4CN(NIC42), IC4FMIN, IC4KIBK
 #endif
-  !
+
 #ifdef W3_IC5
   REAL                    :: IC5MINIG, IC5MINWT,                  &
        IC5MAXKRATIO, IC5MAXKI, IC5MINHW,    &
@@ -970,7 +972,8 @@ MODULE W3GRIDMD
        IC3VISC, IC3ELAS, IC3DENS, IC3HICE
 #endif
 #ifdef W3_IC4
-  NAMELIST /SIC4/  IC4METHOD, IC4KI, IC4FC
+  NAMELIST /SIC4/  IC4METHOD, IC4KI, IC4FC, IC4CN, IC4FMIN,  &
+                   IC4KIBK
 #endif
 #ifdef W3_IC5
   NAMELIST /SIC5/  IC5MINIG, IC5MINWT, IC5MAXKRATIO,        &
@@ -2891,6 +2894,9 @@ SUBROUTINE W3GRID()
     IC4METHOD = 1 !switch for methods within IC4
     IC4KI=0.0
     IC4FC=0.0
+    IC4CN=0.0
+    IC4FMIN=0.0
+    IC4KIBK=0.0
 #endif
     !
 #ifdef W3_IC5
@@ -5313,6 +5319,9 @@ SUBROUTINE W3GRID()
     IC4PARS(1)=IC4METHOD
     IC4_KI=IC4KI
     IC4_FC=IC4FC
+    IC4_CN=IC4CN
+    IC4_FMIN=IC4FMIN
+    IC4_KIBK=IC4KIBK
 #endif
     !
 #ifdef W3_IC5

model/src/w3iogrmd.F90

@@ -1292,11 +1292,11 @@ SUBROUTINE W3IOGR ( INXOUT, NDSM, IMOD, FEXT &
 #endif
 #ifdef W3_IC4
       WRITE   (NDSM)                                             &
-           IC4PARS,IC4_KI,IC4_FC
+           IC4PARS,IC4_KI,IC4_FC,IC4_CN,IC4_FMIN,IC4_KIBK
 #ifdef W3_ASCII
       WRITE   (NDSA,*)                                           &
-           'IC4PARS,IC4_KI,IC4_FC:',                             &
-           IC4PARS,IC4_KI,IC4_FC
+           'IC4PARS,IC4_KI,IC4_FC,IC4_CN,IC4_FMIN,IC4_KIBK:',    &
+           IC4PARS,IC4_KI,IC4_FC,IC4_CN,IC4_FMIN,IC4_KIBK
 #endif
 #endif
 #ifdef W3_IC5
@@ -1338,7 +1338,7 @@ SUBROUTINE W3IOGR ( INXOUT, NDSM, IMOD, FEXT &
 #endif
 #ifdef W3_IC4
       READ   (NDSM,END=801,ERR=802,IOSTAT=IERR)                                              &
-           IC4PARS,IC4_KI,IC4_FC
+           IC4PARS,IC4_KI,IC4_FC,IC4_CN,IC4_FMIN,IC4_KIBK
 #endif
 #ifdef W3_IC5
       READ   (NDSM,END=801,ERR=802,IOSTAT=IERR)                                              &