Changeset 1053 for palm/trunk/SOURCE/average_3d_data.f90

Timestamp:

Nov 13, 2012 5:11:03 PM (12 years ago)

Author:

hoffmann

Message:

two-moment cloud physics implemented

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/average_3d_data.f90 (modified) (3 diffs)

palm/trunk/SOURCE/average_3d_data.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! averaging of nr, qr added
 !
 ! Former revisions:
@@ -117 +117 @@
              ENDDO
 
+          CASE ( 'nr' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb, nzt+1
+                      nr_av(k,j,i) = nr_av(k,j,i) / REAL( average_count_3d )
+                   ENDDO
+                ENDDO
+             ENDDO
+
           CASE ( 'p' )
              DO  i = nxlg, nxrg
@@ -207 +216 @@
              ENDDO
 
+          CASE ( 'qr' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb, nzt+1
+                      qr_av(k,j,i) = qr_av(k,j,i) / REAL( average_count_3d )
+                   ENDDO
+                ENDDO
+             ENDDO
+
           CASE ( 'qv' )
              DO  i = nxlg, nxrg