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

Timestamp:

Feb 8, 2018 1:24:35 PM (6 years ago)

Author:

suehring

Message:

Output of ground-heat flux at natural- and urban-type surfaces in one output variable; enable restart data of _av variables that belong to both land- and urban-surface model

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable output of ground-heat flux also at urban surfaces.
+! 
+! 2790 2018-02-06 11:57:19Z suehring
 ! Bugfix in summation of surface sensible and latent heat flux
 ! 
@@ -172 +175 @@
 
     USE averaging,                                                             &
-        ONLY:  diss_av, e_av, kh_av, km_av, lpt_av, lwp_av, nc_av, nr_av,      &
+        ONLY:  diss_av, e_av, ghf_av, kh_av, km_av, lpt_av, lwp_av, nc_av,     &
+               nr_av,                                                          &
                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, &
@@ -255 +259 @@
           SELECT CASE ( trimvar )
 
+             CASE ( 'ghf*' )
+                IF ( .NOT. ALLOCATED( ghf_av ) )  THEN
+                   ALLOCATE( ghf_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                ghf_av = 0.0_wp
+
              CASE ( 'e' )
                 IF ( .NOT. ALLOCATED( e_av ) )  THEN
@@ -538 +548 @@
 !--    Store the array chosen on the temporary array.
        SELECT CASE ( trimvar )
+
+          CASE ( 'ghf*' )
+             DO  m = 1, surf_lsm_h%ns
+                i   = surf_lsm_h%i(m)            
+                j   = surf_lsm_h%j(m)
+                ghf_av(j,i) = ghf_av(j,i) + surf_lsm_h%ghf(m)
+             ENDDO
+
+             DO  m = 1, surf_usm_h%ns
+                i   = surf_usm_h%i(m)            
+                j   = surf_usm_h%j(m)
+                ghf_av(j,i) = ghf_av(j,i) + surf_usm_h%frac(0,m)     *          &
+                                            surf_usm_h%wghf_eb(m)        +      &
+                                            surf_usm_h%frac(1,m)     *          &
+                                            surf_usm_h%wghf_eb_green(m)  +      &
+                                            surf_usm_h%frac(2,m)     *          &
+                                            surf_usm_h%wghf_eb_window(m)
+             ENDDO
 
           CASE ( 'e' )