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

Timestamp:

Jan 12, 2018 4:03:39 PM (7 years ago)

Author:

suehring

Message:

In case of natural- and urban-type surfaces output surfaces fluxes in W/m2

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
+! 
+! 2742 2018-01-12 14:59:47Z suehring
 ! Enable output of surface temperature
 ! 
@@ -157 +160 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  dzw, e, nc, nr, p, pt, precipitation_rate, q, qc, ql, ql_c,     &
-               ql_v, qr, rho_ocean, s, sa, u, v, vpt, w
+        ONLY:  dzw, e, heatflux_output_conversion, nc, nr, p, pt,              &
+               precipitation_rate, q, qc, ql, ql_c,                            &
+               ql_v, qr, rho_ocean, s, sa, u, v, vpt, w,                       &
+               waterflux_output_conversion
 
     USE averaging,                                                             &
@@ -173 +178 @@
 
     USE cloud_parameters,                                                      &
-        ONLY:  l_d_cp, pt_d_t
+        ONLY:  cp, l_d_cp, l_v, pt_d_t
 
     USE control_parameters,                                                    &
@@ -740 +745 @@
 
           CASE ( 'qsws*' )
+!
+!--          In case of default surfaces, clean-up flux by density.
+!--          In case of land- and urban-surfaces, convert fluxes into
+!--          dynamic units.
              DO  m = 1, surf_def_h(0)%ns
                 i = surf_def_h(0)%i(m)
                 j = surf_def_h(0)%j(m)
-                qsws_av(j,i) = qsws_av(j,i) + surf_def_h(0)%qsws(m)
+                qsws_av(j,i) = qsws_av(j,i) + surf_def_h(0)%qsws(m) *          &
+                                              waterflux_output_conversion(k)
              ENDDO
              DO  m = 1, surf_lsm_h%ns
                 i = surf_lsm_h%i(m)
                 j = surf_lsm_h%j(m)
-                qsws_av(j,i) = qsws_av(j,i) + surf_lsm_h%qsws(m)
+                qsws_av(j,i) = qsws_av(j,i) + surf_lsm_h%qsws(m) * l_v
              ENDDO
              DO  m = 1, surf_usm_h%ns
                 i = surf_usm_h%i(m)
                 j = surf_usm_h%j(m)
-                qsws_av(j,i) = qsws_av(j,i) + surf_usm_h%qsws(m)
+                qsws_av(j,i) = qsws_av(j,i) + surf_usm_h%qsws(m) * l_v
              ENDDO
 
@@ -813 +823 @@
 
           CASE ( 'shf*' )
+!
+!--          In case of default surfaces, clean-up flux by density.
+!--          In case of land- and urban-surfaces, convert fluxes into
+!--          dynamic units.
              DO  m = 1, surf_def_h(0)%ns
                 i = surf_def_h(0)%i(m)
                 j = surf_def_h(0)%j(m)
-                shf_av(j,i) = shf_av(j,i) + surf_def_h(0)%shf(m)
+                shf_av(j,i) = shf_av(j,i) + surf_def_h(0)%shf(m)  *            &
+                                            heatflux_output_conversion(k)
              ENDDO
              DO  m = 1, surf_lsm_h%ns
                 i = surf_lsm_h%i(m)
                 j = surf_lsm_h%j(m)
-                shf_av(j,i) = shf_av(j,i) + surf_lsm_h%shf(m)
+                shf_av(j,i) = shf_av(j,i) + surf_lsm_h%shf(m) * cp
              ENDDO
              DO  m = 1, surf_usm_h%ns
                 i = surf_usm_h%i(m)
                 j = surf_usm_h%j(m)
-                shf_av(j,i) = shf_av(j,i) + surf_usm_h%shf(m)
+                shf_av(j,i) = shf_av(j,i) + surf_usm_h%shf(m) * cp
              ENDDO