Ignore:
Timestamp:
Aug 20, 2019 11:55:33 AM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/inflow_turbulence.f90

    r3655 r4172  
    2525! -----------------
    2626! $Id$
     27! added optional recycling of absolute values for pt and q
     28!
     29! 3655 2019-01-07 16:51:22Z knoop
    2730! Corrected "Former revisions" section
    2831!
     
    8790       
    8891    USE control_parameters,                                                    &
    89         ONLY:  humidity, passive_scalar, recycling_plane, recycling_yshift
     92        ONLY:  humidity, passive_scalar, recycling_plane, recycling_yshift, raq
    9093       
    9194    USE cpulog,                                                                &
     
    101104
    102105    IMPLICIT NONE
    103 
     106   
    104107    INTEGER(iwp) ::  i        !< loop index
    105108    INTEGER(iwp) ::  j        !< loop index
     
    119122    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,7,nbgp) ::                         &
    120123       local_inflow_dist  !< auxiliary variable for inflow_dist, used for yshift
    121 
     124   
    122125    CALL cpu_log( log_point(40), 'inflow_turbulence', 'start' )
    123 
     126   
    124127!
    125128!-- Carry out spanwise averaging in the recycling plane
     
    187190!
    188191!-- Calculate the disturbances at the recycling plane
     192!-- for recycling of absolute quantities, the disturbance is defined as the absolute value
     193!-- (and not as the deviation from the mean profile)
    189194    i = recycling_plane
    190195
     
    198203                inflow_dist(k,j,2,l) = v(k,j,i)   - avpr(k,2,l)
    199204                inflow_dist(k,j,3,l) = w(k,j,i)   - avpr(k,3,l)
    200                 inflow_dist(k,j,4,l) = pt(k,j,i)  - avpr(k,4,l)
     205                IF ( raq(4) )  THEN
     206                   inflow_dist(k,j,4,l) = pt(k,j,i)
     207                ELSE
     208                   inflow_dist(k,j,4,l) = pt(k,j,i) - avpr(k,4,l)
     209                ENDIF
    201210                inflow_dist(k,j,5,l) = e(k,j,i)   - avpr(k,5,l)
    202                 IF ( humidity )                                                &
    203                    inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
     211                IF ( humidity ) THEN
     212                   IF ( raq(6) ) THEN
     213                      inflow_dist(k,j,6,l) = q(k,j,i)
     214                   ELSE
     215                      inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
     216                   ENDIF
     217                ENDIF
    204218                IF ( passive_scalar )                                          &
    205219                   inflow_dist(k,j,7,l) = s(k,j,i) - avpr(k,7,l)
     
    214228       DO  j = nysg, nyng
    215229          DO  k = nzb, nzt+1
    216 
    217230             inflow_dist(k,j,1,l) = u(k,j,i+1) - avpr(k,1,l)
    218231             inflow_dist(k,j,2,l) = v(k,j,i)   - avpr(k,2,l)
    219232             inflow_dist(k,j,3,l) = w(k,j,i)   - avpr(k,3,l)
    220              inflow_dist(k,j,4,l) = pt(k,j,i)  - avpr(k,4,l)
     233             IF ( raq(4) )  THEN
     234                inflow_dist(k,j,4,l) = pt(k,j,i)
     235             ELSE
     236                inflow_dist(k,j,4,l) = pt(k,j,i) - avpr(k,4,l)
     237             ENDIF
    221238             inflow_dist(k,j,5,l) = e(k,j,i)   - avpr(k,5,l)
    222              IF ( humidity )                                                   &
    223                 inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
     239             IF ( humidity )  THEN
     240                IF ( raq(6) ) THEN
     241                   inflow_dist(k,j,6,l) = q(k,j,i)
     242                ELSE
     243                   inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
     244                ENDIF
     245             ENDIF
    224246             IF ( passive_scalar )                                             &
    225247                inflow_dist(k,j,7,l) = s(k,j,i) - avpr(k,7,l)
     
    291313             w(k,j,-nbgp:-1)  =                                             &
    292314                        inflow_dist(k,j,3,1:nbgp) * inflow_damping_factor(k)
    293              pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) +                 &
    294                         inflow_dist(k,j,4,1:nbgp) * inflow_damping_factor(k)
    295              e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) +                 &
     315             IF ( raq(4) )  THEN
     316                pt(k,j,-nbgp:-1) = inflow_dist(k,j,4,1:nbgp)
     317             ELSE
     318                pt(k,j,-nbgp:-1) = mean_inflow_profiles(k,4) +               &
     319                inflow_dist(k,j,4,1:nbgp) * inflow_damping_factor(k)
     320             ENDIF
     321             e(k,j,-nbgp:-1)  = mean_inflow_profiles(k,5) +                  &
    296322                        inflow_dist(k,j,5,1:nbgp) * inflow_damping_factor(k)
    297323             e(k,j,-nbgp:-1)  = MAX( e(k,j,-nbgp:-1), 0.0_wp )
    298324
    299              IF ( humidity )                                                &
    300                 q(k,j,-nbgp:-1)  = mean_inflow_profiles(k,6) +              &
     325             IF ( humidity )  THEN
     326                IF ( raq(6) )  THEN
     327                   q(k,j,-nbgp:-1)  = inflow_dist(k,j,6,1:nbgp)
     328                ELSE
     329                   q(k,j,-nbgp:-1)  = mean_inflow_profiles(k,6) +           &
    301330                        inflow_dist(k,j,6,1:nbgp) * inflow_damping_factor(k)
     331                ENDIF
     332             ENDIF
    302333             IF ( passive_scalar )                                          &
    303334                s(k,j,-nbgp:-1)  = mean_inflow_profiles(k,7) +              &
    304335                        inflow_dist(k,j,7,1:nbgp) * inflow_damping_factor(k)
    305 
    306336          ENDDO
    307337       ENDDO
Note: See TracChangeset for help on using the changeset viewer.