Ignore:
Timestamp:
Jun 2, 2020 10:22:25 AM (4 years ago)
Author:
suehring
Message:

Bugfix in summation for statistical regions

File:
1 edited

Legend:

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

    r4521 r4551  
    2525! -----------------
    2626! $Id$
     27! Bugfix in summation for statistical regions
     28!
     29! 4521 2020-05-06 11:39:49Z schwenkel
    2730! Rename variable
    2831!
     
    15011504             IF ( i >= nxl  .AND.  i <= nxr  .AND.                             &
    15021505                  j >= nys  .AND.  j <= nyn )  THEN
    1503                 sums_l(nzb,93,tn)  = sums_l(nzb,93,tn) + surf_lsm_h%ghf(m)
    1504                 sums_l(nzb,94,tn)  = sums_l(nzb,94,tn) + surf_lsm_h%qsws_liq(m)
    1505                 sums_l(nzb,95,tn)  = sums_l(nzb,95,tn) + surf_lsm_h%qsws_soil(m)
    1506                 sums_l(nzb,96,tn)  = sums_l(nzb,96,tn) + surf_lsm_h%qsws_veg(m)
    1507                 sums_l(nzb,97,tn)  = sums_l(nzb,97,tn) + surf_lsm_h%r_a(m)
    1508                 sums_l(nzb,98,tn) = sums_l(nzb,98,tn)+ surf_lsm_h%r_s(m)
     1506                sums_l(nzb,93,tn)  = sums_l(nzb,93,tn) + surf_lsm_h%ghf(m)       * rmask(j,i,sr)
     1507                sums_l(nzb,94,tn)  = sums_l(nzb,94,tn) + surf_lsm_h%qsws_liq(m)  * rmask(j,i,sr)
     1508                sums_l(nzb,95,tn)  = sums_l(nzb,95,tn) + surf_lsm_h%qsws_soil(m) * rmask(j,i,sr)
     1509                sums_l(nzb,96,tn)  = sums_l(nzb,96,tn) + surf_lsm_h%qsws_veg(m)  * rmask(j,i,sr)
     1510                sums_l(nzb,97,tn)  = sums_l(nzb,97,tn) + surf_lsm_h%r_a(m)       * rmask(j,i,sr)
     1511                sums_l(nzb,98,tn)  = sums_l(nzb,98,tn) + surf_lsm_h%r_s(m)       * rmask(j,i,sr)
    15091512             ENDIF
    15101513          ENDDO
     
    16431646                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i) )  &
    16441647                            - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) )**2&
    1645                 + w(k,j,i)**2                                        ) * flag
     1648                + w(k,j,i)**2                                        ) * flag * rmask(j,i,sr)
    16461649
    16471650                   sums_ll(k,2) = sums_ll(k,2) + 0.5_wp * w(k,j,i)             &
    16481651                                       * ( ( p(k,j,i) + p(k+1,j,i) )           &
    16491652                                         / momentumflux_output_conversion(k) ) &
    1650                                        * flag
     1653                                       * flag * rmask(j,i,sr)
    16511654
    16521655                ENDDO
     
    16841687                 - (km(k-1,j,i)+km(k,j,i)) * (e(k,j,i)-e(k-1,j,i)) * ddzu(k)   &
    16851688                                                                ) * ddzw(k)    &
    1686                                                                   * flag
     1689                                                                  * flag * rmask(j,i,sr)
    16871690
    16881691                   sums_l(k,69,tn) = sums_l(k,69,tn) - 0.5_wp * (              &
    16891692                   (km(k,j,i)+km(k+1,j,i)) * (e(k+1,j,i)-e(k,j,i)) * ddzu(k+1) &
    1690                                                                 )  * flag
     1693                                                                )  * flag * rmask(j,i,sr)
    16911694
    16921695                ENDDO
     
    20822085          hom(:,1,96,sr)  = sums(:,96)             ! qsws_veg
    20832086          hom(:,1,97,sr)  = sums(:,97)             ! r_a
    2084           hom(:,1,98,sr) = sums(:,98)              ! r_s
     2087          hom(:,1,98,sr)  = sums(:,98)             ! r_s
    20852088
    20862089       ENDIF
Note: See TracChangeset for help on using the changeset viewer.