Ignore:
Timestamp:
Jun 12, 2018 7:03:02 AM (3 years ago)
Author:
Giersch
Message:

New vertical stretching procedure has been introduced

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

    r3051 r3065  
    2525! -----------------
    2626! $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
    2731! Bugfix in calculation of initial Reynolds-stress tensor.
    2832!
     
    313317
    314318    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
    316321
    317322    USE pmc_interface,                                                         &
     
    365370                              'is not allowed'
    366371          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 )
    367378       ENDIF
    368379
     
    607618
    608619!
    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
    610622          nuy(k) = INT( luy * ddy )
    611           nuz(k) = INT( luz / dz  )
     623          nuz(k) = INT( luz / dz(1)  )
    612624          nvy(k) = INT( lvy * ddy )
    613           nvz(k) = INT( lvz / dz  )
     625          nvz(k) = INT( lvz / dz(1)  )
    614626          nwy(k) = INT( lwy * ddy )
    615           nwz(k) = INT( lwz / dz  )
     627          nwz(k) = INT( lwz / dz(1)  )
    616628!
    617629!--       Workaround, assume isotropic turbulence
Note: See TracChangeset for help on using the changeset viewer.