calc_temp_melt.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  c a l c _ t e m p _ m e l t
!
!> @file
!!
!! Computation of the melting 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 melting temperatures.
!<------------------------------------------------------------------------------
subroutine calc_temp_melt()

use sico_types
use sico_variables
use sico_vars

implicit none
integer(i4b) :: i, j, kc, kt
real(dp) :: atm1(0:kcmax), atm2(0:ktmax)

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

do kc=0, kcmax
   atm1(kc) = beta*(1.0_dp-eaz_c_quotient(kc))
end do

do kt=0, ktmax
   atm2(kt) = 1.0_dp-zeta_t(kt)
end do

!-------- Compute the melting temperatures --------

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

   do kt=0, ktmax
      temp_t_m(kt,j,i) = -beta*(h_c(j,i)+h_t(j,i)*atm2(kt))
   end do

   do kc=0, kcmax
      temp_c_m(kc,j,i) = -atm1(kc)*h_c(j,i)
   end do

end do
end do

end subroutine calc_temp_melt
!