detect_field.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 detect_field(HPROGRAM,PFIELD,OITSHERE)
!     ################################################
!
!!****  *DETECT_FIELD* - generic routine to check is a field is non-zero
!!
!!    PURPOSE
!!    -------
!
!
!!**  METHOD
!!    ------
!!    
!!    
!!
!!    EXTERNAL
!!    --------
!!
!!       
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!    REFERENCE
!!    ---------
!!
!!      
!!
!!    AUTHOR
!!    ------
!!      S.Malardel       * Meteo France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original      02/2003
!-----------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!
!
#ifdef SFX_MNH
USE modi_detect_field_mnh
#endif
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   declarations of arguments
!
 CHARACTER (LEN=6),  INTENT(IN) :: hprogram   ! program
!
REAL, DIMENSION(:,:), INTENT(IN)::pfield ! array containing the data field
!
LOGICAL    , INTENT(OUT)         :: oitshere  ! T --> PFIELD is non zero somewhere
REAL(KIND=JPRB) :: zhook_handle
!
!
!*       0.2   declarations of local variables
!
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('DETECT_FIELD',0,zhook_handle)
IF (hprogram=='MESONH') THEN
!
#ifdef SFX_MNH
  CALL detect_field_mnh(hprogram,SIZE(pfield,1),SIZE(pfield,2),pfield,oitshere)
#endif
!
ELSEIF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. hprogram=='TEXTE ' &
   .OR. hprogram=='FA    ' .OR. hprogram=='BINARY' .OR. hprogram=='LFI   ' &
   .OR. hprogram=='NC    ' ) THEN
!
  oitshere = .true.
!
ELSEIF (hprogram=='AROME ') THEN
!
#ifdef SFX_ARO
  CALL detect_field_aro(hprogram,SIZE(pfield,1),SIZE(pfield,2),pfield,oitshere)
#endif
!
ENDIF
IF (lhook) CALL dr_hook('DETECT_FIELD',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE detect_field