source: palm/trunk/SOURCE/print_1d.f90 @ 3248

Last change on this file since 3248 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: 6.1 KB
Line 
1!> @file print_1d.f90
2!------------------------------------------------------------------------------!
3! This file is part of the PALM model system.
4!
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.
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!
17! Copyright 1997-2018 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: print_1d.f90 3241 2018-09-12 15:02:00Z sward $
27! unused format statement removed
28!
29! 2718 2018-01-02 08:49:38Z maronga
30! Corrected "Former revisions" section
31!
32! 2696 2017-12-14 17:12:51Z kanani
33! Change in file header (GPL part)
34!
35! 2101 2017-01-05 16:42:31Z suehring
36!
37! 2000 2016-08-20 18:09:15Z knoop
38! Forced header and separation lines into 80 columns
39!
40! 1697 2015-10-28 17:14:10Z raasch
41! small E- and F-FORMAT changes to avoid informative compiler messages about
42! insufficient field width
43!
44! 1682 2015-10-07 23:56:08Z knoop
45! Code annotations made doxygen readable
46!
47! 1320 2014-03-20 08:40:49Z raasch
48! ONLY-attribute added to USE-statements,
49! kind-parameters added to all INTEGER and REAL declaration statements,
50! kinds are defined in new module kinds,
51! old module precision_kind is removed,
52! revision history before 2012 removed,
53! comment fields (!:) to be used for variable explanations added to
54! all variable declaration statements
55!
56! 1318 2014-03-17 13:35:16Z raasch
57! barrier argument removed from cpu_log,
58! module interfaces removed
59!
60! 1036 2012-10-22 13:43:42Z raasch
61! code put under GPL (PALM 3.9)
62!
63! RCS Log replace by Id keyword, revision history cleaned up
64!
65! Revision 1.1  1997/09/19 07:45:22  raasch
66! Initial revision
67!
68!
69! Description:
70! ------------
71!> List output of 1D-profiles.
72!------------------------------------------------------------------------------!
73 SUBROUTINE print_1d
74 
75
76    USE arrays_3d,                                                             &
77        ONLY:  zu, zw
78
79    USE control_parameters,                                                    &
80        ONLY:  run_description_header, simulated_time_chr
81
82    USE cpulog,                                                                &
83        ONLY:  cpu_log, log_point
84
85    USE indices,                                                               &
86        ONLY:  nzb, nzt
87
88    USE kinds
89
90    USE pegrid
91
92    USE statistics,                                                            &
93        ONLY:  flow_statistics_called, hom, region, statistic_regions
94
95    IMPLICIT NONE
96
97
98    CHARACTER (LEN=20) ::  period_chr  !<
99
100    INTEGER(iwp) ::  k   !<
101    INTEGER(iwp) ::  sr  !<
102
103
104!
105!-- If required, compute statistics.
106    IF ( .NOT. flow_statistics_called )  CALL flow_statistics
107
108!
109!-- Flow_statistics has its own cpu-time measuring.
110    CALL cpu_log( log_point(18), 'print_1d', 'start' )
111
112    IF ( myid == 0 )  THEN
113!
114!--    Open file for list output of profiles.
115       CALL check_open( 16 )
116
117!
118!--    Prepare header.
119       period_chr = ' no time-average!'
120
121!
122!--    Output for the total domain (and each subregion, if applicable).
123       DO  sr = 0, statistic_regions
124!
125!--       Write header.
126          WRITE ( 16, 112 )
127          WRITE ( 16, 100 )  TRIM( run_description_header ) // '    ' // &
128                             TRIM( region( sr ) ), TRIM( period_chr ), 'uv'
129          WRITE ( 16, 105 )  TRIM( simulated_time_chr )
130          WRITE ( 16, 111 )
131
132!
133!--       Output of values on the scalar levels.
134          WRITE ( 16, 120 )
135          WRITE ( 16, 111 )
136          DO  k = nzt+1, nzb, -1
137             WRITE ( 16, 121)  k, zu(k), hom(k,1,1,sr),           &
138                               hom(k,1,1,sr) - hom(k,1,5,sr),     &
139                               hom(k,1,2,sr),                     &
140                               hom(k,1,2,sr) - hom(k,1,6,sr),     &
141                               hom(k,1,4,sr),                     &
142                               hom(k,1,4,sr) - hom(k,1,7,sr),     &
143                               hom(k,1,8,sr), hom(k,1,9,sr),      &
144                               hom(k,1,10,sr), hom(k,1,11,sr), zu(k), k
145          ENDDO
146          WRITE ( 16, 111 )
147          WRITE ( 16, 120 )
148          WRITE ( 16, 111 )
149
150!
151!--       Output of values on the w-levels.
152          WRITE ( 16, 112 )
153          WRITE ( 16, 100 )  TRIM( run_description_header ) // '    ' // &
154                             TRIM( region( sr ) ), TRIM( period_chr ), 'w'
155          WRITE ( 16, 105 )  TRIM( simulated_time_chr )
156          WRITE ( 16, 111 )
157
158          WRITE ( 16, 130 )
159          WRITE ( 16, 111 )
160          DO  k = nzt+1, nzb, -1
161             WRITE ( 16, 131)  k, zw(k), hom(k,1,16,sr),            &
162                               hom(k,1,18,sr), hom(k,1,12,sr), &
163                               hom(k,1,19,sr), hom(k,1,14,sr), &
164                               hom(k,1,20,sr), zw(k), k
165          ENDDO
166          WRITE ( 16, 111 )
167          WRITE ( 16, 130 )
168          WRITE ( 16, 111 )
169
170       ENDDO
171
172    ENDIF
173
174    CALL cpu_log( log_point(18), 'print_1d', 'stop' )
175
176!
177!-- Formats.
178100 FORMAT (1X,A/1X,10('-')/ &
179            ' Horizontally',A,' averaged profiles on the ',A,'-level')
180105 FORMAT (' Time: ',A)
181111 FORMAT (1X,131('-'))
182112 FORMAT (/)
183120 FORMAT ('   k     zu      u     du     v     dv     pt    dpt    ', &
184            ' e      Km    Kh     l      zu      k')
185121 FORMAT (1X,I4,1X,F7.1,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,2X,F6.2,1X,F5.2, &
186            1X,F7.4,1X,F5.2,1X,F5.2,1X,F6.2,1X,F7.1,2X,I4)
187130 FORMAT ('   k     zw       w''pt''      wpt        w''u''       wu       ',&
188            '  w''v''       wv        zw      k')
189131 FORMAT (1X,I4,1X,F7.1,6(1X,E10.3),1X,F7.1,2X,I4)
190
191
192 END SUBROUTINE print_1d
Note: See TracBrowser for help on using the repository browser.