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

Timestamp:

Sep 4, 2014 12:03:25 PM (10 years ago)

Author:

suehring

Message:

Bugfix: set dirichlet boundary condition for passive_scalar at model domain top

File:

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

palm/trunk/SOURCE/boundary_conds.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: set dirichlet boundary condition for passive_scalar at model domain top 
 ! 
 ! Former revisions:
@@ -118 +118 @@
         ONLY:  bc_pt_t_val, bc_q_t_val, constant_diffusion,                    &
                cloud_physics, dt_3d, humidity,                                 &
-               ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_sa_t, ibc_uv_b, ibc_uv_t,      &
-               icloud_scheme, inflow_l, inflow_n, inflow_r, inflow_s,          &
+               ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_sa_t, ibc_uv_b,       &
+               ibc_uv_t, icloud_scheme, inflow_l, inflow_n, inflow_r, inflow_s,&
                intermediate_timestep_count, large_scale_forcing, ocean,        &
                outflow_l, outflow_n, outflow_r, outflow_s, passive_scalar,     &
@@ -288 +288 @@
 !
 !--    Top boundary
-       q_p(nzt+1,:,:) = q_p(nzt,:,:)   + bc_q_t_val * dzu(nzt+1)
+       IF ( ibc_q_t == 0 ) THEN
+          q_p(nzt+1,:,:) = q(nzt+1,:,:)
+       ELSEIF ( ibc_q_t == 1 ) THEN
+          q_p(nzt+1,:,:) = q_p(nzt,:,:)   + bc_q_t_val * dzu(nzt+1)
+       ENDIF
 
        IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  precipitation )  THEN