Changeset 420 for palm/trunk/SOURCE/disturb_field.f90
- Timestamp:
- Jan 13, 2010 3:10:53 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/disturb_field.f90
r77 r420 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! A loop has been splitted to make runs reproducible on HLRN systems 7 7 ! 8 8 ! Former revisions: … … 98 98 ! 99 99 !-- Exchange of ghost points for the random perturbation 100 100 101 CALL exchange_horiz( dist1 ) 101 102 … … 104 105 !-- Neighboured grid points in all three directions are used for the 105 106 !-- filter operation. 106 DO i = nxl, nxr 107 DO j = nys, nyn 107 !-- Loop has been splitted to make runs reproducible on HLRN systems using 108 !-- compiler option -O3 109 DO i = nxl, nxr 110 DO j = nys, nyn 108 111 DO k = disturbance_level_ind_b-1, disturbance_level_ind_t+1 109 dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1) + dist1(k,j-1,i) & 110 + dist1(k,j+1,i) + dist1(k+1,j,i) + dist1(k-1,j,i) & 112 dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1) & 113 + dist1(k,j+1,i) + dist1(k+1,j,i) & 114 ) / 12.0 115 ENDDO 116 DO k = disturbance_level_ind_b-1, disturbance_level_ind_t+1 117 dist2(k,j,i) = dist2(k,j,i) + ( dist1(k,j-1,i) + dist1(k-1,j,i) & 111 118 + 6.0 * dist1(k,j,i) & 112 119 ) / 12.0 113 120 ENDDO 114 ENDDO115 ENDDO121 ENDDO 122 ENDDO 116 123 117 124 ! … … 119 126 !-- Afterwards, filter operation and exchange of ghost points are repeated. 120 127 CALL exchange_horiz( dist2 ) 128 121 129 DO i = nxl, nxr 122 130 DO j = nys, nyn … … 129 137 ENDDO 130 138 ENDDO 139 131 140 CALL exchange_horiz( dist1 ) 132 141
Note: See TracChangeset
for help on using the changeset viewer.