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/sum_up_3d_data.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
     
    135140
    136141    USE arrays_3d,                                                             &
    137         ONLY:  dzw, e, nr, p, pt, precipitation_rate, q, qc, ql, ql_c, ql_v,   &
    138                qr, rho_ocean, s, sa, u, v, vpt, w
     142        ONLY:  dzw, e, nc, nr, p, pt, precipitation_rate, q, qc, ql, ql_c,     &
     143               ql_v, qr, rho_ocean, s, sa, u, v, vpt, w
    139144
    140145    USE averaging,                                                             &
    141         ONLY:  e_av, lpt_av, lwp_av, nr_av, ol_av, p_av, pc_av, pr_av, prr_av, &
    142                precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av,      &
    143                ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_ocean_av, s_av, sa_av,  &
    144                shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, &
    145                z0h_av, z0q_av
     146        ONLY:  e_av, lpt_av, lwp_av, nc_av, nr_av, ol_av, p_av, pc_av, pr_av, &
     147               prr_av, precipitation_rate_av, pt_av, q_av, qc_av, ql_av,       &
     148               ql_c_av, ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_ocean_av,&
     149               s_av, sa_av, shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, &
     150               w_av, z0_av, z0h_av, z0q_av
    146151
    147152    USE cloud_parameters,                                                      &
     
    230235                lwp_av = 0.0_wp
    231236
     237             CASE ( 'nc' )
     238                IF ( .NOT. ALLOCATED( nc_av ) )  THEN
     239                   ALLOCATE( nc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     240                ENDIF
     241                nc_av = 0.0_wp
     242
    232243             CASE ( 'nr' )
    233244                IF ( .NOT. ALLOCATED( nr_av ) )  THEN
     
    483494                   lwp_av(j,i) = lwp_av(j,i) + SUM( ql(nzb:nzt,j,i)            &
    484495                                               * dzw(1:nzt+1) ) * rho_surface
     496                ENDDO
     497             ENDDO
     498
     499          CASE ( 'nc' )
     500             DO  i = nxlg, nxrg
     501                DO  j = nysg, nyng
     502                   DO  k = nzb, nzt+1
     503                      nc_av(k,j,i) = nc_av(k,j,i) + nc(k,j,i)
     504                   ENDDO
    485505                ENDDO
    486506             ENDDO
Note: See TracChangeset for help on using the changeset viewer.