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

Last change on this file since 4 was 4, checked in by raasch, 15 years ago

Id keyword set as property for all *.f90 files

  • Property svn:keywords set to Id
File size: 4.2 KB
Line 
1 SUBROUTINE print_1d
2
3!------------------------------------------------------------------------------!
4! Actual revisions:
5! -----------------
6!
7!
8! Former revisions:
9! -----------------
10! $Id: print_1d.f90 4 2007-02-13 11:33:16Z raasch $
11! RCS Log replace by Id keyword, revision history cleaned up
12!
13! Revision 1.11  2006/02/23 12:50:43  raasch
14! Preliminary no output of time-averaged data
15!
16! Revision 1.1  1997/09/19 07:45:22  raasch
17! Initial revision
18!
19!
20! Description:
21! ------------
22! List output of 1D-profiles.
23!------------------------------------------------------------------------------!
24
25    USE arrays_3d
26    USE control_parameters
27    USE cpulog
28    USE indices
29    USE interfaces
30    USE pegrid
31    USE statistics
32
33    IMPLICIT NONE
34
35
36    CHARACTER (LEN=20) ::  period_chr
37    INTEGER ::  k, sr
38
39
40!
41!-- If required, compute statistics.
42    IF ( .NOT. flow_statistics_called )  CALL flow_statistics
43
44!
45!-- Flow_statistics has its own cpu-time measuring.
46    CALL cpu_log( log_point(18), 'print_1d', 'start' )
47
48    IF ( myid == 0 )  THEN
49!
50!--    Open file for list output of profiles.
51       CALL check_open( 16 )
52
53!
54!--    Prepare header.
55       period_chr = ' no time-average!'
56
57!
58!--    Output for the total domain (and each subregion, if applicable).
59       DO  sr = 0, statistic_regions
60!
61!--       Write header.
62          WRITE ( 16, 112 )
63          WRITE ( 16, 100 )  TRIM( run_description_header ) // '    ' // &
64                             TRIM( region( sr ) ), TRIM( period_chr ), 'uv'
65          WRITE ( 16, 105 )  TRIM( simulated_time_chr )
66!          ELSE
67!             WRITE ( 16, 106 )  TRIM( simulated_time_chr ),           &
68!                                averaging_interval_pr, average_count_pr
69!          ENDIF
70          WRITE ( 16, 111 )
71
72!
73!--       Output of values on the scalar levels.
74          WRITE ( 16, 120 )
75          WRITE ( 16, 111 )
76          DO  k = nzt+1, nzb, -1
77             WRITE ( 16, 121)  k, zu(k), hom(k,1,1,sr),           &
78                               hom(k,1,1,sr) - hom(k,1,5,sr),     &
79                               hom(k,1,2,sr),                     &
80                               hom(k,1,2,sr) - hom(k,1,6,sr),     &
81                               hom(k,1,4,sr),                     &
82                               hom(k,1,4,sr) - hom(k,1,7,sr),     &
83                               hom(k,1,8,sr), hom(k,1,9,sr),      &
84                               hom(k,1,10,sr), hom(k,1,11,sr), zu(k), k
85          ENDDO
86          WRITE ( 16, 111 )
87          WRITE ( 16, 120 )
88          WRITE ( 16, 111 )
89
90!
91!--       Output of values on the w-levels.
92          WRITE ( 16, 112 )
93          WRITE ( 16, 100 )  TRIM( run_description_header ) // '    ' // &
94                             TRIM( region( sr ) ), TRIM( period_chr ), 'w'
95          WRITE ( 16, 105 )  TRIM( simulated_time_chr )
96!          ELSE
97!             WRITE ( 16, 106 )  TRIM( simulated_time_chr ),           &
98!                                averaging_interval_pr, average_count_pr
99!          ENDIF
100          WRITE ( 16, 111 )
101
102          WRITE ( 16, 130 )
103          WRITE ( 16, 111 )
104          DO  k = nzt+1, nzb, -1
105             WRITE ( 16, 131)  k, zw(k), hom(k,1,16,sr),            &
106                               hom(k,1,18,sr), hom(k,1,12,sr), &
107                               hom(k,1,19,sr), hom(k,1,14,sr), &
108                               hom(k,1,20,sr), zw(k), k
109          ENDDO
110          WRITE ( 16, 111 )
111          WRITE ( 16, 130 )
112          WRITE ( 16, 111 )
113
114       ENDDO
115
116    ENDIF
117
118    CALL cpu_log( log_point(18), 'print_1d','stop', 'nobarrier' )
119
120!
121!-- Formats.
122100 FORMAT (1X,A/1X,10('-')/ &
123            ' Horizontally',A,' averaged profiles on the ',A,'-level')
124105 FORMAT (' Time: ',A)
125106 FORMAT (' Time: ',A,18X,'averaged over',F7.1,' s (',I4, &
126            ' Single times)')
127111 FORMAT (1X,131('-'))
128112 FORMAT (/)
129120 FORMAT ('   k     zu      u     du     v     dv     pt    dpt    ', &
130            'e      Km    Kh     l      zu      k')
131121 FORMAT (1X,I4,1X,F7.1,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,2X,F6.2,1X,F5.2, &
132            1X,F6.4,1X,F5.2,1X,F5.2,1X,F6.2,1X,F7.1,2X,I4)
133130 FORMAT ('   k     zw      w''pt''     wpt       w''u''      wu       ', &
134            ' w''v''      wv        zw      k')
135131 FORMAT (1X,I4,1X,F7.1,6(1X,E9.3),1X,F7.1,2X,I4)
136
137
138 END SUBROUTINE print_1d
Note: See TracBrowser for help on using the repository browser.