pt_by_pt_treatment.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 pt_by_pt_treatment (UG, U, USS, &
                                     KLUOUT,PLAT,PLON,PVALUE,HSUBROUTINE,&
                                     KNBLINES,PNODATA,OMULTITYPE,KFACT)
!     ###################################################################
!
!!**** *PT_BY_PT_TREATMENT* 
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson          Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    12/09/95
!!                27/03/96 (V. Masson) modify the arguments for the call of 
!!                         interpolation subroutine
!!                06/2009 (B. Decharme) call Topographic index statistics calculation
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_sso_n, ONLY : sso_t
!
USE modi_average1_cover
USE modi_average1_orography
USE modi_average1_cti
USE modi_average1_ldb
USE modi_average1_mesh
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
TYPE(surf_atm_t), INTENT(INOUT) :: U
TYPE(sso_t), INTENT(INOUT) :: USS
!
INTEGER,           INTENT(IN) :: KLUOUT
REAL,DIMENSION(:), INTENT(IN) :: PLAT
REAL,DIMENSION(:), INTENT(IN) :: PLON
REAL,DIMENSION(:), INTENT(IN) :: PVALUE
 CHARACTER(LEN=6), INTENT(IN)  :: HSUBROUTINE   ! Name of the subroutine to call
INTEGER, OPTIONAL, INTENT(IN) :: KNBLINES
REAL, OPTIONAL,    INTENT(IN) :: PNODATA
LOGICAL, OPTIONAL, INTENT(IN) :: OMULTITYPE
INTEGER, OPTIONAL, INTENT(IN) :: KFACT
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
LOGICAL :: GMULTITYPE
INTEGER :: ITYPE, IFACT
INTEGER :: INBLINES
!----------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PT_BY_PT_TREATMENT',0,zhook_handle)
!
inblines = 0
IF (PRESENT(knblines)) inblines = knblines
!
gmultitype = .false.
IF (PRESENT(omultitype)) gmultitype = omultitype
!
ifact = 1
IF (PRESENT(kfact)) ifact = kfact
!
SELECT CASE (hsubroutine)

  CASE ('A_COVR')
    IF (PRESENT(pnodata)) THEN
      CALL average1_cover(ug, u, kluout,inblines,plat,plon,pvalue,pnodata)
    ELSE
      CALL average1_cover(ug, u, kluout,inblines,plat,plon,pvalue)
    ENDIF

  CASE ('A_OROG')
    IF (PRESENT(pnodata)) THEN
      CALL average1_orography(ug, uss, kluout,inblines,plat,plon,pvalue,pnodata)
    ELSE
      CALL average1_orography(ug, uss, kluout,inblines,plat,plon,pvalue)
    ENDIF

  CASE ('A_CTI ')
    IF (PRESENT(pnodata)) THEN
      CALL average1_cti(ug, kluout,inblines,plat,plon,pvalue,pnodata)
    ELSE
      CALL average1_cti(ug, kluout,inblines,plat,plon,pvalue)
    ENDIF

  CASE ('A_LDBD')
    IF (PRESENT(pnodata)) THEN
      CALL average1_ldb(ug, kluout,inblines,plat,plon,pvalue,'D',pnodata)
    ELSE
      CALL average1_ldb(ug, kluout,inblines,plat,plon,pvalue,'D')
    ENDIF
        
  CASE ('A_LDBS')
    IF (PRESENT(pnodata)) THEN
      CALL average1_ldb(ug, kluout,inblines,plat,plon,pvalue,'S',pnodata)
    ELSE
      CALL average1_ldb(ug, kluout,inblines,plat,plon,pvalue,'S')
    ENDIF
        
  CASE ('A_MESH')
    IF (PRESENT(pnodata)) THEN
      CALL average1_mesh(ug, kluout,inblines,plat,plon,pvalue,gmultitype,ifact,pnodata)
    ELSE
      CALL average1_mesh(ug,kluout,inblines,plat,plon,pvalue,gmultitype,ifact)
    ENDIF
        
END SELECT
IF (lhook) CALL dr_hook('PT_BY_PT_TREATMENT',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE pt_by_pt_treatment