Ignore:
Timestamp:
Feb 23, 2007 4:53:48 AM (14 years ago)
Author:
raasch
Message:

preliminary version of modified boundary conditions at top

File:
1 edited

Legend:

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

    r4 r19  
    44! Actual revisions:
    55! -----------------
    6 !
     6! fluxes at top modified (tswst, qswst)
    77!
    88! Former revisions:
     
    347347                                 + ( w(k,j,i)   - w(k,j-1,i) ) * ddy           &
    348348                                                               ) * rmask(j,i,sr)
     349             ENDDO
     350
     351             DO  k = nzb_diff_s_outer(j,i)-1, nzt_diff
    349352!
    350353!--             Heat flux w"pt"
     
    417420                   sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + &
    418421                                       qsws(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     422                ENDIF
     423             ENDIF
     424
     425!
     426!--          Subgridscale fluxes at the top surface
     427             IF ( use_top_fluxes )  THEN
     428                sums_l(nzt,16,tn) = sums_l(nzt,16,tn) + &
     429                                    tswst(j,i)  * rmask(j,i,sr)   ! w"pt"
     430                sums_l(nzt,58,tn) = sums_l(nzt,58,tn) + &
     431                                    0.0 * rmask(j,i,sr)           ! u"pt"
     432                sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
     433                                    0.0 * rmask(j,i,sr)           ! v"pt"
     434                IF ( moisture )  THEN
     435                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
     436                                       qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     437                   IF ( cloud_physics )  THEN
     438                      sums_l(nzt,45,tn) = sums_l(nzt,45,tn) + (           &
     439                                          ( 1.0 + 0.61 * q(nzt,j,i) ) *   &
     440                                          tswst(j,i) + 0.61 * pt(nzt,j,i) * &
     441                                                     qsws(j,i)            &
     442                                                              )
     443!
     444!--                   Formula does not work if ql(nzb) /= 0.0
     445                      sums_l(nzt,51,tn) = sums_l(nzt,51,tn) + &   ! w"q" (w"qv")
     446                                          qswst(j,i) * rmask(j,i,sr)
     447                   ENDIF
     448                ENDIF
     449                IF ( passive_scalar )  THEN
     450                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
     451                                       qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
    419452                ENDIF
    420453             ENDIF
Note: See TracChangeset for help on using the changeset viewer.