prognostics.F90

Go to the documentation of this file.

module prognostics_mod
  use grids_mod, only : primal_grid
  use datadefn_mod, only : default_precision
  use mpi, only : mpi_request_null
  implicit none

#ifndef TEST_MODE
  private
#endif

  type, public :: prognostic_field_type
    logical :: active = .false., halo_allocated=.false.
    integer, dimension(3) :: grid
    real(kind=DEFAULT_PRECISION), dimension(:,:,:), allocatable :: data
    real(kind=DEFAULT_PRECISION), dimension(:), allocatable :: flux_previous_y, flux_y_buffer
    real(kind=DEFAULT_PRECISION), dimension(:,:), allocatable :: flux_previous_x
    integer :: async_flux_handle = mpi_request_null
  end type prognostic_field_type

  type, public :: prognostic_field_ptr_type
     type(prognostic_field_type), pointer :: ptr
  end type prognostic_field_ptr_type

  public get_field_interpolation_index
contains

  function get_field_interpolation_index(field)
    type(prognostic_field_type), intent(inout) :: field
    logical, dimension(3) :: get_field_interpolation_index

    integer :: i
    do i=1,3
      get_field_interpolation_index(i) = merge(.true., .false., field%grid(i) == primal_grid)
    end do
  end function get_field_interpolation_index
end module prognostics_mod