Changeset 1788 for palm/trunk/SOURCE/sum_up_3d_data.f90

Timestamp:

Mar 10, 2016 11:01:04 AM (8 years ago)

Author:

maronga

Message:

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

File:

• palm/trunk/SOURCE/sum_up_3d_data.f90 (modified) (5 diffs)

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added z0q and z0q_av
 ! 
 ! Former revisions:
@@ -92 +92 @@
     USE arrays_3d,                                                             &
         ONLY:  dzw, e, nr, ol, p, pt, q, qc, ql, ql_c, ql_v, qr, qsws, rho, sa,&
-               shf, ts, u, us, v, vpt, w, z0, z0h
+               shf, ts, u, us, v, vpt, w, z0, z0h, z0q
 
     USE averaging,                                                             &
@@ -98 +98 @@
                precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av,      &
                ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av,  &
-               shf_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av
+               shf_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av,  &
+               z0q_av
 
     USE cloud_parameters,                                                      &
@@ -494 +495 @@
                 ENDIF
                 z0h_av = 0.0_wp
+
+             CASE ( 'z0q*' )
+                IF ( .NOT. ALLOCATED( z0q_av ) )  THEN
+                   ALLOCATE( z0q_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                z0q_av = 0.0_wp
 
              CASE DEFAULT
@@ -1006 +1013 @@
              ENDDO
 
+          CASE ( 'z0q*' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   z0q_av(j,i) = z0q_av(j,i) + z0q(j,i)
+                ENDDO
+             ENDDO
+
           CASE DEFAULT
 !