teb_irrig.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 teb_irrig(OIRRIG, PTSTEP, KMONTH, PSOLAR_TIME,   &
               pstart_month, pend_month, pstart_hour, pend_hour,&
               p24h_irrig, pirrig                               )
!   ##########################################################################
!
!!****  *TEB_IRRIG*  
!!
!!    PURPOSE
!!    -------
!
! Computes the amount of water for irrigation
!         
!     
!!**  METHOD
!     ------
!
!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!      
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!
!!      C. DeMunck & V. Masson           * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    05/2013
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
LOGICAL,            INTENT(IN)  :: oirrig   ! Flag for irrigation or not
REAL,               INTENT(IN)  :: ptstep   ! time step
INTEGER,            INTENT(IN)  :: kmonth   ! Present Month
REAL, DIMENSION(:), INTENT(IN)  :: psolar_time ! solar time (s since midnight)
REAL, DIMENSION(:), INTENT(IN)  :: pstart_month! First month for irrigation (included)
REAL, DIMENSION(:), INTENT(IN)  :: pend_month  ! Last  month for irrigation (included)
REAL, DIMENSION(:), INTENT(IN)  :: pstart_hour ! First month for irrigation (included)
REAL, DIMENSION(:), INTENT(IN)  :: pend_hour   ! Last  month for irrigation (ecluded)
REAL, DIMENSION(:), INTENT(IN)  :: p24h_irrig  ! 24h total irrigation
REAL, DIMENSION(:), INTENT(OUT) :: pirrig      ! present irrigation
!
REAL(KIND=JPRB) :: zhook_handle
!
!*      0.2    declarations of local variables
!
LOGICAL, DIMENSION(SIZE(PIRRIG)):: gmonth!if irrigation occurs during present month
LOGICAL, DIMENSION(SIZE(PIRRIG)):: ghour !if irrigation occurs at     present time
REAL,    DIMENSION(SIZE(PIRRIG)):: zduration ! duration of irrigation per day (s)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('TEB_IRRIG',0,zhook_handle)
!
!*      1.     Default
!              -------
!
pirrig = 0.
!
IF (.NOT. oirrig) THEN
  IF (lhook) CALL dr_hook('TEB_IRRIG',1,zhook_handle)
  RETURN
END IF
!
!-------------------------------------------------------------------------------
!
!*      2.     Is present month irrigated ?
!              ----------------------------
!
WHERE (pstart_month <= pend_month)   ! for summer in North hemisphere
  gmonth = kmonth>=pstart_month .AND. kmonth<= pend_month
ELSEWHERE                            ! for summer in South hemisphere
                                     ! change of year during irrigation period
  gmonth = kmonth>=pstart_month .OR.  kmonth<= pend_month
END WHERE
!
!-------------------------------------------------------------------------------
!
!*      3.     Is present time irrigated (for each location) ?
!              -----------------------------------------------
!
ghour = .false.
zduration = 0.
!
WHERE (pstart_hour  <= pend_hour)      ! typically irrigation during day
  ghour  = psolar_time+ptstep>=pstart_hour*3600. .AND. psolar_time<=pend_hour*3600.
  zduration = 3600.* (pend_hour-pstart_hour)
ELSE WHERE                             ! typically irrigation during night
                                       ! midnight is during irrigation period
  ghour  = psolar_time+ptstep>=pstart_hour*3600. .OR.  psolar_time< pend_hour*3600.
  zduration = 3600.* (pend_hour + 24.-pstart_hour)
END WHERE
!
!-------------------------------------------------------------------------------
!
!*      4.     Computes irrigation
!              -------------------
!
!* if duration is very short, all water is irrigated during the timestep
zduration = max(zduration, ptstep)
!
!* Concentration of the total 24h irrigation during the period of irrigation
WHERE(gmonth(:) .AND. ghour(:)) pirrig(:) = p24h_irrig(:) / zduration(:)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('TEB_IRRIG',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE teb_irrig