Changeset 3241 for palm/trunk/SOURCE/lpm_advec.f90
- Timestamp:
- Sep 12, 2018 3:02:00 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/lpm_advec.f90
r3207 r3241 25 25 ! ----------------- 26 26 ! $Id$ 27 ! unused variables removed 28 ! 29 ! 3207 2018-08-27 12:55:33Z schwenkel 27 30 ! Minor bugfix for sgs-velocities in case of cloud droplets 28 31 ! … … 175 178 176 179 USE control_parameters, & 177 ONLY: atmos_ocean_sign, cloud_droplets, constant_flux_layer, dt_3d,&178 d t_3d_reached_l, dz, g, kappa, topography, u_gtrans, v_gtrans180 ONLY: cloud_droplets, constant_flux_layer, dt_3d, dt_3d_reached_l, & 181 dz, g, kappa, topography, u_gtrans, v_gtrans 179 182 180 183 USE grid_variables, & 181 ONLY: d dx, dx, ddy, dy184 ONLY: dx, dy 182 185 183 186 USE indices, & … … 187 190 188 191 USE particle_attributes, & 189 ONLY: block_offset, c_0, dt_min_part, grid_particles, 190 iran_part, log_z_z0, number_of_particles, number_of_sublayers,&191 particles, particle_groups, offset_ocean_nzt, sgs_wf_part,&192 ONLY: block_offset, c_0, dt_min_part, grid_particles, iran_part, & 193 log_z_z0, number_of_particles, number_of_sublayers, & 194 particles, particle_groups, sgs_wf_part, & 192 195 use_sgs_for_particles, vertical_particle_advection, z0_av_global 193 196 … … 202 205 LOGICAL :: subbox_at_wall !< flag to see if the current subgridbox is adjacent to a wall 203 206 204 INTEGER(iwp) :: agp !< loop variable205 INTEGER(iwp) :: gp_outside_of_building(1:8) !< number of grid points used for particle interpolation in case of topography206 207 INTEGER(iwp) :: i !< index variable along x 207 208 INTEGER(iwp) :: ip !< index variable along x … … 214 215 INTEGER(iwp) :: n !< loop variable over all particles in a grid box 215 216 INTEGER(iwp) :: nb !< block number particles are sorted in 216 INTEGER(iwp) :: num_gp !< number of adjacent grid points inside topography217 217 INTEGER(iwp) :: surf_start !< Index on surface data-type for current grid box 218 218 … … 223 223 REAL(wp) :: bb !< dummy argument for horizontal particle interpolation 224 224 REAL(wp) :: cc !< dummy argument for horizontal particle interpolation 225 REAL(wp) :: d_sum !< dummy argument for horizontal particle interpolation in case of topography226 225 REAL(wp) :: d_z_p_z0 !< inverse of interpolation length for logarithmic interpolation 227 226 REAL(wp) :: dd !< dummy argument for horizontal particle interpolation … … 246 245 REAL(wp) :: gg !< dummy argument for horizontal particle interpolation 247 246 REAL(wp) :: height_p !< dummy argument for logarithmic interpolation 248 REAL(wp) :: location(1:30,1:3) !< wall locations249 247 REAL(wp) :: log_z_z0_int !< logarithmus used for surface_layer interpolation 250 248 REAL(wp) :: random_gauss !< … … 276 274 REAL(wp), PARAMETER :: d0_rog = 0.745_wp !< separation diameter 277 275 278 REAL(wp), DIMENSION(1:30) :: d_gp_pl !< dummy argument for particle interpolation scheme in case of topography279 REAL(wp), DIMENSION(1:30) :: de_dxi !< horizontal TKE gradient along x at adjacent wall280 REAL(wp), DIMENSION(1:30) :: de_dyi !< horizontal TKE gradient along y at adjacent wall281 REAL(wp), DIMENSION(1:30) :: de_dzi !< horizontal TKE gradient along z at adjacent wall282 REAL(wp), DIMENSION(1:30) :: dissi !< dissipation at adjacent wall283 REAL(wp), DIMENSION(1:30) :: ei !< TKE at adjacent wall284 285 276 REAL(wp), DIMENSION(number_of_particles) :: term_1_2 !< flag to communicate whether a particle is near topography or not 286 277 REAL(wp), DIMENSION(number_of_particles) :: dens_ratio !< … … 1151 1142 REAL(wp) :: fac !< flag to identify adjacent topography 1152 1143 REAL(wp) :: fs_n !< weighting factor to prevent that subgrid-scale particle speed becomes too large 1153 REAL(wp) :: sgs_w !< constant (1/3)1154 1144 REAL(wp) :: rg_n !< random number 1155 1145 REAL(wp) :: term1 !< memory term … … 1183 1173 ! 1184 1174 !-- Random term 1185 term3 = SQRT( MAX( a1, 1E-20 ) ) * ( rg_n - 1.0_wp ) * SQRT( dt_n )1175 term3 = SQRT( MAX( a1, 1E-20_wp ) ) * ( rg_n - 1.0_wp ) * SQRT( dt_n ) 1186 1176 ! 1187 1177 !-- In cese one of the adjacent grid-boxes belongs to topograhy, the previous
Note: See TracChangeset
for help on using the changeset viewer.