Test PR to compare difference between COAWST vs UFS

CMakeLists.txt

@@ -20,13 +20,10 @@ if(hasParent)
   remove_definitions(-DDEBUG)
 endif()
 
-set(valid_caps "MULTI_ESMF" "NUOPC_MESH")
-
-set(UFS_CAP "" CACHE STRING "Valid options are ${valid_caps}")
+set(MULTI_ESMF OFF CACHE BOOL "Build ww3_multi_esmf library")
 set(NETCDF ON CACHE BOOL "Build NetCDF programs (requires NetCDF)")
 set(ENDIAN "BIG" CACHE STRING "Endianness of unformatted output files. Valid values are 'BIG', 'LITTLE', 'NATIVE'.") 
 set(EXCLUDE_FIND "" CACHE STRING "Don't try and search for these libraries (assumd to be handled by the compiler/wrapper)")
-set(BUILD_TOOLS OFF CACHE BOOL "Build WW3 tools")
 
 # make sure all "exclude_find" entries are lower case
 list(TRANSFORM EXCLUDE_FIND TOLOWER)
@@ -47,12 +44,6 @@ else()
   endif()
 endif()
 
-if(UFS_CAP)
-  if(NOT UFS_CAP IN_LIST valid_caps)
-    message(FATAL_ERROR "Invalid UFS_CAP selection. Valids options are ${valid_caps}")
-  endif()
-endif()
-
 message(STATUS "Build with switch: ${switch_file}")
 # Copy switch file to build dir
 configure_file(${switch_file} ${CMAKE_BINARY_DIR}/switch COPYONLY)
@@ -67,9 +58,3 @@ if(NOT CMAKE_BUILD_TYPE MATCHES "^(Debug|Release|RelWithDebInfo|MinSizeRel)$")
 endif()
 
 add_subdirectory(model)
-
-# Turn on unit testing.
-include(CTest)
-if(BUILD_TESTING)
-  add_subdirectory(regtests/unittests)
-endif()

manual/defs.tex

@@ -94,9 +94,6 @@
 \newcommand{\cR}{{\cal R}}
 \newcommand{\cS}{{\cal S}}
 
