write_ecoclimap2_data.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 write_ecoclimap2_data (DGU, U, &
                                        hprogram)
!     #######################
!
!
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modi_write_surf
!
USE modd_data_cover,     ONLY : tdata_seed, tdata_reap, xdata_watsup, xdata_irrig,&
                                  ldata_irrig, xdata_vegtype, lclim_lai  
USE modd_data_cover_par, ONLY : jpcover, nvt_irr
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!* dummy arguments
!  ---------------
!
!
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
 CHARACTER(LEN=6),     INTENT(IN)    :: hprogram  ! program calling surf. schemes
!
!
!* local variables
!  ---------------
!
 CHARACTER(LEN=12) :: yrecfm     ! Name of the article to be read
 CHARACTER(LEN=100):: ycomment   ! Comment
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('WRITE_ECOCLIMAP2_DATA',0,zhook_handle)
yrecfm='DATA_IRRIG'
ycomment='FLAG TO READ USER IRRIGATION DATA FOR ECOCLIMAP2'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,ldata_irrig,iresp,ycomment)
!
yrecfm='LCLIM_LAI'
ycomment='FLAG TO USE CLIMATOLOGICAL LAI'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,lclim_lai,iresp,ycomment)
!
IF (.NOT. ldata_irrig .AND. lhook) CALL dr_hook('WRITE_ECOCLIMAP2_DATA',1,zhook_handle)
IF (.NOT. ldata_irrig) RETURN
!
DO jcover=1,jpcover
  IF (xdata_vegtype(jcover,nvt_irr)==0.) cycle
  WRITE(yrecfm,fmt='(A6,I3.3)') 'IRRIG_',jcover
  WRITE(ycomment,fmt='(A47,I3.3)') &
    'SEED MONTH&DAY, REAP MONTH&DAY, WATSUP, IRRIG  ',jcover  
  zwork(1) = tdata_seed(jcover,nvt_irr)%TDATE%MONTH
  zwork(2) = tdata_seed(jcover,nvt_irr)%TDATE%DAY
  zwork(3) = tdata_reap(jcover,nvt_irr)%TDATE%MONTH
  zwork(4) = tdata_reap(jcover,nvt_irr)%TDATE%DAY
  zwork(5) = xdata_watsup(jcover,nvt_irr)
  zwork(6) = xdata_irrig(jcover,nvt_irr)
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,zwork,iresp,ycomment,'-',"Irrig_parameters")
END DO
IF (lhook) CALL dr_hook('WRITE_ECOCLIMAP2_DATA',1,zhook_handle)
!
!------------------------------------------------------------------------------
!
END SUBROUTINE write_ecoclimap2_data