read_eco2_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 read_eco2_irrig (&
                                   dtco, &
                                  hprogram)
!     #######################
!
!
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
!
USE modi_read_surf
!
USE modd_data_cover,     ONLY : tdata_seed, tdata_reap, xdata_watsup, xdata_irrig,&
                                  ldata_irrig 
USE modd_data_cover_par, ONLY : jpcover, nvt_irr
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!* dummy arguments
!  ---------------
!
!
!
!
TYPE(data_cover_t), INTENT(INOUT) :: dtco
!
 CHARACTER(LEN=6),     INTENT(IN)    :: hprogram  ! program calling surf. schemes
!
!
!* local variables
!  ---------------
!
 CHARACTER(LEN=12) :: yrecfm     ! Name of the article to be read
INTEGER           :: iresp      ! reading return code
!
INTEGER           :: iversion   ! surface version
INTEGER           :: ibugfix    ! surface bugfix
!
INTEGER           :: jcover     ! loop counter
!
REAL, DIMENSION(6) :: zwork
REAL(KIND=JPRB) :: zhook_handle
!------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('READ_ECO2_IRRIG',0,zhook_handle)
yrecfm='VERSION'
 CALL read_surf(&
                hprogram,yrecfm,iversion,iresp)
yrecfm='BUG'
 CALL read_surf(&
                hprogram,yrecfm,ibugfix,iresp)
!
IF (iversion<4 .OR. iversion==4 .AND. ibugfix<2) THEN
  ldata_irrig  = .false.
  IF (lhook) CALL dr_hook('READ_ECO2_IRRIG',1,zhook_handle)
  RETURN
END IF
!
yrecfm='DATA_IRRIG'
 CALL read_surf(&
                hprogram,yrecfm,ldata_irrig,iresp)
!
IF (.NOT. ldata_irrig .AND. lhook) CALL dr_hook('READ_ECO2_IRRIG',1,zhook_handle)
IF (.NOT. ldata_irrig) RETURN
!
DO jcover=1,jpcover
  IF (dtco%XDATA_VEGTYPE(jcover,nvt_irr)==0.) cycle
  WRITE(yrecfm,fmt='(A6,I3.3)') 'IRRIG_',jcover
  CALL read_surf(&
                hprogram,yrecfm,zwork,iresp,hdir='-')
  tdata_seed(jcover,nvt_irr)%TDATE%MONTH = nint(zwork(1))
  tdata_seed(jcover,nvt_irr)%TDATE%DAY   = nint(zwork(2))
  tdata_reap(jcover,nvt_irr)%TDATE%MONTH = nint(zwork(3))
  tdata_reap(jcover,nvt_irr)%TDATE%DAY   = nint(zwork(4))
  xdata_watsup(jcover,nvt_irr)             = zwork(5)
  xdata_irrig(jcover,nvt_irr)             = zwork(6)
  tdata_seed(jcover,nvt_irr)%TDATE%YEAR  = 9999
  tdata_seed(jcover,nvt_irr)%TIME        = 0.
  tdata_reap(jcover,nvt_irr)%TDATE%YEAR  = 9999
  tdata_reap(jcover,nvt_irr)%TIME        = 0.
END DO
IF (lhook) CALL dr_hook('READ_ECO2_IRRIG',1,zhook_handle)
!
!------------------------------------------------------------------------------
!
END SUBROUTINE read_eco2_irrig