Changeset 1933 for palm/trunk/SOURCE/boundary_conds.f90
- Timestamp:
- Jun 13, 2016 7:12:51 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/boundary_conds.f90
r1823 r1933 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 22 ! 21 ! 22 ! 23 23 ! Former revisions: 24 24 ! ----------------- 25 25 ! $Id$ 26 ! 27 ! 1823 2016-04-07 08:57:52Z hoffmann 28 ! Initial version of purely vertical nesting introduced. 26 29 ! 27 30 ! 1822 2016-04-07 07:49:42Z hoffmann … … 159 162 USE pegrid 160 163 164 USE pmc_interface, & 165 ONLY : nesting_mode 166 161 167 162 168 IMPLICIT NONE … … 359 365 ! 360 366 !-- The same restoration for u at i=nxl and v at j=nys as above must be made 361 !-- in case of nest boundaries. Note however, that the above ELSEIF-structure is 362 !-- not appropriate here as there may be more than one nest boundary on a 363 !-- PE-domain. Furthermore Neumann conditions for SGS-TKE are not required here. 364 IF ( nest_bound_s ) THEN 365 v_p(:,nys,:) = v_p(:,nys-1,:) 366 ENDIF 367 IF ( nest_bound_l ) THEN 368 u_p(:,:,nxl) = u_p(:,:,nxl-1) 367 !-- in case of nest boundaries. This must not be done in case of vertical nesting 368 !-- mode as in that case the lateral boundaries are actually cyclic. 369 IF ( nesting_mode /= 'vertical' ) THEN 370 IF ( nest_bound_s ) THEN 371 v_p(:,nys,:) = v_p(:,nys-1,:) 372 ENDIF 373 IF ( nest_bound_l ) THEN 374 u_p(:,:,nxl) = u_p(:,:,nxl-1) 375 ENDIF 369 376 ENDIF 370 377
Note: See TracChangeset
for help on using the changeset viewer.