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

Last change on this file since 1 was 1, checked in by raasch, 17 years ago

Initial repository layout and content

File size: 5.0 KB
Line 
1 SUBROUTINE print_1d
2
3!------------------------------------------------------------------------------!
4! Actual revisions:
5! -----------------
6!
7!
8! Former revisions:
9! -----------------
10! $Log: print_1d.f90,v $
11! Revision 1.11  2006/02/23 12:50:43  raasch
12! Preliminary no output of time-averaged data
13!
14! Revision 1.10  2001/03/30 07:48:24  raasch
15! Translation of remaining German identifiers (variables, subroutines, etc.)
16!
17! Revision 1.9  2001/01/22 20:51:40  letzel
18! All comments translated into English.
19!
20! Revision 1.8  2001/01/22 07:55:14  raasch
21! Module test_variables removed
22!
23! Revision 1.7  1998/07/06 12:31:17  raasch
24! Mittelungszeitraum wird direkt durch average_period_pr1d angegeben,
25! + USE test_variables
26!
27! Revision 1.6  1998/03/30 11:38:16  raasch
28! Ausgaben fuer die (Teil-) Gebiete
29!
30! Revision 1.5  1998/03/25 13:56:31  raasch
31! Erweiterungen fuer Gebietsstatistiken
32!
33! Revision 1.4  1998/03/18 20:18:33  raasch
34! Zusaetzliche Header werden ausgegeben
35!
36! Revision 1.3  1998/03/03 08:01:28  raasch
37! Aufruf von flow_statistics erfolgt jetzt gegebenenfalls hier
38! Ausgabe von zeitlich gemittelten Profilen
39!
40! Revision 1.2  1998/02/10 15:09:16  raasch
41! Ausdruck von Werten im w-Niveau
42!
43! Revision 1.1  1997/09/19 07:45:22  raasch
44! Initial revision
45!
46!
47! Description:
48! ------------
49! List output of 1D-profiles.
50!------------------------------------------------------------------------------!
51
52    USE arrays_3d
53    USE control_parameters
54    USE cpulog
55    USE indices
56    USE interfaces
57    USE pegrid
58    USE statistics
59
60    IMPLICIT NONE
61
62
63    CHARACTER (LEN=20) ::  period_chr
64    INTEGER ::  k, sr
65
66
67!
68!-- If required, compute statistics.
69    IF ( .NOT. flow_statistics_called )  CALL flow_statistics
70
71!
72!-- Flow_statistics has its own cpu-time measuring.
73    CALL cpu_log( log_point(18), 'print_1d', 'start' )
74
75    IF ( myid == 0 )  THEN
76!
77!--    Open file for list output of profiles.
78       CALL check_open( 16 )
79
80!
81!--    Prepare header.
82       period_chr = ' no time-average!'
83
84!
85!--    Output for the total domain (and each subregion, if applicable).
86       DO  sr = 0, statistic_regions
87!
88!--       Write header.
89          WRITE ( 16, 112 )
90          WRITE ( 16, 100 )  TRIM( run_description_header ) // '    ' // &
91                             TRIM( region( sr ) ), TRIM( period_chr ), 'uv'
92          WRITE ( 16, 105 )  TRIM( simulated_time_chr )
93!          ELSE
94!             WRITE ( 16, 106 )  TRIM( simulated_time_chr ),           &
95!                                averaging_interval_pr, average_count_pr
96!          ENDIF
97          WRITE ( 16, 111 )
98
99!
100!--       Output of values on the scalar levels.
101          WRITE ( 16, 120 )
102          WRITE ( 16, 111 )
103          DO  k = nzt+1, nzb, -1
104             WRITE ( 16, 121)  k, zu(k), hom(k,1,1,sr),           &
105                               hom(k,1,1,sr) - hom(k,1,5,sr),     &
106                               hom(k,1,2,sr),                     &
107                               hom(k,1,2,sr) - hom(k,1,6,sr),     &
108                               hom(k,1,4,sr),                     &
109                               hom(k,1,4,sr) - hom(k,1,7,sr),     &
110                               hom(k,1,8,sr), hom(k,1,9,sr),      &
111                               hom(k,1,10,sr), hom(k,1,11,sr), zu(k), k
112          ENDDO
113          WRITE ( 16, 111 )
114          WRITE ( 16, 120 )
115          WRITE ( 16, 111 )
116
117!
118!--       Output of values on the w-levels.
119          WRITE ( 16, 112 )
120          WRITE ( 16, 100 )  TRIM( run_description_header ) // '    ' // &
121                             TRIM( region( sr ) ), TRIM( period_chr ), 'w'
122          WRITE ( 16, 105 )  TRIM( simulated_time_chr )
123!          ELSE
124!             WRITE ( 16, 106 )  TRIM( simulated_time_chr ),           &
125!                                averaging_interval_pr, average_count_pr
126!          ENDIF
127          WRITE ( 16, 111 )
128
129          WRITE ( 16, 130 )
130          WRITE ( 16, 111 )
131          DO  k = nzt+1, nzb, -1
132             WRITE ( 16, 131)  k, zw(k), hom(k,1,16,sr),            &
133                               hom(k,1,18,sr), hom(k,1,12,sr), &
134                               hom(k,1,19,sr), hom(k,1,14,sr), &
135                               hom(k,1,20,sr), zw(k), k
136          ENDDO
137          WRITE ( 16, 111 )
138          WRITE ( 16, 130 )
139          WRITE ( 16, 111 )
140
141       ENDDO
142
143    ENDIF
144
145    CALL cpu_log( log_point(18), 'print_1d','stop', 'nobarrier' )
146
147!
148!-- Formats.
149100 FORMAT (1X,A/1X,10('-')/ &
150            ' Horizontally',A,' averaged profiles on the ',A,'-level')
151105 FORMAT (' Time: ',A)
152106 FORMAT (' Time: ',A,18X,'averaged over',F7.1,' s (',I4, &
153            ' Single times)')
154111 FORMAT (1X,131('-'))
155112 FORMAT (/)
156120 FORMAT ('   k     zu      u     du     v     dv     pt    dpt    ', &
157            'e      Km    Kh     l      zu      k')
158121 FORMAT (1X,I4,1X,F7.1,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,2X,F6.2,1X,F5.2, &
159            1X,F6.4,1X,F5.2,1X,F5.2,1X,F6.2,1X,F7.1,2X,I4)
160130 FORMAT ('   k     zw      w''pt''     wpt       w''u''      wu       ', &
161            ' w''v''      wv        zw      k')
162131 FORMAT (1X,I4,1X,F7.1,6(1X,E9.3),1X,F7.1,2X,I4)
163
164
165 END SUBROUTINE print_1d
Note: See TracBrowser for help on using the repository browser.