Changeset 1646 for palm

Timestamp:

Sep 2, 2015 4:00:10 PM (9 years ago)

Author:

hoffmann

Message:

bugfix: wrong calculation of d_mean in microphysics

File:

• palm/trunk/SOURCE/microphysics.f90 (modified) (5 diffs)

palm/trunk/SOURCE/microphysics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: Wrong computation of d_mean.
 ! 
 ! Former revisions:
@@ -877 +877 @@
 
                 DO k = nzb_s_inner(j,i)+1, nzt
-                   d_mean = 0.5_wp * ( qr(k+1,j,i) + qr(k-1,j,i) )
+                   d_mean = 0.5_wp * ( qr(k+1,j,i) - qr(k-1,j,i) )
                    d_min  = qr(k,j,i) - MIN( qr(k+1,j,i), qr(k,j,i), qr(k-1,j,i) )
                    d_max  = MAX( qr(k+1,j,i), qr(k,j,i), qr(k-1,j,i) ) - qr(k,j,i)
@@ -885 +885 @@
                                                               ABS( d_mean ) )
 
-                   d_mean = 0.5_wp * ( nr(k+1,j,i) + nr(k-1,j,i) )
+                   d_mean = 0.5_wp * ( nr(k+1,j,i) - nr(k-1,j,i) )
                    d_min  = nr(k,j,i) - MIN( nr(k+1,j,i), nr(k,j,i), nr(k-1,j,i) )
                    d_max  = MAX( nr(k+1,j,i), nr(k,j,i), nr(k-1,j,i) ) - nr(k,j,i)
@@ -1685 +1685 @@
 
           DO k = nzb_s_inner(j,i)+1, nzt
-             d_mean = 0.5_wp * ( qr_1d(k+1) + qr_1d(k-1) )
+             d_mean = 0.5_wp * ( qr_1d(k+1) - qr_1d(k-1) )
              d_min  = qr_1d(k) - MIN( qr_1d(k+1), qr_1d(k), qr_1d(k-1) )
              d_max  = MAX( qr_1d(k+1), qr_1d(k), qr_1d(k-1) ) - qr_1d(k)
@@ -1693 +1693 @@
                                                         ABS( d_mean ) )
 
-             d_mean = 0.5_wp * ( nr_1d(k+1) + nr_1d(k-1) )
+             d_mean = 0.5_wp * ( nr_1d(k+1) - nr_1d(k-1) )
              d_min  = nr_1d(k) - MIN( nr_1d(k+1), nr_1d(k), nr_1d(k-1) )
              d_max  = MAX( nr_1d(k+1), nr_1d(k), nr_1d(k-1) ) - nr_1d(k)