Changeset 1361 for palm/trunk/SOURCE/advec_s_bc.f90

Timestamp:

Apr 16, 2014 3:17:48 PM (10 years ago)

Author:

hoffmann

Message:

improved version of two-moment cloud physics

File:

• palm/trunk/SOURCE/advec_s_bc.f90 (modified) (2 diffs)

palm/trunk/SOURCE/advec_s_bc.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! nr and qr added
 ! 
 ! Former revisions:
@@ -926 +926 @@
        ENDIF
 
+    ELSEIF ( sk_char == 'qr' )  THEN
+
+!
+!--    Rain water content boundary condition at the bottom boundary:
+!--    Dirichlet (fixed surface rain water content).
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
+             sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
+             sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
+          ENDDO
+       ENDDO
+
+!
+!--    Rain water content boundary condition at the top boundary: Dirichlet
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             sk_p(nzt+2,j,i)   = sk_p(nzt+1,j,i)
+             sk_p(nzt+3,j,i)   = sk_p(nzt+1,j,i)
+          ENDDO
+       ENDDO
+
+    ELSEIF ( sk_char == 'nr' )  THEN
+
+!
+!--    Rain drop concentration boundary condition at the bottom boundary:
+!--    Dirichlet (fixed surface rain drop concentration).
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             sk_p(nzb,j,i)   = sk_p(nzb+1,j,i)
+             sk_p(nzb-1,j,i) = sk_p(nzb,j,i)
+             sk_p(nzb-2,j,i) = sk_p(nzb,j,i)
+          ENDDO
+       ENDDO
+
+!
+!--    Rain drop concentration boundary condition at the top boundary: Dirichlet
+       DO  i = nxl, nxr
+          DO  j = nys, nyn
+             sk_p(nzt+2,j,i)   = sk_p(nzt+1,j,i)
+             sk_p(nzt+3,j,i)   = sk_p(nzt+1,j,i)
+          ENDDO
+       ENDDO
+
     ELSEIF ( sk_char == 'e' )  THEN