hor_interpol_cartesian.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_cartesian(KLUOUT,PFIELDIN,PFIELDOUT)
!     #################################################################################
!
!
USE modd_surfex_mpi, ONLY : nrank, npio
USE modd_prep,           ONLY : xlat_out, xlon_out, xx_out, xy_out, linterp
USE modd_grid_cartesian, ONLY : xx, xy, nx, ny, xcx, xcy, ncij
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
!
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

!-------------------------------------------------------------------------------------
!
!*      1.    Allocations
!
IF (lhook) CALL dr_hook('HOR_INTERPOL_CARTESIAN',0,zhook_handle)
!
!*      4.    Interpolation with bilinear
!
 CALL bilin_value(kluout,nx,ny,pfieldin,xcx,xcy,ncij(:,1),ncij(:,2),pfieldout)
 CALL bilin_extrap(kluout,nx,ny,ncij,xx,xy,pfieldin,xx_out,xy_out,pfieldout,linterp)
!
IF (lhook) CALL dr_hook('HOR_INTERPOL_CARTESIAN',1,zhook_handle)
!
!-------------------------------------------------------------------------------------
END SUBROUTINE hor_interpol_cartesian