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

Timestamp:

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

Author:

raasch

Message:

preliminary changes for precipitation output

File:

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

palm/trunk/SOURCE/data_output_2d.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Output of precipitation amount/rate and roughness length
 !
 ! Former revisions:
@@ -257 +257 @@
                 ENDIF
 
+             CASE ( 'pra*_xy' )        ! 2d-array / integral quantity => no av
+                CALL exchange_horiz_2d( precipitation_amount )
+                DO  i = nxl-1, nxr+1
+                   DO  j = nys-1, nyn+1
+                      local_pf(i,j,nzb+1) =  precipitation_amount(j,i)
+                   ENDDO
+                ENDDO
+                precipitation_amount = 0.0   ! reset for next integ. interval
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
+
+             CASE ( 'prr*_xy' )        ! 2d-array
+                IF ( av == 0 )  THEN
+                   CALL exchange_horiz_2d( precipitation_rate )
+                   DO  i = nxl-1, nxr+1
+                      DO  j = nys-1, nyn+1 
+                         local_pf(i,j,nzb+1) =  precipitation_rate(j,i)
+                      ENDDO
+                   ENDDO
+                ELSE
+                   CALL exchange_horiz_2d( precipitation_rate_av )
+                   DO  i = nxl-1, nxr+1
+                      DO  j = nys-1, nyn+1 
+                         local_pf(i,j,nzb+1) =  precipitation_rate_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
+
              CASE ( 'pt_xy', 'pt_xz', 'pt_yz' )
                 IF ( av == 0 )  THEN
@@ -419 +451 @@
                 ENDIF
                 IF ( mode == 'xy' )  level_z = zw
+
+             CASE ( 'z0*_xy' )        ! 2d-array
+                IF ( av == 0 ) THEN
+                   DO  i = nxl-1, nxr+1
+                      DO  j = nys-1, nyn+1 
+                         local_pf(i,j,nzb+1) =  z0(j,i)
+                      ENDDO
+                   ENDDO
+                ELSE
+                   DO  i = nxl-1, nxr+1
+                      DO  j = nys-1, nyn+1 
+                         local_pf(i,j,nzb+1) =  z0_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
 
              CASE DEFAULT