Changeset 3909
- Timestamp:
- Apr 17, 2019 9:13:25 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
r3900 r3909 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Minor bugfix for last commit 28 ! 29 ! 3900 2019-04-16 15:17:43Z suehring 27 30 ! Missing re-calculation of perturbation seeds in case of restarts 28 31 ! … … 923 926 924 927 IF ( a22(k) .NE. 0._wp ) THEN 925 fv_xz(k,i) = ( v(k,j,i) / mc_factor - a21(k) * fu_ yz(k,i) - &928 fv_xz(k,i) = ( v(k,j,i) / mc_factor - a21(k) * fu_xz(k,i) - & 926 929 v_init(k) ) / a22(k) 927 930 ELSE … … 930 933 931 934 IF ( a33(k) .NE. 0._wp ) THEN 932 fw_xz(k,i) = ( w(k,j,i) / mc_factor - a31(k) * fu_ yz(k,i) - &933 a32(k) * fv_ yz(k,i) ) / a33(k)935 fw_xz(k,i) = ( w(k,j,i) / mc_factor - a31(k) * fu_xz(k,i) - & 936 a32(k) * fv_xz(k,i) ) / a33(k) 934 937 ELSE 935 938 fw_xz = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.