Ignore:
Timestamp:
Mar 9, 2007 12:05:41 PM (15 years ago)
Author:
raasch
Message:

preliminary update of further changes, advec_particles is not running!

File:
1 edited

Legend:

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

    r39 r57  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Reference temperature pt_reference can be used in buoyancy term
    77!
    88! Former revisions:
     
    8282                dpt_dz = ( theta(k+1,j,i) - theta(k-1,j,i) ) * dd2zu(k)
    8383                IF ( dpt_dz > 0.0 ) THEN
    84                    l_stable = 0.76 * SQRT( e(k,j,i) ) / &
    85                                      SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
     84                   IF ( use_pt_reference )  THEN
     85                      l_stable = 0.76 * SQRT( e(k,j,i) ) / &
     86                                        SQRT( g / pt_reference * dpt_dz ) + 1E-5
     87                   ELSE
     88                      l_stable = 0.76 * SQRT( e(k,j,i) ) / &
     89                                        SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
     90                   ENDIF
    8691                ELSE
    8792                   l_stable = l_grid(k)
     
    195200          dpt_dz = ( theta(k+1,j,i) - theta(k-1,j,i) ) * dd2zu(k)
    196201          IF ( dpt_dz > 0.0 ) THEN
    197              l_stable = 0.76 * SQRT( e(k,j,i) ) / &
    198                                SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
     202             IF ( use_pt_reference )  THEN
     203                l_stable = 0.76 * SQRT( e(k,j,i) ) / &
     204                                  SQRT( g / pt_reference * dpt_dz ) + 1E-5
     205             ELSE
     206                l_stable = 0.76 * SQRT( e(k,j,i) ) / &
     207                                  SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
     208             ENDIF
    199209          ELSE
    200210             l_stable = l_grid(k)
Note: See TracChangeset for help on using the changeset viewer.