Ignore:
Timestamp:
Jun 29, 2020 8:49:58 AM (12 months ago)
Author:
suehring
Message:

mesoscale nesting: omit explicit pressure forcing via geostrophic wind components; chemistry: enable profile output of vertical fluxes; urban-surface: bugfix in initialization in case of cyclic_fill

File:
1 edited

Legend:

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

    r4509 r4581  
    2424! -----------------
    2525! $Id$
     26! Enable output of resolved-scale vertical fluxes of chemical species.
     27!
     28! 4509 2020-04-26 15:57:55Z raasch
    2629! file re-formatted to follow the PALM coding standard
    2730!
     
    229232
    230233    USE statistics,                                                                                &
    231         ONLY:  hom, weight_substep,                                                                &
    232                sums_salsa_ws_l, sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l,                       &
    233                sums_wsncs_ws_l, sums_wsnrs_ws_l,                                                   &
    234                sums_wspts_ws_l, sums_wsqs_ws_l, sums_wsss_ws_l,                                    &
    235                sums_wsqcs_ws_l, sums_wsqrs_ws_l,                                                   &
    236                sums_wsqis_ws_l, sums_wsnis_ws_l,                                                   &
    237                sums_wssas_ws_l, sums_wsus_ws_l, sums_wsvs_ws_l
     234        ONLY:  hom,                                                                                &
     235               sums_salsa_ws_l,                                                                    &
     236               sums_us2_ws_l,                                                                      &
     237               sums_vs2_ws_l,                                                                      &
     238               sums_ws2_ws_l,                                                                      &
     239               sums_wschs_ws_l,                                                                    &
     240               sums_wsncs_ws_l,                                                                    &
     241               sums_wsnrs_ws_l,                                                                    &
     242               sums_wspts_ws_l,                                                                    &
     243               sums_wsqs_ws_l,                                                                     &
     244               sums_wsss_ws_l,                                                                     &
     245               sums_wsqcs_ws_l,                                                                    &
     246               sums_wsqrs_ws_l,                                                                    &
     247               sums_wsqis_ws_l,                                                                    &
     248               sums_wsnis_ws_l,                                                                    &
     249               sums_wssas_ws_l,                                                                    &
     250               sums_wsus_ws_l,                                                                     &
     251               sums_wsvs_ws_l,                                                                     &
     252               weight_substep
    238253
    239254
    240255
    241256    IMPLICIT NONE
    242    
     257
    243258    REAL(wp) ::  adv_mom_1            !< 1/4 - constant used in 5th-order advection scheme for momentum advection (1st-order part)
    244259    REAL(wp) ::  adv_mom_3            !< 1/24 - constant used in 5th-order advection scheme for momentum advection (3rd-order part)
     
    18751890          ENDDO
    18761891
    1877       CASE ( 'aerosol_mass', 'aerosol_number', 'salsa_gas' )
     1892       CASE ( 'aerosol_mass', 'aerosol_number', 'salsa_gas' )
    18781893
    18791894          DO  k = nzb, nzt
     
    18861901          ENDDO
    18871902
    1888 !       CASE ( 'kc' )
    1889         !kk Has to be implemented for kpp chemistry
     1903       CASE ( 'kc' )
     1904          DO  k = nzb, nzt
     1905             sums_wschs_ws_l(k,tn)  = sums_wschs_ws_l(k,tn) +                                      &
     1906                                      ( flux_t(k) / ( w(k,j,i) + SIGN( 1.0E-20_wp, w(k,j,i) ) )    &
     1907                                                  * ( w(k,j,i) - hom(k,1,3,0)                 )    &
     1908                                      + diss_t(k) / ( ABS(w(k,j,i)) + 1.0E-20_wp              )    &
     1909                                                  *   ABS( w(k,j,i) - hom(k,1,3,0)            )    &
     1910                                      ) * weight_substep(intermediate_timestep_count)
     1911          ENDDO
    18901912
    18911913    END SELECT
Note: See TracChangeset for help on using the changeset viewer.