Ignore:
Timestamp:
Feb 23, 2007 4:53:48 AM (15 years ago)
Author:
raasch
Message:

preliminary version of modified boundary conditions at top

File:
1 edited

Legend:

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

    r4 r19  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Calculation extended for gridpoint nzt
    77!
    88! Former revisions:
     
    5656       REAL, DIMENSION(:,:), POINTER   ::  rif
    5757       REAL, DIMENSION(:,:,:), POINTER ::  e, km, theta
    58        REAL, DIMENSION(nzb+1:nzt-1,nys:nyn) ::  dissipation, l, ll
     58       REAL, DIMENSION(nzb+1:nzt,nys:nyn) ::  dissipation, l, ll
    5959 
    6060
     
    7272             ENDIF
    7373
    74              DO  k = nzb_s_inner(j,i)+1, nzt-1
     74             DO  k = nzb_s_inner(j,i)+1, nzt
    7575!
    7676!--             Calculate the mixing length (for dissipation)
     
    102102!--       Calculate the tendency terms
    103103          DO  j = nys, nyn
    104              DO  k = nzb_s_inner(j,i)+1, nzt-1
     104             DO  k = nzb_s_inner(j,i)+1, nzt
    105105
    106106                 dissipation(k,j) = ( 0.19 + 0.74 * l(k,j) / ll(k,j) ) * &
     
    130130          IF ( use_sgs_for_particles )  THEN
    131131             DO  j = nys, nyn
    132                 DO  k = nzb_s_inner(j,i)+1, nzt-1
     132                DO  k = nzb_s_inner(j,i)+1, nzt
    133133                   diss(k,j,i) = dissipation(k,j)
    134134                ENDDO
     
    171171       REAL, DIMENSION(:,:), POINTER   ::  rif
    172172       REAL, DIMENSION(:,:,:), POINTER ::  e, km, theta
    173        REAL, DIMENSION(nzb+1:nzt-1)    ::  dissipation, l, ll
     173       REAL, DIMENSION(nzb+1:nzt)    ::  dissipation, l, ll
    174174
    175175
     
    187187!
    188188!--    Calculate the mixing length (for dissipation)
    189        DO  k = nzb_s_inner(j,i)+1, nzt-1
     189       DO  k = nzb_s_inner(j,i)+1, nzt
    190190          dpt_dz = ( theta(k+1,j,i) - theta(k-1,j,i) ) * dd2zu(k)
    191191          IF ( dpt_dz > 0.0 ) THEN
     
    234234!--    Store dissipation if needed for calculating the sgs particle velocities
    235235       IF ( use_sgs_for_particles )  THEN
    236           DO  k = nzb_s_inner(j,i)+1, nzt-1
     236          DO  k = nzb_s_inner(j,i)+1, nzt
    237237             diss(k,j,i) = dissipation(k)
    238238          ENDDO
Note: See TracChangeset for help on using the changeset viewer.