From b64b61b8374a5358b49e174ba2a5060f0bc28de1 Mon Sep 17 00:00:00 2001 From: Shivaprakash Muruganandham Date: Thu, 15 Feb 2024 09:06:20 -0800 Subject: [PATCH] Update draft dependence function for basalMassBalAdjustment Add slopeGen and interceptGen coefficients for the draft dependence --- .../src/mode_forward/mpas_li_iceshelf_melt.F | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_iceshelf_melt.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_iceshelf_melt.F index 4c769b6f14b7..0cab6eb3d937 100644 --- a/components/mpas-albany-landice/src/mode_forward/mpas_li_iceshelf_melt.F +++ b/components/mpas-albany-landice/src/mode_forward/mpas_li_iceshelf_melt.F @@ -681,8 +681,8 @@ subroutine calculate_aislens_melt_variability_adjustment( & ! slopeSer and interceptSer variables: should they be defined as ! input or local variables for this subroutine? - real (kind=RKIND) :: slopeSer ! slope of relation between depth and melt rate - real (kind=RKIND) :: interceptSer ! depth at which melting goes to 0 + real (kind=RKIND), pointer :: config_basal_mass_bal_adj_aislens_slopeGen ! slope of relation between depth and melt rate + real (kind=RKIND), pointer :: config_basal_mass_bal_adj_aislens_interceptGen ! intercept for draft dependence parameterization of basal melt melting !real (kind=RKIND) :: maxMeltSer ! maximum allowable melt rate !real (kind=RKIND) :: sillDepth ! depth below which melt rate no longer increases !real (kind=RKIND), pointer :: config_basal_mass_bal_seroussi_amplitude @@ -698,11 +698,16 @@ subroutine calculate_aislens_melt_variability_adjustment( & !Input field values for slopeSer(iCell) and interceptSer(iCell) !from file + call mpas_pool_get_config(liConfigs, 'config_basal_mass_bal_adj_aislens_slopeGen') !, & + !config_basal_mass_bal_seroussi_amplitude) ! meters + call mpas_pool_get_config(liConfigs, 'config_basal_mass_bal_adj_aislens_interceptGen') !, & + ! config_basal_mass_bal_seroussi_period + do iCell = 1, nCellsSolve zDraft(iCell) = lowerSurface(iCell) - config_sea_level - floatingBasalMassBalAdjustment(iCell) = interceptSer(iCell) + & - zDraft(iCell) * slopeSer(iCell) + floatingBasalMassBalAdjustment(iCell) = config_basal_mass_bal_adj_aislens_interceptGen(iCell) + & + zDraft(iCell) * config_basal_mass_bal_adj_aislens_slopeGen(iCell) enddo ! iCell ! Optional if below is required within subroutine: