writesurf_pgd_teb_gardenn.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_garden_n (HSELECT, TOP, T, KTIME, IO, K, P, HPROGRAM)
!     ###############################################
!
!!****  *WRITE_PGD_TEB_GARDEN_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_teb_option_n, ONLY : teb_options_t
USE modd_teb_n, ONLY : teb_t
!
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY : isba_k_t, isba_p_t
!
USE modd_surf_par,          ONLY : xundef, nundef
!
USE modi_write_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
 CHARACTER(LEN=*), DIMENSION(:), INTENT(IN) :: HSELECT
!
TYPE(teb_options_t), INTENT(INOUT) :: TOP
TYPE(teb_t), INTENT(INOUT) :: T
!
INTEGER, INTENT(IN) :: KTIME
!
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(isba_k_t), INTENT(INOUT) :: K
TYPE(isba_p_t), INTENT(INOUT) :: P
!
 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
!
REAL, DIMENSION(:), ALLOCATABLE :: ZWORK
INTEGER :: JJ, JL
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_TEB_GARDEN_N',0,zhook_handle)
!
!* soil scheme option
!
yrecfm='GD_ISBA'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%CISBA,iresp,hcomment=ycomment)
!
!* Reference grid for DIF
!
IF(io%CISBA=='DIF') THEN
  DO jl=1,io%NGROUND_LAYER
     WRITE(ylvl,'(I4)') jl     
     yrecfm='GD_SGRID'//adjustl(ylvl(:len_trim(ylvl)))
     ycomment='Depth of TEB Garden soilgrid layer '//adjustl(ylvl(:len_trim(ylvl)))
     CALL write_surf(hselect, hprogram,yrecfm,io%XSOILGRID(jl),iresp,hcomment=ycomment)
  END DO 
ENDIF
!
!* number of soil layers
!
yrecfm='GD_LAYER'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%NGROUND_LAYER,iresp,hcomment=ycomment)
!
!* number of time data for vegetation characteristics (VEG, LAI, EMIS, Z0) 
!
yrecfm='GD_NTIME'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,ktime,iresp,hcomment=ycomment)
!
! * clay fraction
!
yrecfm='GD_CLAY'
ycomment='X_Y_GD_CLAY'
 CALL write_surf(hselect,hprogram,yrecfm,k%XCLAY(:,1),iresp,hcomment=ycomment)
!        
! * sand fraction
!
yrecfm='GD_SAND'
ycomment='X_Y_GD_SAND'
 CALL write_surf(hselect,hprogram,yrecfm,k%XSAND(:,1),iresp,hcomment=ycomment)
!        
! * orographic runoff coefficient
!
yrecfm='GD_RUNOFFB'
ycomment='X_Y_GD_RUNOFFB'
 CALL write_surf(hselect,hprogram,yrecfm,k%XRUNOFFB,iresp,hcomment=ycomment)
!        
! * subgrid drainage coefficient
!
yrecfm='GD_WDRAIN'
ycomment='X_Y_GD_WDRAIN'
 CALL write_surf(hselect,hprogram,yrecfm,k%XWDRAIN,iresp,hcomment=ycomment)
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_TEB_GARDEN_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE writesurf_pgd_teb_garden_n