Changeset 1936 for palm/trunk/SOURCE/lpm_advec.f90

Timestamp:

Jun 13, 2016 1:37:44 PM (8 years ago)

Author:

suehring

Message:

deallocation of unused particle memory, formatting adjustments

File:

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

palm/trunk/SOURCE/lpm_advec.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-!  
+! Formatting adjustments
 ! 
 ! Former revisions:
@@ -264 +264 @@
           jlog = ( particles(n)%y + 0.5_wp * dy ) * ddy
 
-          IF ( constant_flux_layer .AND. zv(n) - zw(nzb_s_inner(jlog,ilog)) < z_p )  THEN
+          IF ( constant_flux_layer  .AND.                                      &
+               zv(n) - zw(nzb_s_inner(jlog,ilog)) < z_p )  THEN
 !
 !--          Resolved-scale horizontal particle velocity is zero below z0. 
@@ -272 +273 @@
 !
 !--             Determine the sublayer. Further used as index.
-                height_p     = ( zv(n) - zw(nzb_s_inner(jlog,ilog)) - z0_av_global )            &
-                                     * REAL( number_of_sublayers, KIND=wp ) &
+                height_p = ( zv(n) - zw(nzb_s_inner(jlog,ilog)) - z0_av_global ) &
+                                     * REAL( number_of_sublayers, KIND=wp )    &
                                      * d_z_p_z0 
 !
@@ -340 +341 @@
           ilog = ( particles(n)%x + 0.5_wp * dx ) * ddx
           jlog = ( particles(n)%y + 0.5_wp * dy ) * ddy
-          IF ( constant_flux_layer .AND. zv(n) - zw(nzb_s_inner(jlog,ilog)) < z_p )  THEN
+          IF ( constant_flux_layer  .AND.                                      &
+               zv(n) - zw(nzb_s_inner(jlog,ilog)) < z_p )  THEN
 
              IF ( zv(n) - zw(nzb_s_inner(jlog,ilog)) < z0_av_global )  THEN
@@ -352 +354 @@
 !--             topography particle on u-grid can be above surface-layer height,
 !--             whereas it can be below on v-grid. 
-                height_p     = ( zv(n) - zw(nzb_s_inner(jlog,ilog)) - z0_av_global ) &
-                                  * REAL( number_of_sublayers, KIND=wp )    &
+                height_p = ( zv(n) - zw(nzb_s_inner(jlog,ilog)) - z0_av_global ) &
+                                  * REAL( number_of_sublayers, KIND=wp )       &
                                   * d_z_p_z0 
 !
@@ -367 +369 @@
 !--             friction velocity can become very small, resulting in a too
 !--             large particle speed.
-                us_int   = MAX( 0.5_wp * ( us(jlog,ilog) + us(jlog-1,ilog) ), &
+                us_int   = MAX( 0.5_wp * ( us(jlog,ilog) + us(jlog-1,ilog) ),  &
                                 0.01_wp )    
 !
@@ -376 +378 @@
 !--             as sensitivity studies revealed no significant effect of 
 !--             using the neutral solution also for un/stable situations.
-                v_int(n) = -vsws(jlog,ilog) / ( us_int * kappa + 1E-10_wp )          &
+                v_int(n) = -vsws(jlog,ilog) / ( us_int * kappa + 1E-10_wp )    &
                          * log_z_z0_int - v_gtrans