calc_thk_water_bas_m.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Module :  c a l c _ t h k _ w a t e r _ b a s _ m
!
!> @file
!!
!! Computation of the thickness of the water column under the ice base.
!!
!! @section Copyright
!!
!! Copyright 2009-2017 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/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Computation of the thickness of the water column under the ice base.
!<------------------------------------------------------------------------------
module calc_thk_water_bas_m

  use sico_types_m
  use sico_variables_m
  use sico_vars_m

#if (BASAL_HYDROLOGY==1)
  use hydro_m
#endif

  implicit none

  private
  public :: calc_thk_water_bas

contains

!-------------------------------------------------------------------------------
!> Main subroutine of calc_thk_water_bas_m:
!! Computation of the thickness of the water column under the ice base.
!<------------------------------------------------------------------------------
  subroutine calc_thk_water_bas(z_sl)

  implicit none

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

  logical, save :: firstcall = .true.

#if (BASAL_HYDROLOGY==1)
  real(dp), save                     :: rho_rho_w_ratio
  integer , dimension(0:IMAX,0:JMAX) :: hydro_icemask
  real(dp), dimension(0:IMAX,0:JMAX) :: hydro_topg, hydro_thk, &
                                        hydro_temppabase, hydro_supply, &
                                        hydro_sflux, hydro_bwat
  type(hydro_t), save                :: hydro
                                        !!! Does this need a save attribute?
#endif

!-------- Water column --------

#if (BASAL_HYDROLOGY==1)

  if (firstcall) then

     rho_rho_w_ratio = rho/rho_w

     call hydro_init(hydro, xi, eta)
     call hydro_gen_conf(hydro, &
          & method='quinn', &
          & avoid_frz=.false., &
          & filter_len=0.0_dp, &
          & rho_seawater=rho_sw, &
          & rho_freshwater=rho_w, &
          & rho_ice=rho)

  end if

  hydro_topg       = transpose(zl)-z_sl
  hydro_temppabase = transpose(temph_b)

  where (transpose(maske)==0_i2b)   ! grounded ice
     hydro_icemask = 1
     hydro_thk     = transpose(h_c+h_t)
     hydro_supply  = rho_rho_w_ratio*transpose(q_b_tot)
  elsewhere
     hydro_icemask = 0
     hydro_thk     = 0.0_dp
     hydro_supply  = 0.0_dp
  end where

  call hydro_set_topg(hydro, hydro_topg)
  call hydro_set_thk(hydro, hydro_thk)
  call hydro_set_temppabase(hydro, hydro_temppabase)
  call hydro_set_supply(hydro, hydro_supply)
  call hydro_set_mask(hydro, hydro_icemask)

  call hydro_update(hydro)

  call hydro_get_sflux(hydro, hydro_sflux)
  call hydro_get_bwat(hydro, hydro_bwat)

  h_w = transpose(hydro_bwat)

#else
  where (maske==0_i2b) h_w = 0.0_dp   ! grounded ice
#endif

  where (maske==2_i2b)   ! ocean
     h_w = z_sl-zl
  elsewhere (maske==3_i2b)   ! floating ice
     h_w = zb-zl
  elsewhere (maske==1_i2b)   ! ice-free land
     h_w = 0.0_dp
  end where

  if (firstcall) firstcall = .false.

  end subroutine calc_thk_water_bas

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

end module calc_thk_water_bas_m
!