calc_temp_bas.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  c a l c _ t e m p _ b a s
!
!> @file
!!
!! Computation of the basal temperatures.
!!
!! @section Copyright
!!
!! Copyright 2009-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 basal temperatures.
!<------------------------------------------------------------------------------
subroutine calc_temp_bas()

use sico_types
use sico_variables
use sico_vars

implicit none
integer(i4b) :: i, j

!-------- Computation of the basal temperatures --------

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

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

      if (n_cts(j,i) == -1) then   ! cold ice base

         temp_b(j,i)  = temp_c(0,j,i)
         temph_b(j,i) = temp_c(0,j,i) - temp_c_m(0,j,i)
                        ! relative to the pressure melting point

      else   ! n_cts(j,i) == 0 or 1, temperate ice base

         temp_b(j,i)  = temp_t_m(0,j,i)
         temph_b(j,i) = 0.0_dp
                        ! relative to the pressure melting point

      end if

   else   ! maske(j,i) == 1 or 2, ice-free land or sea

      temp_b(j,i)  = temp_c(0,j,i)
      temph_b(j,i) = temp_c(0,j,i) - temp_c_m(0,j,i)
                     ! relative to the pressure melting point

   end if

end do
end do

end subroutine calc_temp_bas
!