Changeset 402 for palm

Timestamp:

Oct 21, 2009 11:59:41 AM (15 years ago)

Author:

maronga

Message:

bugfix: calculation of time-averaged surface-heatfluxes

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• sum_up_3d_data.f90 (modified) (6 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -14 +14 @@
 Bugfix: exchange of ghost points for prho included (time_integration)
 
-time_integration
+Bugfix: calculation of time-averaged surface heatfluxes (sum_up_3d_data)
+
+sum_up_3d_data, time_integration
 
 

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix in calculation of shf*_av, qsws*_av 
 !
 ! Former revisions:
@@ -321 +321 @@
                 ENDDO
              ENDDO
-             
+
           CASE ( 'ql' )
              DO  i = nxl-1, nxr+1
@@ -358 +358 @@
              ENDDO
 
+          CASE ( 'qsws*' )
+             DO  i = nxl-1, nxr+1
+                DO  j = nys-1, nyn+1
+                   qsws_av(j,i) = qsws_av(j,i) + qsws(j,i)
+                ENDDO
+             ENDDO
+
           CASE ( 'qv' )
              DO  i = nxl-1, nxr+1
@@ -375 +382 @@
                 ENDDO
              ENDDO
-             
+
           CASE ( 's' )
              DO  i = nxl-1, nxr+1
@@ -384 +391 @@
                 ENDDO
              ENDDO
-             
+
           CASE ( 'sa' )
              DO  i = nxl-1, nxr+1
@@ -393 +400 @@
                 ENDDO
              ENDDO
-             
+
+          CASE ( 'shf*' )
+             DO  i = nxl-1, nxr+1
+                DO  j = nys-1, nyn+1
+                   shf_av(j,i) = shf_av(j,i) + shf(j,i)
+                ENDDO
+             ENDDO
+
           CASE ( 't*' )
              DO  i = nxl-1, nxr+1