Changeset 1557 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Mar 5, 2015 4:43:04 PM (10 years ago)

Author:

suehring

Message:

Enable monotone advection for scalars in combination with fifth-order scheme using monotonic limiter.

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (13 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
+! Adjustments for monotonic limiter
 ! 
 ! Former revisions:
@@ -143 +143 @@
         ONLY:   average_count_pr, cloud_droplets, cloud_physics, do_sum,       &
                 dt_3d, g, humidity, icloud_scheme, kappa, large_scale_forcing, &
-                large_scale_subsidence, max_pr_user, message_string, ocean,    &
+                large_scale_subsidence, max_pr_user, message_string,           &
+                monotonic_adjustment, ocean,                                   &
                 passive_scalar, precipitation, simulated_time,                 &
                 use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
@@ -281 +282 @@
        ENDIF
 
-       IF ( ws_scheme_sca .AND. sr == 0 )  THEN
+       IF ( ws_scheme_sca .AND. .NOT. monotonic_adjustment                    &
+           .AND. sr == 0 )  THEN
 
           DO  i = 0, threads_per_task-1
@@ -802 +804 @@
 !--             but so far there is no other suitable place to calculate)
                 IF ( ocean )  THEN
-                   IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
+                   IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.      &
+                             sr /= 0 )  THEN
                       pts = 0.5_wp * ( sa(k,j,i)   - hom(k,1,23,sr) +          &
                                        sa(k+1,j,i) - hom(k+1,1,23,sr) )
@@ -853 +856 @@
                       ENDIF
                    ELSE
-                      IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
+                      IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.   &
+                                sr /= 0 )  THEN
                          pts = 0.5_wp * ( vpt(k,j,i)   - hom(k,1,44,sr) +      &
                                           vpt(k+1,j,i) - hom(k+1,1,44,sr) )
                          sums_l(k,46,tn) = sums_l(k,46,tn) + pts * w(k,j,i) *  &
                                                              rmask(j,i,sr)
-                      ELSE IF ( ws_scheme_sca .AND. sr == 0 )  THEN
+                      ELSE IF ( ws_scheme_sca .AND. .NOT. monotonic_adjustment &
+                                .AND. sr == 0 )  THEN
                          sums_l(k,46,tn) = ( 1.0_wp + 0.61_wp *                & 
                                              hom(k,1,41,sr) ) *                &
@@ -870 +875 @@
 !--             Passive scalar flux
                 IF ( passive_scalar .AND. ( .NOT. ws_scheme_sca                &
-                     .OR. sr /= 0 ) )  THEN
+                     .OR. monotonic_adjustment .OR. sr /= 0 ) )  THEN
                    pts = 0.5_wp * ( q(k,j,i)   - hom(k,1,41,sr) +              &
                                     q(k+1,j,i) - hom(k+1,1,41,sr) )
@@ -914 +919 @@
 
        ENDIF
-       IF ( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
+       IF ( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR. sr /= 0 )  THEN
          !$OMP DO
          DO  i = nxl, nxr
@@ -1525 +1530 @@
         ONLY :  average_count_pr, cloud_droplets, cloud_physics, do_sum,       &
                 dt_3d, g, humidity, icloud_scheme, kappa, large_scale_forcing, &
-                large_scale_subsidence, max_pr_user, message_string, ocean,    &
+                large_scale_subsidence, max_pr_user, message_string,           &
+                monotonic_adjustment, ocean,                                   &
                 passive_scalar, precipitation, simulated_time,                 &
                 use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
@@ -1663 +1669 @@
        ENDIF
 
-       IF ( ws_scheme_sca .AND. sr == 0 )  THEN
+       IF ( ws_scheme_sca .AND. .NOT. monotonic_adjustment .AND. sr == 0 )  THEN
 
           DO  i = 0, threads_per_task-1
@@ -2486 +2492 @@
        IF ( ocean )  THEN
 
-          IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
+          IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR. sr /= 0 )  THEN
 
              !$acc parallel loop gang present( hom, rflags_invers, rmask, sa, sums_l, w ) create( s1 )
@@ -2628 +2634 @@
           ELSE
 
-             IF( .NOT. ws_scheme_sca  .OR.  sr /= 0 )  THEN
+             IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.  sr /= 0 )  THEN
 
                 !$acc parallel loop gang present( hom, rflags_invers, rmask, sums_l, vpt, w ) create( s1 )
@@ -2660 +2666 @@
 !
 !--    Passive scalar flux
-       IF ( passive_scalar  .AND.  ( .NOT. ws_scheme_sca  .OR.  sr /= 0 ) )  THEN
+       IF ( passive_scalar  .AND.  ( .NOT. ws_scheme_sca .OR. monotonic_adjustment &
+                                     .OR.  sr /= 0 ) )  THEN
 
           !$acc parallel loop gang present( hom, q, rflags_invers, rmask, sums_l, w ) create( s1 )
@@ -2714 +2721 @@
        ENDIF
 
-       IF ( .NOT. ws_scheme_sca  .OR.  sr /= 0 )  THEN
+       IF ( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.  sr /= 0 )  THEN
 
           !$OMP DO