Changeset 19 for palm/trunk/SOURCE/flow_statistics.f90
 Timestamp:
 Feb 23, 2007 4:53:48 AM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/flow_statistics.f90
r4 r19 4 4 ! Actual revisions: 5 5 !  6 ! 6 ! fluxes at top modified (tswst, qswst) 7 7 ! 8 8 ! Former revisions: … … 347 347 + ( w(k,j,i)  w(k,j1,i) ) * ddy & 348 348 ) * rmask(j,i,sr) 349 ENDDO 350 351 DO k = nzb_diff_s_outer(j,i)1, nzt_diff 349 352 ! 350 353 ! Heat flux w"pt" … … 417 420 sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + & 418 421 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") 419 452 ENDIF 420 453 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.