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/diffusivities.f90

    r4 r57  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Reference temperature pt_reference can be used in buoyancy term
    77!
    88! Former revisions:
     
    9090             dpt_dz = ( theta(k+1,j,i) - theta(k-1,j,i) ) * dd2zu(k)
    9191             IF ( dpt_dz > 0.0 ) THEN
    92                 l_stable = 0.76 * sqrt_e(k) / &
    93                                   SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
     92                IF ( use_pt_reference )  THEN
     93                   l_stable = 0.76 * sqrt_e(k) / &
     94                                     SQRT( g / pt_reference * dpt_dz ) + 1E-5
     95                ELSE
     96                   l_stable = 0.76 * sqrt_e(k) / &
     97                                     SQRT( g / theta(k,j,i) * dpt_dz ) + 1E-5
     98                ENDIF
    9499             ELSE
    95100                l_stable = l_grid(k)
Note: See TracChangeset for help on using the changeset viewer.