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

Timestamp:

Jul 26, 2018 5:12:48 PM (6 years ago)

Author:

suehring

Message:

Major bugfix in calculation of ol and ts at building roofs; bugfix in restart data for surface elements; in 2D data output, mask latent heat flux at urban-type surfaces

File:

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

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove output of latent heat flux at urban-surfaces and set fill values 
+! instead 
+! 
+! 3052 2018-05-31 06:11:20Z raasch
 ! Do not open FORTRAN binary files in case of parallel netCDF I/O
 ! 
@@ -933 +937 @@
              CASE ( 'qsws*_xy' )        ! 2d-array
                 IF ( av == 0 ) THEN
+                   local_pf(:,:,nzb+1) = REAL( fill_value, KIND = wp )
 !
 !--                In case of default surfaces, clean-up flux by density.
-!--                In case of land- and urban-surfaces, convert fluxes into
+!--                In case of land-surfaces, convert fluxes into
 !--                dynamic units
                    DO  m = 1, surf_def_h(0)%ns
@@ -949 +954 @@
                       k = surf_lsm_h%k(m)
                       local_pf(i,j,nzb+1) = 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)
-                      k = surf_usm_h%k(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%qsws(m) * l_v
                    ENDDO
                 ELSE