Ignore:
Timestamp:
Mar 10, 2016 11:01:04 AM (6 years ago)
Author:
maronga
Message:

added support for water and paved surfaced in land surface model / minor changes

File:
1 edited

Legend:

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

    r1694 r1788  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Added z0q and z0q_av
    2222!
    2323! Former revisions:
     
    9292    USE arrays_3d,                                                             &
    9393        ONLY:  dzw, e, nr, ol, p, pt, q, qc, ql, ql_c, ql_v, qr, qsws, rho, sa,&
    94                shf, ts, u, us, v, vpt, w, z0, z0h
     94               shf, ts, u, us, v, vpt, w, z0, z0h, z0q
    9595
    9696    USE averaging,                                                             &
     
    9898               precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av,      &
    9999               ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av,  &
    100                shf_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av
     100               shf_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av,  &
     101               z0q_av
    101102
    102103    USE cloud_parameters,                                                      &
     
    494495                ENDIF
    495496                z0h_av = 0.0_wp
     497
     498             CASE ( 'z0q*' )
     499                IF ( .NOT. ALLOCATED( z0q_av ) )  THEN
     500                   ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
     501                ENDIF
     502                z0q_av = 0.0_wp
    496503
    497504             CASE DEFAULT
     
    10061013             ENDDO
    10071014
     1015          CASE ( 'z0q*' )
     1016             DO  i = nxlg, nxrg
     1017                DO  j = nysg, nyng
     1018                   z0q_av(j,i) = z0q_av(j,i) + z0q(j,i)
     1019                ENDDO
     1020             ENDDO
     1021
    10081022          CASE DEFAULT
    10091023!
Note: See TracChangeset for help on using the changeset viewer.