Changeset 4328 for palm/trunk


Ignore:
Timestamp:
Dec 9, 2019 6:53:04 PM (5 years ago)
Author:
suehring
Message:

Minor change in nzb_max computation; commentation added; minor formatting adjustments in advec_ws

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4327 r4328  
    2525! -----------------
    2626! $Id$
     27! Minor formatting adjustments
     28!
     29! 4327 2019-12-06 14:48:31Z Giersch
    2730! Setting of advection flags for vertical fluxes of w revised, air density for
    2831! vertical flux calculation of w at k=1 is considered now
     
    20882091                                           )                                 
    20892092                                                                             
    2090               diss_l_u(k,j,tn) = - ABS( u_comp_l ) * (                         &
     2093             diss_l_u(k,j,tn) = - ABS( u_comp_l ) * (                          &
    20912094                              ( 10.0_wp * ibit2 * adv_mom_5                    &
    20922095                           +     3.0_wp * ibit1 * adv_mom_3                    &
     
    21012104                              ) *                                              &
    21022105                                        ( u(k,j,i+2) - u(k,j,i-3) )            &
    2103                                                      )
     2106                                                    )
    21042107
    21052108          ENDDO
     
    23872390                - ( u(k,j,i)   + u(k,j,i-1)   )                                &
    23882391                                    * (                                        &
    2389                      REAL( IBITS(advc_flags_m(k,j,i-1),0,1),  KIND = wp )      &
     2392                     REAL( IBITS(advc_flags_m(k,j,i-1),0,1), KIND = wp )       &
    23902393                   + REAL( IBITS(advc_flags_m(k,j,i-1),1,1), KIND = wp )       &
    23912394                   + REAL( IBITS(advc_flags_m(k,j,i-1),2,1), KIND = wp )       &
     
    25912594                                            )
    25922595
    2593               diss_l_v(k,j,tn) = - ABS( u_comp(k) ) * (                       &
     2596             diss_l_v(k,j,tn) = - ABS( u_comp(k) ) * (                        &
    25942597                              ( 10.0_wp * ibit11 * adv_mom_5                  &
    25952598                           +     3.0_wp * ibit10 * adv_mom_3                  &
     
    26042607                              ) *                                             &
    26052608                                        ( v(k,j,i+2) - v(k,j,i-3) )           &
    2606                                                       )
     2609                                                     )
    26072610
    26082611          ENDDO
     
    29352938          ibit11 = REAL( IBITS(advc_flags_m(k,j,i),11,1), KIND = wp )
    29362939          ibit10 = REAL( IBITS(advc_flags_m(k,j,i),10,1), KIND = wp )
    2937           ibit9  = REAL( IBITS(advc_flags_m(k,j,i),9,1), KIND = wp )
     2940          ibit9  = REAL( IBITS(advc_flags_m(k,j,i),9,1),  KIND = wp )
    29382941         
    29392942          ibit14 = REAL( IBITS(advc_flags_m(k,j,i),14,1), KIND = wp )
     
    32133216                                            )
    32143217
    3215                diss_l_w(k,j,tn) = - ABS( u_comp(k) ) * (                      &
     3218             diss_l_w(k,j,tn) = - ABS( u_comp(k) ) * (                        &
    32163219                             ( 10.0_wp * ibit20 * adv_mom_5                   &
    32173220                          +     3.0_wp * ibit19 * adv_mom_3                   &
     
    32263229                             ) *                                              &
    32273230                                        ( w(k,j,i+2) - w(k,j,i-3) )           &
    3228                                                        )
     3231                                                     )
    32293232
    32303233          ENDDO
  • palm/trunk/SOURCE/init_grid.f90

    r4314 r4328  
    2525! -----------------
    2626! $Id$
     27! Minor change in nzb_max computation. Commentation added.
     28!
     29! 4314 2019-11-29 10:29:20Z suehring
    2730! Set additional topography flag 4 to mark topography grid points emerged
    2831! from the filtering process.
     
    612615!-- Determine the maximum level of topography. It is used for
    613616!-- steering the degradation of order of the applied advection scheme,
    614 !-- as well in the lpm. 
     617!-- as well in the lpm.
    615618    k_top = 0
    616619    DO  i = nxl, nxr
     
    622625    ENDDO
    623626#if defined( __parallel )
    624     CALL MPI_ALLREDUCE( k_top + 1, nzb_max, 1, MPI_INTEGER,                    & !is +1 really necessary here?
     627    CALL MPI_ALLREDUCE( k_top, nzb_max, 1, MPI_INTEGER,                        &
    625628                        MPI_MAX, comm2d, ierr )
    626629#else
    627     nzb_max = k_top + 1
     630    nzb_max = k_top
    628631#endif
    629 !   
    630 !-- If topography extents up to the model top, limit nzb_max to nzt.
    631     nzb_max = MIN( nzb_max, nzt )
     632!
     633!-- Increment nzb_max by 1 in order to allow for proper diverengence correction.
     634!-- Further, in case topography extents up to the model top, limit to nzt.
     635    nzb_max = MIN( nzb_max+1, nzt )
    632636!
    633637!-- Determine minimum index of topography. Usually, this will be nzb. In case
Note: See TracChangeset for help on using the changeset viewer.