Changeset 4144 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
 Timestamp:
 Aug 6, 2019 9:11:47 AM (23 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
r4071 r4144 25 25 !  26 26 ! $Id$ 27 ! relational operators .EQ., .NE., etc. replaced by ==, /=, etc. 28 ! 29 ! 4071 20190703 20:02:00Z suehring 27 30 ! Bugfix, initialize mean_inflow_profiles in case turbulence and inflow 28 31 ! information is not read from file. … … 996 999 DO j = nysg, nyng 997 1000 DO k = nzb, nzt+1 998 IF ( a11(k) .NE. 0._wp )THEN1001 IF ( a11(k) /= 0.0_wp ) THEN 999 1002 fu_yz(k,j) = ( u(k,j,i)  u_init(k) ) / a11(k) 1000 1003 ELSE 1001 fu_yz(k,j) = 0. _wp1004 fu_yz(k,j) = 0.0_wp 1002 1005 ENDIF 1003 1006 1004 IF ( a22(k) .NE. 0._wp )THEN1007 IF ( a22(k) /= 0.0_wp ) THEN 1005 1008 fv_yz(k,j) = ( v(k,j,i)  & 1006 1009 a21(k) * fu_yz(k,j)  v_init(k) ) / a22(k) 1007 1010 ELSE 1008 fv_yz(k,j) = 0. _wp1011 fv_yz(k,j) = 0.0_wp 1009 1012 ENDIF 1010 1013 1011 IF ( a33(k) .NE. 0._wp )THEN1014 IF ( a33(k) /= 0.0_wp ) THEN 1012 1015 fw_yz(k,j) = ( w(k,j,i)  & 1013 1016 a31(k) * fu_yz(k,j)  a32(k) * & … … 1029 1032 DO k = nzb, nzt+1 1030 1033 1031 IF ( a11(k) .NE. 0._wp )THEN1034 IF ( a11(k) /= 0.0_wp ) THEN 1032 1035 fu_xz(k,i) = ( u(k,j,i)  u_init(k) ) / a11(k) 1033 1036 ELSE 1034 fu_xz(k,i) = 0. _wp1037 fu_xz(k,i) = 0.0_wp 1035 1038 ENDIF 1036 1039 1037 IF ( a22(k) .NE. 0._wp )THEN1040 IF ( a22(k) /= 0.0_wp ) THEN 1038 1041 fv_xz(k,i) = ( v(k,j,i)  & 1039 1042 a21(k) * fu_xz(k,i)  v_init(k) ) / a22(k) 1040 1043 ELSE 1041 fv_xz(k,i) = 0. _wp1044 fv_xz(k,i) = 0.0_wp 1042 1045 ENDIF 1043 1046 1044 IF ( a33(k) .NE. 0._wp )THEN1047 IF ( a33(k) /= 0.0_wp ) THEN 1045 1048 fw_xz(k,i) = ( w(k,j,i)  & 1046 1049 a31(k) * fu_xz(k,i)  &
Note: See TracChangeset
for help on using the changeset viewer.