Ignore:
Timestamp:
Jul 12, 2016 4:34:24 PM (5 years ago)
Author:
suehring
Message:

Separate balance equations for humidity and passive_scalar

File:
1 edited

Legend:

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

    r1950 r1960  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Scalar surface flux added
    2222!
    2323! Former revisions:
     
    101101    USE arrays_3d,                                                             &
    102102        ONLY:  dzw, e, nr, ol, p, pt, precipitation_rate, q, qc, ql, ql_c,     &
    103                ql_v, qr, qsws, rho, sa, shf, ts, u, us, v, vpt, w, z0, z0h, z0q
     103               ql_v, qr, qsws, rho, sa, shf, ssws, ts, u, us, v, vpt, w, z0,   &
     104               z0h, z0q
    104105
    105106    USE averaging,                                                             &
     
    107108               precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av,      &
    108109               ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av,  &
    109                shf_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av, &
    110                z0q_av
     110               shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, &
     111               z0h_av, z0q_av
    111112
    112113    USE cloud_parameters,                                                      &
     
    949950             ENDDO
    950951
     952          CASE ( 'ssws*' )
     953             DO  i = nxlg, nxrg
     954                DO  j = nysg, nyng
     955                   ssws_av(j,i) = ssws_av(j,i) + ssws(j,i)
     956                ENDDO
     957             ENDDO
     958
    951959          CASE ( 't*' )
    952960             DO  i = nxlg, nxrg
Note: See TracChangeset for help on using the changeset viewer.