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

Timestamp:

Mar 4, 2015 5:44:27 PM (9 years ago)

Author:

maronga

Message:

LSM output of r_a and r_s added

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added output of r_a and r_s
 ! 
 ! Former revisions:
@@ -102 +102 @@
                m_soil_av, nzb_soil, nzt_soil, qsws_eb, qsws_eb_av,             &
                qsws_liq_eb, qsws_liq_eb_av, qsws_soil_eb, qsws_soil_eb_av,     &
-               qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, t_soil,         &
-               t_soil_av
+               qsws_veg_eb, qsws_veg_eb_av, shf_eb, shf_eb_av, r_a, r_a_av,    &
+               r_s, r_s_av, t_soil, t_soil_av
 
     USE particle_attributes,                                                   &
@@ -334 +334 @@
                 ENDIF
                 rho_av = 0.0_wp
+
+             CASE ( 'r_a*' )
+                IF ( .NOT. ALLOCATED( r_a_av ) )  THEN
+                   ALLOCATE( r_a_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                r_a_av = 0.0_wp
+
+             CASE ( 'r_s*' )
+                IF ( .NOT. ALLOCATED( r_s_av ) )  THEN
+                   ALLOCATE( r_s_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                r_s_av = 0.0_wp
 
              CASE ( 's' )
@@ -724 +736 @@
              ENDDO
 
+          CASE ( 'r_a*' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   r_a_av(j,i) = r_a_av(j,i) + r_a(j,i)
+                ENDDO
+             ENDDO
+
+          CASE ( 'r_s*' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   r_s_av(j,i) = r_s_av(j,i) + r_s(j,i)
+                ENDDO
+             ENDDO
+
           CASE ( 'rho' )
              DO  i = nxlg, nxrg