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