sources_c_accel.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: sources_c_accel.f90                                               #
!#                                                                           #
!# Copyright (C) 2009-2019                                                   #
!# Björn Sperling   <sperling@astrophysik.uni-kiel.de>                       #
!# Tobias Illenseer <tillense@astrophysik.uni-kiel.de>                       #
!# Jannes Klee      <tillense@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 sources_c_accel_mod
  USE sources_base_mod
  USE fluxes_base_mod
  USE physics_base_mod
  USE mesh_base_mod
  USE marray_base_mod
  USE marray_compound_mod
  USE common_dict
  IMPLICIT NONE
  !--------------------------------------------------------------------------!
  PRIVATE
    CHARACTER(LEN=32) :: source_name = "constant acceleration"

  TYPE, EXTENDS(sources_base) :: sources_c_accel
    TYPE(marray_base)                 :: accel
  CONTAINS
    PROCEDURE :: initsources_c_accel
    PROCEDURE :: externalsources_single
    PROCEDURE :: calctimestep_single
    PROCEDURE :: infosources
    PROCEDURE :: finalize
  END TYPE

  PUBLIC :: &
       ! classes
       sources_c_accel

CONTAINS

  SUBROUTINE initsources_c_accel(this,Mesh,Physics,Fluxes,config,IO)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(sources_c_accel), INTENT(INOUT) :: this
    CLASS(mesh_base),       INTENT(IN)    :: Mesh
    CLASS(fluxes_base),     INTENT(IN)    :: Fluxes
    CLASS(physics_base),    INTENT(IN)    :: Physics
    TYPE(Dict_TYP),           POINTER     :: config, IO
    !------------------------------------------------------------------------!
    INTEGER :: k,stype
    REAL    :: accel(3)
    !------------------------------------------------------------------------!
    CALL getattr(config,"stype",stype)
    CALL this%InitLogging(stype,source_name)
    CALL this%InitSources(mesh,fluxes,physics,config,io)

    this%accel = marray_base(physics%VDIM)
    this%accel%data1d(:) = 0.0

    ! initialize constant acceleration
    CALL getattr(config, "xaccel", accel(1), 0.0)
    CALL getattr(config, "yaccel", accel(2), 0.0)
    CALL getattr(config, "zaccel", accel(3), 0.0)
    DO k=1,physics%VDIM
      this%accel%data2d(:,k) = accel(k)
    END DO

  END SUBROUTINE initsources_c_accel

  SUBROUTINE externalsources_single(this,Mesh,Physics,Fluxes,Sources,time,dt,pvar,cvar,sterm)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(sources_c_accel),INTENT(INOUT):: this
    CLASS(mesh_base),INTENT(IN)         :: Mesh
    CLASS(physics_base),INTENT(INOUT)   :: Physics
    CLASS(fluxes_base),INTENT(IN)       :: Fluxes
    CLASS(sources_base), INTENT(INOUT)  :: Sources
    REAL,INTENT(IN)                     :: time, dt
    CLASS(marray_compound),INTENT(INOUT):: pvar,cvar,sterm
    !------------------------------------------------------------------------!
    ! compute source terms due to constant acceleration
    CALL physics%ExternalSources(this%accel,pvar,cvar,sterm)
  END SUBROUTINE

  SUBROUTINE calctimestep_single(this,Mesh,Physics,Fluxes,pvar,cvar,time,dt)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(sources_c_accel), INTENT(INOUT) :: this
    CLASS(mesh_base),       INTENT(IN)    :: Mesh
    CLASS(physics_base),    INTENT(INOUT) :: Physics
    CLASS(fluxes_base),     INTENT(IN)    :: Fluxes
    CLASS(marray_compound), INTENT(INOUT) :: pvar,cvar
    REAL,                   INTENT(IN)    :: time
    REAL,                   INTENT(OUT)   :: dt
    !------------------------------------------------------------------------!
    dt = huge(dt)
  END SUBROUTINE calctimestep_single

  SUBROUTINE infosources(this,Mesh)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(sources_c_accel), INTENT(IN) :: this
    CLASS(mesh_base),       INTENT(IN) :: Mesh
    !------------------------------------------------------------------------!
  END SUBROUTINE infosources

  SUBROUTINE finalize(this)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(sources_c_accel), INTENT(INOUT) :: this
    !------------------------------------------------------------------------!
    CALL this%accel%Destroy()
  END SUBROUTINE finalize

END MODULE sources_c_accel_mod