Changeset 4500 for palm/trunk/SOURCE/data_output_2d.f90

Timestamp:

Apr 17, 2020 10:12:45 AM (4 years ago)

Author:

suehring

Message:

Surface output: correct output of ground/wall-heat flux at USM surfaces; add conversion factor to heat- and momentum-flux outputs; data_output_2d: Unify output conversion of sensible and latent heat flux; data-output module: avoid uninitialized variables; restart_data_mpi_io: fix overlong lines

File:

• palm/trunk/SOURCE/data_output_2d.f90 (modified) (6 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Unify output conversion of sensible and latent heat flux
+! 
+! 4457 2020-03-11 14:20:43Z raasch
 ! use statement for exchange horiz added
 ! 
@@ -94 +97 @@
 
     USE basic_constants_and_equations_mod,                                     &
-        ONLY:  c_p, lv_d_cp, l_v
+        ONLY:  lv_d_cp
 
     USE bulk_cloud_model_mod,                                                  &
@@ -716 +719 @@
                       j = surf_lsm_h%j(m)
                       k = surf_lsm_h%k(m)
-                      local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m) * l_v
+                      local_pf(i,j,nzb+1) = surf_lsm_h%qsws(m) * waterflux_output_conversion(k)
                    ENDDO
                    DO  m = 1, surf_usm_h%ns
@@ -722 +725 @@
                       j = surf_usm_h%j(m)
                       k = surf_usm_h%k(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * l_v
+                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * waterflux_output_conversion(k)
                    ENDDO
                 ELSE
@@ -820 +823 @@
                       j = surf_lsm_h%j(m)
                       k = surf_lsm_h%k(m)
-                      local_pf(i,j,nzb+1) = surf_lsm_h%shf(m) * c_p
+                      local_pf(i,j,nzb+1) = surf_lsm_h%shf(m) * heatflux_output_conversion(k)
                    ENDDO
                    DO  m = 1, surf_usm_h%ns
@@ -826 +829 @@
                       j = surf_usm_h%j(m)
                       k = surf_usm_h%k(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%shf(m) * c_p
+                      local_pf(i,j,nzb+1) = surf_usm_h%shf(m) * heatflux_output_conversion(k)
                    ENDDO
                 ELSE