Ignore:
Timestamp:
Jun 18, 2019 10:32:41 AM (5 years ago)
Author:
suehring
Message:

diagnostic output: Modularize diagnostic output, rename subroutines; formatting adjustments; allocate arrays only when required; add output of uu, vv, ww to enable variance calculation via temporal EC method; radiation: bugfix in masked data output; flow_statistics: Correct conversion to kinematic vertical scalar fluxes in case of pw-scheme and statistic regions

File:
1 edited

Legend:

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

    r3828 r4039  
    2525! -----------------
    2626! $Id$
     27! Correct conversion to kinematic scalar fluxes in case of pw-scheme and
     28! statistic regions
     29!
     30! 3828 2019-03-27 19:36:23Z raasch
    2731! unused variables removed
    2832!
     
    14971501                                    vpt(k+1,j,i) - hom(k+1,1,44,sr) )
    14981502                      sums_l(k,46,tn) = sums_l(k,46,tn) + pts * w(k,j,i) *     &
     1503                                               rho_air_zw(k) *                 &
    14991504                                               heatflux_output_conversion(k) * &
    15001505                                                          rmask(j,i,sr) * flag
     
    15091514                              hom(k+1,1,42,sr) )
    15101515                         sums_l(k,52,tn) = sums_l(k,52,tn) + pts * w(k,j,i) *  &
     1516                                             rho_air_zw(k) *                   &
    15111517                                             waterflux_output_conversion(k) *  &
    15121518                                                             rmask(j,i,sr)  *  &
     
    15381544                                          vpt(k+1,j,i) - hom(k+1,1,44,sr) )
    15391545                         sums_l(k,46,tn) = sums_l(k,46,tn) + pts * w(k,j,i) *  &
     1546                                              rho_air_zw(k) *                  &
    15401547                                              heatflux_output_conversion(k) *  &
    15411548                                                             rmask(j,i,sr)  *  &
     
    16851692                           ( pt(k,j,i)   - hom(k,1,4,sr) +                     &
    16861693                             pt(k+1,j,i) - hom(k+1,1,4,sr) )                   &
     1694                           * rho_air_zw(k)                                     &
    16871695                           * heatflux_output_conversion(k)                     &
    16881696                           * w(k,j,i) * rmask(j,i,sr) * flag
     
    16911699                                      q(k+1,j,i) - hom(k+1,1,41,sr) )
    16921700                      sums_l(k,49,tn) = sums_l(k,49,tn) + pts * w(k,j,i) *     &
     1701                                       rho_air_zw(k) *                         &
    16931702                                       waterflux_output_conversion(k) *        &
    16941703                                       rmask(j,i,sr) * flag
Note: See TracChangeset for help on using the changeset viewer.