hor_interpol_conf_proj.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
SUBROUTINE hor_interpol_conf_proj(GCP,KLUOUT,PFIELDIN,PFIELDOUT)
!     #################################################################################
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    MODIFICATION
!!    ------------
!!
!!    02/04/12 M. Tomasini  Add an index in the second dimension of the ISBA 
!!                          array variables for BILIN interpolation routine to 
!!                          not bug in case 2D (this is not the more beautiful
!!                          method; the BILIN routine should better be adapted)
!!                          Search  ! Ajout MT
!!    10/02/15 M.Moge  using SIZE(PFIELDOUT,1) instead of SIZE(XLAT_OUT)
!-------------------------------------------------------------------------------
!
USE modd_grid_conf_proj_n, ONLY : grid_conf_proj_t, xy, xx, xcx, xcy, ncij
!
USE modd_surfex_mpi, ONLY : nrank,npio
USE modd_prep,           ONLY : xx_out, xy_out, linterp
USE modd_surf_par,   ONLY : xundef
!
USE modi_bilin_value
USE modi_bilin_extrap
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
TYPE(grid_conf_proj_t),INTENT(INOUT) :: GCP
!
INTEGER,            INTENT(IN)  :: KLUOUT    ! logical unit of output listing
REAL, DIMENSION(:,:), INTENT(IN)  :: PFIELDIN  ! field to interpolate horizontally
REAL, DIMENSION(:,:), INTENT(OUT) :: PFIELDOUT ! interpolated field
!
!*      0.2    declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
LOGICAL, DIMENSION(:), ALLOCATABLE :: GINTERP_DUPLIQUE ! .true. where physical value is needed ! Ajout MT
!-------------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('HOR_INTERPOL_CONF_PROJ',0,zhook_handle)
!
!*      4.    Interpolation with bilinear
!
 CALL bilin_value(kluout,gcp%NX,gcp%NY,pfieldin,xcx,xcy,ncij(:,1),ncij(:,2),pfieldout)
 CALL bilin_extrap(kluout,gcp%NX,gcp%NY,ncij,xx,xy,pfieldin,xx_out,xy_out,pfieldout,linterp)
!
IF (lhook) CALL dr_hook('HOR_INTERPOL_CONF_PROJ',1,zhook_handle)
!
!-------------------------------------------------------------------------------------
END SUBROUTINE hor_interpol_conf_proj