Ignore:
Timestamp:
Feb 7, 2018 2:48:48 PM (7 years ago)
Author:
hellstea
Message:

Bugfix in nesting anterpolation relaxation functions

File:
1 edited

Legend:

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

    r2773 r2795  
    2525! -----------------
    2626! $Id$
     27! Bugfix in computation of the anterpolation under-relaxation functions.
     28!
     29! 2773 2018-01-30 14:12:54Z suehring
    2730! - Nesting for chemical species
    2831! - Bugfix in setting boundary condition at downward-facing walls for passive
     
    26322635       IF ( nesting_mode /= 'vertical' )  THEN
    26332636          DO  ii = icl, icr
    2634              IF ( nest_bound_l )  THEN
     2637             IF ( ifuu(ii) < ( nx + 1 ) / 2 )  THEN   
    26352638                xi = ( MAX( 0.0_wp, ( cg%coord_x(ii) -                          &
    2636                        lower_left_coord_x ) ) / anterp_relax_length_l )**4
    2637              ELSEIF ( nest_bound_r )  THEN
     2639                     lower_left_coord_x ) ) / anterp_relax_length_l )**4
     2640                frax(ii) = xi / ( 1.0_wp + xi )
     2641             ELSE
    26382642                xi = ( MAX( 0.0_wp, ( lower_left_coord_x + ( nx + 1 ) * dx -    &
    26392643                                      cg%coord_x(ii) ) ) /                      &
    26402644                       anterp_relax_length_r )**4
     2645                frax(ii) = xi / ( 1.0_wp + xi )               
     2646             ENDIF
     2647          ENDDO
     2648
     2649
     2650          DO  jj = jcs, jcn
     2651             IF ( jfuv(jj) < ( ny + 1 ) / 2 )  THEN
     2652                eta = ( MAX( 0.0_wp, ( cg%coord_y(jj) -                         &
     2653                     lower_left_coord_y ) ) / anterp_relax_length_s )**4
     2654                fray(jj) = eta / ( 1.0_wp + eta )
    26412655             ELSE
    2642                 xi = 999999.9_wp
    2643              ENDIF
    2644              frax(ii) = xi / ( 1.0_wp + xi )
    2645           ENDDO
    2646 
    2647 
    2648           DO  jj = jcs, jcn
    2649              IF ( nest_bound_s )  THEN
    2650                 eta = ( MAX( 0.0_wp, ( cg%coord_y(jj) -                         &
    2651                         lower_left_coord_y ) ) / anterp_relax_length_s )**4
    2652              ELSEIF ( nest_bound_n )  THEN
    26532656                eta = ( MAX( 0.0_wp, ( lower_left_coord_y + ( ny + 1 ) * dy -   &
    26542657                                       cg%coord_y(jj)) ) /                      &
    26552658                        anterp_relax_length_n )**4
    2656              ELSE
    2657                 eta = 999999.9_wp
     2659                fray(jj) = eta / ( 1.0_wp + eta )
    26582660             ENDIF
    2659              fray(jj) = eta / ( 1.0_wp + eta )
    26602661          ENDDO
    26612662       ENDIF
Note: See TracChangeset for help on using the changeset viewer.