calc_vz_static.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  c a l c _ v z _ s t a t i c
!
!> @file
!!
!! Computation of the vertical velocity vz for static ice.
!!
!! @section Copyright
!!
!! Copyright 2014-2016 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 vertical velocity vz for static ice.
!<------------------------------------------------------------------------------
subroutine calc_vz_static()

use sico_types
use sico_variables
use sico_vars

implicit none

integer(i4b) :: i, j, kc, kt

do i=0, imax
do j=0, jmax

   if ((maske(j,i) == 0).or.(maske(j,i) == 3)) then   ! grounded or floating ice

      vz_b(j,i) = dzb_dtau(j,i)-q_b_tot(j,i)
                  ! kinematic boundary condition at the ice base

      do kt=0, ktmax-1
         vz_t(kt,j,i) = vz_b(j,i)
      end do

      vz_m(j,i) = vz_b(j,i)

      do kc=0, kcmax-1
         vz_c(kc,j,i) = vz_b(j,i)
      end do

      vz_s(j,i) = vz_b(j,i)

   else   ! maske(j,i) == (1 or 2)

      vz_b(j,i) = 0.0_dp

      do kt=0, ktmax-1
         vz_t(kt,j,i) = 0.0_dp
      end do

      vz_m(j,i) = 0.0_dp

      do kc=0, kcmax-1
         vz_c(kc,j,i) = 0.0_dp
      end do

      vz_s(j,i) = 0.0_dp

   end if

end do
end do

end subroutine calc_vz_static
!