psmile_trs_give_neighbors_gauss.F90

Go to the documentation of this file.

!------------------------------------------------------------------------
! Copyright 2006-2010, CERFACS, Toulouse, France.
! Copyright 2006-2010, NEC Europe Ltd., London, UK.
! All rights reserved. Use is subject to OASIS4 license terms.
!-----------------------------------------------------------------------
!BOP
!
! !ROUTINE: PSMILe_Trs_give_neighbors_gauss
!
! !INTERFACE:

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

  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

! indicator whether the transformer has to operate on pole cells
  INTEGER, DIMENSION(id_epio_tgt_size), INTENT (In) :: 
                                              ida_neigh_bascule
!
! !OUTPUT PARAMETERS:
!
  INTEGER, INTENT (Out)                    :: id_err   ! error value

! !DESCRIPTION:
!
! Subroutine "PSMILe_Trs_give_neighbors_gauss" insures the exchange of information,
! when an interpolation is required, for Gauss-reduced source grids
!
! !REVISION HISTORY:
!   Date      Programmer   Description
! ----------  ----------   -----------
! 19/02/2009  R. Redler    Creation
!
!EOP
!----------------------------------------------------------------------
! $Id: psmile_trs_give_neighbors_gauss.F90,v 1.19.2.4 2008/10/31 12:32:35 redler Exp $
! $Author: redler $
!----------------------------------------------------------------------
!
! 0. Local declarations
!
  CHARACTER(LEN=len_cvs_string), SAVE  :: mycvs = 
     '$Id: psmile_trs_give_neighbors_gauss.F90,v 1.19.2.4 2008/10/31 12:32:35 redler Exp $'

  INTEGER, DIMENSION(PSMILe_trans_Header_length) :: ila_args 
!
!----------------------------------------------------------------------
!
#ifdef VERBOSE
      print *, trim(ch_id), ': PSMILe_Trs_give_neighbors_gauss: 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
  ila_args(11) = PRISM_Gaussreduced_regvrt
!
! 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)
!
! 4. Send the information about special of the target points
!    Point-based interpolations from a Gauss-reduced source grid
!    require additional information:
!
  CALL MPI_Send(ida_neigh_bascule, id_epio_tgt_size, &
     MPI_Integer, id_trans_rank, 5, comm_trans, id_err)
!
#ifdef VERBOSE
      print *, trim(ch_id), &
         ': PSMILe_Trs_give_neighbors_gauss: eof ierror =', id_err
      call psmile_flushstd
#endif /* VERBOSE */
  
END SUBROUTINE PSMILe_Trs_give_neighbors_gauss