Changeset 1888

Timestamp:

Apr 21, 2016 12:20:49 PM (9 years ago)

Author:

suehring

Message:

Bugfix concerning logarithmic interpolation of particle speed

File:

• palm/trunk/SOURCE/lpm_advec.f90 (modified) (3 diffs)

palm/trunk/SOURCE/lpm_advec.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-!
+! Bugfix concerning logarithmic interpolation of particle speed
 ! 
 ! Former revisions:
@@ -273 +273 @@
 !--             using the neutral solution also for un/stable situations.
 !--             Calculated left and bottom index on u grid.
-                us_int = 0.5_wp * ( us(j,i) + us(j,i-1) )  
-
-                u_int  = -usws(j,i) / ( us_int * kappa + 1E-10_wp )           & 
+                us_int   = 0.5_wp * ( us(j,i) + us(j,i-1) )  
+
+                u_int(n) = -usws(j,i) / ( us_int * kappa + 1E-10_wp )          & 
                          * log_z_z0_int(n) - u_gtrans
 
@@ -331 +331 @@
 !--             using the neutral solution also for un/stable situations.
 !--             Calculated left and bottom index on v grid.
-                us_int = 0.5_wp * ( us(j,i) + us(j-1,i) )    
-
-                v_int  = -vsws(j,i) / ( us_int * kappa + 1E-10_wp )             &
+                us_int   = 0.5_wp * ( us(j,i) + us(j-1,i) )    
+
+                v_int(n) = -vsws(j,i) / ( us_int * kappa + 1E-10_wp )          &
                          * log_z_z0_int(n) - v_gtrans
              ENDIF