Ignore:
Timestamp:
Dec 14, 2017 5:12:51 PM (4 years ago)
Author:
kanani
Message:

Merge of branch palm4u into trunk

Location:
palm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk

  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/sor.f90

    r2101 r2696  
    11!> @file sor.f90
    22!------------------------------------------------------------------------------!
    3 ! This file is part of PALM.
     3! This file is part of the PALM model system.
    44!
    55! PALM is free software: you can redistribute it and/or modify it under the
     
    2525! -----------------
    2626! $Id$
     27! - Large-scale forcing implemented (MS)
     28!
     29! 2101 2017-01-05 16:42:31Z suehring
    2730!
    2831! 2037 2016-10-26 11:15:40Z knoop
     
    7578
    7679    USE control_parameters,                                                    &
    77         ONLY:  bc_lr_cyc, bc_ns_cyc, ibc_p_b, ibc_p_t, inflow_l, inflow_n,     &
    78                inflow_r, inflow_s, nest_bound_l, nest_bound_n, nest_bound_r,   &
    79                nest_bound_s, n_sor, omega_sor, outflow_l, outflow_n,           &
    80                outflow_r, outflow_s
     80        ONLY:  bc_lr_cyc, bc_ns_cyc, force_bound_l, force_bound_n,             &
     81               force_bound_r, force_bound_s, ibc_p_b, ibc_p_t, inflow_l,       &
     82               inflow_n, inflow_r, inflow_s, nest_bound_l, nest_bound_n,       &
     83               nest_bound_r, nest_bound_s, n_sor, omega_sor, outflow_l,        &
     84               outflow_n, outflow_r, outflow_s
    8185
    8286    IMPLICIT NONE
     
    167171!--    Horizontal (Neumann) boundary conditions in case of non-cyclic boundaries
    168172       IF ( .NOT. bc_lr_cyc )  THEN
    169           IF ( inflow_l .OR. outflow_l .OR. nest_bound_l )  p(:,:,nxl-1) = p(:,:,nxl)
    170           IF ( inflow_r .OR. outflow_r .OR. nest_bound_r )  p(:,:,nxr+1) = p(:,:,nxr)
     173          IF ( inflow_l      .OR.  outflow_l  .OR.                                 &
     174               nest_bound_l  .OR.  force_bound_l )  p(:,:,nxl-1) = p(:,:,nxl)
     175          IF ( inflow_r      .OR.  outflow_r  .OR.                                 &
     176               nest_bound_r  .OR.  force_bound_r )  p(:,:,nxr+1) = p(:,:,nxr)
    171177       ENDIF
    172178       IF ( .NOT. bc_ns_cyc )  THEN
    173           IF ( inflow_n .OR. outflow_n .OR. nest_bound_n )  p(:,nyn+1,:) = p(:,nyn,:)
    174           IF ( inflow_s .OR. outflow_s .OR. nest_bound_s )  p(:,nys-1,:) = p(:,nys,:)
     179          IF ( inflow_n      .OR.  outflow_n  .OR.                                 &
     180               nest_bound_n  .OR.  force_bound_n )  p(:,nyn+1,:) = p(:,nyn,:)
     181          IF ( inflow_s      .OR.  outflow_s  .OR.                                 &
     182               nest_bound_s  .OR.  force_bound_s )  p(:,nys-1,:) = p(:,nys,:)
    175183       ENDIF
    176184
     
    229237!--    Horizontal (Neumann) boundary conditions in case of non-cyclic boundaries
    230238       IF ( .NOT. bc_lr_cyc )  THEN
    231           IF ( inflow_l .OR. outflow_l .OR. nest_bound_l )  p(:,:,nxl-1) = p(:,:,nxl)
    232           IF ( inflow_r .OR. outflow_r .OR. nest_bound_r )  p(:,:,nxr+1) = p(:,:,nxr)
     239          IF ( inflow_l      .OR.  outflow_l  .OR.                             &
     240               nest_bound_l  .OR.  force_bound_l )  p(:,:,nxl-1) = p(:,:,nxl)
     241          IF ( inflow_r      .OR.  outflow_r  .OR.                             &
     242               nest_bound_r  .OR.  force_bound_r )  p(:,:,nxr+1) = p(:,:,nxr)
    233243       ENDIF
    234244       IF ( .NOT. bc_ns_cyc )  THEN
    235           IF ( inflow_n .OR. outflow_n .OR. nest_bound_n )  p(:,nyn+1,:) = p(:,nyn,:)
    236           IF ( inflow_s .OR. outflow_s .OR. nest_bound_s )  p(:,nys-1,:) = p(:,nys,:)
     245          IF ( inflow_n      .OR.  outflow_n  .OR.                             &
     246               nest_bound_n  .OR.  force_bound_n )  p(:,nyn+1,:) = p(:,nyn,:)
     247          IF ( inflow_s      .OR.  outflow_s  .OR.                             &
     248               nest_bound_s  .OR.  force_bound_s )  p(:,nys-1,:) = p(:,nys,:)
    237249       ENDIF
    238250
Note: See TracChangeset for help on using the changeset viewer.