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

Timestamp:

Mar 19, 2007 8:20:46 AM (17 years ago)

Author:

raasch

Message:

preliminary changes for precipitation output

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! +sum-up of precipitation rate and roughness length (prr*, z0*)
 !
 ! Former revisions:
@@ -80 +80 @@
                 pr_av = 0.0
 
+             CASE ( 'prr*' )
+                IF ( .NOT. ALLOCATED( precipitation_rate_av ) )  THEN
+                   ALLOCATE( precipitation_rate_av(nys-1:nyn+1,nxl-1:nxr+1) )
+                ENDIF
+                precipitation_rate_av = 0.0
+
              CASE ( 'pt' )
                 IF ( .NOT. ALLOCATED( pt_av ) )  THEN
@@ -163 +169 @@
                 ENDIF
                 w_av = 0.0
+
+             CASE ( 'z0*' )
+                IF ( .NOT. ALLOCATED( z0_av ) )  THEN
+                   ALLOCATE( z0_av(nys-1:nyn+1,nxl-1:nxr+1) )
+                ENDIF
+                z0_av = 0.0
 
              CASE DEFAULT
@@ -239 +251 @@
              ENDDO
 
+          CASE ( 'pr*' )
+             DO  i = nxl-1, nxr+1
+                DO  j = nys-1, nyn+1
+                   precipitation_rate_av(j,i) = precipitation_rate_av(j,i) + &
+                                                precipitation_rate(j,i)
+                ENDDO
+             ENDDO
+
           CASE ( 'pt' )
              IF ( .NOT. cloud_physics ) THEN
@@ -372 +392 @@
              ENDDO
 
+          CASE ( 'z0*' )
+             DO  i = nxl-1, nxr+1
+                DO  j = nys-1, nyn+1
+                   z0_av(j,i) = z0_av(j,i) + z0(j,i)
+                ENDDO
+             ENDDO
+
           CASE DEFAULT
 !