writesurf_pgd_teb_vegn.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_veg_n (DGU, U, &
                                           dtgd, tgdo, tgdp, tvg, &
                                          hprogram)
!     ###############################################
!
!!****  *WRITE_PGD_TEB_VEG_n* - writes ISBA fields describing urban gardens
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      A. Lemonsu & C. de Munck   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    06/2011 
!!
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_data_teb_garden_n, ONLY : data_teb_garden_t
USE modd_teb_garden_option_n, ONLY : teb_garden_options_t
USE modd_teb_garden_pgd_n, ONLY : teb_garden_pgd_t
USE modd_teb_veg_n, ONLY : teb_veg_options_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(data_teb_garden_t), INTENT(INOUT) :: dtgd
TYPE(teb_garden_options_t), INTENT(INOUT) :: tgdo
TYPE(teb_garden_pgd_t), INTENT(INOUT) :: tgdp
TYPE(teb_veg_options_t), INTENT(INOUT) :: tvg
!
 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
 CHARACTER(LEN=4 ) :: ylvl
!
INTEGER :: jj, jlayer
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_TEB_VEG_N',0,zhook_handle)
!
!* soil scheme option
!
yrecfm='GD_ISBA'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tvg%CISBA,iresp,hcomment=ycomment)
!
!* Reference grid for DIF
!
IF(tvg%CISBA=='DIF') THEN
  DO jlayer=1,tgdo%NGROUND_LAYER
     WRITE(ylvl,'(I4)') jlayer     
     yrecfm='GD_SGRID'//adjustl(ylvl(:len_trim(ylvl)))
     ycomment='Depth of TEB Garden soilgrid layer '//adjustl(ylvl(:len_trim(ylvl)))
     CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tgdo%XSOILGRID(jlayer),iresp,hcomment=ycomment)
  END DO 
ENDIF
!
!* number of soil layers
!
yrecfm='GD_LAYER'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tgdo%NGROUND_LAYER,iresp,hcomment=ycomment)
!
!* number of time data for vegetation characteristics (VEG, LAI, EMIS, Z0) 
!
yrecfm='GD_NTIME'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,dtgd%NTIME,iresp,hcomment=ycomment)
!
! * clay fraction
!
yrecfm='GD_CLAY'
ycomment='X_Y_GD_CLAY'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tgdp%XCLAY(:,1),iresp,hcomment=ycomment)
!        
! * sand fraction
!
yrecfm='GD_SAND'
ycomment='X_Y_GD_SAND'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tgdp%XSAND(:,1),iresp,hcomment=ycomment)
!        
! * orographic runoff coefficient
!
yrecfm='GD_RUNOFFB'
ycomment='X_Y_GD_RUNOFFB'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tgdp%XRUNOFFB,iresp,hcomment=ycomment)
!        
! * subgrid drainage coefficient
!
yrecfm='GD_WDRAIN'
ycomment='X_Y_GD_WDRAIN'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,tgdp%XWDRAIN,iresp,hcomment=ycomment)
!
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_TEB_VEG_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE writesurf_pgd_teb_veg_n