make_mask_isba_to_topd.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 make_mask_isba_to_topd(KI)
!     #######################
!
!!****  *MAKE_MASK_ISBA_TO_TOPD*  
!!
!!    PURPOSE
!!    -------
!
!     Create a mask for each Surfex mesh and each catchment. 
!         
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!                          
!!    REFERENCE
!!    ---------
     
!!    AUTHOR
!!    ------
!!
!!      K. Chancibault  * CNRM * 
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      Original   16/03/2005
!!                 11/2011 : Loops simplified (Vincendon)
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_topodyn,       ONLY : nncat, nnmc
USE modd_coupling_topd, ONLY : nmaskt, nmaski, nnpix
USE modd_surf_par,        ONLY : nundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
INTEGER, INTENT(IN) :: ki    ! Grid dimensions
!
!*      0.2    declarations of local variables
!
INTEGER, DIMENSION(KI)  :: inbpix_in_mesh ! number of pixel in each ISBA mesh
INTEGER                 :: jcat, jpix, inumpix
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MAKE_MASK_ISBA_TO_TOPD',0,zhook_handle)
!
inumpix=maxval(nnpix)
!
ALLOCATE(nmaski(ki,nncat,inumpix))
nmaski(:,:,:) = nundef
!
inbpix_in_mesh(:) = 0
!
DO jcat = 1,nncat
  !
  DO jpix = 1,nnmc(jcat)
    !si le point du bassin versant est dans une maille isba
    IF ((nmaskt(jcat,jpix)/=0).AND.(nmaskt(jcat,jpix)/=nundef)) THEN
      !indice du point du bassin versant dans la maille isba
      inbpix_in_mesh(nmaskt(jcat,jpix)) = inbpix_in_mesh(nmaskt(jcat,jpix)) + 1
      ! nmaski associe à la maille isba, au bassin versant et au numéro du point 
      ! du bassin versant dans la maille isba, l'indice du point dans le bassin
      ! versant
      nmaski(nmaskt(jcat,jpix),jcat,inbpix_in_mesh(nmaskt(jcat,jpix))) = jpix
    ENDIF
    !
  ENDDO
  !
ENDDO
! write(*,*) 'NMASKT min et max',MINVAL(NMASKT(1,:)),MAXVAL(NMASKT(1,:))
! write(*,*) 'NMASKT min et max',MINVAL(NMASKT(2,:)),MAXVAL(NMASKT(2,:))
! write(*,*) 'NMASKT min et max',MINVAL(NMASKT(3,:)),MAXVAL(NMASKT(3,:))
! write(*,*) 'NMASKT min et max',MINVAL(NMASKT(4,:)),MAXVAL(NMASKT(4,:))
! write(*,*) 'NMASKI 3132 min et max',MINVAL(NMASKI(MINVAL(NMASKT(1,:)),1,:)),MAXVAL(NMASKI(MINVAL(NMASKT(1,:)),1,:))
! write(*,*) 'NMASKI 6662 min et max',MINVAL(NMASKI(MAXVAL(NMASKT(1,:)),1,:)),MAXVAL(NMASKI(MAXVAL(NMASKT(1,:)),1,:))
! stop
!
IF (lhook) CALL dr_hook('MAKE_MASK_ISBA_TO_TOPD',1,zhook_handle)
!
END SUBROUTINE make_mask_isba_to_topd