minor changes

UMEP-dev · Aug 20, 2023 · 4f2d46a · 4f2d46a
1 parent efdec82
commit 4f2d46a
Show file tree

Hide file tree

Showing 2 changed files with 25 additions and 5 deletions.
diff --git a/src/suews/src/suews_ctrl_driver.f95 b/src/suews/src/suews_ctrl_driver.f95
@@ -14,7 +14,8 @@ MODULE SUEWS_Driver
                             SITE_PRM, LUMPS_PRM, EHC_PRM, LC_PAVED_PRM, LC_BLDG_PRM, LC_DECTR_PRM, LC_EVETR_PRM, &
                             LC_GRASS_PRM, LC_BSOIL_PRM, LC_WATER_PRM, anthroHEAT_STATE, &
                             OHM_STATE, PHENOLOGY_STATE, SNOW_STATE, SUEWS_FORCING, SUEWS_TIMER, &
-                            HYDRO_STATE, HEAT_STATE
+                            HYDRO_STATE, HEAT_STATE,&
+                            ROUGHNESS_STATE
    USE meteo, ONLY: qsatf, RH2qa, qa2RH
    USE AtmMoistStab_module, ONLY: cal_AtmMoist, cal_Stab, stab_psi_heat, stab_psi_mom
    USE NARP_MODULE, ONLY: NARP_cal_SunPosition, NARP_cal_SunPosition_DTS
@@ -2029,21 +2030,21 @@ SUBROUTINE SUEWS_cal_Main_DTS( &
       ! aggregated heat storage of all wall facets
       REAL(KIND(1D0)), DIMENSION(nlayer) :: QS_wall ! heat storage flux for wall component [W m-2]
       !interface temperature between depth layers
-      REAL(KIND(1D0)), DIMENSION(nlayer, ndepth) :: temp_out_wall !interface temperature between depth layers [degC]
+      ! REAL(KIND(1D0)), DIMENSION(nlayer, ndepth) :: temp_out_wall !interface temperature between depth layers [degC]
 
       ! energy fluxes of individual surfaces
-      REAL(KIND(1D0)), DIMENSION(nlayer) :: QG_wall ! heat flux used in ESTM_ehc as forcing of wall surface [W m-2]
+      ! REAL(KIND(1D0)), DIMENSION(nlayer) :: QG_wall ! heat flux used in ESTM_ehc as forcing of wall surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nlayer) :: QN_wall ! net all-wave radiation of wall surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nlayer) :: qe_wall ! latent heat flux of wall surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nlayer) :: qh_wall ! sensible heat flux of wall surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nlayer) :: qh_resist_wall ! resistance based sensible heat flux of wall surface [W m-2]
 
       ! standard suews surfaces
       !interface temperature between depth layers
-      REAL(KIND(1D0)), DIMENSION(nsurf, ndepth) :: temp_out_surf !interface temperature between depth layers[degC]
+      ! REAL(KIND(1D0)), DIMENSION(nsurf, ndepth) :: temp_out_surf !interface temperature between depth layers[degC]
 
       ! energy fluxes of individual surfaces
-      REAL(KIND(1D0)), DIMENSION(nsurf) :: QG_surf ! heat flux used in ESTM_ehc as forcing of individual surface [W m-2]
+      ! REAL(KIND(1D0)), DIMENSION(nsurf) :: QG_surf ! heat flux used in ESTM_ehc as forcing of individual surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nsurf) :: QN_surf ! net all-wave radiation of individual surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nsurf) :: qs_surf ! aggregated heat storage of of individual surface [W m-2]
       REAL(KIND(1D0)), DIMENSION(nsurf) :: qe0_surf ! latent heat flux from PM of individual surface [W m-2]

diff --git a/src/suews/src/suews_ctrl_type.f95 b/src/suews/src/suews_ctrl_type.f95
@@ -488,6 +488,25 @@ MODULE SUEWS_DEF_DTS
       PROCEDURE :: DEALLOCATE => deallocHeatState_c
    END TYPE HEAT_STATE
 
+   TYPE :: ROUGHNESS_STATE
+   ! this type is used to collect the intermediate results in the SUEWS model
+
+   ! calculated values of FAI
+      REAL(KIND(1D0)) :: FAIBldg_use
+      REAL(KIND(1D0)) :: FAIEveTree_use
+      REAL(KIND(1D0)) :: FAIDecTree_use
+
+      REAL(KIND(1D0)) :: FAI
+      REAL(KIND(1D0)) :: PAI
+      REAL(KIND(1D0)) :: Zh ! effective height of bluff bodies
+      REAL(KIND(1D0)) :: z0m ! aerodynamic roughness length
+      REAL(KIND(1D0)) :: zdm ! zero-plance displacement
+      REAL(KIND(1D0)) :: ZZD ! z-zdm
+
+
+   end type ROUGHNESS_STATE
+
+
 CONTAINS
    SUBROUTINE allocate_hydro_state(self, nlayer)
       IMPLICIT NONE