Changeset 3065 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
- Timestamp:
- Jun 12, 2018 7:03:02 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
r3051 r3065 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Error message related to vertical stretching has been added, dz was replaced 28 ! by dz(1) 29 ! 30 ! 3051 2018-05-30 17:43:55Z suehring 27 31 ! Bugfix in calculation of initial Reynolds-stress tensor. 28 32 ! … … 313 317 314 318 USE control_parameters, & 315 ONLY: bc_lr, bc_ns, forcing, nest_domain, rans_mode, turbulent_inflow 319 ONLY: bc_lr, bc_ns, forcing, nest_domain, number_stretch_level_start, & 320 rans_mode, turbulent_inflow 316 321 317 322 USE pmc_interface, & … … 365 370 'is not allowed' 366 371 CALL message( 'stg_check_parameters', 'PA0039', 1, 2, 0, 6, 0 ) 372 ENDIF 373 374 IF ( number_stretch_level_start > 0 ) THEN 375 message_string = 'Using synthetic turbulence generator ' // & 376 'in combination with stretching is not allowed' 377 CALL message( 'stg_check_parameters', 'PA0420', 1, 2, 0, 6, 0 ) 367 378 ENDIF 368 379 … … 607 618 608 619 ! 609 !-- Convert length scales from meter to number of grid points 620 !-- Convert length scales from meter to number of grid points. Attention: 621 !-- Does not work if grid stretching is used 610 622 nuy(k) = INT( luy * ddy ) 611 nuz(k) = INT( luz / dz )623 nuz(k) = INT( luz / dz(1) ) 612 624 nvy(k) = INT( lvy * ddy ) 613 nvz(k) = INT( lvz / dz )625 nvz(k) = INT( lvz / dz(1) ) 614 626 nwy(k) = INT( lwy * ddy ) 615 nwz(k) = INT( lwz / dz )627 nwz(k) = INT( lwz / dz(1) ) 616 628 ! 617 629 !-- Workaround, assume isotropic turbulence
Note: See TracChangeset
for help on using the changeset viewer.