Changeset 3182 for palm/trunk/SOURCE/turbulence_closure_mod.f90
 Timestamp:
 Jul 27, 2018 1:36:03 PM (3 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/turbulence_closure_mod.f90
r3145 r3182 20 20 ! Current revisions: 21 21 !  22 ! 22 ! Rename variables and remove unused variable from USE statement 23 23 ! 24 24 ! Former revisions: … … 149 149 150 150 USE control_parameters, & 151 ONLY: constant_diffusion, dt_3d, e_init, humidity, inflow_l,&151 ONLY: constant_diffusion, dt_3d, e_init, humidity, & 152 152 initializing_actions, intermediate_timestep_count, & 153 153 intermediate_timestep_count_max, kappa, km_constant, & … … 212 212 213 213 REAL(wp), DIMENSION(:), ALLOCATABLE :: l_black !< mixing length according to Blackadar 214 214 215 REAL(wp), DIMENSION(:), ALLOCATABLE :: l_grid !< geometric mean of grid sizes dx, dy, dz 215 216 … … 378 379 379 380 USE control_parameters, & 380 ONLY: message_string, nest_domain, neutral, turbulent_inflow, & 381 turbulent_outflow 381 ONLY: message_string, neutral, turbulent_inflow, turbulent_outflow 382 382 383 383 IMPLICIT NONE … … 1049 1049 1050 1050 USE control_parameters, & 1051 ONLY: complex_terrain, dissipation_1d, topography1051 ONLY: bc_dirichlet_l, complex_terrain, dissipation_1d, topography 1052 1052 1053 1053 USE model_1d_mod, & … … 1252 1252 ! Use these mean profiles at the inflow (provided that Dirichlet 1253 1253 ! conditions are used) 1254 IF ( inflow_l ) THEN1254 IF ( bc_dirichlet_l ) THEN 1255 1255 DO j = nysg, nyng 1256 1256 DO k = nzb, nzt+1 … … 4613 4613 4614 4614 USE control_parameters, & 4615 ONLY: e_min, outflow_l, outflow_n, outflow_r, outflow_s 4615 ONLY: bc_radiation_l, bc_radiation_n, bc_radiation_r, bc_radiation_s, & 4616 e_min 4616 4617 4617 4618 USE surface_layer_fluxes_mod, & … … 4772 4773 ! Set Neumann boundary conditions at the outflow boundaries in case of 4773 4774 ! noncyclic lateral boundaries 4774 IF ( outflow_l ) THEN4775 IF ( bc_radiation_l ) THEN 4775 4776 km(:,:,nxl1) = km(:,:,nxl) 4776 4777 kh(:,:,nxl1) = kh(:,:,nxl) 4777 4778 ENDIF 4778 IF ( outflow_r ) THEN4779 IF ( bc_radiation_r ) THEN 4779 4780 km(:,:,nxr+1) = km(:,:,nxr) 4780 4781 kh(:,:,nxr+1) = kh(:,:,nxr) 4781 4782 ENDIF 4782 IF ( outflow_s ) THEN4783 IF ( bc_radiation_s ) THEN 4783 4784 km(:,nys1,:) = km(:,nys,:) 4784 4785 kh(:,nys1,:) = kh(:,nys,:) 4785 4786 ENDIF 4786 IF ( outflow_n ) THEN4787 IF ( bc_radiation_n ) THEN 4787 4788 km(:,nyn+1,:) = km(:,nyn,:) 4788 4789 kh(:,nyn+1,:) = kh(:,nyn,:)
Note: See TracChangeset
for help on using the changeset viewer.