Ignore:
Timestamp:
Jul 12, 2016 4:34:24 PM (5 years ago)
Author:
suehring
Message:

Separate balance equations for humidity and passive_scalar

File:
1 edited

Legend:

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

    r1818 r1960  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Separate humidity and passive scalar
    2222!
    2323! Former revisions:
     
    7171
    7272    USE arrays_3d,                                                             &
    73         ONLY:  e, inflow_damping_factor, mean_inflow_profiles, pt, q, u, v, w
     73        ONLY:  e, inflow_damping_factor, mean_inflow_profiles, pt, q, s, u, v, w
    7474       
    7575    USE control_parameters,                                                    &
     
    9898    INTEGER(iwp) ::  prev     !< ID of sending PE for y-shift
    9999
    100     REAL(wp), DIMENSION(nzb:nzt+1,6,nbgp)           ::                         &
     100    REAL(wp), DIMENSION(nzb:nzt+1,7,nbgp)           ::                         &
    101101       avpr               !< stores averaged profiles at recycling plane
    102     REAL(wp), DIMENSION(nzb:nzt+1,6,nbgp)           ::                         &
     102    REAL(wp), DIMENSION(nzb:nzt+1,7,nbgp)           ::                         &
    103103       avpr_l             !< auxiliary variable to calculate avpr
    104     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,6,nbgp) ::                         &
     104    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,7,nbgp) ::                         &
    105105       inflow_dist        !< turbulence signal of vars, added at inflow boundary
    106     REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,6,nbgp) ::                         &
     106    REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,7,nbgp) ::                         &
    107107       local_inflow_dist  !< auxiliary variable for inflow_dist, used for yshift
    108108
     
    112112!-- Carry out spanwise averaging in the recycling plane
    113113    avpr_l = 0.0_wp
    114     ngp_pr = ( nzt - nzb + 2 ) * 6 * nbgp
     114    ngp_pr = ( nzt - nzb + 2 ) * 7 * nbgp
    115115    ngp_ifd = ngp_pr * ( nyn - nys + 1 + 2 * nbgp )
    116116
     
    131131                avpr_l(k,4,l) = avpr_l(k,4,l) + pt(k,j,i)
    132132                avpr_l(k,5,l) = avpr_l(k,5,l) + e(k,j,i)
    133                 IF ( humidity  .OR.  passive_scalar )                          &
     133                IF ( humidity )                                                &
    134134                   avpr_l(k,6,l) = avpr_l(k,6,l) + q(k,j,i)
     135                IF ( passive_scalar )                                          &
     136                   avpr_l(k,7,l) = avpr_l(k,7,l) + s(k,j,i)
    135137
    136138             ENDDO
     
    156158             avpr_l(k,4,l) = avpr_l(k,4,l) + pt(k,j,i)
    157159             avpr_l(k,5,l) = avpr_l(k,5,l) + e(k,j,i)
    158              IF ( humidity  .OR.  passive_scalar )                             &
     160             IF ( humidity )                                                   &
    159161                avpr_l(k,6,l) = avpr_l(k,6,l) + q(k,j,i)
     162             IF ( passive_scalar )                                             &
     163                avpr_l(k,7,l) = avpr_l(k,7,l) + s(k,j,i)
    160164
    161165          ENDDO
     
    183187                inflow_dist(k,j,4,l) = pt(k,j,i)  - avpr(k,4,l)
    184188                inflow_dist(k,j,5,l) = e(k,j,i)   - avpr(k,5,l)
    185                 IF ( humidity  .OR.  passive_scalar )                          &
     189                IF ( humidity )                                                &
    186190                   inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
     191                IF ( passive_scalar )                                          &
     192                   inflow_dist(k,j,7,l) = s(k,j,i) - avpr(k,7,l)
    187193            ENDDO
    188194          ENDDO
     
    201207             inflow_dist(k,j,4,l) = pt(k,j,i)  - avpr(k,4,l)
    202208             inflow_dist(k,j,5,l) = e(k,j,i)   - avpr(k,5,l)
    203              IF ( humidity  .OR.  passive_scalar )                             &
     209             IF ( humidity )                                                   &
    204210                inflow_dist(k,j,6,l) = q(k,j,i) - avpr(k,6,l)
     211             IF ( passive_scalar )                                             &
     212                inflow_dist(k,j,7,l) = s(k,j,i) - avpr(k,7,l)
    205213             
    206214          ENDDO
     
    276284             e(k,j,-nbgp:-1)  = MAX( e(k,j,-nbgp:-1), 0.0_wp )
    277285
    278              IF ( humidity  .OR.  passive_scalar )                          &
     286             IF ( humidity )                                                &
    279287                q(k,j,-nbgp:-1)  = mean_inflow_profiles(k,6) +              &
    280288                        inflow_dist(k,j,6,1:nbgp) * inflow_damping_factor(k)
     289             IF ( passive_scalar )                                          &
     290                s(k,j,-nbgp:-1)  = mean_inflow_profiles(k,7) +              &
     291                        inflow_dist(k,j,7,1:nbgp) * inflow_damping_factor(k)
    281292
    282293          ENDDO
Note: See TracChangeset for help on using the changeset viewer.