source: palm/trunk/SOURCE/disturb_heatflux.f90 @ 3421

Last change on this file since 3421 was 3241, checked in by raasch, 6 years ago

various changes to avoid compiler warnings (mainly removal of unused variables)

  • Property svn:keywords set to Id
File size: 5.5 KB
RevLine 
[1682]1!> @file disturb_heatflux.f90
[2000]2!------------------------------------------------------------------------------!
[2696]3! This file is part of the PALM model system.
[1036]4!
[2000]5! PALM is free software: you can redistribute it and/or modify it under the
6! terms of the GNU General Public License as published by the Free Software
7! Foundation, either version 3 of the License, or (at your option) any later
8! version.
[1036]9!
10! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
13!
14! You should have received a copy of the GNU General Public License along with
15! PALM. If not, see <http://www.gnu.org/licenses/>.
16!
[2718]17! Copyright 1997-2018 Leibniz Universitaet Hannover
[1320]18!------------------------------------------------------------------------------!
[1036]19!
[484]20! Current revisions:
[1322]21! ------------------
[1354]22!
[2233]23!
[1321]24! Former revisions:
25! -----------------
26! $Id: disturb_heatflux.f90 3241 2018-09-12 15:02:00Z gronemeier $
[3241]27! unused variable removed
28!
29! 2718 2018-01-02 08:49:38Z maronga
[2716]30! Corrected "Former revisions" section
31!
32! 2696 2017-12-14 17:12:51Z kanani
33! Change in file header (GPL part)
[1321]34!
[2716]35! 2233 2017-05-30 18:08:54Z suehring
36!
[2233]37! 2232 2017-05-30 17:47:52Z suehring
38! Adjustment according new surface data type.
39! Implemented parallel random number generator to obtain always the same
40! random number distribution regardless of the processor distribution.
41!
[2038]42! 2037 2016-10-26 11:15:40Z knoop
43! Anelastic approximation implemented
44!
[2001]45! 2000 2016-08-20 18:09:15Z knoop
46! Forced header and separation lines into 80 columns
47!
[1683]48! 1682 2015-10-07 23:56:08Z knoop
49! Code annotations made doxygen readable
50!
[1354]51! 1353 2014-04-08 15:21:23Z heinze
52! REAL constants provided with KIND-attribute
53!
[1323]54! 1322 2014-03-20 16:38:49Z raasch
55! REAL constants defined as wp-kind
56!
[1321]57! 1320 2014-03-20 08:40:49Z raasch
[1320]58! ONLY-attribute added to USE-statements,
59! kind-parameters added to all INTEGER and REAL declaration statements,
60! kinds are defined in new module kinds,
61! revision history before 2012 removed,
62! comment fields (!:) to be used for variable explanations added to
63! all variable declaration statements
[1]64!
[1037]65! 1036 2012-10-22 13:43:42Z raasch
66! code put under GPL (PALM 3.9)
67!
[1]68! Revision 1.1  1998/03/25 20:03:47  raasch
69! Initial revision
70!
71!
72! Description:
73! ------------
[1682]74!> Generate random, normally distributed heatflux values and store them as the
75!> near-surface heatflux.
[1]76!------------------------------------------------------------------------------!
[2232]77 SUBROUTINE disturb_heatflux( surf )
[1682]78 
[1]79
[1320]80    USE arrays_3d,                                                             &
[2232]81        ONLY:  heatflux_input_conversion
[1320]82       
83    USE control_parameters,                                                    &
[2232]84        ONLY:  iran, surface_heatflux, random_generator, wall_heatflux
[1320]85       
86    USE cpulog,                                                                &
87        ONLY:  cpu_log, log_point
88       
89    USE kinds
90   
91    USE indices,                                                               &
[2232]92        ONLY:  nzb
[1]93
[2232]94    USE random_generator_parallel,                                             &
95        ONLY:  random_number_parallel, random_seed_parallel, random_dummy,     &
[3241]96               seq_random_array
[2232]97
98    USE surface_mod,                                                           &
99        ONLY:  surf_type
100
[1]101    IMPLICIT NONE
102
[2232]103    INTEGER(iwp) ::  i  !< grid index, x direction
104    INTEGER(iwp) ::  j  !< grid index, y direction
105    INTEGER(iwp) ::  k  !< grid index, z direction
106    INTEGER(iwp) ::  m  !< loop variables over surface elements
[1320]107   
[1682]108    REAL(wp) ::  random_gauss  !<
109    REAL(wp) ::  randomnumber  !<
[1]110
[2232]111    TYPE(surf_type) ::  surf   !< surface-type variable
[1]112
[2232]113
[1]114    CALL cpu_log( log_point(23), 'disturb_heatflux', 'start' )
115
116!
[2232]117!-- Generate random disturbances and store them. Note, if
118!-- random_generator /= 'random-parallel' it is not guaranteed to obtain 
119!-- the same random distribution if the number of processors is changed.
120    IF ( random_generator /= 'random-parallel' )  THEN
121
122       DO  m = 1, surf%ns
123
124          k = surf%k(m)
125
126          randomnumber = random_gauss( iran, 5.0_wp )     
[1]127!
[2232]128!--       k-1 is topography top index. If this is 0, set surface heatflux. Over
129!--       topography surface_heatflux is replaced by wall_heatflux(0).
130          IF ( k-1 == 0 )  THEN
131             surf%shf(m) = randomnumber * surface_heatflux                     &
132                              * heatflux_input_conversion(nzb)
133          ELSE
134             surf%shf(m) = randomnumber * wall_heatflux(0)                     &
135                              * heatflux_input_conversion(k-1)
[1]136          ENDIF
137       ENDDO
[2232]138    ELSE
[1]139
[2232]140       DO  m = 1, surf%ns
141
142          i = surf%i(m)
143          j = surf%j(m)
144          k = surf%k(m)
145
146          CALL random_seed_parallel( put=seq_random_array(:, j, i) )
147          CALL random_number_parallel( random_dummy )   
[1]148!
[2232]149!--       k-1 is topography top index. If this is 0, set surface heatflux. Over
150!--       topography surface_heatflux is replaced by wall_heatflux(0).
151          IF ( k-1 == 0 )  THEN
152             surf%shf(m) = ( random_dummy - 0.5_wp ) * surface_heatflux        &
153                              * heatflux_input_conversion(nzb)
154          ELSE
155             surf%shf(m) = ( random_dummy - 0.5_wp ) * wall_heatflux(0)        &
156                              * heatflux_input_conversion(k-1)
157          ENDIF
[1]158
[2232]159          CALL random_seed_parallel( get=seq_random_array(:, j, i) )
160
161       ENDDO
162
163    ENDIF
164
[1]165    CALL cpu_log( log_point(23), 'disturb_heatflux', 'stop' )
166
167
168 END SUBROUTINE disturb_heatflux
Note: See TracBrowser for help on using the repository browser.