Changeset 366 for palm/trunk/SOURCE/diffusion_w.f90

Timestamp:

Aug 25, 2009 8:06:27 AM (15 years ago)

Author:

raasch

Message:

speed optomizations +bugfix in init_ocean

File:

• palm/trunk/SOURCE/diffusion_w.f90 (modified) (5 diffs)

palm/trunk/SOURCE/diffusion_w.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! bc_lr/bc_ns replaced by bc_lr_cyc/bc_ns_cyc
 !
 ! Former revisions:
@@ -100 +100 @@
 !--             velocity components parallel to the outflow boundary in
 !--             the direction normal to the outflow boundary.
-                IF ( bc_lr /= 'cyclic' )  THEN
+                IF ( .NOT. bc_lr_cyc )  THEN
                    kmxp_x = MAX( kmxp_x, km_damp_x(i) )
                    kmxm_x = MAX( kmxm_x, km_damp_x(i) )
                 ENDIF
-                IF ( bc_ns /= 'cyclic' )  THEN
+                IF ( .NOT. bc_ns_cyc )  THEN
                    kmyp_y = MAX( kmyp_y, km_damp_y(j) )
                    kmym_y = MAX( kmym_y, km_damp_y(j) )
@@ -150 +150 @@
 !--                velocity components parallel to the outflow boundary in
 !--                the direction normal to the outflow boundary.
-                   IF ( bc_lr /= 'cyclic' )  THEN
+                   IF ( .NOT. bc_lr_cyc )  THEN
                       kmxp_x = MAX( kmxp_x, km_damp_x(i) )
                       kmxm_x = MAX( kmxm_x, km_damp_x(i) )
                    ENDIF
-                   IF ( bc_ns /= 'cyclic' )  THEN
+                   IF ( .NOT. bc_ns_cyc )  THEN
                       kmyp_y = MAX( kmyp_y, km_damp_y(j) )
                       kmym_y = MAX( kmym_y, km_damp_y(j) )
@@ -231 +231 @@
 !--       parallel to the outflow boundary in the direction normal to the
 !--       outflow boundary.
-          IF ( bc_lr /= 'cyclic' )  THEN
+          IF ( .NOT. bc_lr_cyc )  THEN
              kmxp_x = MAX( kmxp_x, km_damp_x(i) )
              kmxm_x = MAX( kmxm_x, km_damp_x(i) )
           ENDIF
-          IF ( bc_ns /= 'cyclic' )  THEN
+          IF ( .NOT. bc_ns_cyc )  THEN
              kmyp_y = MAX( kmyp_y, km_damp_y(j) )
              kmym_y = MAX( kmym_y, km_damp_y(j) )
@@ -293 +293 @@
 !--          velocity components parallel to the outflow boundary in
 !--          the direction normal to the outflow boundary.
-             IF ( bc_lr /= 'cyclic' )  THEN
+             IF ( .NOT. bc_lr_cyc )  THEN
                 kmxp_x = MAX( kmxp_x, km_damp_x(i) )
                 kmxm_x = MAX( kmxm_x, km_damp_x(i) )
              ENDIF
-             IF ( bc_ns /= 'cyclic' )  THEN
+             IF ( .NOT. bc_ns_cyc )  THEN
                 kmyp_y = MAX( kmyp_y, km_damp_y(j) )
                 kmym_y = MAX( kmym_y, km_damp_y(j) )