get_ign_maskall.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 get_ign_maskall (UG, KNI, PX, PY, OTOT)
!     #######################################################
!!****  *GET_IGN_MASKALL* - 
!!
!!    PURPOSE
!!    -------
!!
      !!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      S. Faroux   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    06/2010 
!!      07/2011     add specific computation for IGN grid (B. Decharme)
!-------------------------------------------------------------------------------                         
!
!
!
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
!
#ifdef SFX_OL
USE modd_io_surf_ol, ONLY: nmask_ign
#endif
!
USE mode_gridtype_ign
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
!
INTEGER, INTENT(IN)              :: KNI
REAL,DIMENSION(:), INTENT(OUT)   :: PX, PY
LOGICAL, INTENT(IN), OPTIONAL :: OTOT
!
REAL, DIMENSION(KNI)             :: ZXX, ZYY
INTEGER                          :: JI, JJ, JK, JL
LOGICAL :: GTOT
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('GET_IGN_MASKALL',0,zhook_handle)
!
gtot = .true.
IF (PRESENT(otot)) gtot = otot
!
IF (ASSOCIATED(ug%XGRID_FULL_PAR).AND.gtot) THEN
  CALL get_gridtype_ign(ug%XGRID_FULL_PAR,px=zxx,py=zyy,pxall=px,pyall=py)
ENDIF
IF (.NOT.gtot .AND. ASSOCIATED(ug%G%XGRID_PAR)) THEN
  IF (ASSOCIATED(ug%XGRID_FULL_PAR)) THEN
    CALL get_gridtype_ign(ug%XGRID_FULL_PAR,pxall=px,pyall=py)
    CALL get_gridtype_ign(ug%G%XGRID_PAR,px=zxx,py=zyy)
  ELSE
    CALL get_gridtype_ign(ug%G%XGRID_PAR,px=zxx,py=zyy,pxall=px,pyall=py)
  ENDIF
ENDIF
!
#ifdef SFX_OL
IF (.NOT.ALLOCATED(nmask_ign))THEN
  ALLOCATE(nmask_ign(kni))
  jl=0
  DO jj=1,SIZE(py)  
    DO ji=1,SIZE(px)
      jl=jl+1
      DO jk=1,kni
        IF((zxx(jk)==px(ji)).AND.(zyy(jk)==py(jj)))THEN
          nmask_ign(jk) = jl
        ENDIF
      ENDDO
    ENDDO
  ENDDO
ENDIF
#endif
!
IF (lhook) CALL dr_hook('GET_IGN_MASKALL',1,zhook_handle)
!
END SUBROUTINE get_ign_maskall