Changeset 3152 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

Timestamp:

Jul 19, 2018 1:26:52 PM (6 years ago)

Author:

suehring

Message:

Further adjustments for surface structure

File:

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

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! q_surface is now part of surface structure
+! 
+! 3149 2018-07-19 05:48:09Z maronga
 ! Correct accidental last commit
 ! 
@@ -338 +341 @@
                 CALL calc_pt_surface
                 IF ( humidity )  THEN
+                   CALL calc_q_surface
                    CALL calc_vpt_surface
                 ENDIF
@@ -1595 +1599 @@
 
 !
-!-- Calculate potential temperature at first grid level
+!-- Set potential temperature at surface grid level.
 !-- ( only for upward-facing surfs )
     SUBROUTINE calc_pt_surface
@@ -1619 +1623 @@
 
 !
-!-- Calculate virtual potential temperature at first grid level
+!-- Set mixing ratio at surface grid level. ( Only for upward-facing surfs. )
+    SUBROUTINE calc_q_surface
+
+       IMPLICIT NONE
+
+       INTEGER(iwp) ::  k_off   !< index offset between surface and atmosphere grid point (-1 for upward-, +1 for downward-facing walls) 
+       INTEGER(iwp) ::  m       !< loop variable over all horizontal surf elements 
+       
+       k_off = surf%koff
+       !$OMP PARALLEL DO PRIVATE( i, j, k )
+       DO  m = 1, surf%ns 
+
+          i   = surf%i(m)            
+          j   = surf%j(m)
+          k   = surf%k(m)
+
+          surf%q_surface(m) = q(k+k_off,j,i)
+
+       ENDDO
+
+    END SUBROUTINE calc_q_surface
+    
+!
+!-- Set virtual potential temperature at surface grid level.
 !-- ( only for upward-facing surfs )
     SUBROUTINE calc_vpt_surface
@@ -1781 +1808 @@
                    z_mo = surf%z_mo(m)
 
-                   surf%qs(m) = kappa * ( surf%qv1(m) - q(k-1,j,i) )           &
+                   surf%qs(m) = kappa * ( surf%qv1(m) - surf%q_surface(m) )    &
                                         / ( LOG( z_mo / surf%z0q(m) )          &
                                             - psi_h( z_mo / surf%ol(m) )       &