Changeset 4144 for palm/trunk/SOURCE/biometeorology_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/biometeorology_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
+! 
+! 4127 2019-07-30 14:47:10Z suehring
 ! Output for bio_mrt added (merge from branch resler)
 ! 
@@ -4319 +4322 @@
 !
 !-- calculation of solar azimuth angle
-    IF (woz .LE. 12.0_wp) alpha = pi - acos( ( sin(thetasr) * sin( lat * dtor ) -                                     &
+    IF (woz <= 12.0_wp) alpha = pi - acos( ( sin(thetasr) * sin( lat * dtor ) -                                     &
     sin( declination * dtor ) ) / ( cos(thetasr) * cos( lat * dtor ) ) )    
-    IF (woz .GT. 12.0_wp) alpha = pi + acos( ( sin(thetasr) * sin( lat * dtor ) -                                     &
+    IF (woz > 12.0_wp)  alpha = pi + acos( ( sin(thetasr) * sin( lat * dtor ) -                                     &
     sin( declination * dtor ) ) / ( cos(thetasr) * cos( lat * dtor ) ) )    
     saa = alpha / dtor
@@ -4347 +4350 @@
     CALL uvem_solar_position
      
-    IF (sza  .GE.  90)  THEN
+    IF (sza  >=  90)  THEN
        vitd3_exposure(:,:) = 0.0_wp
     ELSE
@@ -4446 +4449 @@
              IF (consider_obstructions )  THEN
                 obstruction_temp1 = building_obstruction_f%var_3d(:,j,i)
-                IF (obstruction_temp1(0)  .NE.  9)  THEN 
+                IF ( obstruction_temp1(0)  /=  9 )  THEN
                    DO  pobi = 0, 44 
                       DO  bi = 0, 7 
@@ -4468 +4471 @@
          
              obstruction_direct_beam = obstruction( nint(startpos_saa_float), nint( sza / 10.0_wp ) ) 
-             IF (sza  .GE.  89.99_wp)  THEN
+             IF (sza  >=  89.99_wp)  THEN
                 sza = 89.99999_wp
              ENDIF