Changeset 856 for palm/trunk/SOURCE

Timestamp:

Mar 20, 2012 6:29:17 PM (13 years ago)

Author:

suehring

Message:

Bug concerning diss_2nd. Return value is set to zero.

File:

• palm/trunk/SOURCE/advec_ws.f90 (modified) (2 diffs)

palm/trunk/SOURCE/advec_ws.f90

@@ -4 +4 @@
 ! Current revisions:
 ! ------------------
+! Bug concerning numerical dissipation at the first grid level. Function
+! diss_2nd always returns zero.  Because diss_2nd will be replaced by numerical
+! dissipation which is more consistent to the WS-schemes, the return value
+! is set to zero for now.
 !
 ! Former revisions:
@@ -3224 +3228 @@
        REAL, INTENT(IN)  :: v2, v1, v0, vm1, vm2, vel_comp, factor,            &
                             grid
-       REAL :: value_min_m, value_max_m, value_min, value_max,                 &
-               value_min_p, value_max_p, diss_m, diss_0, diss_p
-             
-       value_min_m = MIN( v0, vm1, vm2 )
-       value_max_m = MAX( v0, vm1, vm2 )
-       value_min   = MIN( v1, v0,  vm2 )
-       value_max   = MAX( v1, v0,  vm2 )
-       value_min_p = MIN( v2, v1,  v0  )
-       value_max_p = MAX( v2, v1,  v0  )
-       
-       diss_m = MAX( 0.0, vm1 - value_max_m ) + MIN( 0.0, vm1 - value_min_m )
-       diss_0 = MAX( 0.0, v0  - value_max   ) + MIN( 0.0, v0  - value_min   )
-       diss_p = MAX( 0.0, v1  - value_max_p ) + MIN( 0.0, v1  - value_min_p )
-       
-       dissip = ABS( vel_comp ) * ( diss_p - 2.0 * diss_0 + diss_m )          &
-                                * grid**2 * factor
+
+       dissip = 0.0
               
     END FUNCTION diss_2nd