Changeset 4519 for palm/trunk

Timestamp:

May 5, 2020 5:33:30 PM (5 years ago)

Author:

suehring

Message:

Add missing computation of passive scalar scaling parameter

File:

• palm/trunk/SOURCE/surface_layer_fluxes_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Add missing computation of passive scalar scaling parameter
+! 
+! 4370 2020-01-10 14:00:44Z raasch
 ! bugfix: openacc porting for vector version of OL calculation added
 ! 
@@ -1586 +1589 @@
                 surf%ss(m) = -surf%ssws(m) / ( surf%us(m) + 1E-30_wp )
              ENDDO
+          ELSEIF ( .NOT. surf_vertical )  THEN
+
+             !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo )
+             DO  m = 1, surf%ns   
+                i    = surf%i(m)
+                j    = surf%j(m)
+                k    = surf%k(m)
+                z_mo = surf%z_mo(m)
+
+                surf%ss(m) = kappa * ( s(k,j,i) - s(k-1,j,i) )              &
+                                  / ( LOG( z_mo / surf%z0h(m) )             &
+                                      - psi_h( z_mo / surf%ol(m) )          &
+                                      + psi_h( surf%z0h(m) / surf%ol(m) ) )
+
+             ENDDO
           ENDIF
 !