<!-- ******************************* -->
<!-- Standard ocean model parameters -->
<!--                                 -->
<!-- ******************************* -->

<!-- Parameters that are repeated in coupledmodel_params.xml will be
     overwritten when using run_coupled. -->
<ParameterList name="Ocean parameters">

  <!-- Specify whether the ocean model should initialize with -->
  <!-- an old state (true) or with a zero solution (false).                  -->
  <!-- In a coupled situation this will be partly overruled by the parameters-->
  <!-- specified in coupledmodel_params.xml                                  -->
  <Parameter name="Load state" type="bool" value="false"/>
  <Parameter name="Save state" type="bool" value="true"/>
  <Parameter name="Input file"  type="string" value="ocean_input.h5" />
  <Parameter name="Output file" type="string" value="ocean_output.h5" />

  <!-- Parameters that affect THCM { -->
  <ParameterList name="THCM">
    <!-- a descriptive name to identify the settings -->
    <Parameter name="Problem Description" type="string" value="Test"/>
    <!-- western domain bound -->
    <Parameter name="Global Bound xmin" type="double" value="0"/>
    <!-- eastern domain bound -->
    <Parameter name="Global Bound xmax" type="double" value="359.99"/>
    <!-- southern domain bound -->
    <Parameter name="Global Bound ymin" type="double" value="-85.5"/>
    <!-- northern domain bound -->
    <Parameter name="Global Bound ymax" type="double" value="85.5"/>
    <!-- periodic boundary conditions in x-direction -->
    <Parameter name="Periodic" type="bool" value="1"/>

    <!-- Depth  -->
    <Parameter name="Depth hdim" type="double" value="5000"/>

    <!-- grid stretching in z-direction (defaults to 1.0)          -->
    <!-- I think the preconditioner does not handle this correctly,-->
    <!-- but in some cases it may work                             -->
    <Parameter name="Grid Stretching qz" type="double" value="2"/>
    <!-- global number of grid points in x-direction -->
    <Parameter name="Global Grid-Size n" type="int" value="96"/>
    <!-- global number of grid points in y-direction -->
    <Parameter name="Global Grid-Size m" type="int" value="38"/>
    <!-- global number of grid points in z-direction -->
    <Parameter name="Global Grid-Size l" type="int" value="12"/>

    <!-- Supply 0-based surface coordinates for the integral condition
         row, default is n-1, m-1. Make sure that the coordinate is
         located on the last processors. -->
    <!-- <Parameter name="Integral row coordinate i" type="int" value="47"/> -->
    <!-- <Parameter name="Integral row coordinate j" type="int" value="16"/> -->
    
    <!-- remove bottom profile? (flat in THCM)      -->
    <Parameter name="Flat Bottom" type="bool" value="1"/>
    <!-- read mask from a file in mkmask The filename is hard-coded in   -->
    <!-- topo.F90.                                                       -->
    <Parameter name="Read Land Mask" type="bool" value="1"/>
    <!-- name of land mask file (assumed to be in topdir/data/mkmask/)   -->
    <Parameter name="Land Mask" type="string" value="mask_global_96x38x12"/>

    <!-- ==================== Coupling ================================= -->
    <!-- Flag enabling the coupling with an atmosphere                   -->
    <!--   0: standalone thcm, 1: accepting external atmosphere          -->
    <Parameter name="Coupled Temperature"  type="int" value="1"/>
    <Parameter name="Coupled Salinity"     type="int" value="0"/>
    <!-- =============================================================== -->

    <!-- Starting Parameters                                          -->
    <!-- coupledmodel_params can overwrite these when running coupled -->
    <ParameterList name="Starting Parameters">
      <Parameter name="Combined Forcing"    type="double" value="0.0"/>
      <Parameter name="Solar Forcing"       type="double" value="1.0"/>
      <Parameter name="Salinity Forcing"    type="double" value="0.0"/>
      <Parameter name="Wind Forcing"        type="double" value="1.0"/>
      <Parameter name="Temperature Forcing" type="double" value="1.0"/>

      <!-- <Parameter name="Rayleigh-Number"     type="double" value="0.0"/> -->
      <Parameter name="Horizontal Ekman-Number" type="double" value="0.0027037"/>

      <Parameter name="SPL1"                type="double" value="1.0"/>
      <Parameter name="SPL2"                type="double" value="0.01"/>

    </ParameterList>

    <!-- Mixing parameter (vmix_GLB in THCM)                                   -->
    <Parameter name="Mixing" type="int" value="1"/>

    <!-- Fred's convective adjustment (rho_mixing in THCM) -->
    <Parameter name="Rho Mixing" type="bool" value="0"/>
    <!-- Taper (tap in usr.F90) -->
    <Parameter name="Taper" type="int" value="1"/>
    <!-- Restoring Temperature-Profile (TRES in THCM) -->
    <!-- 0: non-restoring                             -->
    <!-- 1: restoring                                 -->
    <Parameter name="Restoring Temperature Profile" type="int" value="1"/>
    
    <!-- Restoring Salt-Profile (SRES in THCM)        -->
    <!-- 0: non-restoring                             -->
    <!-- 1: restoring                                 -->
    <Parameter name="Restoring Salinity Profile" type="int" value="1"/>
    
    <!-- Temperature forcing                        -->
    <!-- Levitus Temperature (ite in THCM)          -->
    <!-- 0: data, 1: idealized                      -->
    <Parameter name="Levitus T" type="int" value="1"/>
    
    <!-- Levitus Salt (its in THCM)  -->
    <!-- 0: data, 1: idealized       -->
    <Parameter name="Levitus S" type="int" value="1"/>
    
    <!-- Wind Forcing (iza in THCM)       -->
    <!-- 0: data,                         -->
    <!-- 1: zonally averaged              -->
    <!-- 2: idealized                     -->
    <Parameter name="Wind Forcing Type" type="int" value="2"/>
    
    <!-- Type of scaling applied to the linear systems.          -->
    <!-- We currently support "None" and "THCM"                  -->
    <!-- "None" is not really recomended.                        -->
    <Parameter name="Scaling" type="string" value="THCM"/>

  </ParameterList> <!-- } THCM -->

</ParameterList> <!-- } -->
