Add MALI buildnml support for a TF-coupling parameters file

When OCN/GLC coupling is active, MALI requires a separate input file
with parameters for the ice-shelf basal melt parameterization.  These
parameters are specific to each MALI mesh.  This commit adds MALI
buildnml support for this additional input file if TF coupling is active
AND the params file has been set up for a given MALI mesh.  For MALI
meshes that do not have a param file set up, ice-shelf basal melting
will be zero.

Note that in the future, these parameters may also be a function of OCN
mesh and compset.  The needed level of complexity is not yet known, so
this logic will be adjusted to accommodate that when the time comes.

Author: matthewhoffman

components/mpas-albany-landice/cime_config/buildnml

@@ -61,6 +61,7 @@ def buildnml(case, caseroot, compname):
     decomp_prefix = ''
     datamode_date = ''
     region_mask_date = ''
+    tf_params_date = ''
 
     if glc_grid == 'mpas.aisgis20km':
         grid_date += '20190326'
@@ -79,12 +80,14 @@ def buildnml(case, caseroot, compname):
         decomp_date += '240507'
         decomp_prefix += 'mpasli.graph.info.'
         region_mask_date += '20231222'
+        tf_params_date += '20250724'
     elif glc_grid == 'mpas.ais4to20km':
         grid_date += '20250625'
         grid_prefix += 'ais_4to20km'
         decomp_date += '240507'
         decomp_prefix += 'mpasli.graph.info.'
         region_mask_date += '20230105'
+        tf_params_date += '20250724'
     elif glc_grid == 'mpas.gis20km':
         grid_date += '20210824'
         grid_prefix += 'gis_20km_r01'
@@ -129,6 +132,20 @@ def buildnml(case, caseroot, compname):
     if region_mask_date != '':
         region_mask_file = f"{din_loc_root}/glc/mpasli/{glc_grid}/{grid_prefix}_region_mask.{region_mask_date}.nc"
 
+    #------------------------
+    # Set TF param file if needed
+    # It is only activated if TF-based ocean coupling is on AND
+    # a param. file has been set up for MALI mesh being used.
+    # Note: in the future, we will likely need to make the params file
+    # a function of also the OCN mesh and compset being used, as they will
+    # require different deltaT bias-correction fields.
+    # It is unclear at this point what level of complexity is required,
+    # so nothing is implemented yet.
+    #------------------------
+    tf_params_file = ''
+    if ocn_glc_ismf_coupling == 'tf' and tf_params_date != '':
+        tf_params_file = f"{din_loc_root}/glc/mpasli/{glc_grid}/{grid_prefix}_tf_params.{tf_params_date}.nc"
+
     #------------------------
     # Set the datamode file if needed
     #------------------------
@@ -145,6 +162,9 @@ def buildnml(case, caseroot, compname):
         if region_mask_date != '':
             input_list.write(f"region_mask = {region_mask_file}\n")
 
+        if tf_params_file != '':
+            input_list.write(f"tf_params = {tf_params_file}\n")
+
         if mali_prognostic_mode == 'DATA':
             input_list.write(f"datamode = {datamode_file}\n")
 
@@ -246,11 +266,16 @@ def buildnml(case, caseroot, compname):
         lines.append('                      input_interval="initial_only"')
         lines.append('                      output_interval="00-00-00_00:00:01"/>')
         lines.append('')
-        lines.append('    <!--')
-        lines.append('    output is the main history output stream. You can add auxiliary streams to')
-        lines.append('    this stream to include more fields.')
-        lines.append('    -->')
-        lines.append('')
+        if tf_params_file != '':
+            lines.append('    <stream name="tfParamsInput"')
+            lines.append('            type="input"')
+            lines.append('            filename_template="{}"'.format(tf_params_file))
+            lines.append('            input_interval="initial_only">')
+            lines.append('        <var name="ismip6shelfMelt_basin"/>')
+            lines.append('        <var name="ismip6shelfMelt_gamma0"/>')
+            lines.append('        <var name="ismip6shelfMelt_deltaT"/>')
+            lines.append('    </stream>')
+            lines.append('')
         if region_mask_date != '':
             lines.append('    <stream name="regionsInput"')
             lines.append('            type="input"')
@@ -259,6 +284,11 @@ def buildnml(case, caseroot, compname):
             lines.append('        <var name="regionCellMasks"/>')
             lines.append('    </stream>')
             lines.append('')
+        lines.append('    <!--')
+        lines.append('    output is the main history output stream. You can add auxiliary streams to')
+        lines.append('    this stream to include more fields.')
+        lines.append('    -->')
+        lines.append('')
         lines.append('    <stream name="output"')
         lines.append('            type="output"')
         lines.append('            filename_template="{}.mali.hist.$Y-$M-$D_$S.nc"'.format(casename))