Changeset 3981 for palm/trunk/SOURCE
 Timestamp:
 May 15, 2019 2:43:01 PM (5 years ago)
 Location:
 palm/trunk/SOURCE
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/lpm_init.f90
r3655 r3981 25 25 !  26 26 ! $Id$ 27 ! Bugfix in particle nesting, initialize TKE gradient with zero 28 ! 29 ! 3655 20190107 16:51:22Z knoop 27 30 ! Changed number concentration unit to SI standard 28 31 ! … … 431 434 de_dy(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 432 435 de_dz(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 436 de_dx = 0.0_wp 437 de_dy = 0.0_wp 438 de_dz = 0.0_wp 433 439 434 440 sgs_wf_part = 1.0_wp / 3.0_wp 
palm/trunk/SOURCE/lpm_init_sgs_tke.f90
r3655 r3981 25 25 !  26 26 ! $Id$ 27 ! Bugfix in particle nesting, set boundary conditions at noncyclic boundaries 28 ! 29 ! 3655 20190107 16:51:22Z knoop 27 30 ! unused variables removed 28 31 ! … … 81 84 ONLY: de_dx, de_dy, de_dz, diss, e, u, v, w, zu 82 85 86 USE control_parameters, & 87 ONLY: bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r, bc_dirichlet_s 88 83 89 USE grid_variables, & 84 90 ONLY: ddx, ddy … … 206 212 CALL exchange_horiz( de_dz, nbgp ) 207 213 CALL exchange_horiz( diss, nbgp ) 208 209 214 ! 215 ! Set boundary conditions at nonperiodic boundaries. Note, at nonperiod 216 ! boundaries zerogradient boundary conditions are set, so that the gradient 217 ! normal to the boundary is zero. 218 IF ( bc_dirichlet_l ) THEN 219 de_dx(:,:,1) = 0.0_wp 220 de_dy(:,:,1) = de_dy(:,:,0) 221 de_dz(:,:,1) = de_dz(:,:,0) 222 ENDIF 223 IF ( bc_dirichlet_r ) THEN 224 de_dx(:,:,nxr+1) = 0.0_wp 225 de_dy(:,:,nxr+1) = de_dy(:,:,nxr) 226 de_dz(:,:,nxr+1) = de_dz(:,:,nxr) 227 ENDIF 228 IF ( bc_dirichlet_n ) THEN 229 de_dx(:,nyn+1,:) = de_dx(:,nyn,:) 230 de_dy(:,nyn+1,:) = 0.0_wp 231 de_dz(:,nyn+1,:) = de_dz(:,nyn,:) 232 ENDIF 233 IF ( bc_dirichlet_s ) THEN 234 de_dx(:,nys1,:) = de_dx(:,nys,:) 235 de_dy(:,nys1,:) = 0.0_wp 236 de_dz(:,nys1,:) = de_dz(:,nys,:) 237 ENDIF 210 238 ! 211 239 ! Calculate the horizontally averaged profiles of SGS TKE and resolved
Note: See TracChangeset
for help on using the changeset viewer.