read_phys_para_value.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  r e a d _ p h y s _ p a r a _ v a l u e
!
!> @file
!!
!! Reading of a value of a physical parameter from the phys_para file.
!!
!! @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 a value of a physical parameter from the phys_para file.
!<------------------------------------------------------------------------------
subroutine read_phys_para_value(para, n_unit)

use sico_types

implicit none

integer(i4b), intent(in)  :: n_unit
real(dp),     intent(out) :: para

character                   :: check, ch_dummy
character(len=*), parameter :: fmt1='(a)', fmt3='(15x)'
logical                     :: first_read

first_read = .true.

do
   if (first_read) then
      read(unit=n_unit, fmt=trim(fmt1), advance='no') check
      first_read=.false.
   else
      read(unit=n_unit, fmt=trim(fmt1)) ch_dummy
      read(unit=n_unit, fmt=trim(fmt1), advance='no') check
   end if

   if (check /= '%') then   ! no comment line
      read(unit=n_unit, fmt=trim(fmt3), advance='no')
      read(unit=n_unit, fmt=*) para
      exit
   end if
end do

end subroutine read_phys_para_value
!