Ignore:
Timestamp:
Feb 25, 2016 12:31:13 PM (5 years ago)
Author:
hellstea
Message:

Introduction of nested domain system

File:
1 edited

Legend:

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

    r1683 r1762  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Introduction of nested domain feature
    2222!
    2323! Former revisions:
     
    6464    USE control_parameters,                                                    &
    6565        ONLY:  bc_lr_cyc, bc_ns_cyc, ibc_p_b, ibc_p_t, inflow_l, inflow_n,     &
    66                inflow_r, inflow_s, n_sor, omega_sor, outflow_l, outflow_n,     &
     66               inflow_r, inflow_s, nest_bound_l, nest_bound_n, nest_bound_r,   &
     67               nest_bound_s, n_sor, omega_sor, outflow_l, outflow_n,           &
    6768               outflow_r, outflow_s
    6869
     
    154155!--    Horizontal (Neumann) boundary conditions in case of non-cyclic boundaries
    155156       IF ( .NOT. bc_lr_cyc )  THEN
    156           IF ( inflow_l .OR. outflow_l )  p(:,:,nxl-1) = p(:,:,nxl)
    157           IF ( inflow_r .OR. outflow_r )  p(:,:,nxr+1) = p(:,:,nxr)
     157          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l )  p(:,:,nxl-1) = p(:,:,nxl)
     158          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r )  p(:,:,nxr+1) = p(:,:,nxr)
    158159       ENDIF
    159160       IF ( .NOT. bc_ns_cyc )  THEN
    160           IF ( inflow_n .OR. outflow_n )  p(:,nyn+1,:) = p(:,nyn,:)
    161           IF ( inflow_s .OR. outflow_s )  p(:,nys-1,:) = p(:,nys,:)
     161          IF ( inflow_n .OR. outflow_n .OR. nest_bound_n )  p(:,nyn+1,:) = p(:,nyn,:)
     162          IF ( inflow_s .OR. outflow_s .OR. nest_bound_s )  p(:,nys-1,:) = p(:,nys,:)
    162163       ENDIF
    163164
     
    216217!--    Horizontal (Neumann) boundary conditions in case of non-cyclic boundaries
    217218       IF ( .NOT. bc_lr_cyc )  THEN
    218           IF ( inflow_l .OR. outflow_l )  p(:,:,nxl-1) = p(:,:,nxl)
    219           IF ( inflow_r .OR. outflow_r )  p(:,:,nxr+1) = p(:,:,nxr)
     219          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l )  p(:,:,nxl-1) = p(:,:,nxl)
     220          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r )  p(:,:,nxr+1) = p(:,:,nxr)
    220221       ENDIF
    221222       IF ( .NOT. bc_ns_cyc )  THEN
    222           IF ( inflow_n .OR. outflow_n )  p(:,nyn+1,:) = p(:,nyn,:)
    223           IF ( inflow_s .OR. outflow_s )  p(:,nys-1,:) = p(:,nys,:)
     223          IF ( inflow_n .OR. outflow_n .OR. nest_bound_n )  p(:,nyn+1,:) = p(:,nyn,:)
     224          IF ( inflow_s .OR. outflow_s .OR. nest_bound_s )  p(:,nys-1,:) = p(:,nys,:)
    224225       ENDIF
    225226
Note: See TracChangeset for help on using the changeset viewer.