Ignore:
Timestamp:
Aug 6, 2019 9:11:47 AM (23 months ago)
Author:
raasch
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

    r4071 r4144  
    2525! -----------------
    2626! $Id$
     27! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
     28!
     29! 4071 2019-07-03 20:02:00Z suehring
    2730! Bugfix, initialize mean_inflow_profiles in case turbulence and inflow
    2831! information is not read from file.
     
    996999          DO  j = nysg, nyng
    9971000             DO  k = nzb, nzt+1
    998                 IF  ( a11(k) .NE. 0._wp ) THEN
     1001                IF  ( a11(k) /= 0.0_wp ) THEN
    9991002                   fu_yz(k,j) = ( u(k,j,i) - u_init(k) ) / a11(k)
    10001003                ELSE
    1001                    fu_yz(k,j) = 0._wp
     1004                   fu_yz(k,j) = 0.0_wp
    10021005                ENDIF
    10031006
    1004                 IF  ( a22(k) .NE. 0._wp ) THEN
     1007                IF  ( a22(k) /= 0.0_wp ) THEN
    10051008                   fv_yz(k,j) = ( v(k,j,i) -                                  &
    10061009                                  a21(k) * fu_yz(k,j) - v_init(k) ) / a22(k)
    10071010                ELSE
    1008                    fv_yz(k,j) = 0._wp
     1011                   fv_yz(k,j) = 0.0_wp
    10091012                ENDIF
    10101013
    1011                 IF  ( a33(k) .NE. 0._wp ) THEN
     1014                IF  ( a33(k) /= 0.0_wp ) THEN
    10121015                   fw_yz(k,j) = ( w(k,j,i) -                                   &
    10131016                                  a31(k) * fu_yz(k,j) - a32(k) *               &
     
    10291032             DO  k = nzb, nzt+1
    10301033
    1031                 IF  ( a11(k) .NE. 0._wp ) THEN
     1034                IF  ( a11(k) /= 0.0_wp ) THEN
    10321035                   fu_xz(k,i) = ( u(k,j,i) - u_init(k) ) / a11(k)
    10331036                ELSE
    1034                    fu_xz(k,i) = 0._wp
     1037                   fu_xz(k,i) = 0.0_wp
    10351038                ENDIF
    10361039
    1037                 IF  ( a22(k) .NE. 0._wp ) THEN
     1040                IF  ( a22(k) /= 0.0_wp ) THEN
    10381041                   fv_xz(k,i) = ( v(k,j,i) -                                   &
    10391042                                  a21(k) * fu_xz(k,i) - v_init(k) ) / a22(k)
    10401043                ELSE
    1041                    fv_xz(k,i) = 0._wp
     1044                   fv_xz(k,i) = 0.0_wp
    10421045                ENDIF
    10431046
    1044                 IF  ( a33(k) .NE. 0._wp ) THEN
     1047                IF  ( a33(k) /= 0.0_wp ) THEN
    10451048                   fw_xz(k,i) = ( w(k,j,i) -                                   &
    10461049                                  a31(k) * fu_xz(k,i) -                        &
Note: See TracChangeset for help on using the changeset viewer.