add FATES hydro solver selection switch to namelist

This commit facilitates moving this switch from the FATES parameter file
to the namelist

Author: glemieux

bld/CLMBuildNamelist.pm

@@ -811,7 +811,7 @@ sub setup_cmdl_fates_mode {
                       "flandusepftdat","use_fates_potentialveg","use_fates_lupft","fates_history_dimlevel",
                       "use_fates_daylength_factor", "use_fates_photosynth_acclimation", "fates_stomatal_model",
                       "fates_stomatal_assimilation", "fates_leafresp_model", "fates_cstarvation_model",
-                      "fates_regeneration_model"
+                      "fates_regeneration_model", "fates_hydro_solver"
                    );
 
        # dis-allow fates specific namelist items with non-fates runs
@@ -4709,7 +4709,7 @@ sub setup_logic_fates {
                        "fates_harvest_mode","fates_parteh_mode", "use_fates_cohort_age_tracking","use_fates_tree_damage",
                        "use_fates_daylength_factor", "use_fates_photosynth_acclimation", "fates_stomatal_model",
                        "fates_stomatal_assimilation", "fates_leafresp_model", "fates_cstarvation_model",
-                       "fates_regeneration_model"
+                       "fates_regeneration_model", "fates_hydro_solver"
                     );
 
         foreach my $var ( @list ) {

bld/namelist_files/namelist_defaults_ctsm.xml

@@ -2394,6 +2394,7 @@ lnd/clm2/surfdata_esmf/NEON/ctsm5.3.0/surfdata_1x1_NEON_TOOL_hist_2000_78pfts_c2
 <fates_leafresp_model          use_fates=".true.">ryan1991</fates_leafresp_model>
 <fates_cstarvation_model       use_fates=".true.">linear</fates_cstarvation_model>
 <fates_regeneration_model      use_fates=".true.">default</fates_regeneration_model>
+<fates_hydro_solver            use_fates=".true.">1D_Taylor</fates_hydro_solver>
 <use_fates_planthydro          use_fates=".true.">.false.</use_fates_planthydro>
 <use_fates_tree_damage         use_fates=".true.">.false.</use_fates_tree_damage>
 <use_fates_cohort_age_tracking use_fates=".true.">.false.</use_fates_cohort_age_tracking>

bld/namelist_files/namelist_definition_ctsm.xml

@@ -775,6 +775,11 @@ Allowed values are:
 This option is older than the luhdata options and may be depricated at some point in the future.
 </entry>
 
+<entry id="fates_hydro_solver" type="char*256" category="physics"
+        group="clm_inparm" valid_values="1D_Taylor, 2D_Picard, 2D_Newton" value="1D_Taylor">
+Set the FATES hydro solver method
+</entry>
+
 <entry id="fates_regeneration_model" type="char*256" category="physics"
         group="clm_inparm" valid_values="default, trs, trs_no_seed_dyn" value="default">
 Set the FATES seed regeneration model

src/main/clm_varctl.F90

@@ -336,6 +336,7 @@ module clm_varctl
   character(len=256), public :: fates_leafresp_model = ''               ! Leaf maintenance respiration model, Ryan or Atkin
   character(len=256), public :: fates_cstarvation_model = ''            ! linear or exponential function
   character(len=256), public :: fates_regeneration_model = ''           ! default, TRS, or TRS without seed dynamics
+  character(len=256), public :: fates_hydro_solver = ''                 ! 1D Taylor, 2D Picard, 2D Newton
   logical, public            :: use_fates_planthydro = .false.          ! true => turn on fates hydro
   logical, public            :: use_fates_cohort_age_tracking = .false. ! true => turn on cohort age tracking
   logical, public            :: use_fates_ed_st3   = .false.            ! true => static stand structure

src/main/controlMod.F90

@@ -247,6 +247,7 @@ subroutine control_init(dtime)
           fates_leafresp_model,                         &
           fates_cstarvation_model,                      &
           fates_regeneration_model,                     &
+          fates_hydro_solver,                           &
           fates_parteh_mode,                            &
           fates_seeddisp_cadence,                       &
           use_fates_tree_damage,                        &
@@ -811,6 +812,7 @@ subroutine control_spmd()
     call mpi_bcast (fates_leafresp_model, len(fates_leafresp_model) , MPI_CHARACTER, 0, mpicom, ier)
     call mpi_bcast (fates_cstarvation_model, len(fates_cstarvation_model) , MPI_CHARACTER, 0, mpicom, ier)
     call mpi_bcast (fates_regeneration_model, len(fates_regeneration_model) , MPI_CHARACTER, 0, mpicom, ier)
+    call mpi_bcast (fates_hydro_solver, len(fates_hydro_solver) , MPI_CHARACTER, 0, mpicom, ier)
     call mpi_bcast (use_fates_planthydro, 1, MPI_LOGICAL, 0, mpicom, ier)
     call mpi_bcast (use_fates_tree_damage, 1, MPI_LOGICAL, 0, mpicom, ier)
     call mpi_bcast (use_fates_cohort_age_tracking, 1, MPI_LOGICAL, 0, mpicom, ier)
@@ -1219,6 +1221,7 @@ subroutine control_print ()
        write(iulog, *) '    fates_leafresp_model = ', fates_leafresp_model
        write(iulog, *) '    fates_cstarvation_model = ', fates_cstarvation_model
        write(iulog, *) '    fates_regeneration_model = ', fates_regeneration_model
+       write(iulog, *) '    fates_hydro_solver = ', fates_hydro_solver
        write(iulog, *) '    fates_paramfile = ', fates_paramfile
        write(iulog, *) '    fates_parteh_mode = ', fates_parteh_mode
        write(iulog, *) '    use_fates_planthydro = ', use_fates_planthydro

src/utils/clmfates_interfaceMod.F90

@@ -62,6 +62,7 @@ module CLMFatesInterfaceMod
    use clm_varctl        , only : fates_leafresp_model
    use clm_varctl        , only : fates_cstarvation_model
    use clm_varctl        , only : fates_regeneration_model
+   use clm_varctl        , only : fates_hydro_solver
    use clm_varctl        , only : use_fates_inventory_init
    use clm_varctl        , only : use_fates_fixed_biogeog
    use clm_varctl        , only : use_fates_nocomp
@@ -413,6 +414,7 @@ subroutine CLMFatesGlobals2()
      integer                                        :: pass_leafresp_model
      integer                                        :: pass_cstarvation_model
      integer                                        :: pass_regeneration_model
+     integer                                        :: pass_hydro_solver
 
      call t_startf('fates_globals2')
 
@@ -580,6 +582,15 @@ subroutine CLMFatesGlobals2()
         end if
         call set_fates_ctrlparms('regeneration_model',ival=pass_regeneration_model)
 
+        if (trim(fates_hydro_solver) == '1D_Taylor') then
+           pass_hydro_solver = 1
+        else if (trim(fates_hydro_solver) == '2D_Picard') then
+           pass_hydro_solver = 2
+        else if (trim(fates_hydro_solver) == '2D_Taylor') then
+           pass_hydro_solver = 3
+        end if
+        call set_fates_ctrlparms('hydr_solver',ival=pass_hydro_solver)
+
         ! FATES logging and harvest modes
         pass_logging = 0
         pass_lu_harvest = 0