prep_teb.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 prep_teb (DTCO, UG, U, USS, GCP, TOP, BOP, NB, TG, SB, NT, GDM, GRM, &
                     HPROGRAM,HATMFILE,HATMFILETYPE,HPGDFILE,HPGDFILETYPE,YDCTL)
!     #################################################################################
!
!!****  *PREP_TEB* - prepares TEB fields
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     V. Masson 
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      S. Riette   06/2009 PREP_TEB_CANOPY has no more argument
!!      P. Marguinaud10/2014, Support for a 2-part PREP
!!------------------------------------------------------------------
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_sso_n, ONLY : sso_t
USE modd_grid_conf_proj_n, ONLY : grid_conf_proj_t
!
USE modd_teb_option_n, ONLY : teb_options_t
USE modd_bem_option_n, ONLY : bem_options_t
USE modd_bem_n, ONLY : bem_np_t
USE modd_sfx_grid_n, ONLY : grid_t
USE modd_canopy_n, ONLY : canopy_t
USE modd_teb_n, ONLY : teb_np_t
USE modd_surfex_n, ONLY : teb_garden_model_t, teb_greenroof_model_t
!
USE mode_prep_ctl, ONLY : prep_ctl
!
USE modi_prep_hor_teb_field
USE modi_prep_ver_teb
USE modi_prep_output_grid
USE modi_get_luout
USE modi_prep_sbl
USE modi_prep_teb_garden
USE modi_prep_teb_greenroof
!
USE modn_prep_teb
!
USE modd_read_namelist, ONLY : lnam_read
!
USE modd_prep,       ONLY : xzs_ls
!
USE modd_prep_teb_garden, ONLY : xwsnow_gd, xrsnow_gd, xtsnow_gd, xlwcsnow_gd, &
                                 xagesnow_gd
!
USE modd_prep_teb_greenroof, ONLY : xwsnow_gr, xrsnow_gr, xtsnow_gr, xlwcsnow_gr, &
                                    xagesnow_gr
!
USE modd_surf_atm,   ONLY : lvertshift
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_clean_prep_output_grid
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
TYPE(surf_atm_t), INTENT(INOUT) :: U
TYPE(sso_t), INTENT(INOUT) :: USS
TYPE(grid_conf_proj_t),INTENT(INOUT) :: GCP
!
TYPE(teb_options_t), INTENT(INOUT) :: TOP
TYPE(bem_options_t), INTENT(INOUT) :: BOP
TYPE(bem_np_t), INTENT(INOUT) :: NB
TYPE(grid_t), INTENT(INOUT) :: TG
TYPE(canopy_t), INTENT(INOUT) :: SB
TYPE(teb_np_t), INTENT(INOUT) :: NT
TYPE(teb_garden_model_t), INTENT(INOUT) :: GDM
TYPE(teb_greenroof_model_t), INTENT(INOUT) :: GRM
!
type(prep_ctl),    INTENT(INOUT) :: ydctl
!
 CHARACTER(LEN=6),   INTENT(IN)  :: HPROGRAM  ! program calling surf. schemes
 CHARACTER(LEN=28),  INTENT(IN)  :: HATMFILE    ! name of the Atmospheric file
 CHARACTER(LEN=6),   INTENT(IN)  :: HATMFILETYPE! type of the Atmospheric file
 CHARACTER(LEN=28),  INTENT(IN)  :: HPGDFILE    ! name of the Atmospheric file
 CHARACTER(LEN=6),   INTENT(IN)  :: HPGDFILETYPE! type of the Atmospheric file
!
!*      0.2    declarations of local variables
!
INTEGER :: ILUOUT
INTEGER :: JP         ! TEB patch number
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------------
!
!*      1.     Default of configuration
!
!
IF (lhook) CALL dr_hook('PREP_TEB',0,zhook_handle)
 CALL get_luout(hprogram,iluout)
!
 CALL prep_output_grid(ug%G, tg, u%NSIZE_FULL, iluout)
!
!-------------------------------------------------------------------------------------
!
!*      2.     Reading and horizontal interpolations
!
!
!
!* option for roads
!
top%CROAD_DIR = croad_dir
top%CWALL_OPT = cwall_opt
!
DO jp=1,top%NTEB_PATCH
  !
  !*      2.0    Large scale orography
  !
  IF (jp==1) &
    CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                            hprogram,'ZS     ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype,1,ydctl)


  !*      2.1    Water reservoirs
  !
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'WS_ROOF',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'WS_ROAD',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  !
  !*      2.2    Building temperature
  !
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'TI_BLD ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  !
  !*      2.3    Road deep temperature
  !
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'TI_ROAD',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  !
  !*      2.4    Temperature profiles
  !
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_ROAD ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_WALLA',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_WALLB',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_ROOF ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  !
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_WIN1 ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  IF (top%CBEM == 'BEM') THEN
    CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'QI_BLD ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
    CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_WIN2 ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
    CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_FLOOR',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
    CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_MASS ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  ENDIF  
  !*      2.5    Snow variables
  !
  nt%AL(jp)%TSNOW_ROOF%SCHEME='1-L'
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'SN_ROOF',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  nt%AL(jp)%TSNOW_ROAD%SCHEME='1-L'
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'SN_ROAD',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  !
  !*      2.6    Canyon air variables
  !
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'T_CAN  ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  CALL prep_hor_teb_field(nb%AL(jp), bop, dtco, u, gcp, tg, nt%AL(jp), top, &
                         hprogram,'Q_CAN  ',hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  !
  !-------------------------------------------------------------------------------------
  !
  !*      3.     Vertical interpolations of all variables
  !
  IF(lvertshift)THEN
    CALL prep_ver_teb(nb%AL(jp), nt%AL(jp), top%XZS, top%CBEM)
  ENDIF
  !
  !-------------------------------------------------------------------------------------
  !
  !*      4.     Urban green areas
  !
  
  IF (top%LGARDEN)    CALL prep_teb_garden(dtco, ug, u, uss,  gcp, tg, top, &
                                           gdm%O, gdm%S, gdm%K, gdm%P, gdm%NPE%AL(jp),  &
                                              hprogram,hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
  IF (top%LGREENROOF) THEN
     CALL prep_teb_greenroof(dtco, ug, u, uss,  gcp, tg, top, grm%O, grm%S, grm%K, grm%P, grm%NPE%AL(jp), &
                             hprogram,hatmfile,hatmfiletype,hpgdfile,hpgdfiletype, jp,ydctl)
    !
    ! Initializing deep GR temp. with that of the outer layer of the structural roof 
    !
    grm%K%XTDEEP(:) = nt%AL(jp)%XT_ROOF(:,1)  
    !
  ENDIF
  !
ENDDO
!
DEALLOCATE(xwsnow_gd,xrsnow_gd,xtsnow_gd,xlwcsnow_gd,xagesnow_gd)
DEALLOCATE(xwsnow_gr,xrsnow_gr,xtsnow_gr,xlwcsnow_gr,xagesnow_gr)
!
!-------------------------------------------------------------------------------------
!
!*      5.     Preparation of canopy air variables
!
top%LCANOPY = lteb_canopy
IF (top%LCANOPY) CALL prep_sbl(tg%NDIM, sb)
!
IF (ydctl%LPART5) THEN
  DEALLOCATE(xzs_ls)
ENDIF
!
!-------------------------------------------------------------------------------------
 CALL clean_prep_output_grid
IF (lhook) CALL dr_hook('PREP_TEB',1,zhook_handle)
!
!
!-------------------------------------------------------------------------------------
!
END SUBROUTINE prep_teb