Fix dCalvingThickness dimension in remove_small_islands

Fix dCalvingThickness dimension in remove_small_islands. Must be nCells
rather than nCellsSolve in size.

Author: trhille

components/mpas-albany-landice/src/mode_forward/mpas_li_calving.F

@@ -625,6 +625,7 @@ subroutine li_restore_calving_front(domain, err)
 
          ! Now check for marine regions that have advanced
          allocate(dCalvingThickness(nCells+1))
+         dCalvingThickness(:) = 0.0_RKIND
 
          ! loop over locally owned cells
          do iCell = 1, nCells
@@ -1084,7 +1085,7 @@ subroutine remove_small_islands(meshPool, geometryPool)
                                         floatingMaskForMassBudget
       integer, dimension(:,:), pointer :: cellsOnCell ! list of cells that neighbor each cell
       integer, dimension(:), pointer :: nEdgesOnCell ! number of cells that border each cell
-      integer, pointer :: nCellsSolve
+      integer, pointer :: nCellsSolve, nCells
       integer :: iCell, jCell, n, nIceNeighbors, nIceNeighbors2, neighborWithIce
       integer :: nOpenOceanNeighbors, nOpenOceanNeighbors2
 
@@ -1095,6 +1096,7 @@ subroutine remove_small_islands(meshPool, geometryPool)
 
       call mpas_pool_get_config(liConfigs, 'config_sea_level', config_sea_level)
       call mpas_pool_get_dimension(meshPool, 'nCellsSolve', nCellsSolve)
+      call mpas_pool_get_dimension(meshPool, 'nCells', nCells)
       call mpas_pool_get_array(meshPool, 'cellsOnCell', cellsOnCell)
       call mpas_pool_get_array(meshPool, 'nEdgesOnCell', nEdgesOnCell)
       call mpas_pool_get_array(geometryPool, 'calvingThickness', calvingThickness)
@@ -1107,7 +1109,7 @@ subroutine remove_small_islands(meshPool, geometryPool)
       call mpas_pool_get_array(geometryPool, 'cellMask', cellMask)
       call mpas_pool_get_array(geometryPool, 'bedTopography', bedTopography)
 
-      allocate(dCalvingThickness(nCellsSolve+1))
+      allocate(dCalvingThickness(nCells+1))
       dCalvingThickness(:) = 0.0_RKIND
 
       do iCell = 1, nCellsSolve