constants_SI.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: constants_SI.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_si_mod
  USE constants_base_mod
  IMPLICIT NONE
  !--------------------------------------------------------------------------!
  PRIVATE
  TYPE, EXTENDS(constants_base) :: constants_si
  CONTAINS
    PROCEDURE :: initconstants_si
  END TYPE
  CHARACTER(LEN=32), PARAMETER :: units_name = 'SI'
  !--------------------------------------------------------------------------!
  ! basic numerical constants
  !!#### Physical constants in SI units
  REAL, PARAMETER :: &
     c  = 2.99792458e+08, &   
     gn = 6.6742e-11, &       
     kb = 1.3806505e-23, &    
     na = 6.022e+23, &        
     sb = 5.6704e-8, &        
     ke = 3.48e-02, &         
     au = 1.49597870691d+11, &
     msun     = 1.9885d+30, & 
     mjupiter = 1.89819d+27, &
     mearth   = 5.9723d+24, & 
     rsun     = 6.957e+8, &   
     rjupiter = 6.9911e+7, &  
     rearth   = 6.371e+6, &   
     day      = 8.64e+4, &    
     pi = 3.1415926535897932384626433832795028842
  !--------------------------------------------------------------------------!
  PUBLIC :: &
       ! types
       constants_si, &
       c, gn, kb, na, sb, ke, au, &
       msun, mjupiter, mearth, rsun, rjupiter, rearth, day, pi
  !--------------------------------------------------------------------------!
 
CONTAINS
 
  SUBROUTINE initconstants_si(this)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(constants_si), INTENT(INOUT) :: this
    !------------------------------------------------------------------------!
    ! assign numerical values of physical constants in SI units;
    ! (C, GN, etc. are defined in constants_common)
    this%C  = c
    this%GN = gn
    this%KB = kb
    this%NA = na
    this%SB = sb
    this%KE = ke
    this%AU = au
    this%MSUN     = msun
    this%MJUPITER = mjupiter
    this%MEARTH   = mearth
    this%RSUN     = rsun
    this%RJUPITER = rjupiter
    this%REARTH   = rearth
    this%DAY      = day
    this%PI       = pi
    ! conversion factors to SI units are unity
    this%cf_time        = 1.0
    this%cf_mass        = 1.0
    this%cf_momentum    = 1.0
    this%cf_energy      = 1.0
    this%cf_power       = 1.0
    this%cf_temperature = 1.0
    this%cf_density     = 1.0
    this%cf_opacity     = 1.0
 
    CALL this%InitConstants(si,units_name)
  END SUBROUTINE initconstants_si
 
END MODULE constants_si_mod