mask_update.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Function :  m a s k _ u p d a t e
!
!> @file
!!
!! Update of the topography mask due to changes of the sea level.
!!
!! @section Copyright
!!
!! Copyright 2009-2013 Ralf Greve
!!
!! @section License
!!
!! This file is part of SICOPOLIS.
!!
!! SICOPOLIS is free software: you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation, either version 3 of the License, or
!! (at your option) any later version.
!!
!! SICOPOLIS is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with SICOPOLIS.  If not, see <http://www.gnu.org/licenses/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Update of the topography mask due to changes of the sea level.
!<------------------------------------------------------------------------------
function mask_update(z_sl, i, j)

use sico_types
use sico_variables

implicit none
integer(i2b)             :: mask_update
integer(i4b), intent(in) :: i, j
real(dp), intent(in)     :: z_sl

real(dp) :: rhosw_rho_ratio, h_ice, h_sea

rhosw_rho_ratio = rho_sw/rho

!-------- Previously ice-free land point or sea point --------

if ( (maske(j,i) == 1_i2b).or.(maske(j,i) == 2_i2b) ) then

   if (zl(j,i) > z_sl) then
      mask_update = 1_i2b   ! now ice-free land
      return
   else
      mask_update = 2_i2b   ! now sea point
      return
   end if

!-------- Previously glaciated land point or floating-ice point!

else   ! (maske(j,i) == 0_i2b, 3_i2b)

   if (zl(j,i) > z_sl) then

      mask_update = 0_i2b   ! now grounded ice
      return

   else

      h_ice = zs(j,i)-zb(j,i)   ! ice thickness
      h_sea = z_sl   -zl(j,i)   ! sea depth

      if ( h_ice < (rhosw_rho_ratio*h_sea) ) then

#if ( MARGIN==1 || ( MARGIN==2 && MARINE_ICE_FORMATION==1 ) )
         mask_update = 2_i2b     ! ice becomes floating, therefore
                                 ! now sea point (ice cut off)
#elif ( MARGIN==2 && MARINE_ICE_FORMATION==2 )
         mask_update = 0_i2b     ! now "underwater ice"
#elif ( MARGIN==3 )
         mask_update = 3_i2b     ! now floating ice
#endif
         return

      else

         mask_update = 0_i2b     ! now grounded ice
         return

      end if

   end if

end if

end function mask_update
!