Changeset 4332 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
 Timestamp:
 Dec 10, 2019 7:44:12 PM (3 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
r4329 r4332 25 25 !  26 26 ! $Id$ 27 ! Limit initial velocity seeds in restart runs, if not the seed calculation 28 ! may become unstable. Further, minor bugfix in initial velocity seed 29 ! calculation. 30 ! 31 ! 4329 20191210 15:46:36Z motisi 27 32 ! Renamed wall_flags_0 to wall_flags_static_0 28 33 ! … … 903 908 DO j = nysg, nyng 904 909 DO k = nzb, nzt+1 905 IF ( a11(k) /= 0.0_wp ) THEN910 IF ( a11(k) > 10E8_wp ) THEN 906 911 fu_yz(k,j) = ( u(k,j,i)  u_init(k) ) / a11(k) 907 912 ELSE 908 fu_yz(k,j) = 0.0_wp913 fu_yz(k,j) = 10E8_wp 909 914 ENDIF 910 915 911 IF ( a22(k) /= 0.0_wp ) THEN916 IF ( a22(k) > 10E8_wp ) THEN 912 917 fv_yz(k,j) = ( v(k,j,i)  & 913 918 a21(k) * fu_yz(k,j)  v_init(k) ) / a22(k) 914 919 ELSE 915 fv_yz(k,j) = 0.0_wp920 fv_yz(k,j) = 10E8_wp 916 921 ENDIF 917 922 918 IF ( a33(k) /= 0.0_wp ) THEN923 IF ( a33(k) > 10E8_wp ) THEN 919 924 fw_yz(k,j) = ( w(k,j,i)  & 920 925 a31(k) * fu_yz(k,j)  a32(k) * & 921 926 fv_yz(k,j) ) / a33(k) 922 927 ELSE 923 fw_yz = 0._wp928 fw_yz(k,j) = 10E8_wp 924 929 ENDIF 925 926 930 ENDDO 927 931 ENDDO … … 936 940 DO k = nzb, nzt+1 937 941 938 IF ( a11(k) /= 0.0_wp ) THEN942 IF ( a11(k) > 10E8_wp ) THEN 939 943 fu_xz(k,i) = ( u(k,j,i)  u_init(k) ) / a11(k) 940 944 ELSE 941 fu_xz(k,i) = 0.0_wp945 fu_xz(k,i) = 10E8_wp 942 946 ENDIF 943 947 944 IF ( a22(k) /= 0.0_wp ) THEN948 IF ( a22(k) > 10E8_wp ) THEN 945 949 fv_xz(k,i) = ( v(k,j,i)  & 946 950 a21(k) * fu_xz(k,i)  v_init(k) ) / a22(k) 947 951 ELSE 948 fv_xz(k,i) = 0.0_wp952 fv_xz(k,i) = 10E8_wp 949 953 ENDIF 950 954 951 IF ( a33(k) /= 0.0_wp ) THEN955 IF ( a33(k) > 10E8_wp ) THEN 952 956 fw_xz(k,i) = ( w(k,j,i)  & 953 957 a31(k) * fu_xz(k,i)  & 954 958 a32(k) * fv_xz(k,i) ) / a33(k) 955 959 ELSE 956 fw_xz = 0._wp960 fw_xz(k,i) = 10E8_wp 957 961 ENDIF 958 962
Note: See TracChangeset
for help on using the changeset viewer.