PageRenderTime 60ms CodeModel.GetById 19ms RepoModel.GetById 0ms app.codeStats 1ms

/wrfv2_fire/phys/module_sf_bep_bem.F

http://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 4101 lines | 2560 code | 730 blank | 811 comment | 59 complexity | 5fab623c5deab0dceb0cde135a58a9af MD5 | raw file
Possible License(s): AGPL-1.0

Large files files are truncated, but you can click here to view the full file

    

  1. MODULE module_sf_bep_bem
    2. 

  2. 

  3. !USE module_model_constants
    4. 

  4.  USE module_sf_urban
    5. 

  5.  USE module_sf_bem
    6. 

  6. 

  7. ! SGClarke 09/11/2008
    8. 

  8. ! Access urban_param.tbl values through calling urban_param_init in module_physics_init
    9. 

  9. ! for CASE (BEPSCHEME) select sf_urban_physics
    10. 

  10. !
    11. 

  11. ! -----------------------------------------------------------------------
    12. 

  12. !  Dimension for the array used in the BEP module
    13. 

  13. ! -----------------------------------------------------------------------
    14. 

  14. 

  15.       integer nurbm           ! Maximum number of urban classes    
    16. 

  16.       parameter (nurbm=3)
    17. 

  17. 

  18.       integer ndm             ! Maximum number of street directions 
    19. 

  19.       parameter (ndm=2)
    20. 

  20. 

  21.       integer nz_um           ! Maximum number of vertical levels in the urban grid
    22. 

  22.       parameter(nz_um=13)
    23. 

  23. 

  24.       integer ng_u            ! Number of grid levels in the ground
    25. 

  25.       parameter (ng_u=10)
    26. 

  26.       integer nwr_u            ! Number of grid levels in the walls or roofs
    27. 

  27.       parameter (nwr_u=10)
    28. 

  28. 

  29.       integer nf_u             !Number of grid levels in the floors (BEM)
    30. 

  30.       parameter (nf_u=10)
    31. 

  31. 

  32.       integer ngb_u            !Number of grid levels in the ground below building (BEM)
    33. 

  33.       parameter (ngb_u=10)
    34. 

  34. 

  35.       real dz_u                ! Urban grid resolution
    36. 

  36.       parameter (dz_u=5.) 
    37. 

  37. 

  38.       integer nbui_max         !maximum number of types of buildings in an urban class
    39. 

  39.       parameter (nbui_max=4)   !must be less or equal than nz_um 
    40. 

  40. 

  41. !---------------------------------------------------------------------------------
    42. 

  42. !Parameters of the windows. The glasses of windows are considered without films  -
    43. 

  43. !Read the paper of J.Karlsson and A.Roos(2000):"modelling the angular behaviour  -
    44. 

  44. !of the total solar energy transmittance of windows".Solar Energy Vol.69,No.4,   -
    45. 

  45. !pp. 321-329, for more details.                                                  -
    46. 

  46. !---------------------------------------------------------------------------------
    47. 

  47.       integer p_num            !number of panes in the windows (1,2 or 3)
    48. 

  48.       parameter (p_num=2)
    49. 

  49.       integer q_num            !category number for the windows (q_num= 4, standard glasses)
    50. 

  50.       parameter(q_num=4)       !Possible values 1,2,...,10
    51. 

  51. 

  52. 

  53. ! The change of ng_u, nwr_u should be done in agreement with the block data
    54. 

  54. !  in the routine "surf_temp" 
    55. 

  55. ! -----------------------------------------------------------------------
    56. 

  56. !  Constant used in the BEP module
    57. 

  57. ! -----------------------------------------------------------------------
    58. 

  58.            
    59. 

  59.       real vk                 ! von Karman constant
    60. 

  60.       real g_u                ! Gravity acceleration
    61. 

  61.       real pi                 !
    62. 

  62.       real r                  ! Perfect gas constant
    63. 

  63.       real cp_u               ! Specific heat at constant pressure
    64. 

  64.       real rcp_u              !
    65. 

  65.       real sigma              !
    66. 

  66.       real p0                 ! Reference pressure at the sea level
    67. 

  67.       real cdrag              ! Drag force constant
    68. 

  68.       real latent             ! Latent heat of vaporization [J/kg] (used in BEM)
    69. 

  69.       parameter(vk=0.40,g_u=9.81,pi=3.141592653,r=287.,cp_u=1004.)        
    70. 

  70.       parameter(rcp_u=r/cp_u,sigma=5.67e-08,p0=1.e+5,cdrag=0.4,latent=2.45e+06)
    71. 

  71. 

  72. ! -----------------------------------------------------------------------     
    73. 

  73. 

  74. 

  75. 

  76. 

  77.    CONTAINS
    78. 

  78.  
    79. 

  79.       subroutine BEP_BEM(FRC_URB2D,UTYPE_URB2D,itimestep,dz8w,dt,u_phy,v_phy,      &
    80. 

  80.                       th_phy,rho,p_phy,swdown,glw,                                 &
    81. 

  81.                       gmt,julday,xlong,xlat,                                       &
    82. 

  82.                       declin_urb,cosz_urb2d,omg_urb2d,                             &
    83. 

  83.                       num_urban_layers,                                            &
    84. 

  84.                       trb_urb4d,tw1_urb4d,tw2_urb4d,tgb_urb4d,                     &
    85. 

  85.                       tlev_urb3d,qlev_urb3d,tw1lev_urb3d,tw2lev_urb3d,             &
    86. 

  86.                       tglev_urb3d,tflev_urb3d,sf_ac_urb3d,lf_ac_urb3d,             &
    87. 

  87.                       cm_ac_urb3d,sfvent_urb3d,lfvent_urb3d,                       &
    88. 

  88.                       sfwin1_urb3d,sfwin2_urb3d,                                   &
    89. 

  89.                       sfw1_urb3d,sfw2_urb3d,sfr_urb3d,sfg_urb3d,                   &
    90. 

  90.                       a_u,a_v,a_t,a_e,b_u,b_v,                                     &
    91. 

  91.                       b_t,b_e,b_q,dlg,dl_u,sf,vl,                                  &
    92. 

  92.                       rl_up,rs_abs,emiss,grdflx_urb,qv_phy,                        &
    93. 

  93.                       ids,ide, jds,jde, kds,kde,                                   &
    94. 

  94.                       ims,ime, jms,jme, kms,kme,                                   &
    95. 

  95.                       its,ite, jts,jte, kts,kte)                    
    96. 

  96. 

  97.       implicit none
    98. 

  98. 

  99. !------------------------------------------------------------------------
    100. 

  100. !     Input
    101. 

  101. !------------------------------------------------------------------------
    102. 

  102.    INTEGER ::                       ids,ide, jds,jde, kds,kde,  &
    103. 

  103.                                     ims,ime, jms,jme, kms,kme,  &
    104. 

  104.                                     its,ite, jts,jte, kts,kte,  &
    105. 

  105.                                     itimestep
    106. 

  106.  
    107. 

  107. 

  108.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   DZ8W
    109. 

  109.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   P_PHY
    110. 

  110.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   RHO
    111. 

  111.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   TH_PHY
    112. 

  112.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   T_PHY
    113. 

  113.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   U_PHY
    114. 

  114.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   V_PHY
    115. 

  115.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   U
    116. 

  116.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme )::   V
    117. 

  117.    REAL, DIMENSION( ims:ime , jms:jme )        ::   GLW
    118. 

  118.    REAL, DIMENSION( ims:ime , jms:jme )        ::   swdown
    119. 

  119.    REAL, DIMENSION( ims:ime, jms:jme )         ::   UST
    120. 

  120.    INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   UTYPE_URB2D
    121. 

  121.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   FRC_URB2D
    122. 

  122.    REAL, INTENT(IN  )   ::                                   GMT 
    123. 

  123.    INTEGER, INTENT(IN  ) ::                               JULDAY
    124. 

  124.    REAL, DIMENSION( ims:ime, jms:jme ),                           &
    125. 

  125.          INTENT(IN   )  ::                           XLAT, XLONG
    126. 

  126.    REAL, INTENT(IN) :: DECLIN_URB
    127. 

  127.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: COSZ_URB2D
    128. 

  128.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: OMG_URB2D
    129. 

  129.    INTEGER, INTENT(IN  ) :: num_urban_layers
    130. 

  130.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: trb_urb4d
    131. 

  131.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw1_urb4d
    132. 

  132.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw2_urb4d
    133. 

  133.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tgb_urb4d
    134. 

  134. !New variables used for BEM
    135. 

  135.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ):: qv_phy
    136. 

  136.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tlev_urb3d
    137. 

  137.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: qlev_urb3d
    138. 

  138.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw1lev_urb3d
    139. 

  139.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw2lev_urb3d
    140. 

  140.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tglev_urb3d
    141. 

  141.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tflev_urb3d
    142. 

  142.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: lf_ac_urb3d
    143. 

  143.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: sf_ac_urb3d
    144. 

  144.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: cm_ac_urb3d
    145. 

  145.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: sfvent_urb3d
    146. 

  146.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: lfvent_urb3d
    147. 

  147.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfwin1_urb3d
    148. 

  148.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfwin2_urb3d
    149. 

  149. !End variables
    150. 

  150.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfw1_urb3d
    151. 

  151.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfw2_urb3d
    152. 

  152.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfr_urb3d
    153. 

  153.    REAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfg_urb3d
    154. 

  154. 

  155.    real z(ims:ime,kms:kme,jms:jme)            ! Vertical coordinates
    156. 

  156.    REAL, INTENT(IN )::   DT      ! Time step
    157. 

  157. 

  158. !      
    159. 

  159. !------------------------------------------------------------------------
    160. 

  160. !     Output
    161. 

  161. !------------------------------------------------------------------------ 
    162. 

  162. !
    163. 

  163. !    Implicit and explicit components of the source and sink terms at each levels,
    164. 

  164. !     the fluxes can be computed as follow: FX = A*X + B   example: t_fluxes = a_t * pt + b_t
    165. 

  165.       real a_u(ims:ime,kms:kme,jms:jme)         ! Implicit component for the momemtum in X-direction (center)
    166. 

  166.       real a_v(ims:ime,kms:kme,jms:jme)         ! Implicit component for the momemtum in Y-direction (center)
    167. 

  167.       real a_t(ims:ime,kms:kme,jms:jme)         ! Implicit component for the temperature
    168. 

  168.       real a_e(ims:ime,kms:kme,jms:jme)         ! Implicit component for the TKE
    169. 

  169.       real b_u(ims:ime,kms:kme,jms:jme)         ! Explicit component for the momemtum in X-direction (center)
    170. 

  170.       real b_v(ims:ime,kms:kme,jms:jme)         ! Explicit component for the momemtum in Y-direction (center)
    171. 

  171.       real b_t(ims:ime,kms:kme,jms:jme)         ! Explicit component for the temperature
    172. 

  172.       real b_e(ims:ime,kms:kme,jms:jme)         ! Explicit component for the TKE
    173. 

  173.       real b_q(ims:ime,kms:kme,jms:jme)         ! Explicit component for the Humidity
    174. 

  174.       real dlg(ims:ime,kms:kme,jms:jme)         ! Height above ground (L_ground in formula (24) of the BLM paper). 
    175. 

  175.       real dl_u(ims:ime,kms:kme,jms:jme)        ! Length scale (lb in formula (22) ofthe BLM paper).
    176. 

  176. ! urban surface and volumes        
    177. 

  177.       real sf(ims:ime,kms:kme,jms:jme)           ! surface of the urban grid cells
    178. 

  178.       real vl(ims:ime,kms:kme,jms:jme)             ! volume of the urban grid cells
    179. 

  179. ! urban fluxes
    180. 

  180.       real rl_up(ims:ime,jms:jme) ! upward long wave radiation
    181. 

  181.       real rs_abs(ims:ime,jms:jme) ! absorbed short wave radiation
    182. 

  182.       real emiss(ims:ime,jms:jme)  ! emissivity averaged for urban surfaces
    183. 

  183.       real grdflx_urb(ims:ime,jms:jme)  ! ground heat flux for urban areas
    184. 

  184. !------------------------------------------------------------------------
    185. 

  185. !     Local
    186. 

  186. !------------------------------------------------------------------------
    187. 

  187. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    188. 

  188. !    Building parameters      
    189. 

  189.       real alag_u(nurbm)                      ! Ground thermal diffusivity [m^2 s^-1]
    190. 

  190.       real alaw_u(nurbm)                      ! Wall thermal diffusivity [m^2 s^-1]
    191. 

  191.       real alar_u(nurbm)                      ! Roof thermal diffusivity [m^2 s^-1]
    192. 

  192.       real csg_u(nurbm)                       ! Specific heat of the ground material [J m^3 K^-1]
    193. 

  193.       real csw_u(nurbm)                       ! Specific heat of the wall material [J m^3 K^-1]
    194. 

  194.       real csr_u(nurbm)                       ! Specific heat of the roof material [J m^3 K^-1]
    195. 

  195.       real twini_u(nurbm)                     ! Initial temperature inside the building's wall [K]
    196. 

  196.       real trini_u(nurbm)                     ! Initial temperature inside the building's roof [K]
    197. 

  197.       real tgini_u(nurbm)                     ! Initial road temperature
    198. 

  198. !
    199. 

  199. ! for twini_u, and trini_u the initial value at the deepest level is kept constant during the simulation
    200. 

  200. !
    201. 

  201. !    Radiation paramters
    202. 

  202.       real albg_u(nurbm)                      ! Albedo of the ground
    203. 

  203.       real albw_u(nurbm)                      ! Albedo of the wall
    204. 

  204.       real albr_u(nurbm)                      ! Albedo of the roof
    205. 

  205.       real albwin_u(nurbm)                    ! Albedo of the windows
    206. 

  206.       real emwind_u(nurbm)                    ! Emissivity of windows
    207. 

  207.       real emg_u(nurbm)                       ! Emissivity of ground
    208. 

  208.       real emw_u(nurbm)                       ! Emissivity of wall
    209. 

  209.       real emr_u(nurbm)                       ! Emissivity of roof
    210. 

  210. 

  211. !   fww,fwg,fgw,fsw,fsg are the view factors used to compute the long wave
    212. 

  212. !   and the short wave radation. 
    213. 

  213.       real fww(nz_um,nz_um,ndm,nurbm)         !  from wall to wall
    214. 

  214.       real fwg(nz_um,ndm,nurbm)               !  from wall to ground
    215. 

  215.       real fgw(nz_um,ndm,nurbm)               !  from ground to wall
    216. 

  216.       real fsw(nz_um,ndm,nurbm)               !  from sky to wall
    217. 

  217.       real fws(nz_um,ndm,nurbm)               !  from sky to wall
    218. 

  218.       real fsg(ndm,nurbm)                     !  from sky to ground
    219. 

  219. 

  220. !    Roughness parameters
    221. 

  221.       real z0g_u(nurbm)         ! The ground's roughness length
    222. 

  222.       real z0r_u(nurbm)         ! The roof's roughness length
    223. 

    224. 

  224. !    Street parameters
    225. 

  225.       integer nd_u(nurbm)       ! Number of street direction for each urban class 
    226. 

  226.       real strd_u(ndm,nurbm)    ! Street length (fix to greater value to the horizontal length of the cells)
    227. 

  227.       real drst_u(ndm,nurbm)    ! Street direction
    228. 

  228.       real ws_u(ndm,nurbm)      ! Street width
    229. 

  229.       real bs_u(ndm,nurbm)      ! Building width
    230. 

  230.       real h_b(nz_um,nurbm)     ! Bulding's heights
    231. 

  231.       real d_b(nz_um,nurbm)     ! Probability that a building has an height h_b
    232. 

  232.       real ss_u(nz_um,nurbm)  ! Probability that a building has an height equal to z
    233. 

  233.       real pb_u(nz_um,nurbm)  ! Probability that a building has an height greater or equal to z
    234. 

  234. 

  235. 

  236. !    Grid parameters
    237. 

  237.       integer nz_u(nurbm)       ! Number of layer in the urban grid
    238. 

  238.       
    239. 

  239.       real z_u(nz_um)         ! Height of the urban grid levels
    240. 

  240. 

  241. !   MT
    242. 

  242.       real cop_u(nurbm)
    243. 

  243.       real pwin_u(nurbm)
    244. 

  244.       real beta_u(nurbm)
    245. 

  245.       integer sw_cond_u(nurbm)
    246. 

  246.       real time_on_u(nurbm)
    247. 

  247.       real time_off_u(nurbm)
    248. 

  248.       real targtemp_u(nurbm)
    249. 

  249.       real gaptemp_u(nurbm)
    250. 

  250.       real targhum_u(nurbm)
    251. 

  251.       real gaphum_u(nurbm)
    252. 

  252.       real perflo_u(nurbm)
    253. 

  253.       real hsesf_u(nurbm)
    254. 

  254.       real hsequip(24)
    255. 

  255. 

  256. ! 1D array used for the input and output of the routine "urban"
    257. 

  257. 

  258.       real z1D(kms:kme)               ! vertical coordinates
    259. 

  259.       real ua1D(kms:kme)                ! wind speed in the x directions
    260. 

  260.       real va1D(kms:kme)                ! wind speed in the y directions
    261. 

  261.       real pt1D(kms:kme)                ! potential temperature
    262. 

  262.       real da1D(kms:kme)                ! air density
    263. 

  263.       real pr1D(kms:kme)                ! air pressure
    264. 

  264.       real pt01D(kms:kme)               ! reference potential temperature
    265. 

  265.       real zr1D                    ! zenith angle
    266. 

  266.       real deltar1D                ! declination of the sun
    267. 

  267.       real ah1D                    ! hour angle (it should come from the radiation routine)
    268. 

  268.       real rs1D                    ! solar radiation
    269. 

  269.       real rld1D                   ! downward flux of the longwave radiation
    270. 

  270. 

  271.       real tw1D(2*ndm,nz_um,nwr_u,nbui_max) ! temperature in each layer of the wall
    272. 

  272.       real tg1D(ndm,ng_u)                   ! temperature in each layer of the ground
    273. 

  273.       real tr1D(ndm,nz_um,nwr_u)   ! temperature in each layer of the roof
    274. 

  274. !
    275. 

  275. !New variable for BEM
    276. 

  276. !
    277. 

  277.       real tlev1D(nz_um,nbui_max)            ! temperature in each floor and in each different type of building
    278. 

  278.       real qlev1D(nz_um,nbui_max)            ! specific humidity in each floor and in each different type of building
    279. 

  279.       real twlev1D(2*ndm,nz_um,nbui_max)     ! temperature in each window in each floor in each different type of building
    280. 

  280.       real tglev1D(ndm,ngb_u,nbui_max)       ! temperature in each layer of the ground below of a type of building
    281. 

  281.       real tflev1D(ndm,nf_u,nz_um-1,nbui_max)! temperature in each layer of the floors in each building
    282. 

  282.       real lflev1D(nz_um,nz_um)           ! latent heat flux due to the air conditioning systems
    283. 

  283.       real sflev1D(nz_um,nz_um)           ! sensible heat flux due to the air conditioning systems
    284. 

  284.       real lfvlev1D(nz_um,nz_um)          ! latent heat flux due to ventilation
    285. 

  285.       real sfvlev1D(nz_um,nz_um)          ! sensible heat flux due to ventilation
    286. 

  286.       real sfwin1D(2*ndm,nz_um,nbui_max)     ! sensible heat flux from windows
    287. 

  287.       real consumlev1D(nz_um,nz_um)       ! consumption due to the air conditioning systems
    288. 

  288.       real qv1D(kms:kme)                  ! specific humidity
    289. 

  289.       real meso_urb                       ! constant to link meso and urban scales [m¯2]
    290. 

  290.       real d_urb(nz_um)    
    291. 

  291.       real sf_ac
    292. 

  292.       integer ibui,nbui
    293. 

  293.       integer nlev(nz_um)
    294. 

  294. !
    295. 

  295. !End new variables
    296. 

  296. !
    297. 

  297. 

  298.       real sfw1D(2*ndm,nz_um,nbui_max)      ! sensible heat flux from walls
    299. 

  299.       real sfg1D(ndm)              ! sensible heat flux from ground (road)
    300. 

  300.       real sfr1D(ndm,nz_um)      ! sensible heat flux from roofs
    301. 

  301.       real sf1D(kms:kme)              ! surface of the urban grid cells
    302. 

  302.       real vl1D(kms:kme)                ! volume of the urban grid cells
    303. 

  303.       real a_u1D(kms:kme)               ! Implicit component of the momentum sources or sinks in the X-direction
    304. 

  304.       real a_v1D(kms:kme)               ! Implicit component of the momentum sources or sinks in the Y-direction
    305. 

  305.       real a_t1D(kms:kme)               ! Implicit component of the heat sources or sinks
    306. 

  306.       real a_e1D(kms:kme)               ! Implicit component of the TKE sources or sinks
    307. 

  307.       real b_u1D(kms:kme)               ! Explicit component of the momentum sources or sinks in the X-direction
    308. 

  308.       real b_v1D(kms:kme)               ! Explicit component of the momentum sources or sinks in the Y-direction
    309. 

  309.       real b_t1D(kms:kme)               ! Explicit component of the heat sources or sinks
    310. 

  310.       real b_ac1D(kms:kme)
    311. 

  311.       real b_e1D(kms:kme)               ! Explicit component of the TKE sources or sinks
    312. 

  312.       real b_q1D(kms:kme)               ! Explicit component of the Humidity sources or sinks
    313. 

  313.       real dlg1D(kms:kme)               ! Height above ground (L_ground in formula (24) of the BLM paper). 
    314. 

  314.       real dl_u1D(kms:kme)              ! Length scale (lb in formula (22) ofthe BLM paper)
    315. 

  315.       real time_bep
    316. 

  316. ! arrays used to collapse indexes
    317. 

  317.       integer ind_zwd(nbui_max,nz_um,nwr_u,ndm)
    318. 

  318.       integer ind_gd(ng_u,ndm)
    319. 

  319.       integer ind_zd(nbui_max,nz_um,ndm)
    320. 

  320.       integer ind_zdf(nz_um,ndm)
    321. 

  321.       integer ind_zrd(nz_um,nwr_u,ndm)
    322. 

  322. !
    323. 

  323.       integer ind_bd(nbui_max,nz_um)
    324. 

  324.       integer ind_wd(nbui_max,nz_um,ndm)
    325. 

  325.       integer ind_gbd(nbui_max,ngb_u,ndm)  
    326. 

  326.       integer ind_fbd(nbui_max,nf_u,nz_um-1,ndm)
    327. 

  327. !
    328. 

  328.       integer ix,iy,iz,iurb,id,iz_u,iw,ig,ir,ix1,iy1,k
    329. 

  329.       integer it, nint
    330. 

  330.       integer iii
    331. 

  331.       real tempo
    332. 

  332.       logical first
    333. 

  333.       character(len=80) :: text
    334. 

  334.       data first/.true./
    335. 

  335.       save first,time_bep 
    336. 

  336. 

  337.       save alag_u,alaw_u,alar_u,csg_u,csw_u,csr_u,                      &
    338. 

  338.           albg_u,albw_u,albr_u,emg_u,emw_u,emr_u,fww,fwg,fgw,fsw,fws,fsg, &
    339. 

  339.           z0g_u,z0r_u, nd_u,strd_u,drst_u,ws_u,bs_u,h_b,d_b,ss_u,pb_u,  &
    340. 

  340.           nz_u,z_u,albwin_u,emwind_u ,                                  &
    341. 

  341.           cop_u, pwin_u, beta_u, sw_cond_u, time_on_u, time_off_u, targtemp_u,  &
    342. 

  342.           gaptemp_u, targhum_u, gaphum_u, perflo_u, hsesf_u, hsequip
    343. 

  343. 

  344. !------------------------------------------------------------------------
    345. 

  345. !    Calculation of the momentum, heat and turbulent kinetic fluxes
    346. 

  346. !     produced by buildings
    347. 

  347. !
    348. 

  348. ! Reference:
    349. 

  349. ! Martilli, A., Clappier, A., Rotach, M.W.:2002, 'AN URBAN SURFACE EXCHANGE
    350. 

  350. ! PARAMETERISATION FOR MESOSCALE MODELS', Boundary-Layer Meteorolgy 104:
    351. 

  351. ! 261-304
    352. 

  352. !
    353. 

  353. ! F. Salamanca and A. Martilli, 2009: 'A new Building Energy Model coupled 
    354. 

  354. ! with an Urban Canopy Parameterization for urban climate simulations_part II. 
    355. 

  355. ! Validation with one dimension off-line simulations'. Theor Appl Climatol
    356. 

  356. ! DOI 10.1007/s00704-009-0143-8 
    357. 

  357. !------------------------------------------------------------------------
    358. 

  358. !prepare the arrays to collapse indexes
    359. 

  359.       if(num_urban_layers.lt.nbui_max*nz_um*ndm*max(nwr_u,ng_u))then
    360. 

  360.         write(*,*)'num_urban_layers too small, please increase to at least ', nbui_max*nz_um*ndm*max(nwr_u,ng_u)
    361. 

  361.         stop
    362. 

  362.       endif
    363. 

  363. !
    364. 

  364. !New conditions for BEM
    365. 

  365. !
    366. 

  366.       if(num_urban_layers.lt.nbui_max*nz_um)then !limit for indoor temperature and indoor humidity
    367. 

  367.         write(*,*)'num_urban_layers too small, please increase to at least ', nbui_max*nz_um
    368. 

  368.         stop
    369. 

  369.       endif
    370. 

  370. 

  371.       if(num_urban_layers.lt.nbui_max*nz_um*ndm)then !limit for window temperature
    372. 

  372.         write(*,*)'num_urban_layers too small, please increase to at least ', nbui_max*nz_um*ndm
    373. 

  373.         stop
    374. 

  374.       endif
    375. 

  375. 

  376.       if(num_urban_layers.lt.nbui_max*ndm*ngb_u)then !limit for ground temperature below a building
    377. 

  377.         write(*,*)'num_urban_layers too small, please increase to at least ', nbui_max*ndm*ngb_u
    378. 

  378.         stop
    379. 

  379.       endif
    380. 

  380. 

  381.       if(num_urban_layers.lt.(nz_um-1)*nbui_max*ndm*nf_u)then !limit for floor temperature 
    382. 

  382.         write(*,*)'num_urban_layers too small, please increase to at least ', nbui_max*ndm*nf_u*(nz_um-1),num_urban_layers
    383. 

  383.         stop
    384. 

  384.       endif
    385. 

  385. 

  386.       if (ndm.ne.2)then
    387. 

  387.          write(*,*) 'number of directions is not correct',ndm
    388. 

  388.          stop
    389. 

  389.       endif
    390. 

  390. 

  391. !
    392. 

  392. !End of new conditions
    393. 

  393. !
    394. 

  394. !
    395. 

  395. !Initialize collapse indexes
    396. 

  396. !
    397. 

  397.       ind_zwd=0       
    398. 

  398.       ind_gd=0
    399. 

  399.       ind_zd=0
    400. 

  400.       ind_zdf=0
    401. 

  401.       ind_zrd=0
    402. 

  402.       ind_bd=0
    403. 

  403.       ind_wd=0
    404. 

  404.       ind_gbd=0
    405. 

  405.       ind_fbd=0
    406. 

  406. !
    407. 

  407. !End initialization indexes
    408. 

  408. !
    409. 

  409. 

  410.       iii=0
    411. 

  411.       do ibui=1,nbui_max
    412. 

  412.       do iz_u=1,nz_um
    413. 

  413.       do iw=1,nwr_u
    414. 

  414.       do id=1,ndm
    415. 

  415.        iii=iii+1
    416. 

  416.        ind_zwd(ibui,iz_u,iw,id)=iii
    417. 

  417.       enddo
    418. 

  418.       enddo
    419. 

  419.       enddo
    420. 

  420.       enddo
    421. 

  421.       
    422. 

  422.       iii=0
    423. 

  423.       do ig=1,ng_u
    424. 

  424.       do id=1,ndm
    425. 

  425.        iii=iii+1
    426. 

  426.        ind_gd(ig,id)=iii
    427. 

  427.       enddo
    428. 

  428.       enddo
    429. 

  429.       
    430. 

  430.       iii=0
    431. 

  431.       do ibui=1,nbui_max
    432. 

  432.       do iz_u=1,nz_um
    433. 

  433.       do id=1,ndm
    434. 

  434.        iii=iii+1
    435. 

  435.        ind_zd(ibui,iz_u,id)=iii
    436. 

  436.       enddo
    437. 

  437.       enddo
    438. 

  438.       enddo
    439. 

  439.       
    440. 

  440.       iii=0
    441. 

  441.       do iz_u=1,nz_um
    442. 

  442.       do iw=1,nwr_u
    443. 

  443.       do id=1,ndm
    444. 

  444.        iii=iii+1
    445. 

  445.        ind_zrd(iz_u,iw,id)=iii
    446. 

  446.       enddo
    447. 

  447.       enddo
    448. 

  448.       enddo
    449. 

  449.       
    450. 

  450. !
    451. 

  451. !New indexes for BEM
    452. 

  452. !
    453. 

  453.       iii=0
    454. 

  454.       do iz_u=1,nz_um
    455. 

  455.       do id=1,ndm
    456. 

  456.          iii=iii+1
    457. 

  457.          ind_zdf(iz_u,id)=iii
    458. 

  458.       enddo ! id
    459. 

  459.       enddo ! iz_u
    460. 

  460.       
    461. 

  461.       iii=0
    462. 

  462.       do ibui=1,nbui_max  !Type of building
    463. 

  463.       do iz_u=1,nz_um     !vertical levels
    464. 

  464.          iii=iii+1
    465. 

  465.          ind_bd(ibui,iz_u)=iii
    466. 

  466.       enddo !iz_u
    467. 

  467.       enddo !ibui
    468. 

  468.       
    469. 

  469.       
    470. 

  470.       iii=0
    471. 

  471.       do ibui=1,nbui_max !type of building
    472. 

  472.       do iz_u=1,nz_um !vertical levels
    473. 

  473.       do id=1,ndm !direction
    474. 

  474.          iii=iii+1
    475. 

  475.          ind_wd(ibui,iz_u,id)=iii
    476. 

  476.       enddo !id
    477. 

  477.       enddo !iz_u
    478. 

  478.       enddo !ibui
    479. 

  479.       
    480. 

  480.       iii=0
    481. 

  481.       do ibui=1,nbui_max!type of building
    482. 

  482.       do iw=1,ngb_u !layers in the wall (ground below a building)
    483. 

  483.       do id=1,ndm !direction
    484. 

  484.          iii=iii+1
    485. 

  485.          ind_gbd(ibui,iw,id)=iii  
    486. 

  486.       enddo !id
    487. 

  487.       enddo !iw 
    488. 

  488.       enddo !ibui    
    489. 

  489.       
    490. 

  490.       iii=0
    491. 

  491.       do ibui=1,nbui_max !type of building
    492. 

  492.       do iw=1,nf_u !layers in the wall (floor)
    493. 

  493.       do iz_u=1,nz_um-1 !vertical levels
    494. 

  494.       do id=1,ndm  !direction
    495. 

  495.          iii=iii+1
    496. 

  496.          ind_fbd(ibui,iw,iz_u,id)=iii
    497. 

  497.       enddo !id
    498. 

  498.       enddo !iz_u
    499. 

  499.       enddo !iw
    500. 

  500.       enddo !ibui
    501. 

  501.       
    502. 

  502. !
    503. 

  503. !End of new indexes
    504. 

  504. !      
    505. 

  505. 

  506.       do ix=its,ite
    507. 

  507.       do iy=jts,jte
    508. 

  508.        z(ix,kts,iy)=0.
    509. 

  509.        do iz=kts+1,kte+1
    510. 

  510.         z(ix,iz,iy)=z(ix,iz-1,iy)+dz8w(ix,iz-1,iy)
    511. 

  511.        enddo
    512. 

  512.       enddo
    513. 

  513.       enddo
    514. 

  514. 

  515.       if (first) then                           ! True only on first call
    516. 

  516.          call init_para(alag_u,alaw_u,alar_u,csg_u,csw_u,csr_u,&
    517. 

  517.                 twini_u,trini_u,tgini_u,albg_u,albw_u,albr_u,albwin_u,emg_u,emw_u,&
    518. 

  518.                 emr_u,emwind_u,z0g_u,z0r_u,nd_u,strd_u,drst_u,ws_u,bs_u,h_b,d_b,&
    519. 

  519.                 cop_u, pwin_u, beta_u, sw_cond_u, time_on_u, time_off_u, &
    520. 

  520.                 targtemp_u, gaptemp_u, targhum_u, gaphum_u, perflo_u, hsesf_u, hsequip)
    521. 

  521. 

  522. !Initialisation of the urban parameters and calculation of the view factor
    523. 

  523.        call icBEP(fww,fwg,fgw,fsw,fws,fsg,                         & 
    524. 

  524.                  z0g_u,z0r_u,                                      & 
    525. 

  525.                  nd_u,strd_u,ws_u,bs_u,h_b,d_b,ss_u,pb_u,          & 
    526. 

  526.                  nz_u,z_u)                                  
    527. 

  527. 

  528.        first=.false.
    529. 

  529. 

  530.       endif ! first
    531. 

  531.       
    532. 

  532.       do ix=its,ite
    533. 

  533.       do iy=jts,jte
    534. 

  534.         if (FRC_URB2D(ix,iy).gt.0.) then    ! Calling BEP only for existing urban classes.
    535. 

  535. 	
    536. 

  536.          iurb=UTYPE_URB2D(ix,iy)
    537. 

  537.          ibui=0
    538. 

  538.          nlev=0
    539. 

  539.          nbui=0
    540. 

  540.          d_urb=0.
    541. 

  541.          do iz=1,nz_um		   
    542. 

  542.          if(ss_u(iz,iurb).gt.0) then		
    543. 

  543.            ibui=ibui+1		                
    544. 

  544.            nlev(ibui)=iz-1
    545. 

  545.            d_urb(ibui)=ss_u(iz,iurb)
    546. 

  546.            nbui=ibui
    547. 

  547. 	 endif	  
    548. 

  548.          end do  !iz
    549. 

  549.          if (nbui.gt.nbui_max) then
    550. 

  550.             write (*,*) 'nbui_max must be increased to',nbui
    551. 

  551.             stop
    552. 

  552.          endif
    553. 

  553. 

  554.        do iz= kts,kte
    555. 

  555.           ua1D(iz)=u_phy(ix,iz,iy)
    556. 

  556.           va1D(iz)=v_phy(ix,iz,iy)
    557. 

  557. 	  pt1D(iz)=th_phy(ix,iz,iy)
    558. 

  558. 	  da1D(iz)=rho(ix,iz,iy)
    559. 

  559. 	  pr1D(iz)=p_phy(ix,iz,iy)
    560. 

  560. !	  pt01D(iz)=th_phy(ix,iz,iy)
    561. 

  561. 	  pt01D(iz)=300.
    562. 

  562. 	  z1D(iz)=z(ix,iz,iy)
    563. 

  563.           qv1D(iz)=qv_phy(ix,iz,iy)
    564. 

  564.           a_u1D(iz)=0.
    565. 

  565.           a_v1D(iz)=0.
    566. 

  566.           a_t1D(iz)=0.
    567. 

  567.           a_e1D(iz)=0.
    568. 

  568.           b_u1D(iz)=0.
    569. 

  569.           b_v1D(iz)=0.
    570. 

  570.           b_t1D(iz)=0.
    571. 

  571.           b_ac1D(iz)=0.
    572. 

  572.           b_e1D(iz)=0.           
    573. 

  573.          enddo
    574. 

  574. 	 z1D(kte+1)=z(ix,kte+1,iy)
    575. 

  575. 

  576.          do id=1,ndm
    577. 

  577.          do iz_u=1,nz_um
    578. 

  578.          do iw=1,nwr_u
    579. 

  579.          do ibui=1,nbui_max
    580. 

  580.           tw1D(2*id-1,iz_u,iw,ibui)=tw1_urb4d(ix,ind_zwd(ibui,iz_u,iw,id),iy)
    581. 

  581.           tw1D(2*id,iz_u,iw,ibui)=tw2_urb4d(ix,ind_zwd(ibui,iz_u,iw,id),iy)
    582. 

  582.          enddo
    583. 

  583.          enddo
    584. 

  584.          enddo
    585. 

  585.          enddo
    586. 

  586. 	
    587. 

  587.          do id=1,ndm
    588. 

  588.           do ig=1,ng_u
    589. 

  589.             tg1D(id,ig)=tgb_urb4d(ix,ind_gd(ig,id),iy)
    590. 

  590.           enddo
    591. 

  591.           do iz_u=1,nz_um
    592. 

  592.           do ir=1,nwr_u
    593. 

  593.             tr1D(id,iz_u,ir)=trb_urb4d(ix,ind_zrd(iz_u,ir,id),iy)
    594. 

  594.           enddo
    595. 

  595.           enddo
    596. 

  596.          enddo
    597. 

  597. !
    598. 

  598. !Initialize variables for BEM
    599. 

  599. !
    600. 

  600.          tlev1D=0.  !Indoor temperature
    601. 

  601.          qlev1D=0.  !Indoor humidity
    602. 

  602. 

  603.          twlev1D=0. !Window temperature
    604. 

  604.          tglev1D=0. !Ground temperature
    605. 

  605.          tflev1D=0. !Floor temperature
    606. 

  606. 

  607.          sflev1D=0.    !Sensible heat flux from the a.c.
    608. 

  608.          lflev1D=0.    !latent heat flux from the a.c.
    609. 

  609.          consumlev1D=0.!consumption of the a.c.
    610. 

  610.          sfvlev1D=0.   !Sensible heat flux from natural ventilation
    611. 

  611.          lfvlev1D=0.   !Latent heat flux from natural ventilation
    612. 

  612.          sfwin1D=0.    !Sensible heat flux from windows
    613. 

  613.          sfw1D=0.      !Sensible heat flux from walls         
    614. 

  614. 

  615.          do iz_u=1,nz_um    !vertical levels
    616. 

  616.          do ibui=1,nbui_max !Type of building
    617. 

  617.             tlev1D(iz_u,ibui)= tlev_urb3d(ix,ind_bd(ibui,iz_u),iy)  
    618. 

  618.             qlev1D(iz_u,ibui)= qlev_urb3d(ix,ind_bd(ibui,iz_u),iy)  
    619. 

  619.          enddo !ibui
    620. 

  620.          enddo !iz_u
    621. 

  621. 

  622.          do id=1,ndm  !direction
    623. 

  623.             do iz_u=1,nz_um !vertical levels
    624. 

  624.                do ibui=1,nbui_max !type of building
    625. 

  625.                   twlev1D(2*id-1,iz_u,ibui)=tw1lev_urb3d(ix,ind_wd(ibui,iz_u,id),iy)
    626. 

  626.                   twlev1D(2*id,iz_u,ibui)=tw2lev_urb3d(ix,ind_wd(ibui,iz_u,id),iy)
    627. 

  627.                   sfwin1D(2*id-1,iz_u,ibui)=sfwin1_urb3d(ix,ind_wd(ibui,iz_u,id),iy)
    628. 

  628.                   sfwin1D(2*id,iz_u,ibui)=sfwin2_urb3d(ix,ind_wd(ibui,iz_u,id),iy)
    629. 

  629.                enddo !ibui  
    630. 

  630.             enddo !iz_u
    631. 

  631.          enddo !id
    632. 

  632. 

  633.          do id=1,ndm !direction
    634. 

  634.             do iw=1,ngb_u !layer in the wall
    635. 

  635.                do ibui=1,nbui_max !type of building
    636. 

  636.                   tglev1D(id,iw,ibui)=tglev_urb3d(ix,ind_gbd(ibui,iw,id),iy)
    637. 

  637.                enddo !ibui
    638. 

  638.             enddo !iw
    639. 

  639.          enddo !id
    640. 

  640.        
    641. 

  641.          do id=1,ndm !direction
    642. 

  642.             do iw=1,nf_u !layer in the walls
    643. 

  643.                do iz_u=1,nz_um-1 !verticals levels
    644. 

  644.                   do ibui=1,nbui_max !type of building
    645. 

  645.                      tflev1D(id,iw,iz_u,ibui)=tflev_urb3d(ix,ind_fbd(ibui,iw,iz_u,id),iy)
    646. 

  646.                   enddo !ibui
    647. 

  647.                enddo ! iz_u
    648. 

  648.              enddo !iw
    649. 

  649.          enddo !id
    650. 

  650. 

  651. !
    652. 

  652. !End initialization for BEM
    653. 

  653. !         
    654. 

  654. 

  655.          do id=1,ndm
    656. 

  656. 	 do iz=1,nz_um
    657. 

  657.          do ibui=1,nbui_max !type of building
    658. 

  658. 	  sfw1D(2*id-1,iz,ibui)=sfw1_urb3d(ix,ind_zd(ibui,iz,id),iy)
    659. 

  659. 	  sfw1D(2*id,iz,ibui)=sfw2_urb3d(ix,ind_zd(ibui,iz,id),iy)
    660. 

  660.          enddo
    661. 

  661. 	 enddo
    662. 

  662. 	 enddo
    663. 

  663. 	 
    664. 

  664. 	 do id=1,ndm
    665. 

  665. 	  sfg1D(id)=sfg_urb3d(ix,id,iy)
    666. 

  666. 	 enddo
    667. 

  667. 	 
    668. 

  668. 	 do id=1,ndm
    669. 

  669. 	 do iz=1,nz_um
    670. 

  670. 	  sfr1D(id,iz)=sfr_urb3d(ix,ind_zdf(iz,id),iy)
    671. 

  671. 	 enddo
    672. 

  672. 	 enddo
    673. 

  673.          
    674. 

  674.          rs1D=swdown(ix,iy)
    675. 

  675.          rld1D=glw(ix,iy)
    676. 

  676. 

  677.          zr1D=acos(COSZ_URB2D(ix,iy))
    678. 

  678.          deltar1D=DECLIN_URB
    679. 

  679.          ah1D=OMG_URB2D(ix,iy)       
    680. 

  680. 

  681.          call BEP1D(iurb,kms,kme,kts,kte,z1D,dt,ua1D,va1D,pt1D,da1D,pr1D,pt01D,  &
    682. 

  682.                    zr1D,deltar1D,ah1D,rs1D,rld1D,                   & 
    683. 

  683.                    alag_u,alaw_u,alar_u,csg_u,csw_u,csr_u,          & 
    684. 

  684.                    albg_u,albw_u,albr_u,albwin_u,emg_u,emw_u,emr_u, & 
    685. 

  685.                    emwind_u,fww,fwg,fgw,fsw,fws,fsg,                & 
    686. 

  686.                    z0g_u,z0r_u,                                     & 
    687. 

  687.                    nd_u,strd_u,drst_u,ws_u,bs_u,h_b,d_b,ss_u,pb_u,  & 
    688. 

  688.                    nz_u,z_u,                                        & 
    689. 

  689.                    cop_u,pwin_u,beta_u,sw_cond_u,time_on_u,         &
    690. 

  690.                    time_off_u,targtemp_u,gaptemp_u,targhum_u,       &
    691. 

  691.                    gaphum_u, perflo_u, hsesf_u, hsequip,            &
    692. 

  692.                    tw1D,tg1D,tr1D,sfw1D,sfg1D,sfr1D,                & 
    693. 

  693.                    a_u1D,a_v1D,a_t1D,a_e1D,                         & 
    694. 

  694.                    b_u1D,b_v1D,b_t1D,b_ac1D,b_e1D,b_q1D,            & 
    695. 

  695.                    dlg1D,dl_u1D,sf1D,vl1D,rl_up(ix,iy),             &
    696. 

  696.                    rs_abs(ix,iy),emiss(ix,iy),grdflx_urb(ix,iy),    &
    697. 

  697.                    qv1D,tlev1D,qlev1D,sflev1D,lflev1D,consumlev1D,  &
    698. 

  698.                    sfvlev1D,lfvlev1D,twlev1D,tglev1D,tflev1D,sfwin1D,&
    699. 

  699.                    ix,iy)                            
    700. 

  700. 

  701.          do id=1,ndm ! direction 
    702. 

  702. 	    do iz=1,nz_um   !vertical levels
    703. 

  703.                do ibui=1,nbui_max !type of building
    704. 

  704. 	          sfw1_urb3d(ix,ind_zd(ibui,iz,id),iy)=sfw1D(2*id-1,iz,ibui) 
    705. 

  705. 	          sfw2_urb3d(ix,ind_zd(ibui,iz,id),iy)=sfw1D(2*id,iz,ibui) 
    706. 

  706. 	       enddo
    707. 

  707. 	    enddo
    708. 

  708.          enddo
    709. 

  709. 	 do id=1,ndm
    710. 

  710. 	  sfg_urb3d(ix,id,iy)=sfg1D(id) 
    711. 

  711. 	 enddo
    712. 

  712.          
    713. 

  713. 	 do id=1,ndm
    714. 

  714. 	 do iz=1,nz_um
    715. 

  715. 	  sfr_urb3d(ix,ind_zdf(iz,id),iy)=sfr1D(id,iz)
    716. 

  716. 	 enddo
    717. 

  717. 	 enddo
    718. 

  718.          
    719. 

  719.          do id=1,ndm
    720. 

  720.          do iz_u=1,nz_um
    721. 

  721.          do iw=1,nwr_u
    722. 

  722.          do ibui=1,nbui_max
    723. 

  723.           tw1_urb4d(ix,ind_zwd(ibui,iz_u,iw,id),iy)=tw1D(2*id-1,iz_u,iw,ibui)
    724. 

  724.           tw2_urb4d(ix,ind_zwd(ibui,iz_u,iw,id),iy)=tw1D(2*id,iz_u,iw,ibui)
    725. 

  725.          enddo
    726. 

  726.          enddo
    727. 

  727.          enddo
    728. 

  728.          enddo
    729. 

  729.          
    730. 

  730.          do id=1,ndm
    731. 

  731.             do ig=1,ng_u
    732. 

  732.                tgb_urb4d(ix,ind_gd(ig,id),iy)=tg1D(id,ig)
    733. 

  733.             enddo
    734. 

  734.             do iz_u=1,nz_um
    735. 

  735.                do ir=1,nwr_u
    736. 

  736.                   trb_urb4d(ix,ind_zrd(iz_u,ir,id),iy)=tr1D(id,iz_u,ir)
    737. 

  737.                enddo
    738. 

  738.             enddo
    739. 

  739.          enddo
    740. 

  740. !
    741. 

  741. !Outputs of BEM
    742. 

  742. !
    743. 

  743.          do ibui=1,nbui_max !type of building
    744. 

  744.          do iz_u=1,nz_um !vertical levels
    745. 

  745.             tlev_urb3d(ix,ind_bd(ibui,iz_u),iy)=tlev1D(iz_u,ibui)  
    746. 

  746.             qlev_urb3d(ix,ind_bd(ibui,iz_u),iy)=qlev1D(iz_u,ibui)  
    747. 

  747.          enddo !iz_u
    748. 

  748.          enddo !ibui
    749. 

  749.          do id=1,ndm !direction
    750. 

  750.          do iz_u=1,nz_um !vertical levels
    751. 

  751.          do ibui=1,nbui_max !type of building
    752. 

  752.                tw1lev_urb3d(ix,ind_wd(ibui,iz_u,id),iy)=twlev1D(2*id-1,iz_u,ibui)
    753. 

  753.                tw2lev_urb3d(ix,ind_wd(ibui,iz_u,id),iy)=twlev1D(2*id,iz_u,ibui)
    754. 

  754.                sfwin1_urb3d(ix,ind_wd(ibui,iz_u,id),iy)=sfwin1D(2*id-1,iz_u,ibui)
    755. 

  755.                sfwin2_urb3d(ix,ind_wd(ibui,iz_u,id),iy)=sfwin1D(2*id,iz_u,ibui)
    756. 

  756.             enddo !ibui  
    757. 

  757.          enddo !iz_u
    758. 

  758.          enddo !id
    759. 

  759.         
    760. 

  760.          do id=1,ndm !direction
    761. 

  761.             do iw=1,ngb_u !layers in the walls
    762. 

  762.                do ibui=1,nbui_max  !type of building
    763. 

  763.                   tglev_urb3d(ix,ind_gbd(ibui,iw,id),iy)=tglev1D(id,iw,ibui)
    764. 

  764.                enddo !ibui
    765. 

  765.             enddo !iw
    766. 

  766.          enddo !id
    767. 

  767. 

  768.          do id=1,ndm   !direction 
    769. 

  769.             do iw=1,nf_u !layer in the walls
    770. 

  770.                do iz_u=1,nz_um-1 !verticals levels
    771. 

  771.                   do ibui=1,nbui_max  !type of building
    772. 

  772.                      tflev_urb3d(ix,ind_fbd(ibui,iw,iz_u,id),iy)=tflev1D(id,iw,iz_u,ibui)
    773. 

  773.                   enddo !ibui
    774. 

  774.                enddo !iz_u
    775. 

  775.              enddo !iw
    776. 

  776.          enddo !id
    777. 

  777.          sf_ac_urb3d(ix,iy)=0.
    778. 

  778.          lf_ac_urb3d(ix,iy)=0.
    779. 

  779.          cm_ac_urb3d(ix,iy)=0.
    780. 

  780.          sfvent_urb3d(ix,iy)=0.
    781. 

  781.          lfvent_urb3d(ix,iy)=0.
    782. 

  782. 

  783.          meso_urb=(1./4.)*FRC_URB2D(ix,iy)/((bs_u(1,iurb)+ws_u(1,iurb))*bs_u(2,iurb))+ &
    784. 

  784.                   (1./4.)*FRC_URB2D(ix,iy)/((bs_u(2,iurb)+ws_u(2,iurb))*bs_u(1,iurb))
    785. 

  785. 

  786.                   
    787. 

  787.          ibui=0
    788. 

  788.          nlev=0
    789. 

  789.          nbui=0
    790. 

  790.          d_urb=0.
    791. 

  791.          do iz=1,nz_um		   
    792. 

  792.          if(ss_u(iz,iurb).gt.0) then		
    793. 

  793.            ibui=ibui+1		                
    794. 

  794.            nlev(ibui)=iz-1
    795. 

  795.            d_urb(ibui)=ss_u(iz,iurb)
    796. 

  796.            nbui=ibui
    797. 

  797. 	 endif	  
    798. 

  798.          end do  !iz
    799. 

  799. 

  800.          do ibui=1,nbui       !type of building   
    801. 

  801.          do iz_u=1,nlev(ibui) !vertical levels                    
    802. 

  802.                sf_ac_urb3d(ix,iy)=sf_ac_urb3d(ix,iy)+meso_urb*d_urb(ibui)*sflev1D(iz_u,ibui)
    803. 

  803.                lf_ac_urb3d(ix,iy)=lf_ac_urb3d(ix,iy)+meso_urb*d_urb(ibui)*lflev1D(iz_u,ibui)
    804. 

  804.                cm_ac_urb3d(ix,iy)=cm_ac_urb3d(ix,iy)+meso_urb*d_urb(ibui)*consumlev1D(iz_u,ibui)
    805. 

  805.                sfvent_urb3d(ix,iy)=sfvent_urb3d(ix,iy)+meso_urb*d_urb(ibui)*sfvlev1D(iz_u,ibui)
    806. 

  806.                lfvent_urb3d(ix,iy)=lfvent_urb3d(ix,iy)+meso_urb*d_urb(ibui)*lfvlev1D(iz_u,ibui)            
    807. 

  807.          enddo !iz_u
    808. 

  808.          enddo !ibui
    809. 

  809. !
    810. 

  810. !Add the latent heat exchanged throughout the ventilation in the lf_ac_urb3d output variable. 
    811. 

  811. !it is only a print variable
    812. 

  812. !
    813. 

  813. !        lf_ac_urb3d(ix,iy)=lf_ac_urb3d(ix,iy)+lfvent_urb3d(ix,iy)
    814. 

  814. !
    815. 

  815. 

  816.          lf_ac_urb3d(ix,iy)=lf_ac_urb3d(ix,iy)-lfvent_urb3d(ix,iy)        
    817. 

  817. 

  818. !
    819. 

  819. !End outputs of bem
    820. 

  820. !
    821. 

  821. 

  822.         sf_ac=0.
    823. 

  823.         do iz= kts,kte
    824. 

  824.           sf(ix,iz,iy)=sf1D(iz)
    825. 

  825.           vl(ix,iz,iy)=vl1D(iz)
    826. 

  826.           a_u(ix,iz,iy)=a_u1D(iz)
    827. 

  827.           a_v(ix,iz,iy)=a_v1D(iz)
    828. 

  828.           a_t(ix,iz,iy)=a_t1D(iz)
    829. 

  829.           a_e(ix,iz,iy)=a_e1D(iz)
    830. 

  830.           b_u(ix,iz,iy)=b_u1D(iz)
    831. 

  831.           b_v(ix,iz,iy)=b_v1D(iz)
    832. 

  832.           b_t(ix,iz,iy)=b_t1D(iz)
    833. 

  833.           sf_ac=sf_ac+b_ac1D(iz)*da1D(iz)*cp_u*dz8w(ix,iz,iy)*vl1D(iz)*FRC_URB2D(ix,iy)
    834. 

  834.           b_e(ix,iz,iy)=b_e1D(iz)
    835. 

  835.           b_q(ix,iz,iy)=b_q1D(iz)
    836. 

  836.           dlg(ix,iz,iy)=dlg1D(iz)
    837. 

  837.           dl_u(ix,iz,iy)=dl_u1D(iz)
    838. 

  838.          enddo
    839. 

  839.          sf(ix,kte+1,iy)=sf1D(kte+1)
    840. 

  840. 

  841.          endif ! FRC_URB2D
    842. 

  842.    
    843. 

  843.       enddo  ! iy
    844. 

  844.       enddo  ! ix
    845. 

  845. 

  846.         time_bep=time_bep+dt
    847. 

  847. 

  848.       return
    849. 

  849.       end subroutine BEP_BEM
    850. 

  850.             
    851. 

  851. ! ===6=8===============================================================72
    852. 

  852. 

  853.       subroutine BEP1D(iurb,kms,kme,kts,kte,z,dt,ua,va,pt,da,pr,pt0,   &  
    854. 

  854.                       zr,deltar,ah,rs,rld,                             & 
    855. 

  855.                       alag_u,alaw_u,alar_u,csg_u,csw_u,csr_u,          & 
    856. 

  856.                       albg_u,albw_u,albr_u,albwin_u,emg_u,emw_u,emr_u, & 
    857. 

  857.                       emwind_u,fww,fwg,fgw,fsw,fws,fsg,                & 
    858. 

  858.                       z0g_u,z0r_u,                                     & 
    859. 

  859.                       nd_u,strd_u,drst_u,ws_u,bs_u,h_b,d_b,ss_u,pb_u,  & 
    860. 

  860.                       nz_u,z_u,                                        & 
    861. 

  861.                       cop_u,pwin_u,beta_u,sw_cond_u,time_on_u,         &
    862. 

  862.                       time_off_u,targtemp_u,gaptemp_u,targhum_u,       &
    863. 

  863.                       gaphum_u, perflo_u, hsesf_u, hsequip,            &
    864. 

  864.                       tw,tg,tr,sfw,sfg,sfr,                            & 
    865. 

  865.                       a_u,a_v,a_t,a_e,                                 &
    866. 

  866.                       b_u,b_v,b_t,b_ac,b_e,b_q,                        & 
    867. 

  867.                       dlg,dl_u,sf,vl,rl_up,rs_abs,emiss,grdflx_urb,    &
    868. 

  868.                       qv,tlev,qlev,sflev,lflev,consumlev,              &
    869. 

  869.                       sfvlev,lfvlev,twlev,tglev,tflev,sfwin,ix,iy)                             
    870. 

  870. 

  871. ! ----------------------------------------------------------------------
    872. 

  872. ! This routine computes the effects of buildings on momentum, heat and
    873. 

  873. !  TKE (turbulent kinetic energy) sources or sinks and on the mixing length.
    874. 

  874. ! It provides momentum, heat and TKE sources or sinks at different levels of a
    875. 

  875. !  mesoscale grid defined by the altitude of its cell interfaces "z" and
    876. 

  876. !  its number of levels "nz".
    877. 

  877. ! The meteorological input parameters (wind, temperature, solar radiation)
    878. 

  878. !  are specified on the "mesoscale grid".
    879. 

  879. ! The inputs concerning the building and street charateristics are defined
    880. 

  880. !  on a "urban grid". The "urban grid" is defined with its number of levels
    881. 

  881. !  "nz_u" and its space step "dz_u".
    882. 

  882. ! The input parameters are interpolated on the "urban grid". The sources or sinks
    883. 

  883. !  are calculated on the "urban grid". Finally the sources or sinks are 
    884. 

  884. !  interpolated on the "mesoscale grid".
    885. 

  885.  
    886. 

  886. 

  887. !  Mesoscale grid            Urban grid             Mesoscale grid
    888. 

  888. !  
    889. 

  889. ! z(4)    ---                                               ---
    890. 

  890. !          |                                                 |
    891. 

  891. !          |                                                 |
    892. 

  892. !          |   Interpolation                  Interpolation  |
    893. 

  893. !          |            Sources or sinks calculation         |
    894. 

  894. ! z(3)    ---                                               ---
    895. 

  895. !          | ua               ua_u  ---  uv_a         a_u    |
    896. 

  896. !          | va               va_u   |   uv_b         b_u    |
    897. 

  897. !          | pt               pt_u  ---  uh_b         a_v    |
    898. 

  898. ! z(2)    ---                        |    etc...      etc...---
    899. 

  899. !          |                 z_u(1) ---                      |
    900. 

  900. !          |                         |                       |
    901. 

  901. ! z(1) ------------------------------------------------------------
    902. 

  902. 

  903. !     
    904. 

  904. ! Reference:
    905. 

  905. ! Martilli, A., Clappier, A., Rotach, M.W.:2002, 'AN URBAN SURFACE EXCHANGE
    906. 

  906. ! PARAMETERISATION FOR MESOSCALE MODELS', Boundary-Layer Meteorolgy 104:
    907. 

  907. ! 261-304
    908. 

  908.  
    909. 

  909. ! ----------------------------------------------------------------------
    910. 

  910. 

  911.       implicit none
    912. 

  912. 

  913. ! ----------------------------------------------------------------------
    914. 

  914. ! INPUT:
    915. 

  915. ! ----------------------------------------------------------------------
    916. 

  916. 

  917. ! Data relative to the "mesoscale grid"
    918. 

  918. 

  919.       integer kms,kme,kts,kte
    920. 

  920.       real z(kms:kme)               ! Altitude above the ground of the cell interfaces.
    921. 

  921.       real ua(kms:kme)                ! Wind speed in the x direction
    922. 

  922.       real va(kms:kme)                ! Wind speed in the y direction
    923. 

  923.       real pt(kms:kme)                ! Potential temperature
    924. 

  924.       real da(kms:kme)                ! Air density
    925. 

  925.       real pr(kms:kme)                ! Air pressure
    926. 

  926.       real pt0(kms:kme)               ! Reference potential temperature (could be equal to "pt")
    927. 

  927.       real qv(kms:kme)              ! Specific humidity
    928. 

  928.       real dt                    ! Time step
    929. 

  929.       real zr                    ! Zenith angle
    930. 

  930.       real deltar                ! Declination of the sun
    931. 

  931.       real ah                    ! Hour angle
    932. 

  932.       real rs                    ! Solar radiation
    933. 

  933.       real rld                   ! Downward flux of the longwave radiation
    934. 

  934. 

  935. ! Data relative to the "urban grid"
    936. 

  936. 

  937.       integer iurb               ! Current urban class
    938. 

  938. 

  939. !    Building parameters      
    940. 

  940.       real alag_u(nurbm)         ! Ground thermal diffusivity [m^2 s^-1]
    941. 

  941.       real alaw_u(nurbm)         ! Wall thermal diffusivity [m^2 s^-1]
    942. 

  942.       real alar_u(nurbm)         ! Roof thermal diffusivity [m^2 s^-1]
    943. 

  943.       real csg_u(nurbm)          ! Specific heat of the ground material [J m^3 K^-1]
    944. 

  944.       real csw_u(nurbm)          ! Specific heat of the wall material [J m^3 K^-1]
    945. 

  945.       real csr_u(nurbm)          ! Specific heat of the roof material [J m^3 K^-1]
    946. 

  946. 

  947. !    Radiation parameters
    948. 

  948.       real albg_u(nurbm)         ! Albedo of the ground
    949. 

  949.       real albw_u(nurbm)         ! Albedo of the wall
    950. 

  950.       real albr_u(nurbm)         ! Albedo of the roof
    951. 

  951.       real albwin_u(nurbm)       ! Albedo of the windows
    952. 

  952.       real emwind_u(nurbm)       ! Emissivity of windows
    953. 

  953.       real emg_u(nurbm)          ! Emissivity of ground
    954. 

  954.       real emw_u(nurbm)          ! Emissivity of wall
    955. 

  955.       real emr_u(nurbm)          ! Emissivity of roof
    956. 

  956. 

  957. !    fww,fwg,fgw,fsw,fsg are the view factors used to compute the long and 
    958. 

  958. !    short wave radation. 
    959. 

  959. !    The calculation of these factor is explained in the Appendix A of the BLM paper
    960. 

  960.       real fww(nz_um,nz_um,ndm,nurbm)  !  from wall to wall
    961. 

  961.       real fwg(nz_um,ndm,nurbm)        !  from wall to ground
    962. 

  962.       real fgw(nz_um,ndm,nurbm)        !  from ground to wall
    963. 

  963.       real fsw(nz_um,ndm,nurbm)        !  from sky to wall
    964. 

  964.       real fws(nz_um,ndm,nurbm)        !  from wall to sky
    965. 

  965.       real fsg(ndm,nurbm)              !  from sky to ground
    966. 

  966. 

  967. !    Roughness parameters
    968. 

  968.       real z0g_u(nurbm)            ! The ground's roughness length
    969. 

  969.       real z0r_u(nurbm)            ! The roof's roughness length
    970. 

  970.       
    971. 

  971. !    Street parameters
    972. 

  972.       integer nd_u(nurbm)          ! Number of street direction for each urban class 
    973. 

  973.       real strd_u(ndm,nurbm)       ! Street length (set to a greater value then the horizontal length of the cells)
    974. 

  974.       real drst_u(ndm,nurbm)       ! Street direction
    975. 

  975.       real ws_u(ndm,nurbm)         ! Street width
    976. 

  976.       real bs_u(ndm,nurbm)         ! Building width
    977. 

  977.       real h_b(nz_um,nurbm)        ! Bulding's heights
    978. 

  978.       real d_b(nz_um,nurbm)        ! The probability that a building has an height "h_b"
    979. 

  979.       real ss_u(nz_um,nurbm)     ! The probability that a building has an height equal to "z"
    980. 

  980.       real pb_u(nz_um,nurbm)     ! The probability that a building has an height greater or equal to "z"
    981. 

  981.         
    982. 

  982. !    Grid parameters
    983. 

  983.       integer nz_u(nurbm)     ! Number of layer in the urban grid
    984. 

  984.       real z_u(nz_um)       ! Height of the urban grid levels
    985. 

  985. 

  986. !   MT
    987. 

  987.       real cop_u(nurbm)
    988. 

  988.       real pwin_u(nurbm)
    989. 

  989.       real beta_u(nurbm)
    990. 

  990.       integer sw_cond_u(nurbm)
    991. 

  991.       real time_on_u(nurbm)
    992. 

  992.       real time_off_u(nurbm)
    993. 

  993.       real targtemp_u(nurbm)
    994. 

  994.       real gaptemp_u(nurbm)
    995. 

  995.       real targhum_u(nurbm)
    996. 

  996.       real gaphum_u(nurbm)
    997. 

  997.       real perflo_u(nurbm)
    998. 

  998.       real hsesf_u(nurbm)
    999. 

  999.       real hsequip(24)
    1000. 

  1000. 

  1001. ! ----------------------------------------------------------------------
    1002. 

  1002. ! INPUT-OUTPUT
    1003. 

  1003. ! ----------------------------------------------------------------------
    1004. 

  1004. 

  1005. ! Data relative to the "urban grid" which should be stored from the current time step to the next one
    1006. 

  1006. 

  1007.       real tw(2*ndm,nz_um,nwr_u,nbui_max)  ! Temperature in each layer of the wall [K]
    1008. 

  1008.       real tr(ndm,nz_um,nwr_u)  ! Temperature in each layer of the roof [K]
    1009. 

  1009.       real tg(ndm,ng_u)          ! Temperature in each layer of the ground [K]
    1010. 

  1010.       real sfw(2*ndm,nz_um,nbui_max)      ! Sensible heat flux from walls
    1011. 

  1011.       real sfg(ndm)              ! Sensible heat flux from ground (road)
    1012. 

  1012.       real sfr(ndm,nz_um)      ! Sensible heat flux from roofs
    1013. 

  1013.       real gfg(ndm)             ! Heat flux transferred from the surface of the ground (road) towards the interior
    1014. 

  1014.       real gfr(ndm,nz_um)     ! Heat flux transferred from the surface of the roof towards the interior
    1015. 

  1015.       real gfw(2*ndm,nz_um,nbui_max)     ! Heat flux transfered from the surface of the walls towards the interior
    1016. 

  1016. ! ----------------------------------------------------------------------
    1017. 

  1017. ! OUTPUT:
    1018. 

  1018. ! ----------------------------------------------------------------------
    1019. 

  1019. 

  1020. ! Data relative to the "mesoscale grid"
    1021. 

  1021. 

  1022.       real sf(kms:kme)             ! Surface of the "mesoscale grid" cells taking into account the buildings
    1023. 

  1023.       real vl(kms:kme)               ! Volume of the "mesoscale grid" cells taking into account the buildings
    1024. 

  1024.      
    1025. 

  1025. !    Implicit and explicit components of the source and sink terms at each levels,
    1026. 

  1026. !     the fluxes can be computed as follow: FX = A*X + B   example: Heat fluxes = a_t * pt + b_t
    1027. 

  1027.       real a_u(kms:kme)              ! Implicit component of the momentum sources or sinks in the X-direction
    1028. 

  1028.       real a_v(kms:kme)              ! Implicit component of the momentum sources or sinks in the Y-direction
    1029. 

  1029.       real a_t(kms:kme)              ! Implicit component of the heat sources or sinks
    1030. 

  1030.       real a_e(kms:kme)              ! Implicit component of the TKE sources or sinks
    1031. 

  1031.       real b_u(kms:kme)              ! Explicit component of the momentum sources or sinks in the X-direction
    1032. 

  1032.       real b_v(kms:kme)              ! Explicit component of the momentum sources or sinks in the Y-direction
    1033. 

  1033.       real b_t(kms:kme)              ! Explicit component of the heat sources or sinks
    1034. 

  1034.       real b_ac(kms:kme)
    1035. 

  1035.       real b_e(kms:kme)              ! Explicit component of the TKE sources or sinks
    1036. 

  1036.       real b_q(kms:kme)              ! Explicit component of the humidity sources or sinks
    1037. 

  1037.       real dlg(kms:kme)              ! Height above ground (L_ground in formula (24) of the BLM paper). 
    1038. 

  1038.       real dl_u(kms:kme)             ! Length scale (lb in formula (22) ofthe BLM paper).
    1039. 

  1039.       
    1040. 

  1040. ! ----------------------------------------------------------------------
    1041. 

  1041. ! LOCAL:
    1042. 

  1042. ! ----------------------------------------------------------------------
    1043. 

  1043. 

  1044.       real dz(kms:kme)               ! vertical space steps of the "mesoscale grid"
    1045. 

  1045. 

  1046. ! Data interpolated from the "mesoscale grid" to the "urban grid"
    1047. 

  1047. 

  1048.       real ua_u(nz_um)          ! Wind speed in the x direction
    1049. 

  1049.       real va_u(nz_um)          ! Wind speed in the y direction
    1050. 

  1050.       real pt_u(nz_um)          ! Potential temperature
    1051. 

  1051.       real da_u(nz_um)          ! Air density
    1052. 

  1052.       real pt0_u(nz_um)         ! Reference potential temperature
    1053. 

  1053.       real pr_u(nz_um)          ! Air pressure
    1054. 

  1054.       real qv_u(nz_um)          !Specific humidity
    1055. 

  1055. 

  1056. ! Data defining the building and street charateristics
    1057. 

  1057. 

  1058.       real alag(ng_u)           ! Ground thermal diffusivity for the current urban class [m^2 s^-1] 
    1059. 

  1059.       
    1060. 

  1060.       real csg(ng_u)            ! Specific heat of the ground material of the current urban class [J m^3 K^-1]
    1061. 

  1061.       real csr(nwr_u)            ! Specific heat of the roof material for the current urban class [J m^3 K^-1]
    1062. 

  1062.       real csw(nwr_u)            ! Specific heat of the wall material for the current urban class [J m^3 K^-1]
    1063. 

  1063. 

  1064.       real z0(ndm,nz_um)      ! Roughness lengths "profiles"
    1065. 

  1065.       real ws(ndm)              ! Street widths of the current urban class
    1066. 

  1066.       real bs(ndm)              ! Building widths of the current urban class
    1067. 

  1067.       real strd(ndm)            ! Street lengths for the current urban class
    1068. 

  1068.       real drst(ndm)            ! Street directions for the current urban class
    1069. 

  1069.       real ss(nz_um)          ! Probability to have a building with height h
    1070. 

  1070.       real pb(nz_um)          ! Probability to have a building with an height equal
    1071. 

  1071. 

  1072. ! Solar radiation at each level of the "urban grid"
    1073. 

  1073. 

  1074.       real rsg(ndm)             ! Short wave radiation from the ground
    1075. 

  1075.       real rsw(2*ndm,nz_um)     ! Short wave radiation from the walls
    1076. 

  1076.       real rlg(ndm)             ! Long wave radiation from the ground
    1077. 

  1077.       real rlw(2*ndm,nz_um)     ! Long wave radiation from the walls
    1078. 

  1078. 

  1079. ! Potential temperature of the surfaces at each level of the "urban grid"
    1080. 

  1080. 

  1081.       real ptg(ndm)             ! Ground potential temperatures 
    1082. 

  1082.       real ptr(ndm,nz_um)     ! Roof potential temperatures 
    1083. 

  1083.       real ptw(2*ndm,nz_um,nbui_max)     ! Walls potential temperatures 
    1084. 

  1084. 

  1085.  
    1086. 

  1086. ! Explicit and implicit component of the momentum, temperature and TKE sources or sinks on
    1087. 

  1087. ! vertical surfaces (walls) ans horizontal surfaces (roofs and street)
    1088. 

  1088. ! The fluxes can be computed as follow: Fluxes of X = A*X + B
    1089. 

  1089. ! Example: Momentum fluxes on vertical surfaces = uva_u * ua_u + uvb_u
    1090. 

  1090. 

  1091.       real uhb_u(ndm,nz_um)   ! U (wind component) Horizontal surfaces, B (explicit) term
    1092. 

  1092.       real uva_u(2*ndm,nz_um)   ! U (wind component)   Vertical surfaces, A (implicit) term
    1093. 

  1093.       real uvb_u(2*ndm,nz_um)   ! U (wind component)   Vertical surfaces, B (explicit) term
    1094. 

  1094.       real vhb_u(ndm,nz_um)   ! V (wind component) Horizontal surfaces, B (explicit) term
    1095. 

  1095.       real vva_u(2*ndm,nz_um)   ! V (wind component)   Vertical surfaces, A (implicit) term
    1096. 

  1096.       real vvb_u(2*ndm,nz_um)   ! V (wind component)   Vertical surfaces, B (explicit) term
    1097. 

  1097.       real thb_u(ndm,nz_um)   ! Temperature        Horizontal surfaces, B (explicit) term
    1098. 

  1098.       real tva_u(2*ndm,nz_um)   ! Temperature          Vertical surfaces, A (implicit) term
    1099. 

  1099.       real tvb_u(2*ndm,nz_um)   ! Temperature          Vertical surfaces, B (explicit) term
    1100. 

  1100.       real tvb_ac(2*ndm,nz_um)
    1101. 

  1101.       real ehb_u(ndm,nz_um)   ! Energy (TKE)       Horizontal surfaces, B (explicit) term
    1102. 

  1102.       real evb_u(2*ndm,nz_um)   ! Energy (TKE)         Vertical surfaces, B (explicit) term
    1103. 

  1103.       real qhb_u(ndm,nz_um)     ! Humidity      Horizontal surfaces, B (explicit) term
    1104. 

  1104.       real qvb_u(2*ndm,nz_um)   ! Humidity      Vertical surfaces, B (explicit) term      
    1105. 

  1105. !
    1106. 

  1106.       real rs_abs ! solar radiation absorbed by urban surfaces 
    1107. 

  1107.       real rl_up ! longwave radiation emitted by urban surface to the atmosphere 
    1108. 

  1108.       real emiss ! mean emissivity of the urban surface
    1109. 

  1109.       real grdflx_urb ! ground heat flux
    1110. 

  1110.       real dt_int ! internal time step
    1111. 

  1111.       integer nt_int ! number of internal time step
    1112. 

  1112.       integer iz,id, it_int
    1113. 

  1113.       integer iw,ix,iy
    1114. 

  1114. 

  1115. !---------------------------------------
    1116. 

  1116. !New variables uses in BEM
    1117. 

  1117. !----------------------------------------
    1118. 

  1118.    
    1119. 

  1119.       real tmp_u(nz_um)     !Air Temperature [K]
    1120. 

  1120. 

  1121.       real dzw(nwr_u)       !Layer sizes in the walls
    1122. 

  1122.       real dzr(nwr_u)       !Layer sizes in the roofs
    1123. 

  1123.       real dzf(nf_u)        !Layer sizes in the floors
    1124. 

  1124.       real dzgb(ngb_u)      !Layer sizes in the ground below the buildings
    1125. 

  1125. 

  1126.       real csgb(ngb_u)      !Specific heat of the ground material below the buildings 
    1127. 

  1127.                             !of the current urban class at each ground levels[J m^3 K^-1]
    1128. 

  1128.       real csf(nf_u)        !Specific heat of the floors mat…
    1129.

Large files files are truncated, but you can click here to view the full file