Changeset 1320 for palm/trunk/SOURCE/disturb_heatflux.f90
- Timestamp:
- Mar 20, 2014 8:40:49 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/disturb_heatflux.f90
r1319 r1320 16 16 ! 17 17 ! Copyright 1997-2014 Leibniz Universitaet Hannover 18 !------------------------------------------------------------------------------ --!18 !------------------------------------------------------------------------------! 19 19 ! 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! ONLY-attribute added to USE-statements, 23 ! kind-parameters added to all INTEGER and REAL declaration statements, 24 ! kinds are defined in new module kinds, 25 ! old module precision_kind is removed, 26 ! revision history before 2012 removed, 27 ! comment fields (!:) to be used for variable explanations added to 28 ! all variable declaration statements 23 29 ! 24 30 ! Former revisions: … … 26 32 ! $Id$ 27 33 ! 28 ! 1318 2014-03-17 13:35:16Z raasch29 ! module interfaces removed30 !31 34 ! 1036 2012-10-22 13:43:42Z raasch 32 35 ! code put under GPL (PALM 3.9) 33 !34 ! 555 2010-09-07 07:32:53Z raasch35 ! Bugfix in if statement36 !37 ! RCS Log replace by Id keyword, revision history cleaned up38 !39 ! Revision 1.7 2006/08/04 14:35:07 raasch40 ! Additional parameter in function random_gauss which limits the range of the41 ! created random numbers, izuf renamed iran42 36 ! 43 37 ! Revision 1.1 1998/03/25 20:03:47 raasch … … 54 48 !------------------------------------------------------------------------------! 55 49 56 USE arrays_3d 57 USE control_parameters 58 USE cpulog 59 USE grid_variables 60 USE indices 50 USE arrays_3d, & 51 ONLY: shf 52 53 USE control_parameters, & 54 ONLY: iran, surface_heatflux, wall_heatflux 55 56 USE cpulog, & 57 ONLY: cpu_log, log_point 58 59 USE kinds 60 61 USE indices, & 62 ONLY: nx, nxl, nxr, ny, nyn, nys, nzb_s_inner 61 63 62 64 IMPLICIT NONE 63 65 64 INTEGER :: i, j 65 REAL :: random_gauss, randomnumber 66 INTEGER(iwp) :: j !: 67 INTEGER(iwp) :: i !: 68 69 REAL(wp) :: random_gauss !: 70 REAL(wp) :: randomnumber !: 66 71 67 72 … … 73 78 DO j = 0, ny 74 79 randomnumber = random_gauss( iran, 5.0 ) 75 IF ( nxl <= i .AND. nxr >= i .AND. nys <= j .AND. nyn >= j ) &80 IF ( nxl <= i .AND. nxr >= i .AND. nys <= j .AND. nyn >= j ) & 76 81 THEN 77 82 IF ( nzb_s_inner(j,i) == 0 ) THEN
Note: See TracChangeset
for help on using the changeset viewer.