Changeset 4144 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

Timestamp:

Aug 6, 2019 9:11:47 AM (5 years ago)

Author:

raasch

Message:

relational operators .EQ., .NE., etc. replaced by ==, /=, etc.

File:

• palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
+! 
+! 4071 2019-07-03 20:02:00Z suehring
 ! Bugfix, initialize mean_inflow_profiles in case turbulence and inflow
 ! information is not read from file. 
@@ -996 +999 @@
           DO  j = nysg, nyng
              DO  k = nzb, nzt+1
-                IF  ( a11(k) .NE. 0._wp ) THEN
+                IF  ( a11(k) /= 0.0_wp )  THEN
                    fu_yz(k,j) = ( u(k,j,i) - u_init(k) ) / a11(k)
                 ELSE
-                   fu_yz(k,j) = 0._wp
+                   fu_yz(k,j) = 0.0_wp
                 ENDIF
 
-                IF  ( a22(k) .NE. 0._wp ) THEN
+                IF  ( a22(k) /= 0.0_wp )  THEN
                    fv_yz(k,j) = ( v(k,j,i) -                                  &
                                   a21(k) * fu_yz(k,j) - v_init(k) ) / a22(k)
                 ELSE
-                   fv_yz(k,j) = 0._wp
+                   fv_yz(k,j) = 0.0_wp
                 ENDIF
 
-                IF  ( a33(k) .NE. 0._wp ) THEN
+                IF  ( a33(k) /= 0.0_wp )  THEN
                    fw_yz(k,j) = ( w(k,j,i) -                                   &
                                   a31(k) * fu_yz(k,j) - a32(k) *               &
@@ -1029 +1032 @@
              DO  k = nzb, nzt+1
 
-                IF  ( a11(k) .NE. 0._wp ) THEN
+                IF  ( a11(k) /= 0.0_wp )  THEN
                    fu_xz(k,i) = ( u(k,j,i) - u_init(k) ) / a11(k)
                 ELSE
-                   fu_xz(k,i) = 0._wp
+                   fu_xz(k,i) = 0.0_wp
                 ENDIF
 
-                IF  ( a22(k) .NE. 0._wp ) THEN
+                IF  ( a22(k) /= 0.0_wp )  THEN
                    fv_xz(k,i) = ( v(k,j,i) -                                   &
                                   a21(k) * fu_xz(k,i) - v_init(k) ) / a22(k)
                 ELSE
-                   fv_xz(k,i) = 0._wp
+                   fv_xz(k,i) = 0.0_wp
                 ENDIF
 
-                IF  ( a33(k) .NE. 0._wp ) THEN
+                IF  ( a33(k) /= 0.0_wp )  THEN
                    fw_xz(k,i) = ( w(k,j,i) -                                   &
                                   a31(k) * fu_xz(k,i) -                        &