phys_para.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  p h y s _ p a r a
!
!> @file
!!
!! Reading of physical parameters.
!!
!! @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/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Reading of physical parameters.
!<------------------------------------------------------------------------------
subroutine phys_para()

use sico_types
use sico_variables
use sico_vars

implicit none
integer(i4b), parameter :: n_unit=31
integer(i4b) :: ios
integer(i4b) :: n

open(n_unit, iostat=ios, &
     file=inpath//'/'//trim(ch_domain_short)//'/'//phys_para_file, &
     status='old')
if (ios /= 0) stop ' phys_para: Error when opening the phys_para file!'

call read_phys_para_value(rho_i, n_unit)
call read_phys_para_value(rho_w, n_unit)
call read_phys_para_value(rho_sw, n_unit)
call read_phys_para_value(l, n_unit)
call read_phys_para_value(g, n_unit)
call read_phys_para_value(nue, n_unit)
call read_phys_para_value(beta, n_unit)
call read_phys_para_value(omega_max, n_unit)
call read_phys_para_value(rho_c, n_unit)
call read_phys_para_value(kappa_c, n_unit)
call read_phys_para_value(c_c, n_unit)
call read_phys_para_value(h_r, n_unit)
call read_phys_para_value(rho_c_r, n_unit)
call read_phys_para_value(kappa_r, n_unit)
call read_phys_para_value(rho_a, n_unit)
call read_phys_para_value(r_t, n_unit)

call read_phys_para_value(r, n_unit)
call read_phys_para_value(a, n_unit)
call read_phys_para_value(b, n_unit)
call read_phys_para_value(phi0, n_unit)
call read_phys_para_value(lambda0, n_unit)

do n=10, -190, -1
   call read_phys_para_value(rf(n), n_unit)
end do

do n=10, -190, -1
   call read_phys_para_value(kappa(n), n_unit)
end do

do n=10, -190, -1
   call read_phys_para_value(c(n), n_unit)
end do

close(n_unit, status='keep')

end subroutine phys_para
!