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

implementation of new bulk microphysics scheme

File:
1 edited

Legend:

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

    r2233 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! 2233 2017-05-30 18:08:54Z suehring
    2732!
    2833! 2232 2017-05-30 17:47:52Z suehring
     
    116121
    117122    USE arrays_3d,                                                             &
    118         ONLY:  e, kh, km, p, pt, q, ql, qc, nr, prr, precipitation_amount, qr, &
    119                s, sa, u, u_m_l, u_m_n, u_m_r, u_m_s, v, v_m_l, v_m_n, v_m_r,   &
    120                v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s
     123        ONLY:  e, kh, km, p, pt, q, ql, qc, nc, nr, prr, precipitation_amount, &
     124               qr, s, sa, u, u_m_l, u_m_n, u_m_r, u_m_s, v, v_m_l, v_m_n,      &
     125               v_m_r, v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s
    121126       
    122127    USE averaging
     
    124129    USE control_parameters,                                                    &
    125130        ONLY:  iran, humidity, passive_scalar, cloud_physics, cloud_droplets,  &
    126                microphysics_seifert, ocean, topography
     131               microphysics_morrison, microphysics_seifert, ocean, topography
    127132               
    128133    USE indices,                                                               &
     
    222227             IF ( ALLOCATED( qc_av ) )  THEN
    223228                WRITE ( 14 )  'qc_av               ';  WRITE ( 14 )  qc_av
     229             ENDIF
     230             IF ( microphysics_morrison )  THEN
     231                WRITE ( 14 )  'nc                  ';  WRITE ( 14 )  nc
     232                IF ( ALLOCATED( nc_av ) )  THEN
     233                   WRITE ( 14 )  'nc_av               ';  WRITE ( 14 )  nc_av
     234                ENDIF
    224235             ENDIF
    225236             IF ( microphysics_seifert )  THEN
Note: See TracChangeset for help on using the changeset viewer.