Ignore:
Timestamp:
Mar 26, 2013 6:16:16 PM (8 years ago)
Author:
hoffmann
Message:

optimization of two-moments cloud physics

File:
1 edited

Legend:

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

    r1112 r1115  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! ql is calculated by calc_liquid_water_content
    2323!
    2424! Former revisions:
     
    766766                                                          rmask(j,i,sr)
    767767                      IF ( .NOT. cloud_droplets )  THEN
    768                          pts = 0.5 *                                          &
    769                               ( ( q(k,j,i) - ql(k,j,i) ) -                    &
    770                               hom(k,1,42,sr) +                                &
    771                               ( q(k+1,j,i) - ql(k+1,j,i) ) -                  &
     768                         pts = 0.5 *                                           &
     769                              ( ( q(k,j,i) - ql(k,j,i) ) -                     &
     770                              hom(k,1,42,sr) +                                 &
     771                              ( q(k+1,j,i) - ql(k+1,j,i) ) -                   &
    772772                              hom(k+1,1,42,sr) )
    773                          sums_l(k,52,tn) = sums_l(k,52,tn) + pts * w(k,j,i) * &
     773                         sums_l(k,52,tn) = sums_l(k,52,tn) + pts * w(k,j,i) *  &
    774774                                                             rmask(j,i,sr)
    775775                         IF ( icloud_scheme == 0  )  THEN
    776                             sums_l(k,54,tn) = sums_l(k,54,tn) + ( ql(k,j,i) + &
    777                                                                 qr(k,j,i) ) * &
     776                            sums_l(k,54,tn) = sums_l(k,54,tn) + ql(k,j,i) *    &
    778777                                                                rmask(j,i,sr)
    779                             sums_l(k,73,tn) = sums_l(k,73,tn) + nr(k,j,i) *   &
     778                            sums_l(k,75,tn) = sums_l(k,75,tn) + qc(k,j,i) *    &
    780779                                                                rmask(j,i,sr)
    781                             sums_l(k,74,tn) = sums_l(k,74,tn) + qr(k,j,i) *   &
    782                                                                 rmask(j,i,sr)
    783                             sums_l(k,75,tn) = sums_l(k,75,tn) + ql(k,j,i) *   &
    784                                                                 rmask(j,i,sr)
    785                             sums_l(k,76,tn) = sums_l(k,76,tn) + prr(k,j,i) *  &
    786                                                                 rmask(j,i,sr)
     780                            IF ( precipitation )  THEN
     781                               sums_l(k,73,tn) = sums_l(k,73,tn) + nr(k,j,i) * &
     782                                                                   rmask(j,i,sr)
     783                               sums_l(k,74,tn) = sums_l(k,74,tn) + qr(k,j,i) * &
     784                                                                   rmask(j,i,sr)
     785                               sums_l(k,76,tn) = sums_l(k,76,tn) + prr(k,j,i) *&
     786                                                                   rmask(j,i,sr)
     787                            ENDIF
    787788                         ELSE
    788                             sums_l(k,54,tn) = sums_l(k,54,tn) + ql(k,j,i) *   &
     789                            sums_l(k,54,tn) = sums_l(k,54,tn) + ql(k,j,i) *    &
    789790                                                                rmask(j,i,sr)
    790791                         ENDIF
    791792                      ELSE
    792                          sums_l(k,54,tn) = sums_l(k,54,tn) + ql(k,j,i) *      &
     793                         sums_l(k,54,tn) = sums_l(k,54,tn) + ql(k,j,i) *       &
    793794                                                             rmask(j,i,sr)
    794795                      ENDIF
Note: See TracChangeset for help on using the changeset viewer.