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

Timestamp:

Jan 11, 2018 12:01:27 PM (6 years ago)

Author:

suehring

Message:

Output of resistance also urban-type surfaces

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! output of r_a moved from land-surface to consider also urban-type surfaces
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -158 +161 @@
                ol_av, p_av, pc_av, pr_av, prr_av, precipitation_rate_av, pt_av,&
                q_av, qc_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qr_av, qsws_av, &
-               qv_av, rho_ocean_av, s_av, sa_av, shf_av, ssws_av, ts_av, u_av, &
-               us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
+               qv_av, r_a_av, rho_ocean_av, s_av, sa_av, shf_av, ssws_av,      &
+               ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
 #if defined( __chem )
     USE chemistry_model_mod,                                                   &
@@ -364 +367 @@
                 ENDIF
                 qv_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 ( 'rho_ocean' )
@@ -746 +755 @@
              ENDDO
 
+          CASE ( 'r_a*' )
+             DO  m = 1, surf_lsm_h%ns
+                i   = surf_lsm_h%i(m)            
+                j   = surf_lsm_h%j(m)
+                r_a_av(j,i) = r_a_av(j,i) + surf_lsm_h%r_a(m)
+             ENDDO
+!
+!--          Please note, resistance is also applied at urban-type surfaces, 
+!--          and is output only as a single variable. Here, tile approach is
+!--          already implemented, so for each surface fraction resistance 
+!--          need to be summed-up.
+             DO  m = 1, surf_usm_h%ns
+                i   = surf_usm_h%i(m)            
+                j   = surf_usm_h%j(m)
+                r_a_av(j,i) = r_a_av(j,i) +                                    &
+                           ( surf_usm_h%frac(0,m) * surf_usm_h%r_a(m)       +  & 
+                             surf_usm_h%frac(1,m) * surf_usm_h%r_a_green(m) +  & 
+                             surf_usm_h%frac(2,m) * surf_usm_h%r_a_window(m) )
+             ENDDO
+
           CASE ( 'rho_ocean' )
              DO  i = nxlg, nxrg