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/read_3d_binary.f90

    r2269 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! 2269 2017-06-09 11:57:32Z suehring
    2732! Enable restart runs for urban_surface_mod
    2833!
     
    127132
    128133    USE arrays_3d,                                                             &
    129         ONLY:  e, kh, km, p, pt, q, ql, qc, nr, prr, precipitation_amount, qr, &
    130                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,   &
    131                v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s
     134        ONLY:  e, kh, km, p, pt, q, ql, qc, nc, nr, prr, precipitation_amount, &
     135               qr, s, sa, u, u_m_l, u_m_n, u_m_r, u_m_s, v, v_m_l, v_m_n,      &
     136               v_m_r, v_m_s, vpt, w, w_m_l, w_m_n, w_m_r, w_m_s
    132137
    133138    USE averaging
     
    496501                                  tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    497502
     503                CASE ( 'nc' )
     504                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     505                   nc(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     506                                   tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     507
     508                CASE ( 'nc_av' )
     509                   IF ( .NOT. ALLOCATED( nc_av ) )  THEN
     510                      ALLOCATE( nc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     511                   ENDIF
     512                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     513                   nc_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     514                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     515
     516
    498517                CASE ( 'nr' )
    499518                   IF ( k == 1 )  READ ( 13 )  tmp_3d
Note: See TracChangeset for help on using the changeset viewer.