Changeset 3832 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Mar 28, 2019 1:16:58 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r3824 r3832 28 28 ! ----------------- 29 29 ! $Id$ 30 ! instrumented with openmp directives 31 ! 32 ! 3824 2019-03-27 15:56:16Z pavelkrc 30 33 ! Remove unused imports 31 34 ! … … 5332 5335 LOGICAL :: spinup !< if true, no calculation of window temperatures 5333 5336 5337 !$OMP PARALLEL PRIVATE (m, i, j, k, kw, wtend, wintend, win_absorp, wall_mod) 5334 5338 wall_mod=1.0_wp 5335 5339 IF (usm_wall_mod .AND. spinup) THEN … … 5341 5345 ! 5342 5346 !-- For horizontal surfaces 5347 !$OMP DO SCHEDULE (STATIC) 5343 5348 DO m = 1, surf_usm_h%ns 5344 5349 ! … … 5521 5526 ! 5522 5527 !-- For vertical surfaces 5528 !$OMP DO SCHEDULE (STATIC) 5523 5529 DO l = 0, 3 5524 5530 DO m = 1, surf_usm_v(l)%ns … … 5699 5705 ENDDO 5700 5706 ENDDO 5707 !$OMP END PARALLEL 5701 5708 5702 5709 END SUBROUTINE usm_material_heat_model … … 5734 5741 ! 5735 5742 !-- For horizontal surfaces 5743 !$OMP PARALLEL PRIVATE (m, i, j, k, kw, lambda_h_green_sat, ke, lambda_green_temp, gtend, & 5744 !$OMP& tend, h_vg, gamma_green_temp, m_total, root_extr_green) 5745 !$OMP DO SCHEDULE (STATIC) 5736 5746 DO m = 1, surf_usm_h%ns 5737 5747 … … 5980 5990 5981 5991 ENDDO 5992 !$OMP END PARALLEL 5982 5993 5983 5994 ! … … 7752 7763 ! 7753 7764 !-- First, treat horizontal surface elements 7765 !$OMP PARALLEL PRIVATE (m, i, j, k, lambda_surface, lambda_surface_window, & 7766 !$OMP& lambda_surface_green, qv1, rho_cp, rho_lv, drho_l_lv, f_shf, & 7767 !$OMP& f_shf_window, f_shf_green, m_total, f1, f2, e_s, e, f3, f_qsws_veg,& 7768 !$OMP& q_s, f_qsws_liq, f_qsws, e_s_dt, dq_s_dt, coef_1, coef_window_1, & 7769 !$OMP& coef_green_1, coef_2, coef_window_2, coef_green_2, stend_wall, & 7770 !$OMP& stend_window, stend_green, tend, m_liq_max) 7771 !$OMP DO SCHEDULE (STATIC) 7754 7772 DO m = 1, surf_usm_h%ns 7755 7773 ! … … 8197 8215 ENDDO 8198 8216 ! 8199 !-- Now, treat vertical surface elements 8217 !-- Now, treat vertical surface elements 8218 !$OMP DO SCHEDULE (STATIC) 8200 8219 DO l = 0, 3 8201 8220 DO m = 1, surf_usm_v(l)%ns … … 8574 8593 8575 8594 ENDDO 8595 !$OMP END PARALLEL 8596 8576 8597 ! 8577 8598 !-- Add-up anthropogenic heat, for now only at upward-facing surfaces … … 8586 8607 dtime = mod(simulated_time + time_utc_init, 24.0_wp*3600.0_wp) 8587 8608 dhour = INT(dtime/3600.0_wp) 8609 8610 !-- TO_DO: activate, if testcase is available 8611 !-- !$OMP PARALLEL DO PRIVATE (i, j, k, acoef, rho_cp) 8612 !-- it may also improve performance to move get_topography_top_index_ji before the k-loop 8588 8613 DO i = nxl, nxr 8589 8614 DO j = nys, nyn
Note: See TracChangeset
for help on using the changeset viewer.