topd_to_isba_slope.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 topd_to_isba_slope (USS, &
                                     ki)
!     ####################
!
!!****  *TOPD_TO_ISBA*  
!!
!!    PURPOSE
!!    -------
!    
!     
!         
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!    
!!    
!!
!!      
!!    REFERENCE
!!    ---------
!!
!!    
!!      
!!    AUTHOR
!!    ------
!!
!!      B. Vincendon    * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      Original   12/11/2012
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
USE modd_surf_atm_sso_n, ONLY : surf_atm_sso_t
!
USE modd_topodyn,       ONLY : nncat, nnmc, xtanb
USE modd_coupling_topd, ONLY : nmaskt,nnpix
USE modd_surf_par,        ONLY : nundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
!
TYPE(surf_atm_sso_t), INTENT(INOUT) :: uss
!
INTEGER, INTENT(IN)                 :: ki      ! Grid dimensions
!
!*      0.2    declarations of local variables
!
!
INTEGER                :: jcat,jpix,jj          ! loop control 
REAL, DIMENSION(KI)    :: zcount                ! TOPO pixel number in an ISBA pixel
                                                ! on the full grid
REAL, DIMENSION(KI)    :: zsso_slope
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('TOPD_TO_ISBA_SLOPE',0,zhook_handle)
!
!*        1.0  Compute Mean slope over each ISBA_MESH
!            ----------------------------------------------------------------------
!
!write(*,*) 'pente avt topmodel',MINVAL(XSSO_SLOPE),MAXVAL(XSSO_SLOPE),SUM(XSSO_SLOPE,MASK=XSSO_SLOPE/=XUNDEF)
!
zsso_slope = uss%XSSO_SLOPE
!
zcount(:) = REAL(nnpix(:))

WHERE (zcount /= 0.0)
   zsso_slope = 0.
ENDWHERE
!
DO jcat=1,nncat
  DO jpix=1,nnmc(jcat)
    IF (nmaskt(jcat,jpix) /= nundef) THEN
      zsso_slope(nmaskt(jcat,jpix)) = zsso_slope(nmaskt(jcat,jpix)) + xtanb(jcat,jpix)
    ENDIF
  ENDDO
ENDDO
!
WHERE (zcount /= 0.0)
   zsso_slope = zsso_slope / zcount
ENDWHERE
!
uss%XSSO_SLOPE = zsso_slope
!
!write(*,*) 'pente apres modification', &
!           MINVAL(XSSO_SLOPE),MAXVAL(XSSO_SLOPE),COUNT(ZCOUNT/=0.0),SUM(XSSO_SLOPE,MASK=XSSO_SLOPE/=XUNDEF)
!
IF (lhook) CALL dr_hook('TOPD_TO_ISBA_SLOPE',1,zhook_handle)
!
END SUBROUTINE topd_to_isba_slope