Changeset 19 for palm/trunk/SOURCE/boundary_conds.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/boundary_conds.f90 (modified) (4 diffs)

palm/trunk/SOURCE/boundary_conds.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Boundary conditions for e(nzt), pt(nzt), and q(nzt) removed because these
+! values are now calculated by the prognostic equation,
+! Dirichlet and zero gradient condition for pt established at top boundary
 !
 ! Former revisions:
@@ -100 +102 @@
 !
 !--    Temperature at top boundary
-       IF ( ibc_pt_t == 1 )  THEN
-          pt(nzt,:,:)   = pt(nzt-1,:,:) + bc_pt_t_val * dzu(nzt)
+       IF ( ibc_pt_t == 0 )  THEN
+          IF ( timestep_scheme(1:5) /= 'runge' )  THEN
+             pt(nzt+1,:,:) = pt_m(nzt+1,:,:)
+          ELSE
+             pt(nzt+1,:,:) = pt_p(nzt+1,:,:)  ! pt_m not used for Runge-Kutta
+          ENDIF
+       ELSEIF ( ibc_pt_t == 1 )  THEN
+          pt(nzt+1,:,:) = pt(nzt,:,:)
+       ELSEIF ( ibc_pt_t == 2 )  THEN
           pt(nzt+1,:,:) = pt(nzt,:,:)   + bc_pt_t_val * dzu(nzt+1)
        ENDIF
@@ -114 +123 @@
              ENDDO
           ENDDO
-          e(nzt,:,:)   = e(nzt-1,:,:)
           e(nzt+1,:,:) = e(nzt,:,:)
        ENDIF
@@ -147 +155 @@
 !
 !--       Top boundary
-          q(nzt,:,:)   = q(nzt-1,:,:) + bc_q_t_val * dzu(nzt)
           q(nzt+1,:,:) = q(nzt,:,:)   + bc_q_t_val * dzu(nzt+1)
        ENDIF