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

Timestamp:

Apr 17, 2020 4:14:16 PM (4 years ago)

Author:

schwenkel

Message:

Implementation of ice microphysics

File:

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

palm/trunk/SOURCE/advec_s_bc.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Implementation of ice microphysics
+! 
+! 4488 2020-04-03 11:34:29Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -957 +960 @@
           ENDDO
 
-       ELSEIF ( sk_char == 'qr' )  THEN
-
-!
-!--       Rain water content boundary condition at the bottom boundary: Dirichlet (fixed surface
-!--       rain water content).
+       ELSEIF ( sk_char == 'qi' )  THEN
+
+!
+!--       Ice crystal content boundary condition at the bottom boundary:
+!--       Dirichlet (fixed surface rain water content).
           DO  i = nxl, nxr
              DO  j = nys, nyn
@@ -971 +974 @@
 
 !
-!--       Rain water content boundary condition at the top boundary: Dirichlet
+!--       Ice crystal content boundary condition at the top boundary: Dirichlet
           DO  i = nxl, nxr
              DO  j = nys, nyn
@@ -979 +982 @@
           ENDDO
 
-       ELSEIF ( sk_char == 'nc' )  THEN
-
-!
-!--       Cloud drop concentration boundary condition at the bottom boundary: Dirichlet (fixed
-!--       surface cloud drop concentration).
+       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
@@ -993 +996 @@
 
 !
+!--       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 == 'nc' )  THEN
+
+!
+!--       Cloud drop concentration boundary condition at the bottom boundary: Dirichlet (fixed
+!--       surface cloud 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
+
+!
 !--       Cloud 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 == 'ni' )  THEN
+
+!
+!--       Ice crystal concentration boundary condition at the bottom boundary:
+!--       Dirichlet (fixed surface cloud 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
+
+!
+!--       Ice crystal concentration boundary condition at the top boundary: Dirichlet
           DO  i = nxl, nxr
              DO  j = nys, nyn