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 (USS, &
                                     kluout,plat,plon,pvalue,hsubroutine,knblines,pnodata)
!     ###################################################################
!
!!**** *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_sso_n, ONLY : surf_atm_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_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
REAL(KIND=JPRB) :: zhook_handle
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER :: inblines
!----------------------------------------------------------------------------
!
      IF (lhook) CALL dr_hook('PT_BY_PT_TREATMENT',0,zhook_handle)
      inblines = 1
      IF (present(knblines)) inblines = knblines

      SELECT CASE (hsubroutine)

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

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

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

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