Add iwp_subsidence to column type.

Adds ground subsidence variable to col_pp in column type to track subsidence associated with excess ice melting in polygonal tundra.
lanl · May 15, 2024 · 57717bf · 57717bf
1 parent 84d59fc
commit 57717bf
Showing 1 changed file with 15 additions and 14 deletions.
diff --git a/components/elm/src/data_types/ColumnType.F90 b/components/elm/src/data_types/ColumnType.F90
@@ -49,18 +49,19 @@ module ColumnType
      logical , pointer :: active       (:) => null() ! true=>do computations on this column
 
      ! topography
-     real(r8), pointer :: glc_topo     (:) => null() ! surface elevation (m)
-     real(r8), pointer :: micro_sigma  (:) => null() ! microtopography pdf sigma (m)
-     real(r8), pointer :: n_melt       (:) => null() ! SCA shape parameter
-     real(r8), pointer :: topo_slope   (:) => null() ! gridcell topographic slope
-     real(r8), pointer :: topo_std     (:) => null() ! gridcell elevation standard deviation
-     real(r8), pointer :: hslp_p10     (:,:) => null() ! hillslope slope percentiles (unitless)
-     integer, pointer  :: nlevbed      (:) => null() ! number of layers to bedrock
-     real(r8), pointer :: zibed        (:) => null() ! bedrock depth in model (interface level at nlevbed)
-     real(r8), pointer :: iwp_microrel (:) => null() ! ice wedge polygon microtopographic relief (m)
-     real(r8), pointer :: iwp_exclvol  (:) => null() ! ice wedge polygon excluded volume (m)
-     real(r8), pointer :: iwp_ddep     (:) => null() ! ice wedge polygon depression depth (m)
-     real(r8), pointer :: meangradz    (:) => null() ! mean topographic gradient at the column level
+     real(r8), pointer :: glc_topo      (:) => null() ! surface elevation (m)
+     real(r8), pointer :: micro_sigma   (:) => null() ! microtopography pdf sigma (m)
+     real(r8), pointer :: n_melt        (:) => null() ! SCA shape parameter
+     real(r8), pointer :: topo_slope    (:) => null() ! gridcell topographic slope
+     real(r8), pointer :: topo_std      (:) => null() ! gridcell elevation standard deviation
+     real(r8), pointer :: hslp_p10      (:,:) => null() ! hillslope slope percentiles (unitless)
+     integer, pointer  :: nlevbed       (:) => null() ! number of layers to bedrock
+     real(r8), pointer :: zibed         (:) => null() ! bedrock depth in model (interface level at nlevbed)
+     real(r8), pointer :: iwp_microrel  (:) => null() ! ice wedge polygon microtopographic relief (m)
+     real(r8), pointer :: iwp_exclvol   (:) => null() ! ice wedge polygon excluded volume (m)
+     real(r8), pointer :: iwp_ddep      (:) => null() ! ice wedge polygon depression depth (m)
+     real(r8), pointer :: iwp_subsidence(:) => null() ! ice wedge polygon ground subsidence (m)
+     real(r8), pointer :: meangradz     (:) => null() ! mean topographic gradient at the column level
 
      ! vertical levels
      integer , pointer :: snl          (:)   => null() ! number of snow layers
@@ -140,7 +141,7 @@ subroutine col_pp_init(this, begc, endc)
 
     ! Assume that columns are not fates columns until fates initialization begins
     allocate(this%is_fates(begc:endc)); this%is_fates(:) = .false.
-    
+
   end subroutine col_pp_init
 
   !------------------------------------------------------------------------
@@ -179,7 +180,7 @@ subroutine col_pp_clean(this)
     deallocate(this%zibed      )
     deallocate(this%hydrologically_active)
     deallocate(this%is_fates)
-    
+
   end subroutine col_pp_clean
 
 end module ColumnType