Changeset 2759 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Jan 17, 2018 4:24:59 PM (7 years ago)

Author:

suehring

Message:

Major bugfix in horizontal diffusion of all scalar quantities at vertical surfaces; further bugfix: density is considered for vertical fluxes of passive scalar and salinity

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (11 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Bugfix, consider density in vertical fluxes of passive scalar as well as 
+! chemical species.
+! 
+! 2753 2018-01-16 14:16:49Z suehring
 ! +r_a_green, r_a_window
 ! 
@@ -116 +120 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  zu, zw, heatflux_input_conversion, waterflux_input_conversion,  &
-               momentumflux_input_conversion
+        ONLY:  heatflux_input_conversion, momentumflux_input_conversion,       &
+               rho_air, rho_air_zw, zu, zw, waterflux_input_conversion 
 
 #if defined( __chem )
@@ -1753 +1757 @@
                       IF ( .NOT. random_heatflux )  THEN
                          surf%shf(num_h) = surface_heatflux *                  &
-                                                 heatflux_input_conversion(nzb)
+                                                 heatflux_input_conversion(k-1)
 !
 !--                      Check if surface heat flux might be replaced by 
@@ -1784 +1788 @@
                       IF ( constant_waterflux )  THEN
                          surf%qsws(num_h) = surface_waterflux *                &
-                                                 waterflux_input_conversion(nzb)
+                                                 waterflux_input_conversion(k-1)
                          IF ( k-1 /= 0 )  THEN
                             surf%qsws(num_h) = wall_humidityflux(0) *          &
@@ -1801 +1805 @@
                    IF ( upward_facing )  THEN
                       IF ( constant_scalarflux )  THEN
-                         surf%ssws(num_h) = surface_scalarflux
+                         surf%ssws(num_h) = surface_scalarflux * rho_air_zw(k-1)
 
                          IF ( k-1 /= 0 )                                       &
-                            surf%ssws(num_h) = wall_scalarflux(0)
+                            surf%ssws(num_h) = wall_scalarflux(0) *            &
+                                               rho_air_zw(k-1)
 
                       ELSE
@@ -1810 +1815 @@
                       ENDIF
                    ELSE
-                      surf%ssws(num_h) = wall_scalarflux(5)
+                      surf%ssws(num_h) = wall_scalarflux(5) * rho_air_zw(k)
                    ENDIF
                 ENDIF
@@ -1824 +1829 @@
                             IF ( upward_facing )  THEN
                                IF ( constant_csflux(lsp_pr) )  THEN
-                                  surf%cssws(lsp,num_h) = surface_csflux(lsp_pr)
-
-                                  IF ( k-1 /= 0 )                                    &
-                                     surf%cssws(lsp,num_h) = wall_csflux(lsp,0)
+                                  surf%cssws(lsp,num_h) =                      &
+                                                       surface_csflux(lsp_pr) *&
+                                                       rho_air_zw(k-1)
+
+                                  IF ( k-1 /= 0 )                              &
+                                     surf%cssws(lsp,num_h) =                   &
+                                                       wall_csflux(lsp,0) *    &
+                                                       rho_air_zw(k-1) 
                                ELSE
                                   surf%cssws(lsp,num_h) = 0.0_wp
                                ENDIF
                             ELSE
-                               surf%cssws(lsp,num_h) = wall_csflux(lsp,5)
+                               surf%cssws(lsp,num_h) = wall_csflux(lsp,5) *    &
+                                                       rho_air_zw(k)
                             ENDIF
                          ENDIF
@@ -1843 +1853 @@
                 IF ( ocean )  THEN
                    IF ( upward_facing )  THEN 
-                      surf%sasws(num_h) = bottom_salinityflux
+                      surf%sasws(num_h) = bottom_salinityflux * rho_air_zw(k-1)
                    ELSE
                       surf%sasws(num_h) = 0.0_wp
@@ -1907 +1917 @@
 !--          Prescribe top scalar flux
              IF ( passive_scalar .AND. constant_top_scalarflux )               &
-                surf%ssws(num_h) = top_scalarflux
+                surf%ssws(num_h) = top_scalarflux * rho_air_zw(nzt+1)
 !
 !--          Prescribe top chemical species' flux
@@ -1913 +1923 @@
              DO  lsp = 1, nvar
                 IF ( air_chemistry  .AND.  constant_top_csflux(lsp) )  THEN 
-                   surf%cssws(lsp,num_h) = top_csflux(lsp)
+                   surf%cssws(lsp,num_h) = top_csflux(lsp) * rho_air_zw(nzt+1)
                 ENDIF
              ENDDO
@@ -1920 +1930 @@
 !--          Prescribe top salinity flux
              IF ( ocean .AND. constant_top_salinityflux)                       &
-                surf%sasws(num_h) = top_salinityflux
+                surf%sasws(num_h) = top_salinityflux * rho_air_zw(nzt+1)
 !
 !--          Top momentum fluxes