Changeset 4585
 Timestamp:
 Jun 30, 2020 3:05:20 PM (3 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90
r4546 r4585 25 25 !  26 26 ! $Id$ 27 ! Limit logarithmically interpolated particle speed to the velocity component 28 ! at the first prognostic grid point (since no stability corrected interpolation 29 ! is employed the particle speed could be overestimated in unstable conditions). 30 ! 31 ! 4546 20200524 12:16:41Z raasch 27 32 ! Variables iran and iran_part completely removed, added I/O of parallel random numbers to restart 28 33 ! file … … 3717 3722 ! as sensitivity studies revealed no significant effect of 3718 3723 ! using the neutral solution also for un/stable situations. 3724 ! However, as the stability correction is missing the particle 3725 ! speed may be overestimated compared to the velocity component 3726 ! at the first prognostic grid point. Hence, the particle speed 3727 ! is limited. 3719 3728 u_int(n) = usws_int / ( us_int * kappa + 1E10_wp ) & 3720 3729 * log_z_z0_int  u_gtrans 3730 3731 u_int(n) = SIGN( MIN( ABS( u_int(n) ), ABS( u(k_wall+1,jp,ip) ) ), u_int(n) ) 3732 3721 3733 ENDIF 3722 3734 ! … … 3812 3824 ! as sensitivity studies revealed no significant effect of 3813 3825 ! using the neutral solution also for un/stable situations. 3826 ! However, as the stability correction is missing the particle 3827 ! speed may be overestimated compared to the velocity component 3828 ! at the first prognostic grid point. Hence, the particle speed 3829 ! is limited. 3814 3830 v_int(n) = vsws_int / ( us_int * kappa + 1E10_wp ) & 3815 3831 * log_z_z0_int  v_gtrans 3832 3833 v_int(n) = SIGN( MIN( ABS( v_int(n) ), ABS( v(k_wall+1,jp,ip) ) ), v_int(n) ) 3816 3834 3817 3835 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.