Changeset 1682 for palm/trunk/SOURCE/disturb_field.f90
- Timestamp:
- Oct 7, 2015 11:56:08 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/disturb_field.f90
r1426 r1682 1 SUBROUTINE disturb_field( nzb_uv_inner, dist1, field ) 2 1 !> @file disturb_field.f90 3 2 !--------------------------------------------------------------------------------! 4 3 ! This file is part of PALM. … … 20 19 ! Current revisions: 21 20 ! ------------------ 22 ! 21 ! Code annotations made doxygen readable 23 22 ! 24 23 ! Former revisions: … … 55 54 ! Description: 56 55 ! ------------ 57 ! Imposing a random perturbation on a 3D-array.58 ! On parallel computers, the random number generator is as well called for all59 ! gridpoints of the total domain to ensure, regardless of the number of PEs60 ! used, that the elements of the array have the same values in the same61 ! order in every case. The perturbation range is steered by dist_range.56 !> Imposing a random perturbation on a 3D-array. 57 !> On parallel computers, the random number generator is as well called for all 58 !> gridpoints of the total domain to ensure, regardless of the number of PEs 59 !> used, that the elements of the array have the same values in the same 60 !> order in every case. The perturbation range is steered by dist_range. 62 61 !------------------------------------------------------------------------------! 62 SUBROUTINE disturb_field( nzb_uv_inner, dist1, field ) 63 63 64 64 65 USE control_parameters, & … … 85 86 IMPLICIT NONE 86 87 87 INTEGER(iwp) :: i ! :88 INTEGER(iwp) :: j ! :89 INTEGER(iwp) :: k ! :90 91 INTEGER(iwp) :: nzb_uv_inner(nysg:nyng,nxlg:nxrg) ! :92 93 REAL(wp) :: randomnumber ! :94 95 REAL(wp) :: dist1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ! :96 REAL(wp) :: field(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ! :97 98 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dist2 ! :88 INTEGER(iwp) :: i !< 89 INTEGER(iwp) :: j !< 90 INTEGER(iwp) :: k !< 91 92 INTEGER(iwp) :: nzb_uv_inner(nysg:nyng,nxlg:nxrg) !< 93 94 REAL(wp) :: randomnumber !< 95 96 REAL(wp) :: dist1(nzb:nzt+1,nysg:nyng,nxlg:nxrg) !< 97 REAL(wp) :: field(nzb:nzt+1,nysg:nyng,nxlg:nxrg) !< 98 99 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: dist2 !< 99 100 100 101
Note: See TracChangeset
for help on using the changeset viewer.