psmile_trs_give_neighbors3d.F90

Go to the documentation of this file.

!------------------------------------------------------------------------
! Copyright 2006-2010, CERFACS, Toulouse, France.
! All rights reserved. Use is subject to OASIS4 license terms.
!-----------------------------------------------------------------------
!BOP
!
! !ROUTINE: PSMILe_Trs_give_neighbors3d
!
! !INTERFACE:

subroutine psmile_trs_give_neighbors3d(id_epio_id,           &
                                       id_trans_rank,        &
                                       id_epio_tgt_size,     &
                                       id_nb_neighbors,      &
                                       ida_neighbor_indices,        &
                                       id_err)
!
! !USES:
!
  use PRISM_constants
  USE PSMILe, dummy_interface => PSMILe_Trs_give_neighbors3d

  IMPLICIT NONE
!
! !INPUT PARAMETERS:
!
! associated epio id
  INTEGER, INTENT (In)                     :: id_epio_id
  INTEGER, INTENT (In)                     :: id_trans_rank   

! size of the corresponding epiot
  INTEGER, INTENT (In)                     :: id_epio_tgt_size

! number of neighbors
  INTEGER, INTENT (In)                     :: id_nb_neighbors

! index of the src grid neighbors for all target points, computed by the PSMILe
  INTEGER, DIMENSION(id_nb_neighbors*id_epio_tgt_size), INTENT (In) :: 
     ida_neighbor_indices
!
! !OUTPUT PARAMETERS:
!
  INTEGER, INTENT (Out)                    :: id_err   ! error value

! !DESCRIPTION:
!
! Subroutine "PSMILe_Trs_give_neighbors3d" insures the exchange of information,
! when an interpolation is required.
!
! !REVISION HISTORY:
!   Date      Programmer   Description
! ----------  ----------   -----------
! 24/03/2003  D. Declat    Creation
!
!EOP
!----------------------------------------------------------------------
! $Id: psmile_trs_give_neighbors3d.F90 2325 2010-04-21 15:00:07Z valcke $
! $Author: valcke $
!----------------------------------------------------------------------
!
! 0. Local declarations
!
  CHARACTER(LEN=len_cvs_string), SAVE  :: mycvs = 
     '$Id: psmile_trs_give_neighbors3d.F90 2325 2010-04-21 15:00:07Z valcke $'

  INTEGER, DIMENSION(PSMILe_trans_Header_length) :: ila_args 
!
!----------------------------------------------------------------------
!
#ifdef VERBOSE
      print *, trim(ch_id), ': PSMILe_Trs_give_neighbors3d: start'
      call psmile_flushstd
#endif /* VERBOSE */

#ifdef DEBUG
   print *, trim(ch_id), ': - id_epio_id       ', id_epio_id
   print *, trim(ch_id), ': - global_rank      ', global_rank
   print *, trim(ch_id), ': - id_epio_tgt_size ', id_epio_tgt_size
   print *, trim(ch_id), ': - id_nb_neighbors  ', id_nb_neighbors
#endif

! 1. Set the header message
!
! 1.1. Initialize the contents to 999999

  ila_args = 999999

! 1.2. Content : action, global rank, number of the transformation, field_id,
!                comp_id, size (of the EPIOT)

  ila_args(1) = PSMILe_trans_Set_neighbors_info
  ila_args(2) = global_rank
  ila_args(3) = id_epio_id
  ila_args(4) = id_epio_tgt_size
  ila_args(5) = id_nb_neighbors
  ila_args(6) = PSMILe_3D
  ila_args(7) = PSMILe_Undef

!
! 2. Send the header message to the transformer
!
  call psmile_trs_inform(ila_args, id_trans_rank, id_err)

!
! 3. Send the neighbors of the target points
!
  CALL MPI_Send(ida_neighbor_indices(1), id_nb_neighbors*id_epio_tgt_size, &
     MPI_Integer, id_trans_rank, 4, comm_trans, id_err)

!
#ifdef VERBOSE
      print *, trim(ch_id), &
         ': PSMILe_Trs_give_neighbors3d: eof ierror =', id_err
      call psmile_flushstd
#endif /* VERBOSE */
  
END SUBROUTINE PSMILe_Trs_give_neighbors3d