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

Timestamp:

Oct 27, 2011 10:56:21 AM (13 years ago)

Author:

heinze

Message:

Output of liquid water potential temperature in case of cloud_physics=.T. enabled

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! lpt_av
 !
 ! Former revisions:
@@ -71 +72 @@
                 e_av = 0.0
 
+             CASE ( 'lpt' )
+                IF ( .NOT. ALLOCATED( lpt_av ) )  THEN
+                   ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                lpt_av = 0.0
+
              CASE ( 'lwp*' )
                 IF ( .NOT. ALLOCATED( lwp_av ) )  THEN
@@ -239 +246 @@
                    DO  k = nzb, nzt+1
                       e_av(k,j,i) = e_av(k,j,i) + e(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDDO
+
+          CASE ( 'lpt' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb, nzt+1
+                      lpt_av(k,j,i) = lpt_av(k,j,i) + pt(k,j,i)
                    ENDDO
                 ENDDO