Ignore:
Timestamp:
Jun 20, 2017 9:51:42 AM (4 years ago)
Author:
schwenkel
Message:

implementation of new bulk microphysics scheme

File:
1 edited

Legend:

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

    r2270 r2292  
    2525! -----------------
    2626! $Id$
     27! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
     28! includes two more prognostic equations for cloud drop concentration (nc) 
     29! and cloud water content (qc).
     30!
     31! 2270 2017-06-09 12:18:47Z maronga
    2732! Revised numbering (removed 2 timeseries)
    2833!
     
    233238    USE arrays_3d,                                                             &
    234239        ONLY:  ddzu, ddzw, e, heatflux_output_conversion, hyp, km, kh,         &
    235                momentumflux_output_conversion, nr, p, prho, prr, pt, q,        &
     240               momentumflux_output_conversion, nc, nr, p, prho, prr, pt, q,    &
    236241               qc, ql, qr, rho_air, rho_air_zw, rho_ocean, s,                  &
    237242               sa, td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q, time_vert, u,   &
     
    245250                dt_3d, g, humidity, kappa, land_surface, large_scale_forcing,  &
    246251                large_scale_subsidence, max_pr_user, message_string, neutral,  &
    247                 microphysics_seifert, ocean, passive_scalar, simulated_time,   &
    248                 use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
    249                 ws_scheme_mom, ws_scheme_sca
     252                microphysics_morrison, microphysics_seifert, ocean,            &
     253                passive_scalar, simulated_time, use_subsidence_tendencies,     &
     254                use_surface_fluxes, use_top_fluxes, ws_scheme_mom,             &
     255                ws_scheme_sca
    250256       
    251257    USE cpulog,                                                                &
     
    12171223                                                             rmask(j,i,sr) *   &
    12181224                                                             flag
     1225                         IF ( microphysics_morrison )  THEN
     1226                            sums_l(k,123,tn) = sums_l(k,123,tn) + nc(k,j,i) *  &
     1227                                                                rmask(j,i,sr) *&
     1228                                                                flag
     1229                         ENDIF
    12191230                         IF ( microphysics_seifert )  THEN
    12201231                            sums_l(k,73,tn) = sums_l(k,73,tn) + nr(k,j,i) *    &
     
    17591770       hom(:,1,71,sr) = sums(:,71)     ! prho
    17601771       hom(:,1,72,sr) = hyp * 1E-2_wp  ! hyp in hPa
     1772       hom(:,1,123,sr) = sums(:,123)   ! nc
    17611773       hom(:,1,73,sr) = sums(:,73)     ! nr
    17621774       hom(:,1,74,sr) = sums(:,74)     ! qr
Note: See TracChangeset for help on using the changeset viewer.