Changeset 3787 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
- Timestamp:
- Mar 7, 2019 8:43:54 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_layer_fluxes_mod.f90
r3745 r3787 26 26 ! ----------------- 27 27 ! $Id$ 28 ! unused variables removed 29 ! 30 ! 3745 2019-02-15 18:57:56Z suehring 28 31 ! Bugfix, missing calculation of 10cm temperature at vertical building walls, 29 32 ! required for indoor model … … 297 300 INTEGER(iwp) :: k !< loop index z direction 298 301 INTEGER(iwp) :: l !< loop index for surf type 299 INTEGER(iwp) :: li_bnd = 7500 !< Lookup table index of the last time step300 302 301 303 INTEGER(iwp), PARAMETER :: num_steps = 15000 !< number of steps in the lookup table … … 308 310 LOGICAL :: surf_vertical !< Flag indicating vertical surfaces 309 311 310 REAL(wp), DIMENSION(0:num_steps-1) :: rib_tab, & !< Lookup table bulk Richardson number311 ol_tab !< Lookup table values of L312 313 312 REAL(wp) :: e_s, & !< Saturation water vapor pressure 314 313 ol_max = 1.0E6_wp, & !< Maximum Obukhov length 315 rib_max, & !< Maximum Richardson number in lookup table316 rib_min, & !< Minimum Richardson number in lookup table317 314 z_mo !< Height of the constant flux layer where MOST is assumed 318 315 … … 726 723 IMPLICIT NONE 727 724 728 INTEGER(iwp) :: li, & !< Index for loop to create lookup table729 num_steps_n !< Number of non-stretched zeta steps730 731 LOGICAL :: terminate_run_l = .FALSE. !< Flag to terminate run (global)732 733 REAL(wp), PARAMETER :: zeta_stretch = -10.0_wp !< Start of stretching in the free convection limit734 735 REAL(wp), DIMENSION(:), ALLOCATABLE :: zeta_tmp736 737 738 REAL(wp) :: zeta_step, & !< Increment of zeta739 regr = 1.01_wp, & !< Stretching factor of zeta_step in the free convection limit740 regr_old = 1.0E9_wp, & !< Stretching factor of last iteration step741 z0h_min = 0.0_wp, & !< Minimum value of z0h to create table742 z0_min = 0.0_wp !< Minimum value of z0 to create table743 744 725 745 726 CALL location_message( 'initializing surface layer', .FALSE. ) 746 747 727 748 728 ! … … 1015 995 1016 996 INTEGER(iwp) :: iter !< Newton iteration step 1017 INTEGER(iwp) :: li !< look index1018 997 INTEGER(iwp) :: m !< loop variable over all horizontal wall elements 1019 998 … … 2016 1995 IMPLICIT NONE 2017 1996 2018 CHARACTER (LEN = *), INTENT(IN) :: z_char 2019 INTEGER(iwp) :: i, j, k, l,m !< running indices1997 CHARACTER (LEN = *), INTENT(IN) :: z_char !< string identifier to identify z level 1998 INTEGER(iwp) :: i, j, k, m !< running indices 2020 1999 2021 2000
Note: See TracChangeset
for help on using the changeset viewer.