ver_interp_lin_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.
!     ######################
      MODULE modi_ver_interp_lin_surf
!     ######################
INTERFACE ver_interp_lin_surf
!     ##############################################
      FUNCTION ver_interp_lin3d_surf_bis(PVAR1,KKLIN,PCOEFLIN) RESULT(PVAR2)
!     ##############################################
!
! third dimension of the arrays is vertical
!
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 
END FUNCTION ver_interp_lin3d_surf_bis
!     ##############################################
      FUNCTION ver_interp_lin2d_surf(PVAR1,KKLIN,PCOEFLIN) RESULT(PVAR2)
!     ##############################################
!
! second dimension of the arrays is vertical
!
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))                               &
                                       :: pvar2 ! variable values on target  
!                                             ! grid 
END FUNCTION ver_interp_lin2d_surf
!     ##############################################
      FUNCTION ver_interp_lin1d_surf(PVAR1,KKLIN,PCOEFLIN) RESULT(PVAR2)
!     ##############################################
!
! first dimension of the arrays is vertical
!
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)) :: pvar2 ! variable values on target
!                                         ! grid 
END FUNCTION ver_interp_lin1d_surf
!
!
END INTERFACE
END MODULE modi_ver_interp_lin_surf
!     ##############################################
      FUNCTION ver_interp_lin3d_surf_bis(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
!
USE modi_ver_interp_lin3d_surf
!
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_LIN_SURF:VER_INTERP_LIN3D_SURF_BIS',0,zhook_handle)
pvar2(:,:,:)=ver_interp_lin3d_surf(pvar1(:,:,:),kklin(:,:,:),pcoeflin(:,:,:))
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN_SURF:VER_INTERP_LIN3D_SURF_BIS',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END FUNCTION ver_interp_lin3d_surf_bis
!     ##############################################
      FUNCTION ver_interp_lin2d_surf(PVAR1,KKLIN,PCOEFLIN) RESULT(PVAR2)
!     ##############################################
!
!!****  *VER_INTERP_LIN* - vertical linear interpolation
!!
!!    PURPOSE
!!    -------
!
!!
!!**  METHOD
!!    ------
!!
!!    This routine calls the 3D version of VER_INTERP_LIN after rewritting of
!!    the fields under 3D form.
!!
!!    EXTERNAL
!!    --------
!!
!!    function VER_INTERP_LIN3D
!!    module   MODI_VER_INTERP_LIN3D
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    Book 2
!!
!!    AUTHOR
!!    ------
!!
!     V.Masson  Meteo-France
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    17/07/97
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modi_ver_interp_lin3d_surf
!
!
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)) :: pvar2 ! variable values on
!                                                       ! target grid 
!
!*       0.2   Declaration of local variables
!              ------------------------------
!
REAL,   DIMENSION(1,SIZE(PVAR1,1),SIZE(PVAR1,2)) :: zvar1 ! variable values on the initial
!                                                         ! grid
REAL,   DIMENSION(1,SIZE(KKLIN,1),SIZE(KKLIN,2)) :: zvar2 ! variable values on target
!
INTEGER,DIMENSION(1,SIZE(KKLIN,1),SIZE(KKLIN,2)) :: iklin ! lower interpolating level of
!                                             ! grid 1 for each level of grid 2 
REAL,   DIMENSION(1,SIZE(PCOEFLIN,1),SIZE(PCOEFLIN,2)):: zcoeflin ! coefficient for level KKLIN
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN_SURF:VER_INTERP_LIN2D_SURF',0,zhook_handle)
zvar1(1,:,:)=pvar1(:,:)
iklin(1,:,:)=kklin(:,:)
zcoeflin(1,:,:)=pcoeflin(:,:)
!
zvar2(:,:,:)=ver_interp_lin3d_surf(zvar1(:,:,:),iklin(:,:,:),zcoeflin(:,:,:))
!
pvar2(:,:)  =zvar2(1,:,:)
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN_SURF:VER_INTERP_LIN2D_SURF',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END FUNCTION ver_interp_lin2d_surf
!     ##############################################
      FUNCTION ver_interp_lin1d_surf(PVAR1,KKLIN,PCOEFLIN) RESULT(PVAR2)
!     ##############################################
!
!!****  *VER_INTERP_LIN* - vertical linear interpolation
!!
!!    PURPOSE
!!    -------
!
!!
!!**  METHOD
!!    ------
!!
!!    This routine calls the 3D version of VER_INTERP_LIN after rewritting of
!!    the fields under 3D form.
!!
!!    EXTERNAL
!!    --------
!!
!!    function VER_INTERP_LIN3D
!!    module   MODI_VER_INTERP_LIN3D
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    Book 2
!!
!!    AUTHOR
!!    ------
!!
!     V.Masson  Meteo-France
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    17/07/97
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modi_ver_interp_lin3d_surf
!
!
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)) :: pvar2 ! variable values on target
!                                         ! grid 
!
!*       0.2   Declaration of local variables
!              ------------------------------
!
INTEGER,DIMENSION(1,1,SIZE(KKLIN)) :: iklin ! lower interpolating level of
!                                             ! grid 1 for each level of grid 2 
REAL,   DIMENSION(1,1,SIZE(PCOEFLIN)) :: zcoeflin ! coefficient for level KKLIN
!
REAL,   DIMENSION(1,1,SIZE(PVAR1)) :: zvar1 ! variable values on the initial
!                                           ! grid
REAL,   DIMENSION(1,1,SIZE(KKLIN)) :: zvar2 ! variable values on target
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN_SURF:VER_INTERP_LIN1D_SURF',0,zhook_handle)
zvar1(1,1,:)=pvar1(:)
iklin(1,1,:)=kklin(:)
zcoeflin(1,1,:)=pcoeflin(:)
!
zvar2(:,:,:)=ver_interp_lin3d_surf(zvar1(:,:,:),iklin(:,:,:),zcoeflin(:,:,:))
!
pvar2(:)    =zvar2(1,1,:)
IF (lhook) CALL dr_hook('MODI_VER_INTERP_LIN_SURF:VER_INTERP_LIN1D_SURF',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END FUNCTION ver_interp_lin1d_surf