writesurf_pgd_teb_irrign.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 writesurf_pgd_teb_irrig_n (DGU, U, &
                                             tir, &
                                            hprogram)
!     ################################################
!
!!****  *WRITESURF_PGD_TEB_IRRIG_n* - writes TEB irrigation physiographic fields
!!                        
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    05/2013
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_teb_irrig_n, ONLY : teb_irrig_t
!
USE modi_write_surf
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
!
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
TYPE(teb_irrig_t), INTENT(INOUT) :: tir
!
 CHARACTER(LEN=6),  INTENT(IN)  :: hprogram ! program calling
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: iresp          ! IRESP  : return-code if a problem appears
 CHARACTER(LEN=12) :: yrecfm         ! Name of the article to be read
 CHARACTER(LEN=100):: ycomment       ! Comment string
INTEGER           :: jlayer         ! loop index
INTEGER           :: jtime          ! loop index
REAL, DIMENSION(:), ALLOCATABLE :: zwork
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_TEB_IRRIG_n',0,zhook_handle)
!
! Flag for irrigation of gardens
yrecfm='L_PAR_GD_IRR'
ycomment='FLAG FOR SPECIFIED GARDEN IRRIGATION PARAMETERS'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%LPAR_GD_IRRIG,iresp,hcomment=ycomment)
!
! Parameters describing irrigation
IF (tir%LPAR_GD_IRRIG) THEN
!
  yrecfm='D_GD_SM_IRR'
  ycomment='Start Month for Gardens Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGD_START_MONTH(:),iresp,hcomment=ycomment)
!
  yrecfm='D_GD_EM_IRR'
  ycomment='End   Month for Gardens Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGD_END_MONTH  (:),iresp,hcomment=ycomment)
!
  yrecfm='D_GD_SH_IRR'
  ycomment='Start Hour  for Gardens Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGD_START_HOUR (:),iresp,hcomment=ycomment)
!
  yrecfm='D_GD_EH_IRR'
  ycomment='End   Hour  for Gardens Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGD_END_HOUR   (:),iresp,hcomment=ycomment)
!
  yrecfm='D_GD_IRRIG'
  ycomment='24h mean Irrigation rate for Gardens Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGD_24H_IRRIG  (:),iresp,hcomment=ycomment)
!
ENDIF
!
! Flag for irrigation of greenroofs
yrecfm='L_PAR_GR_IRR'
ycomment='FLAG FOR SPECIFIED GREENROOFS IRRIGATION PARAMETERS'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%LPAR_GR_IRRIG,iresp,hcomment=ycomment)
!
! Parameters describing irrigation
IF (tir%LPAR_GR_IRRIG) THEN
!
  yrecfm='D_GR_SM_IRR'
  ycomment='Start Month for Greenroofs Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGR_START_MONTH(:),iresp,hcomment=ycomment)
!
  yrecfm='D_GR_EM_IRR'
  ycomment='End   Month for Greenroofs Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGR_END_MONTH  (:),iresp,hcomment=ycomment)
!
  yrecfm='D_GR_SH_IRR'
  ycomment='Start Hour  for Greenroofs Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGR_START_HOUR (:),iresp,hcomment=ycomment)
!
  yrecfm='D_GR_EH_IRR'
  ycomment='End   Hour  for Greenroofs Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGR_END_HOUR   (:),iresp,hcomment=ycomment)
!
  yrecfm='D_GR_IRRIG'
  ycomment='24h mean Irrigation rate for Greenroofs Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XGR_24H_IRRIG  (:),iresp,hcomment=ycomment)
!
ENDIF
!
! Flag for watering of greenroofs
yrecfm='L_PAR_RD_IRR'
ycomment='FLAG FOR SPECIFIED ROAD IRRIGATION PARAMETERS'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%LPAR_RD_IRRIG,iresp,hcomment=ycomment)
!
! Parameters describing watering
IF (tir%LPAR_RD_IRRIG) THEN
!
  yrecfm='D_RD_SM_IRR'
  ycomment='Start Month for Roads Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XRD_START_MONTH(:),iresp,hcomment=ycomment)
!
  yrecfm='D_RD_EM_IRR'
  ycomment='End   Month for Roads Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XRD_END_MONTH  (:),iresp,hcomment=ycomment)
!
  yrecfm='D_RD_SH_IRR'
  ycomment='Start Hour  for Roads Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XRD_START_HOUR (:),iresp,hcomment=ycomment)
!
  yrecfm='D_RD_EH_IRR'
  ycomment='End   Hour  for Roads Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XRD_END_HOUR   (:),iresp,hcomment=ycomment)
!
  yrecfm='D_RD_IRRIG'
  ycomment='24h mean Irrigation rate for Roads Irrigation'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tir%XRD_24H_IRRIG  (:),iresp,hcomment=ycomment)
!
ENDIF
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_TEB_IRRIG_n',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE writesurf_pgd_teb_irrig_n