common.f90 Source File

module m_common
  implicit none

  integer, parameter :: dp = kind(0.0d0)
  real(dp), parameter :: pi = 4*atan(1.0_dp)

  integer, parameter :: RDR_X2Y = 12, RDR_X2Z = 13, RDR_Y2X = 21, &
                        RDR_Y2Z = 23, RDR_Z2X = 31, RDR_Z2Y = 32, &
                        RDR_C2X = 41, RDR_C2Y = 42, RDR_C2Z = 43, &
                        RDR_X2C = 14, RDR_Y2C = 24, RDR_Z2C = 34
  integer, parameter :: DIR_X = 1, DIR_Y = 2, DIR_Z = 3, DIR_C = 4
  integer, parameter :: POISSON_SOLVER_FFT = 0, POISSON_SOLVER_CG = 1
  integer, parameter :: VERT = 0000, & ! Vertex centered data
                        CELL = 1110, & ! Cell centered data
                        X_FACE = 1100, & ! Data on faces normal to X
                        Y_FACE = 1010, & ! Data on faces normal to Y
                        Z_FACE = 0110, & ! Data on faces normal to Z
                        X_EDGE = 0010, & ! Data on edges along X
                        Y_EDGE = 0100, & ! Data on edges along Y
                        Z_EDGE = 1000, & ! Data on edges along Z
                        NULL_LOC = -0001 ! The location of data isn't specified
  integer, parameter :: BC_PERIODIC = 0, BC_NEUMANN = 1, BC_DIRICHLET = 2, &
                        BC_HALO = -1
  integer, protected :: &
    rdr_map(4, 4) = reshape([0, RDR_Y2X, RDR_Z2X, RDR_C2X, &
                             RDR_X2Y, 0, RDR_Z2Y, RDR_C2Y, &
                             RDR_X2Z, RDR_Y2Z, 0, RDR_C2Z, &
                             RDR_X2C, RDR_Y2C, RDR_Z2C, 0], shape=[4, 4])

contains

  pure subroutine get_dirs_from_rdr(dir_from, dir_to, rdr_dir)
    integer, intent(out) :: dir_from, dir_to
    integer, intent(in) :: rdr_dir
    integer, dimension(2) :: dirs

    dirs = findloc(rdr_map, rdr_dir)
    dir_from = dirs(1)
    dir_to = dirs(2)

  end subroutine

  pure integer function get_rdr_from_dirs(dir_from, dir_to) result(rdr_dir)
      !! Returns RDR_?2? value based on two direction inputs
    integer, intent(in) :: dir_from, dir_to

    rdr_dir = rdr_map(dir_from, dir_to)
  end function get_rdr_from_dirs

  function get_argument(pos) result(arg)
    integer, intent(in) :: pos
    character(:), allocatable :: arg

    character(len=200) :: temp
    integer :: stat

    call get_command_argument(pos, temp, status=stat)

    if (stat > 0) then
      error stop 'Argument retrieval failed!'
    else if (stat == -1) then
      error stop 'Argument is truncated!'
    end if

    arg = trim(temp)
  end function get_argument

  integer function move_data_loc(in_data_loc, dir, move) result(out_data_loc)
    integer, intent(in) :: in_data_loc, dir, move

    out_data_loc = in_data_loc + move*(10**dir)
  end function move_data_loc

end module m_common