timedisc_generic.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: mesh_generic.f90                                                  #
!#                                                                           #
!# Copyright (C) 2016 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 timedisc_generic_mod
  USE timedisc_base_mod
  USE timedisc_modeuler_mod
  USE timedisc_rkfehlberg_mod
  USE timedisc_cashkarp_mod
  USE timedisc_dormand_prince_mod
  USE timedisc_ssprk_mod
  USE mesh_base_mod
  USE physics_base_mod
  USE common_dict

CONTAINS

  SUBROUTINE new_timedisc(Timedisc,Mesh,Physics,config,IO)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(timedisc_base), ALLOCATABLE   :: Timedisc
    CLASS(mesh_base),     INTENT(INOUT) :: Mesh
    CLASS(physics_base),  INTENT(IN)    :: Physics
    TYPE(DICT_TYP),       POINTER       :: config, IO
    !------------------------------------------------------------------------!
    INTEGER                             :: method
    !------------------------------------------------------------------------!
    CALL getattr(config,"method",method)

    ! allocate data
    SELECT CASE(method)
    CASE(modified_euler)
      ALLOCATE(timedisc_modeuler::timedisc)
    CASE(rk_fehlberg)
      ALLOCATE(timedisc_rkfehlberg::timedisc)
    CASE(cash_karp)
      ALLOCATE(timedisc_cashkarp::timedisc)
    CASE(dormand_prince)
      ALLOCATE(timedisc_dormand_prince::timedisc)
    CASE(ssprk)
      ALLOCATE(timedisc_ssprk::timedisc)
    CASE DEFAULT
      ALLOCATE(timedisc_modeuler::timedisc)
      CALL timedisc%Error("new_timedisc","Unknown timedisc integration scheme")
    END SELECT

    ! call initialization
    SELECT TYPE(obj => timedisc)
    TYPE IS (timedisc_modeuler)
      CALL obj%InitTimedisc_modeuler(mesh,physics,config,io)
    TYPE IS (timedisc_rkfehlberg)
      CALL obj%InitTimedisc_rkfehlberg(mesh,physics,config,io)
    TYPE IS (timedisc_cashkarp)
      CALL obj%InitTimedisc_cashkarp(mesh,physics,config,io)
    TYPE IS (timedisc_dormand_prince)
      CALL obj%InitTimedisc_dormand_prince(mesh,physics,config,io)
    TYPE IS (timedisc_ssprk)
      CALL obj%InitTimedisc_ssprk(mesh,physics,config,io)
    END SELECT
  END SUBROUTINE
END MODULE timedisc_generic_mod