source: palm/trunk/SOURCE/data_output_profiles.f90 @ 1267

Last change on this file since 1267 was 1107, checked in by raasch, 12 years ago

last commit documented

  • Property svn:keywords set to Id
File size: 11.9 KB
RevLine 
[1]1 SUBROUTINE data_output_profiles
2
[1036]3!--------------------------------------------------------------------------------!
4! This file is part of PALM.
5!
6! PALM is free software: you can redistribute it and/or modify it under the terms
7! of the GNU General Public License as published by the Free Software Foundation,
8! either version 3 of the License, or (at your option) any later 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-2012  Leibniz University Hannover
18!--------------------------------------------------------------------------------!
19!
[254]20! Current revisions:
[1]21! -----------------
[1107]22!
[392]23!
24! Former revisions:
25! -----------------
26! $Id: data_output_profiles.f90 1107 2013-03-04 06:23:14Z heinze $
27!
[1107]28! 1106 2013-03-04 05:31:38Z raasch
29! bugfix: initial time for preruns of coupled runs is output as
30! -coupling_start_time
31!
[1093]32! 1092 2013-02-02 11:24:22Z raasch
33! unused variables removed
34!
[1037]35! 1036 2012-10-22 13:43:42Z raasch
36! code put under GPL (PALM 3.9)
37!
[965]38! 964 2012-07-26 09:14:24Z raasch
39! code for profil-output removed
40!
[392]41! 345 2009-07-01 14:37:56Z heinze
[345]42! In case of restart runs without extension, initial profiles are not written
43! to NetCDF-file anymore.
[291]44! simulated_time in NetCDF output replaced by time_since_reference_point.
[263]45! Output of NetCDF messages with aid of message handling routine.
[254]46! Output of messages replaced by message handling routine.
[1]47!
[198]48! 197 2008-09-16 15:29:03Z raasch
49! Time coordinate t=0 stored on netcdf-file only if an output is required for
50! this time for at least one of the profiles
51!
52! February 2007
[3]53! RCS Log replace by Id keyword, revision history cleaned up
54!
[90]55! 87 2007-05-22 15:46:47Z raasch
56! var_hom renamed pr_palm
57!
[1]58! Revision 1.18  2006/08/16 14:27:04  raasch
59! PRINT* statements for testing removed
60!
61! Revision 1.1  1997/09/12 06:28:48  raasch
62! Initial revision
63!
64!
65! Description:
66! ------------
67! Plot output of 1D-profiles for PROFIL
68!------------------------------------------------------------------------------!
69
70    USE control_parameters
71    USE cpulog
72    USE indices
73    USE interfaces
74    USE netcdf_control
75    USE pegrid
76    USE profil_parameter
77    USE statistics
78
79    IMPLICIT NONE
80
81
[1092]82    INTEGER ::  i, sr
[1]83
84!
85!-- If required, compute statistics
86    IF ( .NOT. flow_statistics_called )  CALL flow_statistics
87
88!
89!-- Flow_statistics has its own CPU time measurement
90    CALL cpu_log( log_point(15), 'data_output_profiles', 'start' )
91
92!
93!-- If required, compute temporal average
94    IF ( averaging_interval_pr == 0.0 )  THEN
95       hom_sum(:,:,:) = hom(:,1,:,:)
96    ELSE
97       IF ( average_count_pr > 0 )  THEN
98          hom_sum = hom_sum / REAL( average_count_pr )
99       ELSE
100!
101!--       This case may happen if dt_dopr is changed in the d3par-list of
102!--       a restart run
103          RETURN
104       ENDIF
105    ENDIF
106
107   
108    IF ( myid == 0 )  THEN
109
110!
111!--    Plot-output for each (sub-)region
112
113!
114!--    Open file for profile output in NetCDF format
115       IF ( netcdf_output )  THEN
116          CALL check_open( 104 )
117       ENDIF
118
119!
120!--    Increment the counter for number of output times
121       dopr_time_count = dopr_time_count + 1
122
123!
124!--    Output of initial profiles
125       IF ( dopr_time_count == 1 )  THEN
[345]126       
127          IF ( .NOT. output_for_t0 ) THEN
[1]128
[345]129             IF ( netcdf_output )  THEN
130#if defined( __netcdf )         
[1]131!
[1106]132!--             Store initial time to time axis, but only if an output
133!--             is required for at least one of the profiles. The initial time
134!--             is either 0, or, in case of a prerun for coupled atmosphere-ocean
135!--             runs, has a negative value
[345]136                DO  i = 1, dopr_n
[197]137                IF ( dopr_initial_index(i) /= 0 )  THEN
138                   nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr,  &
[1106]139                                           (/ -coupling_start_time /), &
140                                           start = (/ 1 /), count = (/ 1 /) )
[345]141                      CALL handle_netcdf_error( 'data_output_profiles', 329 )
142                      output_for_t0 = .TRUE.
143                      EXIT
144                   ENDIF
145                ENDDO
[1]146
147!
[345]148!--             Store normalization factors
149                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), & ! wpt0
150                                     (/ hom_sum(nzb,18,normalizing_region) /), &
151                                        start = (/ 1 /), count = (/ 1 /) )
152                CALL handle_netcdf_error( 'data_output_profiles', 330 )
[1]153
[345]154                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2
155                           (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
156                                        start = (/ 1 /), count = (/ 1 /) )
157                CALL handle_netcdf_error( 'data_output_profiles', 331 )
158                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
159                           (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]160                                     start = (/ 1 /), count = (/ 1 /) )
[345]161                CALL handle_netcdf_error( 'data_output_profiles', 332 )
162                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
163                           (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
164                                        start = (/ 1 /), count = (/ 1 /) )
165                CALL handle_netcdf_error( 'data_output_profiles', 333 )
[1]166
[345]167                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), &!ws2tsw
168                           (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 *   &
169                              hom_sum(nzb+3,pr_palm,normalizing_region)    /), &
170                                        start = (/ 1 /), count = (/ 1 /) )
171                CALL handle_netcdf_error( 'data_output_profiles', 334 )
[1]172
[345]173                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), &!wstsw2
174                           (/ hom_sum(nzb+8,pr_palm,normalizing_region) *      &
175                              hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
176                                        start = (/ 1 /), count = (/ 1 /) )
177                CALL handle_netcdf_error( 'data_output_profiles', 335 )
[1]178
[345]179                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
180                              (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
181                                        start = (/ 1 /), count = (/ 1 /) )
182                CALL handle_netcdf_error( 'data_output_profiles', 336 )
183             
[1]184#endif
[345]185             ENDIF
[1]186!
[345]187!--          Loop over all 1D variables
188             DO  i = 1, dopr_n
[1]189
[345]190                IF ( dopr_initial_index(i) /= 0 )  THEN
[1]191
192!
[345]193!--                Output for the individual (sub-)regions
194                   DO  sr = 0, statistic_regions
[1]195
[345]196                      IF ( netcdf_output )  THEN
[1]197#if defined( __netcdf )
198!
[345]199!--                      Write data to netcdf file
200                         nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr),    &
201                                       hom(nzb:nzt+1,1,dopr_initial_index(i),sr), &
202                                                 start = (/ 1, 1 /),              &
203                                                 count = (/ nzt-nzb+2, 1 /) )
204                         CALL handle_netcdf_error( 'data_output_profiles', 337 )
[1]205#endif
[345]206                      ENDIF
[1]207
[345]208                   ENDDO
[1]209
[345]210                ENDIF   ! Initial profile available
[1]211
[345]212             ENDDO   ! Loop over dopr_n for initial profiles
[1]213
[345]214             IF ( netcdf_output  .AND.  output_for_t0 )  THEN
215                dopr_time_count = dopr_time_count + 1
216             ENDIF
[1]217
[345]218          END IF
[1]219       ENDIF   ! Initial profiles
220
221       IF ( netcdf_output )  THEN
222#if defined( __netcdf )
[345]223
[1]224!
225!--       Store time to time axis         
[291]226          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr,        &
227                                  (/ time_since_reference_point /), &
228                                  start = (/ dopr_time_count /),    &
[1]229                                  count = (/ 1 /) )
[263]230          CALL handle_netcdf_error( 'data_output_profiles', 338 )
[1]231
232!
233!--       Store normalization factors
234          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), &  ! wpt0
235                                  (/ hom_sum(nzb,18,normalizing_region) /), &
236                                  start = (/ dopr_time_count /),               &
237                                  count = (/ 1 /) )
[263]238          CALL handle_netcdf_error( 'data_output_profiles', 339 )
[1]239
240          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), &  ! ws2
[87]241                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
[1]242                                  start = (/ dopr_time_count /),               &
243                                  count = (/ 1 /) )
[263]244          CALL handle_netcdf_error( 'data_output_profiles', 340 )
[1]245
246          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), &  ! tsw2
[87]247                        (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]248                                  start = (/ dopr_time_count /),               &
249                                  count = (/ 1 /) )
[263]250          CALL handle_netcdf_error( 'data_output_profiles', 341 )
[1]251
252          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), &  ! ws3
[87]253                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
[1]254                                  start = (/ dopr_time_count /),               &
255                                  count = (/ 1 /) )
[263]256          CALL handle_netcdf_error( 'data_output_profiles', 342 )
[1]257
258          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), &  ! ws2tsw
[87]259                        (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 *   &
260                           hom_sum(nzb+3,pr_palm,normalizing_region)    /), &
[1]261                                  start = (/ dopr_time_count /),               &
262                                  count = (/ 1 /) )
[263]263          CALL handle_netcdf_error( 'data_output_profiles', 343 )
264         
[1]265          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), &  ! wstsw2
[87]266                        (/ hom_sum(nzb+8,pr_palm,normalizing_region) *      &
267                           hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]268                                  start = (/ dopr_time_count /),               &
269                                  count = (/ 1 /) )
[263]270          CALL handle_netcdf_error( 'data_output_profiles', 344 )
[1]271
272          nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), &  ! z_i
[87]273                           (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
[1]274                                  start = (/ dopr_time_count /),               &
275                                  count = (/ 1 /) )
[263]276          CALL handle_netcdf_error( 'data_output_profiles', 345 )
[1]277#endif
278       ENDIF
279
280!
281!--    Output of the individual (non-initial) profiles
282       DO  i = 1, dopr_n
283
284!
285!--       Output for the individual (sub-)domains
286          DO  sr = 0, statistic_regions
287
288             IF ( netcdf_output )  THEN
289#if defined( __netcdf )
290!
291!--             Write data to netcdf file
292                nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr),          &
293                                        hom_sum(nzb:nzt+1,dopr_index(i),sr),&
294                                        start = (/ 1, dopr_time_count /),      &
295                                        count = (/ nzt-nzb+2, 1 /) )
[263]296                CALL handle_netcdf_error( 'data_output_profiles', 346 )
[1]297#endif
298             ENDIF
299
300          ENDDO
301
[964]302       ENDDO
[1]303
304    ENDIF  ! Output on PE0
305
306!
307!-- If averaging has been done above, the summation counter must be re-set.
308    IF ( averaging_interval_pr /= 0.0 )  THEN
309       average_count_pr = 0
310    ENDIF
311
312    CALL cpu_log( log_point(15), 'data_output_profiles','stop', 'nobarrier' )
313
314!
315!-- Formats
316100 FORMAT ('#1 ',A,1X,A)
317101 FORMAT (E15.7,1X,E15.7)
318102 FORMAT ('NEXT')
319
320 END SUBROUTINE data_output_profiles
Note: See TracBrowser for help on using the repository browser.