marray_cellvector.f90

Go to the documentation of this file.

!#############################################################################
!#                                                                           #
!# fosite - 3D hydrodynamical simulation program                             #
!# module: marray_cellvector.f90                                             #
!#                                                                           #
!# Copyright (C) 2018                                                        #
!# Tobias Illenseer <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 marray_cellvector_mod
  USE marray_base_mod
  IMPLICIT NONE
  !--------------------------------------------------------------------------!
  PRIVATE
  TYPE, EXTENDS(marray_base) :: marray_cellvector
    REAL, DIMENSION(:,:,:,:), POINTER   :: center, &     !< geometric center
                                         bcenter

    REAL, DIMENSION(:,:,:,:,:), POINTER :: faces, &      !< cell face centers
                                         corners
    CONTAINS
    PROCEDURE :: assignpointers
  END TYPE
  INTERFACE marray_cellvector
    MODULE PROCEDURE createmarray_cellvector
  END INTERFACE
  !--------------------------------------------------------------------------!
  PUBLIC :: marray_cellvector

CONTAINS

  FUNCTION createmarray_cellvector() RESULT(new_cv)
    IMPLICIT NONE
    !-------------------------------------------------------------------!
    TYPE(marray_cellvector) :: new_cv
    !-------------------------------------------------------------------!
    ! only set rank & dims - allocation is done, when cellvector is assigned (with =)
    new_cv%DIMS(1) = 1+1+6+8
    new_cv%DIMS(2) = 3
    new_cv%RANK = 2
  END FUNCTION createmarray_cellvector
  
  SUBROUTINE assignpointers(this)
    IMPLICIT NONE
    !------------------------------------------------------------------------!
    CLASS(marray_cellvector),INTENT(INOUT) :: this
    !------------------------------------------------------------------------!
    CALL this%marray_base%AssignPointers()
    ! assign array pointers
    this%center  => this%RemapBounds(this%data5d(:,:,:,1,:))
    this%bcenter => this%RemapBounds(this%data5d(:,:,:,2,:))
    this%faces   => this%RemapBounds(this%data5d(:,:,:,3:8,:))
    this%corners => this%RemapBounds(this%data5d(:,:,:,9:16,:))    
  END SUBROUTINE assignpointers
  
END MODULE marray_cellvector_mod