!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!MNH_LIC for details. version 1.
!-----------------------------------------------------------------
!--------------- special set of characters for RCS information
!-----------------------------------------------------------------
! $Source: /srv/cvsroot/MNH-VX-Y-Z/src/MNH/ini_cturb.f90,v $ $Revision: 1.2.4.1.2.2.2.1.2.2 $ $Date: 2014/01/09 15:01:55 $
!-----------------------------------------------------------------
!-----------------------------------------------------------------
!       #########################
        MODULE MODI_IBM_VALUECORN_2D
!       #########################
!
INTERFACE
!
FUNCTION IBM_VALUECORN_2D(PVAR,IINDEX) RESULT(PVALUE)

  REAL, DIMENSION(:,:)         ,INTENT(IN)    :: PVAR
  INTEGER,DIMENSION(:)         ,INTENT(IN)    :: IINDEX
  REAL, DIMENSION(8)                          :: PVALUE

END FUNCTION IBM_VALUECORN_2D
!
END INTERFACE
!
END MODULE MODI_IBM_VALUECORN_2D
!
!
!
!       ##################################################
        FUNCTION IBM_VALUECORN_2D(PVAR,IINDEX) RESULT(PVALUE)
!       ##################################################
!
!!****     *IBM_VALUECORN_2D*  - routine to affect values at cornes cell
!!
!!      PURPOSE
!!      -------
!         The purpose of this routine is to compute (VAR) at corners of cell (U,V,W,P)
!
!!      METHOD
!!      ------
!     Index initial value 
!     1 <-> i  ,j  ,k
!     2 <-> i+1,j  ,k
!     3 <-> i  ,j+1,k
!     4 <-> i+1,j+1,k
!     5 <-> i  ,j  ,k+1
!     6 <-> i+1,j  ,k+1
!     7 <-> i  ,j+1,k+1
!     8 <-> i+1,j+1,k+1
!!
!!      EXTERNAL
!!      --------
!!        NONE
!!
!!      IMPLICIT ARGUMENTS
!!      ------------------
!!
!!      REFERENCE
!!      ---------
!!
!!      AUTHOR
!!      ------
!!        Franck Auguste       * CERFACS(AE) *
!!
!!      MODIFICATIONS
!!      -------------
!!        Original          01/01/2015
!!
!------------------------------------------------------------------------------
!       
!**** 0. DECLARATIONS
!     ---------------
!
! module
!
! declaration
USE MODD_IBM_PARAM_n
!
! interface
!
IMPLICIT NONE
!
!------------------------------------------------------------------------------
!
!       0.1  declarations of arguments
REAL, DIMENSION(:,:)         ,INTENT(IN)    :: PVAR                                                ! variable array
INTEGER, DIMENSION(:)        ,INTENT(IN)    :: IINDEX                                              ! IJK reference
REAL, DIMENSION(8)                          :: PVALUE
!
!------------------------------------------------------------------------------
!
!       0.2  declaration of local variables
INTEGER                                     :: JI,JJ,JK,JL                                         ! loop index
!
!------------------------------------------------------------------------------
!
!       0.3  Allocation
!
!
!------------------------------------------------------------------------------
!       
!     Index initial value 
!     1 <-> i  ,j  ,k
!     2 <-> i+1,j  ,k
!     3 <-> i  ,j+1,k
!     4 <-> i+1,j+1,k
!     5 <-> i  ,j  ,k+1
!     6 <-> i+1,j  ,k+1
!     7 <-> i  ,j+1,k+1
!     8 <-> i+1,j+1,k+1

  JK = 2
  DO JL = 1,8
    !
    ! corners index 
    IF (JL==1) THEN
      JI = IINDEX(1)
      JJ = IINDEX(2)
    ENDIF
    IF (JL==2) THEN
      JI = IINDEX(1)+1
      JJ = IINDEX(2)
    ENDIF
    IF (JL==3) THEN
      JI = IINDEX(1)
      JJ = IINDEX(2)+1
    ENDIF
    IF (JL==4) THEN
      JI = IINDEX(1)+1
      JJ = IINDEX(2)+1
    ENDIF
    IF (JL==5) THEN
      JI = IINDEX(1)
      JJ = IINDEX(2)
    ENDIF
    IF (JL==6) THEN
      JI = IINDEX(1)+1
      JJ = IINDEX(2)
    ENDIF
    IF (JL==7) THEN
      JI = IINDEX(1)
      JJ = IINDEX(2)+1
    ENDIF
    IF (JL==8) THEN
      JI = IINDEX(1)+1
      JJ = IINDEX(2)+1
    ENDIF    

    PVALUE(JL) = PVAR(JI,JJ)

  ENDDO
!
!**** X. DEALLOCATIONS/CLOSES
!     -----------------------
!
RETURN
!
END FUNCTION IBM_VALUECORN_2D