geometry_cartesian.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: geometry_cartesian.f90                                            #
!#                                                                           #
!# Copyright (C) 2007-2010                                                   #
!# Tobias Illenseer <tillense@astrophysik.uni-kiel.de>                       #
!# Jannes Klee <jklee@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 geometry_cartesian_mod
  USE geometry_base_mod
  USE common_dict
  IMPLICIT NONE
  !--------------------------------------------------------------------------!
  TYPE, EXTENDS (geometry_base) ::  geometry_cartesian
  CONTAINS
    PROCEDURE :: initgeometry_cartesian
    PROCEDURE :: scalefactors_1
    PROCEDURE :: scalefactors_2
    PROCEDURE :: scalefactors_3
    PROCEDURE :: scalefactors_4
    PROCEDURE :: radius_1
    PROCEDURE :: radius_2
    PROCEDURE :: radius_3
    PROCEDURE :: radius_4
    PROCEDURE :: positionvector_1
    PROCEDURE :: positionvector_2
    PROCEDURE :: positionvector_3
    PROCEDURE :: positionvector_4
    PROCEDURE :: convert2cartesian_coords_1
    PROCEDURE :: convert2cartesian_coords_2
    PROCEDURE :: convert2cartesian_coords_3
    PROCEDURE :: convert2cartesian_coords_4
    PROCEDURE :: convert2curvilinear_coords_1
    PROCEDURE :: convert2curvilinear_coords_2
    PROCEDURE :: convert2curvilinear_coords_3
    PROCEDURE :: convert2curvilinear_coords_4
    PROCEDURE :: convert2cartesian_vectors_1
    PROCEDURE :: convert2cartesian_vectors_2
    PROCEDURE :: convert2cartesian_vectors_3
    PROCEDURE :: convert2cartesian_vectors_4
    PROCEDURE :: convert2curvilinear_vectors_1
    PROCEDURE :: convert2curvilinear_vectors_2
    PROCEDURE :: convert2curvilinear_vectors_3
    PROCEDURE :: convert2curvilinear_vectors_4
    PROCEDURE :: finalize
  END TYPE
  PRIVATE
  CHARACTER(LEN=32), PARAMETER :: geometry_name = "cartesian"
  !--------------------------------------------------------------------------!
  PUBLIC :: geometry_cartesian
  !--------------------------------------------------------------------------!

