ver_interp_lin3d_surf.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.
!     ##############################################
      FUNCTION ver_interp_lin3d_surf(PVAR1,KKLIN,PCOEFLIN) RESULT(PVAR2)
!     ##############################################
!
!!****  *VER_INTERP_LIN* - vertical linear interpolation
!!
!!    PURPOSE
!!    -------
!     This function interpolates the 3D fields from one grid
!     to another using linear interpolation cofficients stored in module
!     MODD_VER_INTERP_LIN.
!
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    Book 2
!!
!!    AUTHOR
!!    ------
!!
!     V.Masson  Meteo-France
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    18/07/97
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declaration of arguments
!              ------------------------
REAL,   DIMENSION(:,:,:), INTENT(IN) :: pvar1 ! variable values on the initial
!                                             ! grid
INTEGER,DIMENSION(:,:,:), INTENT(IN) :: kklin ! lower interpolating level of
!                                             ! grid 1 for each level of grid 2 
REAL,   DIMENSION(:,:,:), INTENT(IN) :: pcoeflin ! coefficient for level KKLIN
!
REAL,   DIMENSION(SIZE(KKLIN,1),SIZE(KKLIN,2),SIZE(KKLIN,3))                   &
                                       :: pvar2 ! variable values on target  
!                                             ! grid 
!
!*       0.2   Declaration of local variables
!              ------------------------------
!
INTEGER                                               :: ji,jj,jk2
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN3D_SURF:VER_INTERP_LIN3D_SURF',0,zhook_handle)
DO jk2=1,SIZE(kklin,3)
  DO jj=1,SIZE(kklin,2)
    DO ji=1,SIZE(kklin,1)
      pvar2(ji,jj,jk2)=    pcoeflin(ji,jj,jk2) *pvar1(ji,jj,kklin(ji,jj,jk2)  )&
                        +(1.-pcoeflin(ji,jj,jk2))*pvar1(ji,jj,kklin(ji,jj,jk2)+1)  
    END DO
  END DO
END DO
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN3D_SURF:VER_INTERP_LIN3D_SURF',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END FUNCTION ver_interp_lin3d_surf