Changeset 2795 for palm/trunk/SOURCE

Timestamp:

Feb 7, 2018 2:48:48 PM (6 years ago)

Author:

hellstea

Message:

Bugfix in nesting anterpolation relaxation functions

File:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in computation of the anterpolation under-relaxation functions. 
+! 
+! 2773 2018-01-30 14:12:54Z suehring
 ! - Nesting for chemical species
 ! - Bugfix in setting boundary condition at downward-facing walls for passive
@@ -2632 +2635 @@
        IF ( nesting_mode /= 'vertical' )  THEN
           DO  ii = icl, icr
-             IF ( nest_bound_l )  THEN
+             IF ( ifuu(ii) < ( nx + 1 ) / 2 )  THEN   
                 xi = ( MAX( 0.0_wp, ( cg%coord_x(ii) -                          &
-                       lower_left_coord_x ) ) / anterp_relax_length_l )**4
-             ELSEIF ( nest_bound_r )  THEN
+                     lower_left_coord_x ) ) / anterp_relax_length_l )**4
+                frax(ii) = xi / ( 1.0_wp + xi )
+             ELSE
                 xi = ( MAX( 0.0_wp, ( lower_left_coord_x + ( nx + 1 ) * dx -    &
                                       cg%coord_x(ii) ) ) /                      &
                        anterp_relax_length_r )**4
+                frax(ii) = xi / ( 1.0_wp + xi )                
+             ENDIF
+          ENDDO
+
+
+          DO  jj = jcs, jcn
+             IF ( jfuv(jj) < ( ny + 1 ) / 2 )  THEN 
+                eta = ( MAX( 0.0_wp, ( cg%coord_y(jj) -                         &
+                     lower_left_coord_y ) ) / anterp_relax_length_s )**4
+                fray(jj) = eta / ( 1.0_wp + eta )
              ELSE
-                xi = 999999.9_wp
-             ENDIF
-             frax(ii) = xi / ( 1.0_wp + xi )
-          ENDDO
-
-
-          DO  jj = jcs, jcn
-             IF ( nest_bound_s )  THEN
-                eta = ( MAX( 0.0_wp, ( cg%coord_y(jj) -                         &
-                        lower_left_coord_y ) ) / anterp_relax_length_s )**4
-             ELSEIF ( nest_bound_n )  THEN
                 eta = ( MAX( 0.0_wp, ( lower_left_coord_y + ( ny + 1 ) * dy -   &
                                        cg%coord_y(jj)) ) /                      &
                         anterp_relax_length_n )**4
-             ELSE
-                eta = 999999.9_wp
+                fray(jj) = eta / ( 1.0_wp + eta )
              ENDIF
-             fray(jj) = eta / ( 1.0_wp + eta )
           ENDDO
        ENDIF