-\newcommand{\rd}{{\mathrm d}}
-
-
 \newcommand{\marbox}[1]{\marginpar{\fbox{{\small #1}}}}
 
 \newcommand{\proddefH}[3]{

manual/eqs/NL1.tex

@@ -1,92 +1,65 @@
+\vsssub
+\subsubsection{~$S_{nl}$: Discrete Interaction Approximation (\dia)} \label{sec:NL1}
+\vsssub
 
-\noindent
-
-
-Resonant nonlinear interactions occur between four wave components
-(quadruplets) with wavenumber vector $\bk$, $\bk_1$, $\bk_2$ and $\bk_3$ are such that 
-% eq:resonance
-\begin{equation} \left .
-\begin{array}{ccc}
-  \bk + \bk_1 & = & \bk_2 + \bk_3     \\
-  f_r + f_{r,1}& =& f_{r,2} +  f_{r,3} 
-\end{array} \:\:\: \right \rbrace \:\:\: , \label{eq:resonance}
-\end{equation}
+\opthead{NL1}{\wam\ model}{H. L. Tolman}
 
-Nonlinear 4-wave interaction theories were 
+\noindent
+Nonlinear wave-wave interactions can be modeled using the discrete interaction
+approximation \citep[\dia,][]{art:Hea85b}. This parameterization was
 originally developed for the spectrum $F(f_r ,\theta)$. To assure the
 conservative nature of $S_{nl}$ for this spectrum (which can be considered as
 the "final product" of the model), this source term is calculated for
 $F(f_r,\theta)$ instead of $N(k,\theta)$, using the conversion
 (\ref{eq:jac_fr}).
 
-\vsssub
-\subsubsection{~$S_{nl}$: Discrete Interaction Approximation (\dia)} \label{sec:NL1}
-\vsssub
-
-\opthead{NL1}{\wam\ model}{H. L. Tolman}
-
-
-
- In the \dia, for each component $\bk$, only 2 quadruplets configuration are 
-used, while there should be thousands for the full integral, and the interaction caused by these 2 quadruplets 
-is scaled so that it gives the right order of magnitude for the flux of energy towards low frequencies. 
-
-Both quadruplets used the DIA use 
+Resonant nonlinear interactions occur between four wave components
+(quadruplets) with wavenumber vector $\bk_1$ through $\bk_4$. In the \dia, it
+is assumed that $\bk_1 = \bk_2$. Resonance conditions then require that
+%--------------------------%
+% Resonance conditions DIA %
+%--------------------------%
 % eq:resonance
 \begin{equation} \left .
 \begin{array}{ccc}
-  \bk_1 & = & \bk\\
-  f_{r,2}  & = & (1+\lambda)f_{r}    \\
-  f_{r,3}  & = & (1-\lambda)f_{r} 
-\end{array} \:\:\: \right \rbrace \:\:\: , \label{eq:DIAchoice}
-\end{equation}
-where $\lambda$ is a constant, usually 0.25, and they only differ by the choice of the interacting angles 
-taking either a plus sign or a minus sign in the following 
-\begin{equation} \left .
-\begin{array}{ccc}
-  \theta_{2,\pm}  & = & \theta \pm \delta_{\theta,2}   \\
- \theta_{3,\pm}  & = & \theta \mp  \delta_{\theta,3}   \\
- \end{array} \:\:\: \right \rbrace \:\:\: , \label{eq:DIAangles}
+  \bk_1 + \bk_2 & = & \bk_3 + \bk_4          \\
+  \sigma_2  & = & \sigma_1                   \\
+  \sigma_3  & = & (1+\lambda_{nl})\sigma_1   \\
+  \sigma_4  & = & (1-\lambda_{nl})\sigma_1
+\end{array} \:\:\: \right \rbrace \:\:\: , \label{eq:resonance}
 \end{equation}
-where $\delta_{\theta,2}$ and $\delta_{\theta,3}$ are only a function of $\lambda$ given by the geometry of
-the interacting wavenumbers along the "figure of 8", namely  
-\begin{eqnarray}
-\cos(\delta_{\theta,2})&=&(1-\lambda)^4+4-(1+\lambda)^4)/[4(1-\lambda)^2], \\
-\sin(\delta_{\theta,3})&=&\sin(\delta_{\theta,2}) (1-\lambda)^2/(1+\lambda)^2.
-\end{eqnarray}
-
-Hence for any $\bk$ one quadruplet selects $\bk_{2,+}$ and $\bk_{3,+}$, and the other quadruplet selects its mirror image 
-$\bk_{2,-}$, $\bk_{2,-}$. Because there are 3 different components interacting in the two DIA-selected quadruplets, any discrete spectral component $(f_r,\theta)$ is actually involved in 6 quadruplets and directly exchanges energy with 12 other components $(f_r',\theta')$. Because the values of $f'_r$ and $\theta'$ do not fall exacly on other discrete components, the spectral density is interpolated using a bilinear interpolation, so that each source term value
-$S_{nl}(\bk)$ contains the direct exchange of energy with 48 other discrete components. 
-we compute the three contributions that correspond to the situation in which $\bk$ takes the role of $\bk$,$\bk_{2,+}$, $\bk_{2,-}$, $\bk_{3,+}$  and $\bk_{3,-}$ in the quadruplet, namely the full source term is, without making explicit that bilinear interpolation, 
-\begin{eqnarray} 
-S_{\mathrm{nl}}(\bk) &=& -2  \left[\delta S_{\mathrm{nl}}(\bk,\bk_{2,+},\bk_{3,+})+\delta S_{\mathrm{nl}}(\bk,\bk_{2,-},\bk_{3,-})\right] \nonumber  \\
- & & + \delta  S_{\mathrm{nl}}(\bk_4,\bk,\bk_5) + \delta  S_{\mathrm{nl}}(\bk_6,\bk,\bk_7)  \\
-     & &        +   \delta S_{\mathrm{nl}}(\bk_8,\bk_9,\bk) + \delta S_{\mathrm{nl}}(\bk_{10},\bk_{11},\bk) . \label{eq:diasum}
-\end{eqnarray}
-where the geometry of the quadruplet $(\bk_4,\bk_4,\bk,\bk_5)$ is obtained from that of $(\bk,\bk,\bk_{2,+},\bk_{3,+})$ by a dilation by a factor $(1+\lambda)^2$ and rotation by the angle $\delta_{\theta,2}$;  $(\bk_6,\bk_6,\bk,\bk_7)$ has the same dilation but the opposite rotation; $(\bk_8,\bk_8,\bk_9,\bk)$ is dilated by a factor  $(1-\lambda)^2$ and rotated by the angle $-\delta_{\theta,3}$: and $(\bk_{10},\bk_{10},\bk_{11},\bk)$ is dilated by the same factor and rotated by the opposite angle. 
-
-
-The elementary contributions $\delta  S_{\mathrm{nl}}(\bk_l,\bk_m,\bk_n)$   are given by  
+where $\lambda_{nl}$ is a constant. For these quadruplets, the contribution
+$\delta S_{nl}$ to the interaction for each discrete $(f_r,\theta)$
+combination of the spectrum corresponding to $\bk_1$ is calculated as
 %----------------------------%
 % Nonlinear interactions DIA %
 %----------------------------%
 % eq:snl_dia
-
-\begin{equation}
-\delta S_{\mathrm{nl}}(\bk_l,\bk_m,\bk_n) =  \frac{C}{g^4} f_{r,l}^{11} \left [ F_l^2 \left ( \frac{F_m}{(1+\lambda)^4} +
-        \frac{F_n}{(1-\lambda)^4} \right ) - \frac{2 F_l F_m F_n}{(1-\lambda^2)^4} \right] ,
-      \label{eq:snl_dia}
-\end{equation}
-where the spectral densities are  $F_l = F(f_{r,l} ,\theta_l)$, etc. 
- $C$ is a proportionality constant that was tuned to reproduce the inverse energy cascade.  Default values for different source term packages are presented in Table~\ref{tab:snl_par}.
+\begin{eqnarray}
+\left ( \begin{array}{c}
+  \delta S_{nl,1} \\ \delta S_{nl,3} \\ \delta S_{nl,4}
+\end{array} \right ) & = & D
+\left ( \begin{array}{r} -2 \\ 1 \\ 1 \end{array} \right )
+C g^{-4} f_{r,1}^{11} \times \nonumber \\
+& & \left [ F_1^2
+\left ( \frac{F_3}{(1+\lambda_{nl})^4} +
+        \frac{F_4}{(1-\lambda_{nl})^4} \right ) -
+\frac{2 F_1 F_3 F_4}{(1-\lambda_{nl}^2)^4}
+\right ] \: , \label{eq:snl_dia}
+\end{eqnarray}
+where $F_1 = F(f_{r,1} ,\theta_1 )$ etc. and $\delta S_{nl,1} = \delta
+S_{nl}(f_{r,1} ,\theta_1 )$ etc., $C$ is a proportionality constant. The
+nonlinear interactions are calculated by considering a limited number of
+combinations $(\lambda_{nl},C)$. In practice, only one combination is
+used. Default values for different source term packages are presented in
+Table~\ref{tab:snl_par}.
 
 
 % tab:snl_par
 
 \begin{table} \begin{center}
  \begin{tabular}{|l|c|c|} \hline
-                    & $\lambda$ &     $C$      \\ \hline
+                    & $\lambda_{nl}$ &     $C$      \\ \hline
 ST6                 &      0.25      & $3.00 \; 10^7$  \\ \hline
 \wam-3              &      0.25      & $2.78 \; 10^7$  \\ \hline
 ST4 (Ardhuin et al.)&      0.25      & $2.50 \; 10^7$  \\ \hline
@@ -95,7 +68,7 @@ \subsubsection{~$S_{nl}$: Discrete Interaction Approximation (\dia)} \label{sec:
 \caption{Default constants in \dia\ for input-dissipation packages.}
 \label{tab:snl_par} \botline \end{table}
 
-This parameterization was developed for deep water, using the appropriate dispersion
+This source term is developed for deep water, using the appropriate dispersion
 relation in the resonance conditions. For shallow water the expression is
 scaled by the factor $D$ (still using the deep-water dispersion relation,
 however)
@@ -159,37 +132,3 @@ \subsubsection{~$S_{nl}$: Discrete Interaction Approximation (\dia)} \label{sec:
 above constants can be reset by the user in the input files of the
 model (see \para\ref{sub:ww3grid}).
 
-\vsssub
-\subsubsection{~$S_{nl}$: Gaussian Quadrature Method  (\dia)} \label{sec:GQM}
-\vsssub
-
-\opthead{NL1 , but with a negative IQTYPE}{TOMAWAC model, M. Benoit}{adaptation to WW3 by S. Siadatmousavi \& F. Ardhuin}
-
-\noindent
-Changing the namelist parameter IQTYPE to a negative value replaces the
-DIA parameterization with the possibility to perform an exact but fast cal-
-culation of $S_{\mathrm{nl}}$ using the Gaussian Quadrature Method of \cite{Lavrenov2001}.
-More details can be found in \cite{Gagnaire-Renou2009}.
-
-
-The quadruplet configurations that are used correspond to the three integrals over $f_1$, $f_2$ and $\theta_1$, with all other frequencies and directions given by the resonance conditions (\ref{eq:resonance}) with only one ambiguity on the angle $\theta_2$ which can be defined by a sign index $s$, as in the DIA.  Starting from eq. (A4) in \cite{Lavrenov2001} as writen in (2.25) of \cite{Gagnaire-Renou2009}, the source term is 
-\begin{equation}
-S_{\mathrm{nl}}(\sigma,\theta) =  8 \sum_s \int_{\sigma_1=0}^\infty \int_{\theta_1=0}^{2 \pi} \int_{\sigma_2=0}^{(\sigma+\sigma_1)/2}  T  \frac{F_2 F_3 (F \sigma_1^4 + F_1 \sigma^4) - F F_1 (F_2 \sigma_3^4 + F_3 \sigma_2^4)}{\sqrt{B}\sqrt{((\left| \bk+\bk_1 \right|/g- \sigma_3^2)^2-\sigma_2^4} } {\mathrm d}\sigma_1 {\mathrm d}\theta_1 {\mathrm d}\sigma_2 ,
-      \label{eq:snl_gqm}
-\end{equation}
-where $B$ is given by eq. (A5) of Lavrenov (2001) and  
-\begin{equation}
-T(\bk,\bk_1,\bk_2,\bk_3) = \frac{\pi g^2 D^2(\bk,\bk_1,\bk_2,\bk_3) }{4 \sigma \sigma_1 \sigma_2 \sigma_3}
-\end{equation}
-where $ D(\bk,\bk_1,\bk_2,\bk_3)$ is given by \cite{Webb1978} in his eq. (A1). 
-
-This triple integral is performed using quadrature functions to best resolve the effect of the singularities in the denominator. It is thus replaced with weighted sums over the 3 dimensions. 
-
-Compared to the DIA, there is no bilinear interpolation and the nearest neighbor is used in frequency and direction. Also, 
-the source term is computed by a loop over the quadruplet configuration, which allows for filtering based on 
-both the value of the coupling coefficient and the energy level at the frequency corresponding to $\bk$. Within 
-that loop, the source term contribution is computed for all 4 interacting components, so that any filtering still 
-conserves energy, action, momentum ... (One may argue that this multiplies by 4 the number of calculations, but it may have the benefit of properly dealing with the high frequency boundary... this is to be verified. The same question arises for the DIA: why have the wavenumber $\bk$ play the role of the other members of the quadruplets when this will also be computed as we loop on the spectral components?). 
-
-If a very aggressive filtering is performed, the source may need to be rescaled.
-

manual/eqs/output.tex

@@ -12,9 +12,9 @@ \subsection{~Output parameters} \label{sub:outpars}
 in \para\ref{sec:ww3shel}. That input file also provides a list of flags  
 indicating if output parameters are available in different field
 output file types (ASCII, grib, igrads, NetCDF). 
-For any details on how these parameters are computed, the user may read the code of the  {\code w3iogo} routine, in the {\code w3iogomd.F90} module. 
+For any details on how these parameters are computed, the user may read the code of the  {\code w3iogo} routine, in the {\code w3iogomd.ftn} module. 
 
-Selection of field outputs in {\code ww3\_shel.nml} or {\code ww3\_shel.inp} is most easily performed by providing a list of the 
+Selection of field outputs in  {\code ww3\_shel.inp} is most easily performed by providing a list of the 
 requested parameters, for example, {\textbf HS DIR SPR}  will request the calculation of significant wave height, mean direction and directional spread. These will thus be stored in the {\code out\_grd.XX} file and can be post-processed, for example in NetCDF using   {\code ww3\_ouf}. Examples are given in \para\ref{sec:ww3multi} and
 \para\ref{sec:ww3ounf}. The names for these namelists are the bold names below, for 
 example \textbf{HS}. 
@@ -26,9 +26,6 @@ \subsection{~Output parameters} \label{sub:outpars}
 file extensions, NetCDF variable names and namelist-based selection (see 
 also \para\ref{sec:ww3ounf}), and the long parameter name/definition.
 
-When the result is not overly sensitive to the contribution of the unresolved part of the spectrum (for $f<f_{NK}$), the contribution of the tail is parameterized assuming a power law decay of the spectrum, by default $E(f,\theta) = E(f_{NK},\theta) (f_{NK}/f)^{-5}$, for some parameters this is either unnecessary or misleading, and the integrals are computed only up to $f_{NK}$. 
-
-
 Finally we note that in all definitions the frequency is the \emph{relative} frequency. Thus, in the presence of currents, these can only be compared to drifting measurement data or data obtained in wavenumber and converted to frequency. Comparison to fixed instrument data requires the use of the full spectrum and proper conversion to the fixed reference frame.
 
 \begin{list}{\Roman{outgrps})\hfill}
@@ -303,12 +300,8 @@ \subsection{~Output parameters} \label{sub:outpars}
 \item \textbf{MSD} Direction of the maximum slope variance mss$_u$
 \item \textbf{MCD} Spectral tail direction
 \item \textbf{QP}  Peakedness parameter \citep{art:G70}
-      \begin{equation} Q_p = \frac{2}{E^2} \int_0^{f_{NK}} f \left( \int_0^{2\pi} 
-      F(f,\theta) \:\rd \theta \right)^2 \: \rd  f \: \label{eq:qp}
-      \end{equation}
-\item \textbf{QKK}  wavenumber peakedness \citep{art:DC23}
-      \begin{equation} Q_{kk} = \frac{1}{E^2} \int_0^{f_{NK}} \int_0^{2\pi} 
-      0.5 \left[ A(k,\theta)+ A(k,\theta+\pi)\right]^2 \frac{\sigma^2}{k C_g} \:\rd \theta  \: \rd  \sigma \: \label{eq:qkk}
+      \begin{equation} Q_p = \frac{2}{E^2} \int_0^{2\pi} \int_0^\infty
+      \sigma\:F(\sigma,\theta)^2\:d\sigma\:d\theta \: \label{eq:qp}
       \end{equation}
 \end{list}
 

manual/impl/switch.tex

@@ -94,7 +94,7 @@ \subsubsection{~Mandatory switches} \label{sub:man_switch}
 Selection of nonlinear interactions:
 \begin{slist}
 \sit{nl0} {No nonlinear interactions used.}
-\sit{nl1} {Discrete interaction approximation (\dia) or Gaussian Quadrature Method (GQM).}
+\sit{nl1} {Discrete interaction approximation (\dia).}
 \sit{nl2} {Exact interaction approximation (\xnl).}
 \sit{nl3} {Generalized Multiple \dia\ (\gmd).}
 \sit{nl4} {Two-scale approximation (TSA).} 

manual/manual.bib

@@ -3654,44 +3654,3 @@ @article{vanVledder2006
 	volume = {53},
 	year = {2006}
 }
-
-@article{art:DC23,
-	author = {De Carlo, Marine and Fabrice Ardhuin and Annabelle Ollivier and Adrien Nigou },
-	journal = {Journal of Geophysical Research - Oceans},
-	keywords = {wave groups,altimetry},
-	pages = {},
-	title = {Wave groups and small scale variability of wave heights observed by altimeters},
-	volume = {},
-	year = {2023}
-}
-
-@ARTICLE{Webb1978,
-   author = "D. J. Webb",
-   title = "Nonlinear transfer between sea waves",
-   journal = DSR,
-   volume = 25,
-   pages = "279--298",
-   year = 1978,
-   where="paper",
-}
-
-
-@ARTICLE{Lavrenov2001,
-   author = "Igor V. Lavrenov",
-   title = "Effect of wind wave parameter fluctuation on the nonlinear spectrum evolution",
-   journal = JPO,
-   volume = 31,
-   pages = "861--873",
-   year = 2001,
-   url="http://ams.allenpress.com/archive/1520-0485/31/4/pdf/i1520-0485-31-4-861",
-   keywords={4-wave interactions,GQM},
-}
-
-
-@PHDTHESIS{Gagnaire-Renou2009,
-   author = "Elodie Gagnaire-Renou",
-   title = "Amelioration de la modelisation spectrale des etats de mer par un calcul quasi-exact des interactions non-lineaires vague-vague",
-   school = "Universit{\'e} du Sud Toulon Var",
-   year = 2010,
-}
-

manual/sys/files_w3.tex

@@ -506,15 +506,11 @@ \subsubsection{~Wave model modules} \label{sec:wave_mod}
 \end{flist}
 
 \noindent
-Nonlinear interaction module (\dia or GQM) \hfill {\file w3snl1md.ftn}
+Nonlinear interaction module (\dia) \hfill {\file w3snl1md.ftn}
 
 \begin{flisti}
 \fit{w3snl1}{Calculation of $S_{nl}$.}
 \fit{insnl1}{Initialization for $S_{nl}$.}
-\fit{w3snlgqm}{Calculation of $S_{nl}$.}
-\fit{w3scouple}{Calculation of coupling coefficient.}
-\fit{gauleg}{Calculation of Gauss-Legendre quadrature coefficients.}
-\fit{INSNLGQM}{Initialization for $S_{nl}$ with GQ method.}
 \end{flisti}
 
 \noindent

model/bin/build_utils.sh

@@ -855,23 +855,15 @@ create_file_list()
             sourcet="$pdlibcode $pdlibyow $db $bt $setupcode w3triamd $stx $nlx $btx $is w3parall $uostmd"
                  IO="w3iogrmd $oasismd $agcmmd $ogcmmd $igcmmd"
                 aux="constants w3servmd w3timemd w3arrymd w3dispmd w3gsrumd $tidecode" ;;
-      ww3_shel|ww3_shel_esmf)
-               if [ "$prog" = "ww3_shel" ]
-               then
-                 core=''
-               else
-                 core='wav_kind_mod wav_shr_mod wav_shel_inp wav_comp_nuopc wav_import_export w3iogoncdmd'
-               fi
-               core="$core w3fldsmd w3initmd w3wavemd w3wdasmd w3updtmd"
-               core="$core wminitmd wmwavemd wmfinlmd wmgridmd wmupdtmd wminiomd"
+      ww3_shel)
+               core='w3fldsmd w3initmd w3wavemd w3wdasmd w3updtmd'
                data="wmmdatmd $memcode w3gdatmd w3wdatmd w3adatmd w3idatmd w3odatmd"
                prop="$pr $smcm"
             sourcet="$pdlibcode $setupcode w3triamd w3srcemd $dsx $flx $ln $st $nl $bt $ic"
             sourcet="$sourcet $is $db $tr $bs $refcode $igcode w3parall $uostmd"
                  IO="w3iogrmd w3iogomd w3iopomd w3iotrmd w3iorsmd w3iobcmd $oasismd $agcmmd $ogcmmd $igcmmd"
-                 IO="$IO w3iosfmd w3partmd wmiopomd"
+                 IO="$IO w3iosfmd w3partmd"
                 aux="constants w3servmd w3timemd $tidecode w3arrymd w3dispmd w3cspcmd w3gsrumd"
-                aux="$aux wmunitmd w3nmlmultimd"
                 aux="$aux w3nmlshelmd $pdlibyow" ;;
       ww3_multi|ww3_multi_esmf)
                if [ "$prog" = "ww3_multi" ]

model/bin/cmplr.env

@@ -233,7 +233,7 @@ fi
 ###############################
 
 
-if [ "$cmplr" == "gnu" ] || [ "$cmplr" == "gnu_debug" ] || [ "$cmplr" == "gnu_prof" ] || \
+if [ "$cmplr" == "gnu" ] || [ "$cmplr" == "gnu_debug" ] || [ "$cmplr" == "gnu_prof" ] || [ "$cmplr" == "gfortran" ] || \
    [ "$cmplr" == "so_gnu" ] || [ "$cmplr" == "so_gnu_debug" ] || [ "$cmplr" == "so_gnu_prof" ] || \
    [ "$cmplr" == "hera.gnu" ] || [ "$cmplr" == "cheyenne.gnu" ] || \
    [ "$cmplr" == "datarmor_gnu" ] || [ "$cmplr" == "datarmor_gnu_debug" ] || [ "$cmplr" == "datarmor_gnu_prof" ] || \