Changeset 4172 for palm/trunk/SOURCE/inflow_turbulence.f90

Timestamp:

Aug 20, 2019 11:55:33 AM (5 years ago)

Author:

oliver.maas

Message:

Added optional method for recycling of absolute values of pt and q instead of recycling only the deviations from the mean profile at the recycling plane. With the new method two problems are solved: 1. A horizontally homogeneous temperature and humidity field is achieved, even for non-neutral boundary layers and thus no thermal circulation is triggered. 2. No gravity waves build up at the inflow due to cyclic boundary conditions for pt and q.

File:

• palm/trunk/SOURCE/inflow_turbulence.f90 (modified) (8 diffs)

palm/trunk/SOURCE/inflow_turbulence.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! added optional recycling of absolute values for pt and q
+! 
+! 3655 2019-01-07 16:51:22Z knoop
 ! Corrected "Former revisions" section
 ! 
@@ -87 +90 @@
         
     USE control_parameters,                                                    &
-        ONLY:  humidity, passive_scalar, recycling_plane, recycling_yshift
+        ONLY:  humidity, passive_scalar, recycling_plane, recycling_yshift, raq
         
     USE cpulog,                                                                &
@@ -101 +104 @@
 
     IMPLICIT NONE
-
+    
     INTEGER(iwp) ::  i        !< loop index
     INTEGER(iwp) ::  j        !< loop index
@@ -119 +122 @@
     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,7,nbgp) ::                         &
        local_inflow_dist  !< auxiliary variable for inflow_dist, used for yshift
-
+    
     CALL cpu_log( log_point(40), 'inflow_turbulence', 'start' )
-
+    
 !
 !-- Carry out spanwise averaging in the recycling plane
@@ -187 +190 @@
 !
 !-- Calculate the disturbances at the recycling plane
+!-- for recycling of absolute quantities, the disturbance is defined as the absolute value
+!-- (and not as the deviation from the mean profile)
     i = recycling_plane
 
@@ -198 +203 @@
                 inflow_dist(k,j,2,l) = v(k,j,i)   - avpr(k,2,l)
                 inflow_dist(k,j,3,l) = w(k,j,i)   - avpr(k,3,l)
-                inflow_dist(k,j,4,l) = pt(k,j,i)  - avpr(k,4,l)
+                IF ( raq(4) )  THEN
+                   inflow_dist(k,j,4,l) = pt(k,j,i)
+                ELSE
+                   inflow_dist(k,j,4,l) = pt(k,j,i) - avpr(k,4,l)
+                ENDIF
                 inflow_dist(k,j,5,l) = e(k,j,i)   - avpr(k,5,l)
-                IF ( humidity )                                                &
-                   inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
+                IF ( humidity ) THEN
+                   IF ( raq(6) ) THEN
+                      inflow_dist(k,j,6,l) = q(k,j,i)
+                   ELSE
+                      inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
+                   ENDIF
+                ENDIF
                 IF ( passive_scalar )                                          &
                    inflow_dist(k,j,7,l) = s(k,j,i) - avpr(k,7,l)
@@ -214 +228 @@
        DO  j = nysg, nyng
           DO  k = nzb, nzt+1
-
              inflow_dist(k,j,1,l) = u(k,j,i+1) - avpr(k,1,l)
              inflow_dist(k,j,2,l) = v(k,j,i)   - avpr(k,2,l)
              inflow_dist(k,j,3,l) = w(k,j,i)   - avpr(k,3,l)
-             inflow_dist(k,j,4,l) = pt(k,j,i)  - avpr(k,4,l)
+             IF ( raq(4) )  THEN
+                inflow_dist(k,j,4,l) = pt(k,j,i)
+             ELSE
+                inflow_dist(k,j,4,l) = pt(k,j,i) - avpr(k,4,l) 
+             ENDIF
              inflow_dist(k,j,5,l) = e(k,j,i)   - avpr(k,5,l)
-             IF ( humidity )                                                   &
-                inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
+             IF ( humidity )  THEN
+                IF ( raq(6) ) THEN
+                   inflow_dist(k,j,6,l) = q(k,j,i)
+                ELSE
+                   inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
+                ENDIF
+             ENDIF
              IF ( passive_scalar )                                             &
                 inflow_dist(k,j,7,l) = s(k,j,i) - avpr(k,7,l)
@@ -291 +313 @@
              w(k,j,-nbgp:-1)  =                                             &
                         inflow_dist(k,j,3,1:nbgp) * inflow_damping_factor(k)
-             pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) +                 &
-                        inflow_dist(k,j,4,1:nbgp) * inflow_damping_factor(k)
-             e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) +                 &
+             IF ( raq(4) )  THEN
+                pt(k,j,-nbgp:-1) = inflow_dist(k,j,4,1:nbgp)
+             ELSE
+                pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) +               &
+                inflow_dist(k,j,4,1:nbgp) * inflow_damping_factor(k)
+             ENDIF
+             e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) +                  &
                         inflow_dist(k,j,5,1:nbgp) * inflow_damping_factor(k)
              e(k,j,-nbgp:-1)  = MAX( e(k,j,-nbgp:-1), 0.0_wp )
 
-             IF ( humidity )                                                &
-                q(k,j,-nbgp:-1)  = mean_inflow_profiles(k,6) +              &
+             IF ( humidity )  THEN
+                IF ( raq(6) )  THEN
+                   q(k,j,-nbgp:-1)  = inflow_dist(k,j,6,1:nbgp)
+                ELSE 
+                   q(k,j,-nbgp:-1)  = mean_inflow_profiles(k,6) +           & 
                         inflow_dist(k,j,6,1:nbgp) * inflow_damping_factor(k)
+                ENDIF
+             ENDIF
              IF ( passive_scalar )                                          &
                 s(k,j,-nbgp:-1)  = mean_inflow_profiles(k,7) +              &
                         inflow_dist(k,j,7,1:nbgp) * inflow_damping_factor(k)
-
           ENDDO
        ENDDO