CONTAINS

  SUBROUTINE initgeometry_cartesian(this,config)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(INOUT) :: this
    TYPE(DICT_TYP),POINTER                   :: config
    !------------------------------------------------------------------------!
    REAL                                     :: dz
    !------------------------------------------------------------------------!
    CALL this%InitGeometry(cartesian,geometry_name,config)
    CALL getattr(config, "dz", dz, 1.0)
  END SUBROUTINE initgeometry_cartesian

  PURE SUBROUTINE scalefactors_1(this,coords,hx,hy,hz)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, INTENT(IN),  DIMENSION(:,:) :: coords
    REAL, INTENT(OUT), DIMENSION(:)   :: hx,hy,hz
    !------------------------------------------------------------------------!
    CALL scalefactors(coords(:,1),coords(:,2),coords(:,3),hx(:),hy(:),hz(:))
  END SUBROUTINE scalefactors_1

  PURE SUBROUTINE scalefactors_2(this,coords,hx,hy,hz)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, INTENT(IN),  DIMENSION(:,:,:) :: coords
    REAL, INTENT(OUT), DIMENSION(:,:)   :: hx,hy,hz
    !------------------------------------------------------------------------!
    CALL scalefactors(coords(:,:,1),coords(:,:,2),coords(:,:,3), &
                      hx(:,:),hy(:,:),hz(:,:))
  END SUBROUTINE scalefactors_2

  PURE SUBROUTINE scalefactors_3(this,coords,hx,hy,hz)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, INTENT(IN),  DIMENSION(:,:,:,:) :: coords
    REAL, INTENT(OUT), DIMENSION(:,:,:)   :: hx,hy,hz
    !------------------------------------------------------------------------!
    CALL scalefactors(coords(:,:,:,1),coords(:,:,:,2),coords(:,:,:,3), &
                      hx(:,:,:),hy(:,:,:),hz(:,:,:))
  END SUBROUTINE scalefactors_3

  PURE SUBROUTINE scalefactors_4(this,coords,hx,hy,hz)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, INTENT(IN),  DIMENSION(:,:,:,:,:) :: coords
    REAL, INTENT(OUT), DIMENSION(:,:,:,:)   :: hx,hy,hz
    !------------------------------------------------------------------------!
    CALL scalefactors(coords(:,:,:,:,1),coords(:,:,:,:,2),coords(:,:,:,:,3), &
                      hx(:,:,:,:),hy(:,:,:,:),hz(:,:,:,:))
  END SUBROUTINE scalefactors_4

  PURE SUBROUTINE radius_1(this,coords,r)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:), INTENT(IN) :: coords
    REAL, DIMENSION(:),  INTENT(OUT) :: r
    !------------------------------------------------------------------------!
    r = radius(coords(:,1),coords(:,2),coords(:,3))
  END SUBROUTINE radius_1

  PURE SUBROUTINE radius_2(this,coords,r)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:), INTENT(IN) :: coords
    REAL, DIMENSION(:,:),  INTENT(OUT) :: r
    !------------------------------------------------------------------------!
    r = radius(coords(:,:,1),coords(:,:,2),coords(:,:,3))
  END SUBROUTINE radius_2

  PURE SUBROUTINE radius_3(this,coords,r)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:), INTENT(IN) :: coords
    REAL, DIMENSION(:,:,:),  INTENT(OUT) :: r
    !------------------------------------------------------------------------!
    r = radius(coords(:,:,:,1),coords(:,:,:,2),coords(:,:,:,3))
  END SUBROUTINE radius_3

  PURE SUBROUTINE radius_4(this,coords,r)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN) :: coords
    REAL, DIMENSION(:,:,:,:),  INTENT(OUT) :: r
    !------------------------------------------------------------------------!
    r = radius(coords(:,:,:,:,1),coords(:,:,:,:,2),coords(:,:,:,:,3))
  END SUBROUTINE radius_4

  PURE SUBROUTINE positionvector_1(this,coords,posvec)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN)  :: this
    REAL, DIMENSION(:,:), INTENT(IN)  :: coords
    REAL, DIMENSION(:,:), INTENT(OUT) :: posvec
    !------------------------------------------------------------------------!
    posvec = coords
  END SUBROUTINE positionvector_1

  PURE SUBROUTINE positionvector_2(this,coords,posvec)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN)  :: this
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: coords
    REAL, DIMENSION(:,:,:), INTENT(OUT) :: posvec
    !------------------------------------------------------------------------!
    posvec = coords
  END SUBROUTINE positionvector_2

  PURE SUBROUTINE positionvector_3(this,coords,posvec)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN)  :: this
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: coords
    REAL, DIMENSION(:,:,:,:), INTENT(OUT) :: posvec
    !------------------------------------------------------------------------!
    posvec = coords
  END SUBROUTINE positionvector_3

  PURE SUBROUTINE positionvector_4(this,coords,posvec)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN)  :: this
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: coords
    REAL, DIMENSION(:,:,:,:,:), INTENT(OUT) :: posvec
    !------------------------------------------------------------------------!
    posvec = coords
  END SUBROUTINE positionvector_4

  PURE SUBROUTINE convert2cartesian_coords_1(this,curv,cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:), INTENT(OUT) :: cart
    !------------------------------------------------------------------------!
    cart = curv
  END SUBROUTINE convert2cartesian_coords_1

  PURE SUBROUTINE convert2cartesian_coords_2(this,curv,cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:), INTENT(OUT) :: cart
    !------------------------------------------------------------------------!
    cart = curv
  END SUBROUTINE convert2cartesian_coords_2

  PURE SUBROUTINE convert2cartesian_coords_3(this,curv,cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:,:), INTENT(OUT) :: cart
    !------------------------------------------------------------------------!
    cart = curv
  END SUBROUTINE convert2cartesian_coords_3

  PURE SUBROUTINE convert2cartesian_coords_4(this,curv,cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:,:,:), INTENT(OUT) :: cart
    !------------------------------------------------------------------------!
    cart = curv
  END SUBROUTINE convert2cartesian_coords_4

  PURE SUBROUTINE convert2curvilinear_coords_1(this,cart,curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:), INTENT(IN)  :: cart
    REAL, DIMENSION(:,:), INTENT(OUT) :: curv
    !------------------------------------------------------------------------!
    curv = cart
  END SUBROUTINE convert2curvilinear_coords_1

  PURE SUBROUTINE convert2curvilinear_coords_2(this,cart,curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: cart
    REAL, DIMENSION(:,:,:), INTENT(OUT) :: curv
    !------------------------------------------------------------------------!
    curv = cart
  END SUBROUTINE convert2curvilinear_coords_2

  PURE SUBROUTINE convert2curvilinear_coords_3(this,cart,curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: cart
    REAL, DIMENSION(:,:,:,:), INTENT(OUT) :: curv
    !------------------------------------------------------------------------!
    curv = cart
  END SUBROUTINE convert2curvilinear_coords_3

  PURE SUBROUTINE convert2curvilinear_coords_4(this,cart,curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: cart
    REAL, DIMENSION(:,:,:,:,:), INTENT(OUT) :: curv
    !------------------------------------------------------------------------!
    curv = cart
  END SUBROUTINE convert2curvilinear_coords_4

  ! vector transformations
  PURE SUBROUTINE convert2cartesian_vectors_1(this,curv,v_curv,v_cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:), INTENT(IN)  :: v_curv
    REAL, DIMENSION(:,:), INTENT(OUT) :: v_cart
    !------------------------------------------------------------------------!
    v_cart = v_curv
  END SUBROUTINE convert2cartesian_vectors_1

  PURE SUBROUTINE convert2cartesian_vectors_2(this,curv,v_curv,v_cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: v_curv
    REAL, DIMENSION(:,:,:), INTENT(OUT) :: v_cart
    !------------------------------------------------------------------------!
    v_cart = v_curv
  END SUBROUTINE convert2cartesian_vectors_2

  PURE SUBROUTINE convert2cartesian_vectors_3(this,curv,v_curv,v_cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: v_curv
    REAL, DIMENSION(:,:,:,:), INTENT(OUT) :: v_cart
    !------------------------------------------------------------------------!
    v_cart = v_curv
  END SUBROUTINE convert2cartesian_vectors_3

  PURE SUBROUTINE convert2cartesian_vectors_4(this,curv,v_curv,v_cart)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: v_curv
    REAL, DIMENSION(:,:,:,:,:), INTENT(OUT) :: v_cart
    !------------------------------------------------------------------------!
    v_cart = v_curv
  END SUBROUTINE convert2cartesian_vectors_4

  PURE SUBROUTINE convert2curvilinear_vectors_1(this,curv,v_cart,v_curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:), INTENT(IN)  :: v_cart
    REAL, DIMENSION(:,:), INTENT(OUT) :: v_curv
    !------------------------------------------------------------------------!
    v_curv = v_cart
  END SUBROUTINE convert2curvilinear_vectors_1

  PURE SUBROUTINE convert2curvilinear_vectors_2(this,curv,v_cart,v_curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:), INTENT(IN)  :: v_cart
    REAL, DIMENSION(:,:,:), INTENT(OUT) :: v_curv
    !------------------------------------------------------------------------!
    v_curv = v_cart
  END SUBROUTINE convert2curvilinear_vectors_2

  PURE SUBROUTINE convert2curvilinear_vectors_3(this,curv,v_cart,v_curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:,:), INTENT(IN)  :: v_cart
    REAL, DIMENSION(:,:,:,:), INTENT(OUT) :: v_curv
    !------------------------------------------------------------------------!
    v_curv = v_cart
  END SUBROUTINE convert2curvilinear_vectors_3

  PURE SUBROUTINE convert2curvilinear_vectors_4(this,curv,v_cart,v_curv)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(IN) :: this
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: curv
    REAL, DIMENSION(:,:,:,:,:), INTENT(IN)  :: v_cart
    REAL, DIMENSION(:,:,:,:,:), INTENT(OUT) :: v_curv
    !------------------------------------------------------------------------!
    v_curv = v_cart
  END SUBROUTINE convert2curvilinear_vectors_4

  SUBROUTINE finalize(this)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(geometry_cartesian), INTENT(INOUT) :: this
    !------------------------------------------------------------------------!
    CALL this%Finalize_base()
  END SUBROUTINE finalize

  ELEMENTAL SUBROUTINE scalefactors(x,y,z,hx,hy,hz)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    REAL, INTENT(IN)  :: x,y,z
    REAL, INTENT(OUT) :: hx,hy,hz
    !------------------------------------------------------------------------!
    hx = 1.
    hy = 1.
    hz = 1.
  END SUBROUTINE scalefactors

  ELEMENTAL FUNCTION radius(x,y,z)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    REAL, INTENT(IN)  :: x,y,z
    REAL :: radius
    !------------------------------------------------------------------------!
    radius = sqrt(x*x+y*y+z*z)
  END FUNCTION radius

END MODULE geometry_cartesian_mod