Changeset 3241 for palm/trunk/SOURCE/turbulence_closure_mod.f90
- Timestamp:
- Sep 12, 2018 3:02:00 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/turbulence_closure_mod.f90
r3183 r3241 25 25 ! ----------------- 26 26 ! $Id$ 27 ! unused variables removed, shortest_distance has wp now 28 ! 29 ! 3183 2018-07-27 14:25:55Z suehring 27 30 ! Rename variables and remove unused variable from USE statement 28 31 ! … … 200 203 REAL(wp) :: c_1 !< model constant for RANS mode 201 204 REAL(wp) :: c_2 !< model constant for RANS mode 202 REAL(wp) :: c_3 !< model constant for RANS mode205 ! REAL(wp) :: c_3 !< model constant for RANS mode 203 206 REAL(wp) :: c_4 !< model constant for RANS mode 204 207 REAL(wp) :: l_max !< maximum length scale for Blackadar mixing length 205 208 REAL(wp) :: dsig_e = 1.0_wp !< factor to calculate Ke from Km (1/sigma_e) 206 209 REAL(wp) :: dsig_diss = 1.0_wp !< factor to calculate K_diss from Km (1/sigma_diss) 207 INTEGER(iwp) :: surf_e !< end index of surface elements at given i-j position208 INTEGER(iwp) :: surf_s !< start index of surface elements at given i-j position209 210 210 211 REAL(wp), DIMENSION(0:4) :: rans_const_c = & !< model constants for RANS mode (namelist param) … … 382 383 383 384 USE control_parameters, & 384 ONLY: message_string, neutral,turbulent_inflow, turbulent_outflow385 ONLY: message_string, turbulent_inflow, turbulent_outflow 385 386 386 387 IMPLICIT NONE … … 459 460 !> Check data output. 460 461 !------------------------------------------------------------------------------! 461 SUBROUTINE tcm_check_data_output( var, unit , i, ilen, k)462 SUBROUTINE tcm_check_data_output( var, unit ) 462 463 463 USE control_parameters, &464 ONLY: data_output, message_string465 466 464 IMPLICIT NONE 467 465 … … 469 467 CHARACTER (LEN=*) :: var !< name of output variable 470 468 471 INTEGER(iwp) :: i !< index of var in data_output472 INTEGER(iwp) :: ilen !< length of var string473 INTEGER(iwp) :: k !< flag if var contains one of '_xy', '_xz' or '_yz'474 469 475 470 SELECT CASE ( TRIM( var ) ) … … 697 692 !------------------------------------------------------------------------------! 698 693 SUBROUTINE tcm_data_output_2d( av, variable, found, grid, mode, local_pf, & 699 two_d,nzb_do, nzt_do )694 nzb_do, nzt_do ) 700 695 701 696 USE averaging, & … … 717 712 718 713 LOGICAL :: found !< flag if output variable is found 719 LOGICAL :: two_d !< flag parameter that indicates 2D variables (horizontal cross sections)720 714 LOGICAL :: resorted !< flag if output is already resorted 721 715 … … 1052 1046 1053 1047 USE control_parameters, & 1054 ONLY: bc_dirichlet_l, complex_terrain, dissipation_1d,topography1048 ONLY: bc_dirichlet_l, complex_terrain, topography 1055 1049 1056 1050 USE model_1d_mod, & 1057 ONLY: diss1d, e1d, kh1d, km1d, l1d1051 ONLY: e1d, kh1d, km1d 1058 1052 1059 1053 USE surface_mod, & … … 1339 1333 INTEGER(iwp) :: ii !< index variable along x 1340 1334 INTEGER(iwp) :: j !< index variable along y 1341 INTEGER(iwp) :: jj !< index variable along y1342 1335 INTEGER(iwp) :: k !< index variable along z 1343 1336 INTEGER(iwp) :: k_max_topo = 0 !< index of maximum topography height … … 1804 1797 !> and debug options for the compiler. This should be fixed 1805 1798 !------------------------------------------------------------------------------! 1806 REAL FUNCTION shortest_distance( array, orientation, pos_i )1799 REAL(wp) FUNCTION shortest_distance( array, orientation, pos_i ) 1807 1800 1808 1801 IMPLICIT NONE … … 2089 2082 2090 2083 USE surface_mod, & 2091 ONLY : surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h, & 2092 surf_usm_v 2084 ONLY : surf_def_h 2093 2085 2094 2086 IMPLICIT NONE … … 2371 2363 2372 2364 USE arrays_3d, & 2373 ONLY: ddzu, diss_l_diss, diss_l_e, diss_s_diss, diss_s_e, & 2374 flux_l_diss, flux_l_e, flux_s_diss, flux_s_e,& 2375 u_p,v_p,w_p 2365 ONLY: diss_l_diss, diss_l_e, diss_s_diss, diss_s_e, flux_l_diss, & 2366 flux_l_e, flux_s_diss, flux_s_e 2376 2367 2377 2368 USE control_parameters, & 2378 ONLY: f, tsc 2379 2380 USE grid_variables, & 2381 ONLY: dx, dy 2382 2383 USE surface_mod, & 2384 ONLY : surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h, & 2385 surf_usm_v 2386 2387 use indices, only: nx, ny 2369 ONLY: tsc 2388 2370 2389 2371 IMPLICIT NONE … … 2393 2375 INTEGER(iwp) :: j !< loop index y direction 2394 2376 INTEGER(iwp) :: k !< loop index z direction 2395 INTEGER(iwp) :: l !< loop index2396 INTEGER(iwp) :: m !< loop index2397 INTEGER(iwp) :: surf_e !< end index of surface elements at given i-j position2398 INTEGER(iwp) :: surf_s !< start index of surface elements at given i-j position2399 2377 INTEGER(iwp) :: tn !< task number of openmp task 2400 2378 2401 INTEGER(iwp) :: pis = 32 !< debug variable, print from i=pis !> @todo remove later2402 INTEGER(iwp) :: pie = 32 !< debug variable, print until i=pie !> @todo remove later2403 INTEGER(iwp) :: pjs = 26 !< debug variable, print from j=pjs !> @todo remove later2404 INTEGER(iwp) :: pje = 26 !< debug variable, print until j=pje !> @todo remove later2405 INTEGER(iwp) :: pkb = 1 !< debug variable, print from k=pkb !> @todo remove later2406 INTEGER(iwp) :: pkt = 7 !< debug variable, print until k=pkt !> @todo remove later2379 ! INTEGER(iwp) :: pis = 32 !< debug variable, print from i=pis !> @todo remove later 2380 ! INTEGER(iwp) :: pie = 32 !< debug variable, print until i=pie !> @todo remove later 2381 ! INTEGER(iwp) :: pjs = 26 !< debug variable, print from j=pjs !> @todo remove later 2382 ! INTEGER(iwp) :: pje = 26 !< debug variable, print until j=pje !> @todo remove later 2383 ! INTEGER(iwp) :: pkb = 1 !< debug variable, print from k=pkb !> @todo remove later 2384 ! INTEGER(iwp) :: pkt = 7 !< debug variable, print until k=pkt !> @todo remove later 2407 2385 2408 2386 REAL(wp), DIMENSION(nzb:nzt+1) :: dum_adv !< debug variable !> @todo remove later … … 2410 2388 REAL(wp), DIMENSION(nzb:nzt+1) :: dum_dif !< debug variable !> @todo remove later 2411 2389 2412 5555 FORMAT(A,7(1X,E12.5)) !> @todo remove later2390 !5555 FORMAT(A,7(1X,E12.5)) !> @todo remove later 2413 2391 2414 2392 ! … … 4812 4790 INTEGER(iwp) :: j !< loop index 4813 4791 INTEGER(iwp) :: k !< loop index 4814 4792 !$ INTEGER(iwp) :: omp_get_thread_num !< opemmp function to get thread number 4815 4793 INTEGER(iwp) :: sr !< statistic region 4816 4794 INTEGER(iwp) :: tn !< thread number … … 4970 4948 REAL(wp) :: vc(-1:1,-1:1) !< v on grid center 4971 4949 REAL(wp) :: wc(-1:1,-1:1) !< w on grid center 4972 REAL(wp) :: u2(-1:1,-1:1) !< u2 on grid center4973 REAL(wp) :: v2(-1:1,-1:1) !< v2 on grid center4974 REAL(wp) :: w2(-1:1,-1:1) !< w2 on grid center4975 REAL(wp) :: uv(-1:1,-1:1) !< u*v on grid center4976 REAL(wp) :: uw(-1:1,-1:1) !< u*w on grid center4977 REAL(wp) :: vw(-1:1,-1:1) !< v*w on grid center4950 ! REAL(wp) :: u2(-1:1,-1:1) !< u2 on grid center 4951 ! REAL(wp) :: v2(-1:1,-1:1) !< v2 on grid center 4952 ! REAL(wp) :: w2(-1:1,-1:1) !< w2 on grid center 4953 ! REAL(wp) :: uv(-1:1,-1:1) !< u*v on grid center 4954 ! REAL(wp) :: uw(-1:1,-1:1) !< u*w on grid center 4955 ! REAL(wp) :: vw(-1:1,-1:1) !< v*w on grid center 4978 4956 4979 4957 REAL(wp) :: ut(nzb:nzt+1,nysg:nyng,nxlg:nxrg) !< test filtered u … … 5229 5207 IMPLICIT NONE 5230 5208 5209 5210 #if defined( __nopointer ) 5231 5211 INTEGER(iwp) :: i !< loop index x direction 5232 5212 INTEGER(iwp) :: j !< loop index y direction 5233 5213 INTEGER(iwp) :: k !< loop index z direction 5214 #endif 5234 5215 INTEGER, INTENT(IN) :: mod_count !< flag defining where pointers point to 5235 5216
Note: See TracChangeset
for help on using the changeset viewer.