constants_geometrical.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: constants_geometrical.f90                                         #
!#                                                                           #
!# Copyright (C) 2007-2016                                                   #
!# Tobias Illenseer <tillense@astrophysik.uni-kiel.de>                       #
!# Manuel Jung <mjung@astrophysik.uni-kiel.de>                               #
!#                                                                           #
!# This program 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 2 of the License, or (at     #
!# your option) any later version.                                           #
!#                                                                           #
!# This program 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, GOOD TITLE or        #
!# NON INFRINGEMENT.  See the GNU General Public License for more            #
!# details.                                                                  #
!#                                                                           #
!# You should have received a copy of the GNU General Public License         #
!# along with this program; if not, write to the Free Software               #
!# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.                 #
!#                                                                           #
!#############################################################################
 
!----------------------------------------------------------------------------!
!----------------------------------------------------------------------------!
MODULE constants_geometrical_mod
  USE constants_base_mod
  USE constants_si_mod
  IMPLICIT NONE
  !--------------------------------------------------------------------------!
  PRIVATE
  TYPE, EXTENDS(constants_si) :: constants_geometrical
  CONTAINS
    PROCEDURE :: initconstants_geometrical
  END TYPE
  CHARACTER(LEN=32), PARAMETER :: units_name = 'geometrical'
  !--------------------------------------------------------------------------!
  PUBLIC :: &
       ! types
       constants_geometrical
  !--------------------------------------------------------------------------!
 
CONTAINS
 
  SUBROUTINE initconstants_geometrical(this)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(constants_geometrical), INTENT(INOUT) :: this
    !------------------------------------------------------------------------!
    REAL                                        :: C, GN, KB, SB, KE
    !------------------------------------------------------------------------!
    ! numerical values of physical constants in geometrical units (c = G = 1)
    this%C  = 1.0                                   ! lightspeed             !
    this%GN = 1.0                                   ! Newtons grav. constant !
    this%KB = 1.0                                   ! Boltzmann constant     !
    this%NA = na  ![1/mol]                          ! Avogadro constant      !
    ! factors for conversion to SI units
    ! time and mass have unit of length i.e. metre
 
    c  = this%constants_SI%C
    gn = this%constants_SI%GN
    kb = this%constants_SI%KB
    sb = this%constants_SI%SB
    ke = this%constants_SI%KE
    this%cf_time = c                 ! i.e. 1 sec [SI] = 3e8 m [geometrical] !
    this%cf_mass = (gn/c)/c
    this%cf_momentum = this%cf_mass/c
    this%cf_energy = this%cf_momentum/c
    this%cf_power = this%cf_energy/c
    this%cf_temperature = kb * this%cf_energy
    this%cf_density = this%cf_mass
    ! some scaled constants
    this%KE = ke*c/(gn/c)               !  [m]        electron scat. opacity !
 
    CALL this%InitConstants(geometrical,units_name)
  END SUBROUTINE initconstants_geometrical
 
END MODULE constants_geometrical_mod