Changeset 2795 for palm/trunk


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

Bugfix in nesting anterpolation relaxation functions

Location:
palm/trunk
Files:
1 edited
1 copied

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
  • palm/trunk/UTIL/Makefile

    r2692 r2795  
    4949CFLAGS = -O
    5050
    51 F90 =
    52 F90_SER =
    53 COPT =
    54 F90FLAGS =
    55 LDFLAGS =
     51F90 = mpif90
     52F90_SER = gfortran
     53COPT = -cpp -DMPI_REAL=MPI_DOUBLE_PRECISION -DMPI_2REAL=MPI_2DOUBLE_PRECISION -D__netcdf
     54F90FLAGS = -fcheck=all -fdefault-real-8 -fconvert=little-endian
     55LDFLAGS  = -fcheck=all -fdefault-real-8 -fconvert=little-endian
     56
     57PKG_CONFIG      := PKG_CONFIG_ALLOW_SYSTEM_CFLAGS=1 pkg-config
     58NETCDF_F90FLAGS != $(PKG_CONFIG) --cflags netcdf-fortran
     59NETCDF_LDFLAGS  != $(PKG_CONFIG) --libs netcdf-fortran
     60F90FLAGS        += $(NETCDF_F90FLAGS)
     61LDFLAGS         += $(NETCDF_LDFLAGS)
    5662
    5763.SUFFIXES: $(SUFFIXES) .f90
Note: See TracChangeset for help on using the changeset viewer.