calving_underwater_ice_m.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Module :  c a l v i n g _ u n d e r w a t e r _ i c e _ m
!
!> @file
!!
!! Calving of "underwater ice".
!!
!! @section Copyright
!!
!! Copyright 2009-2017 Ralf Greve, Thorben Dunse
!!
!! @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/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Calving of "underwater ice".
!<------------------------------------------------------------------------------
module calving_underwater_ice_m

  use sico_types_m
  use sico_variables_m
  use sico_vars_m

  implicit none

  real(dp), dimension(0:JMAX,0:IMAX) :: calv_uw_ice

  public

contains

!-------------------------------------------------------------------------------
!> Main routine: Calving of "underwater ice".
!<------------------------------------------------------------------------------
  subroutine calving_underwater_ice(z_sl)

  implicit none

  real(dp), intent(in) :: z_sl

  real(dp)                           :: year_sec_inv
  real(dp)                           :: rhosw_rho_ratio
  real(dp)                           :: calv_uw_coeff, r1_calv_uw, r2_calv_uw
  real(dp), dimension(0:JMAX,0:IMAX) :: H, H_sea

!-------- Term abbreviations --------

  year_sec_inv = 1.0_dp/year_sec

  rhosw_rho_ratio = rho_sw/rho

!-------- Setting of parameters --------

#if (defined(CALV_UW_COEFF))
  calv_uw_coeff = calv_uw_coeff * year_sec_inv
#else
  stop ' >>> calving_underwater_ice: CALV_UW_COEFF undefined!'
#endif

#if (defined(R1_CALV_UW))
  r1_calv_uw = r1_calv_uw
#else
  stop ' >>> calving_underwater_ice: R1_CALV_UW undefined!'
#endif

#if (defined(R2_CALV_UW))
  r2_calv_uw = r2_calv_uw
#else
  stop ' >>> calving_underwater_ice: R2_CALV_UW undefined!'
#endif

!-------- Ice thickness and sea depth --------

  h     = max(h_c + h_t, 0.0_dp)   ! ice thickness
  h_sea = max(z_sl - zl, 0.0_dp)   ! sea depth

!-------- Calving of "underwater ice" --------

  where ( (maske == 0_i2b).and.(h < rhosw_rho_ratio*h_sea) )
     calv_uw_ice = calv_uw_coeff * h**r1_calv_uw * h_sea**r2_calv_uw
  elsewhere
     calv_uw_ice = 0.0_dp
  end where

  end subroutine calving_underwater_ice

!-------------------------------------------------------------------------------

end module calving_underwater_ice_m
!