Changeset 19 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Feb 23, 2007 4:53:48 AM (17 years ago)

Author:

raasch

Message:

preliminary version of modified boundary conditions at top

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! fluxes at top modified (tswst, qswst)
 !
 ! Former revisions:
@@ -347 +347 @@
                                  + ( w(k,j,i)   - w(k,j-1,i) ) * ddy           &
                                                                ) * rmask(j,i,sr)
+             ENDDO
+
+             DO  k = nzb_diff_s_outer(j,i)-1, nzt_diff
 !
 !--             Heat flux w"pt"
@@ -417 +420 @@
                    sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + &
                                        qsws(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
+                ENDIF
+             ENDIF
+
+!
+!--          Subgridscale fluxes at the top surface
+             IF ( use_top_fluxes )  THEN
+                sums_l(nzt,16,tn) = sums_l(nzt,16,tn) + &
+                                    tswst(j,i)  * rmask(j,i,sr)   ! w"pt"
+                sums_l(nzt,58,tn) = sums_l(nzt,58,tn) + &
+                                    0.0 * rmask(j,i,sr)           ! u"pt"
+                sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
+                                    0.0 * rmask(j,i,sr)           ! v"pt"
+                IF ( moisture )  THEN
+                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
+                                       qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
+                   IF ( cloud_physics )  THEN
+                      sums_l(nzt,45,tn) = sums_l(nzt,45,tn) + (           &
+                                          ( 1.0 + 0.61 * q(nzt,j,i) ) *   &
+                                          tswst(j,i) + 0.61 * pt(nzt,j,i) * &
+                                                     qsws(j,i)            &
+                                                              )
+!
+!--                   Formula does not work if ql(nzb) /= 0.0
+                      sums_l(nzt,51,tn) = sums_l(nzt,51,tn) + &   ! w"q" (w"qv")
+                                          qswst(j,i) * rmask(j,i,sr)
+                   ENDIF
+                ENDIF
+                IF ( passive_scalar )  THEN
+                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
+                                       qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
                 ENDIF
              ENDIF