Changeset 1462 for palm


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/boundary_conds.f90

    r1410 r1462  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix: set dirichlet boundary condition for passive_scalar at model domain top
    2323!
    2424! Former revisions:
     
    118118        ONLY:  bc_pt_t_val, bc_q_t_val, constant_diffusion,                    &
    119119               cloud_physics, dt_3d, humidity,                                 &
    120                ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_sa_t, ibc_uv_b, ibc_uv_t,      &
    121                icloud_scheme, inflow_l, inflow_n, inflow_r, inflow_s,          &
     120               ibc_pt_b, ibc_pt_t, ibc_q_b, ibc_q_t, ibc_sa_t, ibc_uv_b,       &
     121               ibc_uv_t, icloud_scheme, inflow_l, inflow_n, inflow_r, inflow_s,&
    122122               intermediate_timestep_count, large_scale_forcing, ocean,        &
    123123               outflow_l, outflow_n, outflow_r, outflow_s, passive_scalar,     &
     
    288288!
    289289!--    Top boundary
    290        q_p(nzt+1,:,:) = q_p(nzt,:,:)   + bc_q_t_val * dzu(nzt+1)
     290       IF ( ibc_q_t == 0 ) THEN
     291          q_p(nzt+1,:,:) = q(nzt+1,:,:)
     292       ELSEIF ( ibc_q_t == 1 ) THEN
     293          q_p(nzt+1,:,:) = q_p(nzt,:,:)   + bc_q_t_val * dzu(nzt+1)
     294       ENDIF
    291295
    292296       IF ( cloud_physics  .AND.  icloud_scheme == 0  .AND.  precipitation )  THEN
Note: See TracChangeset for help on using the changeset viewer.