Changeset 75 for palm/trunk/SOURCE/disturb_field.f90
- Timestamp:
- Mar 22, 2007 9:54:05 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/disturb_field.f90
r4 r75 1 SUBROUTINE disturb_field( nzb_uv_inner, dist1, field , xrp, ynp)1 SUBROUTINE disturb_field( nzb_uv_inner, dist1, field ) 2 2 3 3 !------------------------------------------------------------------------------! 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! xrp, ynp eliminated, 2nd+3rd argument removed from exchange horiz 7 7 ! 8 8 ! Former revisions: … … 36 36 IMPLICIT NONE 37 37 38 INTEGER :: i, j, k , xrp, ynp38 INTEGER :: i, j, k 39 39 INTEGER :: nzb_uv_inner(nys-1:nyn+1,nxl-1:nxr+1) 40 40 41 41 REAL :: randomnumber, & 42 42 dist1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 43 field(nzb:nzt+1,nys-1:nyn+ ynp+1,nxl-1:nxr+xrp+1)43 field(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) 44 44 REAL, DIMENSION(:,:,:), ALLOCATABLE :: dist2 45 45 … … 94 94 ! 95 95 !-- Exchange of ghost points for the random perturbation 96 CALL exchange_horiz( dist1 , 0, 0)96 CALL exchange_horiz( dist1 ) 97 97 98 98 ! … … 114 114 !-- Exchange of ghost points for the filtered perturbation. 115 115 !-- Afterwards, filter operation and exchange of ghost points are repeated. 116 CALL exchange_horiz( dist2 , 0, 0)116 CALL exchange_horiz( dist2 ) 117 117 DO i = nxl, nxr 118 118 DO j = nys, nyn … … 125 125 ENDDO 126 126 ENDDO 127 CALL exchange_horiz( dist1 , 0, 0)127 CALL exchange_horiz( dist1 ) 128 128 129 129 !
Note: See TracChangeset
for help on using the changeset viewer.