Ignore:
Timestamp:
Jul 27, 2007 9:09:17 AM (14 years ago)
Author:
raasch
Message:

preliminary version for coupled runs

File:
1 edited

Legend:

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

    r98 r102  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Prescribed momentum fluxes at the top surface are used
    77!
    88! Former revisions:
     
    381381!--          ----  strictly speaking the following k-loop would have to be
    382382!--                split up according to the staggered grid.
    383              DO  k = nzb_diff_s_outer(j,i)-1, nzt
     383             DO  k = nzb_diff_s_outer(j,i)-1, nzt_diff
    384384!
    385385!--             Momentum flux w"u"
     
    398398                                 + ( w(k,j,i)   - w(k,j-1,i) ) * ddy           &
    399399                                                               ) * rmask(j,i,sr)
    400              ENDDO
    401 
    402              DO  k = nzb_diff_s_outer(j,i)-1, nzt_diff
    403400!
    404401!--             Heat flux w"pt"
     
    490487!--          Subgridscale fluxes at the top surface
    491488             IF ( use_top_fluxes )  THEN
     489                sums_l(nzt,12,tn) = sums_l(nzt,12,tn) + &
     490                                    uswst(j,i) * rmask(j,i,sr)    ! w"u"
     491                sums_l(nzt,14,tn) = sums_l(nzt,14,tn) + &
     492                                    vswst(j,i) * rmask(j,i,sr)    ! w"v"
    492493                sums_l(nzt,16,tn) = sums_l(nzt,16,tn) + &
    493494                                    tswst(j,i)  * rmask(j,i,sr)   ! w"pt"
Note: See TracChangeset for help on using the changeset viewer.