Changeset 1115 for palm/trunk/SOURCE/boundary_conds.f90

Timestamp:

Mar 26, 2013 6:16:16 PM (11 years ago)

Author:

hoffmann

Message:

optimization of two-moments cloud physics

File:

• palm/trunk/SOURCE/boundary_conds.f90 (modified) (6 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! boundary conditions of two-moment cloud scheme are restricted to Neumann-
+! boundary-conditions
 !  
 ! Former revisions:
@@ -246 +247 @@
        q_p(nzt+1,:,:) = q_p(nzt,:,:)   + bc_q_t_val * dzu(nzt+1)
 
-       IF ( cloud_physics .AND. icloud_scheme == 0 )  THEN
+       IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  &
+            precipitation )  THEN
 !             
-!--       Surface conditions for constant_humidity_flux
-          IF ( ibc_qr_b == 0 ) THEN
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   qr_p(nzb_s_inner(j,i),j,i) = qr(nzb_s_inner(j,i),j,i)
-                ENDDO
-             ENDDO
-          ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   qr_p(nzb_s_inner(j,i),j,i) = qr_p(nzb_s_inner(j,i)+1,j,i)
-                ENDDO
-             ENDDO
-          ENDIF
-!
-!--       Top boundary
-          qr_p(nzt+1,:,:) = qr_p(nzt,:,:) + bc_qr_t_val * dzu(nzt+1)
-!             
-!--       Surface conditions for constant_humidity_flux
-          IF ( ibc_nr_b == 0 ) THEN
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   nr_p(nzb_s_inner(j,i),j,i) = nr(nzb_s_inner(j,i),j,i)
-                ENDDO
-             ENDDO
-          ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   nr_p(nzb_s_inner(j,i),j,i) = nr_p(nzb_s_inner(j,i)+1,j,i)
-                ENDDO
-             ENDDO
-          ENDIF
-!
-!--       Top boundary
-          nr_p(nzt+1,:,:) = nr_p(nzt,:,:) + bc_nr_t_val * dzu(nzt+1)
+!--       Surface conditions rain water (Neumann)
+          DO  i = nxlg, nxrg
+             DO  j = nysg, nyng
+                qr_p(nzb_s_inner(j,i),j,i) = qr_p(nzb_s_inner(j,i)+1,j,i)
+                nr_p(nzb_s_inner(j,i),j,i) = nr_p(nzb_s_inner(j,i)+1,j,i)
+             ENDDO
+          ENDDO
+!
+!--       Top boundary condition for rain water (Neumann)
+          qr_p(nzt+1,:,:) = qr_p(nzt,:,:)
+          nr_p(nzt+1,:,:) = nr_p(nzt,:,:)
+           
        ENDIF
-
 !
 !--    In case of inflow at the south boundary the boundary for v is at nys
@@ -308 +286 @@
           pt_p(:,nys-1,:)     = pt_p(:,nys,:)
           IF ( .NOT. constant_diffusion     )  e_p(:,nys-1,:) = e_p(:,nys,:)
-          IF ( humidity .OR. passive_scalar )  THEN
+          IF ( humidity  .OR.  passive_scalar )  THEN
              q_p(:,nys-1,:) = q_p(:,nys,:)
-             IF ( cloud_physics .AND. icloud_scheme == 0 )  THEN
+             IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  &
+                  precipitation)  THEN
                 qr_p(:,nys-1,:) = qr_p(:,nys,:)
                 nr_p(:,nys-1,:) = nr_p(:,nys,:)
@@ -318 +297 @@
           pt_p(:,nyn+1,:)     = pt_p(:,nyn,:)
           IF ( .NOT. constant_diffusion     )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
-          IF ( humidity .OR. passive_scalar )  THEN
+          IF ( humidity  .OR.  passive_scalar )  THEN
              q_p(:,nyn+1,:) = q_p(:,nyn,:)
-             IF ( cloud_physics .AND. icloud_scheme == 0 )  THEN
+             IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  &
+                  precipitation )  THEN
                 qr_p(:,nyn+1,:) = qr_p(:,nyn,:)
                 nr_p(:,nyn+1,:) = nr_p(:,nyn,:)
@@ -328 +308 @@
           pt_p(:,:,nxl-1)     = pt_p(:,:,nxl)
           IF ( .NOT. constant_diffusion     )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
-          IF ( humidity .OR. passive_scalar )  THEN
+          IF ( humidity  .OR.  passive_scalar )  THEN
              q_p(:,:,nxl-1) = q_p(:,:,nxl)
-             IF ( cloud_physics .AND. icloud_scheme == 0 )  THEN
+             IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  &
+                  precipitation )  THEN
                 qr_p(:,:,nxl-1) = qr_p(:,:,nxl)
                 nr_p(:,:,nxl-1) = nr_p(:,:,nxl)
@@ -340 +321 @@
           IF ( humidity .OR. passive_scalar )  THEN
              q_p(:,:,nxr+1) = q_p(:,:,nxr)
-             IF ( cloud_physics .AND. icloud_scheme == 0 )  THEN
+             IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  precipitation )  THEN
                 qr_p(:,:,nxr+1) = qr_p(:,:,nxr)
                 nr_p(:,:,nxr+1) = nr_p(:,:,nxr)