Changeset 856


Ignore:
Timestamp:
Mar 20, 2012 6:29:17 PM (10 years ago)
Author:
suehring
Message:

Bug concerning diss_2nd. Return value is set to zero.

File:
1 edited

Legend:

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

    r802 r856  
    44! Current revisions:
    55! ------------------
     6! Bug concerning numerical dissipation at the first grid level. Function
     7! diss_2nd always returns zero.  Because diss_2nd will be replaced by numerical
     8! dissipation which is more consistent to the WS-schemes, the return value
     9! is set to zero for now.
    610!
    711! Former revisions:
     
    32243228       REAL, INTENT(IN)  :: v2, v1, v0, vm1, vm2, vel_comp, factor,            &
    32253229                            grid
    3226        REAL :: value_min_m, value_max_m, value_min, value_max,                 &
    3227                value_min_p, value_max_p, diss_m, diss_0, diss_p
    3228              
    3229        value_min_m = MIN( v0, vm1, vm2 )
    3230        value_max_m = MAX( v0, vm1, vm2 )
    3231        value_min   = MIN( v1, v0,  vm2 )
    3232        value_max   = MAX( v1, v0,  vm2 )
    3233        value_min_p = MIN( v2, v1,  v0  )
    3234        value_max_p = MAX( v2, v1,  v0  )
    3235        
    3236        diss_m = MAX( 0.0, vm1 - value_max_m ) + MIN( 0.0, vm1 - value_min_m )
    3237        diss_0 = MAX( 0.0, v0  - value_max   ) + MIN( 0.0, v0  - value_min   )
    3238        diss_p = MAX( 0.0, v1  - value_max_p ) + MIN( 0.0, v1  - value_min_p )
    3239        
    3240        dissip = ABS( vel_comp ) * ( diss_p - 2.0 * diss_0 + diss_m )          &
    3241                                 * grid**2 * factor
     3230
     3231       dissip = 0.0
    32423232             
    32433233    END FUNCTION diss_2nd
Note: See TracChangeset for help on using the changeset viewer.