Ignore:
Timestamp:
Mar 10, 2015 5:57:55 PM (6 years ago)
Author:
suehring
Message:

Bugfixes in monotonic limter.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r1558 r1567  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! Reverse modifications made for monotonic limiter.
    2424!
    2525! Former revisions:
     
    146146        ONLY:   average_count_pr, cloud_droplets, cloud_physics, do_sum,       &
    147147                dt_3d, g, humidity, icloud_scheme, kappa, large_scale_forcing, &
    148                 large_scale_subsidence, max_pr_user, message_string,           &
    149                 monotonic_adjustment, ocean,                                   &
     148                large_scale_subsidence, max_pr_user, message_string, ocean,    &
    150149                passive_scalar, precipitation, simulated_time,                 &
    151150                use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
     
    285284       ENDIF
    286285
    287        IF ( ws_scheme_sca .AND. .NOT. monotonic_adjustment                    &
    288            .AND. sr == 0 )  THEN
     286       IF ( ws_scheme_sca .AND. sr == 0 )  THEN
    289287
    290288          DO  i = 0, threads_per_task-1
     
    807805!--             but so far there is no other suitable place to calculate)
    808806                IF ( ocean )  THEN
    809                    IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.      &
    810                              sr /= 0 )  THEN
     807                   IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
    811808                      pts = 0.5_wp * ( sa(k,j,i)   - hom(k,1,23,sr) +          &
    812809                                       sa(k+1,j,i) - hom(k+1,1,23,sr) )
     
    859856                      ENDIF
    860857                   ELSE
    861                       IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.   &
    862                                 sr /= 0 )  THEN
     858                      IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
    863859                         pts = 0.5_wp * ( vpt(k,j,i)   - hom(k,1,44,sr) +      &
    864860                                          vpt(k+1,j,i) - hom(k+1,1,44,sr) )
    865861                         sums_l(k,46,tn) = sums_l(k,46,tn) + pts * w(k,j,i) *  &
    866862                                                             rmask(j,i,sr)
    867                       ELSE IF ( ws_scheme_sca .AND. .NOT. monotonic_adjustment &
    868                                 .AND. sr == 0 )  THEN
     863                      ELSE IF ( ws_scheme_sca .AND. sr == 0 )  THEN
    869864                         sums_l(k,46,tn) = ( 1.0_wp + 0.61_wp *                &
    870865                                             hom(k,1,41,sr) ) *                &
     
    878873!--             Passive scalar flux
    879874                IF ( passive_scalar .AND. ( .NOT. ws_scheme_sca                &
    880                      .OR. monotonic_adjustment .OR. sr /= 0 ) )  THEN
     875                     .OR. sr /= 0 ) )  THEN
    881876                   pts = 0.5_wp * ( q(k,j,i)   - hom(k,1,41,sr) +              &
    882877                                    q(k+1,j,i) - hom(k+1,1,41,sr) )
     
    922917
    923918       ENDIF
    924        IF ( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR. sr /= 0 )  THEN
     919       IF ( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
    925920         !$OMP DO
    926921         DO  i = nxl, nxr
     
    15331528        ONLY :  average_count_pr, cloud_droplets, cloud_physics, do_sum,       &
    15341529                dt_3d, g, humidity, icloud_scheme, kappa, large_scale_forcing, &
    1535                 large_scale_subsidence, max_pr_user, message_string,           &
    1536                 monotonic_adjustment, ocean,                                   &
     1530                large_scale_subsidence, max_pr_user, message_string, ocean,    &
    15371531                passive_scalar, precipitation, simulated_time,                 &
    15381532                use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
     
    16721666       ENDIF
    16731667
    1674        IF ( ws_scheme_sca .AND. .NOT. monotonic_adjustment .AND. sr == 0 )  THEN
     1668       IF ( ws_scheme_sca .AND. sr == 0 )  THEN
    16751669
    16761670          DO  i = 0, threads_per_task-1
     
    24952489       IF ( ocean )  THEN
    24962490
    2497           IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR. sr /= 0 )  THEN
     2491          IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
    24982492
    24992493             !$acc parallel loop gang present( hom, rflags_invers, rmask, sa, sums_l, w ) create( s1 )
     
    26372631          ELSE
    26382632
    2639              IF( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.  sr /= 0 )  THEN
     2633             IF( .NOT. ws_scheme_sca .OR.  sr /= 0 )  THEN
    26402634
    26412635                !$acc parallel loop gang present( hom, rflags_invers, rmask, sums_l, vpt, w ) create( s1 )
     
    26692663!
    26702664!--    Passive scalar flux
    2671        IF ( passive_scalar  .AND.  ( .NOT. ws_scheme_sca .OR. monotonic_adjustment &
    2672                                      .OR.  sr /= 0 ) )  THEN
     2665       IF ( passive_scalar  .AND.  ( .NOT. ws_scheme_sca  .OR.  sr /= 0 ) )  THEN
    26732666
    26742667          !$acc parallel loop gang present( hom, q, rflags_invers, rmask, sums_l, w ) create( s1 )
     
    27242717       ENDIF
    27252718
    2726        IF ( .NOT. ws_scheme_sca .OR. monotonic_adjustment .OR.  sr /= 0 )  THEN
     2719       IF ( .NOT. ws_scheme_sca .OR.  sr /= 0 )  THEN
    27272720
    27282721          !$OMP DO
Note: See TracChangeset for help on using the changeset viewer.