Ignore:
Timestamp:
Aug 24, 2010 9:17:48 AM (11 years ago)
Author:
maronga
Message:

bugfix: calculation of subgrid-scale flux at nzt for use_top_fluxes=.T.

File:
1 edited

Legend:

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

    r531 r550  
    519519!--          Subgridscale fluxes at the top surface
    520520             IF ( use_top_fluxes )  THEN
    521                 sums_l(nzt,12,tn) = sums_l(nzt,12,tn) + &
     521                sums_l(nzt:nzt+1,12,tn) = sums_l(nzt:nzt+1,12,tn) + &
    522522                                    uswst(j,i) * rmask(j,i,sr)    ! w"u"
    523                 sums_l(nzt,14,tn) = sums_l(nzt,14,tn) + &
     523                sums_l(nzt:nzt+1,14,tn) = sums_l(nzt:nzt+1,14,tn) + &
    524524                                    vswst(j,i) * rmask(j,i,sr)    ! w"v"
    525                 sums_l(nzt,16,tn) = sums_l(nzt,16,tn) + &
     525                sums_l(nzt:nzt+1,16,tn) = sums_l(nzt:nzt+1,16,tn) + &
    526526                                    tswst(j,i)  * rmask(j,i,sr)   ! w"pt"
    527                 sums_l(nzt,58,tn) = sums_l(nzt,58,tn) + &
     527                sums_l(nzt:nzt+1,58,tn) = sums_l(nzt:nzt+1,58,tn) + &
    528528                                    0.0 * rmask(j,i,sr)           ! u"pt"
    529                 sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
    530                                     0.0 * rmask(j,i,sr)           ! v"pt"
     529                sums_l(nzt:nzt+1,61,tn) = sums_l(nzt:nzt+1,61,tn) + &
     530                                    0.0 * rmask(j,i,sr)           ! v"pt"
     531
    531532                IF ( ocean )  THEN
    532533                   sums_l(nzt,65,tn) = sums_l(nzt,65,tn) + &
Note: See TracChangeset for help on using the changeset viewer.