source: palm/trunk/SOURCE/header.f90 @ 767

Last change on this file since 767 was 767, checked in by raasch, 11 years ago

New:
---

Flow field initialization with given (e.g. measured) profiles. Profile data
for u-,v-velocity components + respective heights are given with new
inipar-parameters u_profile, v_profile, and uv_heights. Final profiles are
calculated from these given profiles by linear interpolation.
(check_parameters, header, init_3d_model, modules, parin)

Changed:


ug,vg replaced by u_init,v_init as the Dirichlet top boundary condition
(boundary_conds)

dirichlet_0 conditions moved from init_3d_model to
check_parameters (check_parameters, init_3d_model)

Errors:


bugfix: dirichlet_0 conditions moved from init_3d_model to
check_parameters (check_parameters, init_3d_model)

  • Property svn:keywords set to Id
File size: 77.5 KB
Line 
1 SUBROUTINE header
2
3!------------------------------------------------------------------------------!
4! Current revisions:
5! -----------------
6! Output of given initial u,v-profiles
7!
8! Former revisions:
9! -----------------
10! $Id: header.f90 767 2011-10-14 06:39:12Z raasch $
11!
12! 759 2011-09-15 13:58:31Z raasch
13! output of maximum number of parallel io streams
14!
15! 707 2011-03-29 11:39:40Z raasch
16! bc_lr/ns replaced by bc_lr/ns_cyc
17!
18! 667 2010-12-23 12:06:00Z suehring/gryschka
19! Output of advection scheme.
20! Modified output of Prandtl-layer height.
21!
22! 580 2010-10-05 13:59:11Z heinze
23! Renaming of ws_vertical_gradient to subs_vertical_gradient,
24! ws_vertical_gradient_level to subs_vertical_gradient_level and
25! ws_vertical_gradient_level_ind to subs_vertical_gradient_level_i
26!
27! 493 2010-03-01 08:30:24Z raasch
28! NetCDF data output format extendend for NetCDF4/HDF5
29!
30! 449 2010-02-02 11:23:59Z raasch
31! +large scale vertical motion (subsidence/ascent)
32! Bugfix: index problem concerning gradient_level indices removed
33!
34! 410 2009-12-04 17:05:40Z letzel
35! Masked data output: + dt_domask, mask_01~20_x|y|z, mask_01~20_x|y|z_loop,
36! mask_scale|_x|y|z, masks, skip_time_domask
37!
38! 346 2009-07-06 10:13:41Z raasch
39! initializing_actions='read_data_for_recycling' renamed to 'cyclic_fill'
40! Coupling with independent precursor runs.
41! Output of messages replaced by message handling routine.
42! Output of several additional dvr parameters
43! +canyon_height, canyon_width_x, canyon_width_y, canyon_wall_left,
44! canyon_wall_south, conserve_volume_flow_mode, dp_external, dp_level_b,
45! dp_smooth, dpdxy, u_bulk, v_bulk
46! topography_grid_convention moved from user_header
47! small bugfix concerning 3d 64bit netcdf output format
48!
49! 206 2008-10-13 14:59:11Z raasch
50! Bugfix: error in zu index in case of section_xy = -1
51!
52! 198 2008-09-17 08:55:28Z raasch
53! Format adjustments allowing output of larger revision numbers
54!
55! 197 2008-09-16 15:29:03Z raasch
56! allow 100 spectra levels instead of 10 for consistency with
57! define_netcdf_header,
58! bugfix in the output of the characteristic levels of potential temperature,
59! geostrophic wind, scalar concentration, humidity and leaf area density,
60! output of turbulence recycling informations
61!
62! 138 2007-11-28 10:03:58Z letzel
63! Allow new case bc_uv_t = 'dirichlet_0' for channel flow.
64! Allow two instead of one digit to specify isosurface and slicer variables.
65! Output of sorting frequency of particles
66!
67! 108 2007-08-24 15:10:38Z letzel
68! Output of informations for coupled model runs (boundary conditions etc.)
69! + output of momentumfluxes at the top boundary
70! Rayleigh damping for ocean, e_init
71!
72! 97 2007-06-21 08:23:15Z raasch
73! Adjustments for the ocean version.
74! use_pt_reference renamed use_reference
75!
76! 87 2007-05-22 15:46:47Z raasch
77! Bugfix: output of use_upstream_for_tke
78!
79! 82 2007-04-16 15:40:52Z raasch
80! Preprocessor strings for different linux clusters changed to "lc",
81! routine local_flush is used for buffer flushing
82!
83! 76 2007-03-29 00:58:32Z raasch
84! Output of netcdf_64bit_3d, particles-package is now part of the default code,
85! output of the loop optimization method, moisture renamed humidity,
86! output of subversion revision number
87!
88! 19 2007-02-23 04:53:48Z raasch
89! Output of scalar flux applied at top boundary
90!
91! RCS Log replace by Id keyword, revision history cleaned up
92!
93! Revision 1.63  2006/08/22 13:53:13  raasch
94! Output of dz_max
95!
96! Revision 1.1  1997/08/11 06:17:20  raasch
97! Initial revision
98!
99!
100! Description:
101! ------------
102! Writing a header with all important informations about the actual run.
103! This subroutine is called three times, two times at the beginning
104! (writing information on files RUN_CONTROL and HEADER) and one time at the
105! end of the run, then writing additional information about CPU-usage on file
106! header.
107!-----------------------------------------------------------------------------!
108
109    USE arrays_3d
110    USE control_parameters
111    USE cloud_parameters
112    USE cpulog
113    USE dvrp_variables
114    USE grid_variables
115    USE indices
116    USE model_1d
117    USE particle_attributes
118    USE pegrid
119    USE subsidence_mod
120    USE spectrum
121
122    IMPLICIT NONE
123
124    CHARACTER (LEN=1)  ::  prec
125    CHARACTER (LEN=2)  ::  do2d_mode
126    CHARACTER (LEN=5)  ::  section_chr
127    CHARACTER (LEN=9)  ::  time_to_string
128    CHARACTER (LEN=10) ::  coor_chr, host_chr
129    CHARACTER (LEN=16) ::  begin_chr
130    CHARACTER (LEN=23) ::  ver_rev
131    CHARACTER (LEN=40) ::  output_format
132    CHARACTER (LEN=70) ::  char1, char2, dopr_chr, &
133                           do2d_xy, do2d_xz, do2d_yz, do3d_chr, &
134                           domask_chr, run_classification
135    CHARACTER (LEN=86) ::  coordinates, gradients, learde, slices,  &
136                           temperatures, ugcomponent, vgcomponent
137    CHARACTER (LEN=85) ::  roben, runten
138
139    CHARACTER (LEN=1), DIMENSION(1:3) ::  dir = (/ 'x', 'y', 'z' /)
140
141    INTEGER ::  av, bh, blx, bly, bxl, bxr, byn, bys, ch, count, cwx, cwy,  &
142         cxl, cxr, cyn, cys, dim, i, ihost, io, j, l, ll, m, mpi_type
143    REAL    ::  cpuseconds_per_simulated_second
144
145!
146!-- Open the output file. At the end of the simulation, output is directed
147!-- to unit 19.
148    IF ( ( runnr == 0 .OR. force_print_header )  .AND. &
149         .NOT. simulated_time_at_begin /= simulated_time )  THEN
150       io = 15   !  header output on file RUN_CONTROL
151    ELSE
152       io = 19   !  header output on file HEADER
153    ENDIF
154    CALL check_open( io )
155
156!
157!-- At the end of the run, output file (HEADER) will be rewritten with
158!-- new informations
159    IF ( io == 19 .AND. simulated_time_at_begin /= simulated_time ) REWIND( 19 )
160
161!
162!-- Determine kind of model run
163    IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
164       run_classification = '3D - restart run'
165    ELSEIF ( TRIM( initializing_actions ) == 'cyclic_fill' )  THEN
166       run_classification = '3D - run with cyclic fill of 3D - prerun data'
167    ELSEIF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 )  THEN
168       run_classification = '3D - run without 1D - prerun'
169    ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
170       run_classification = '3D - run with 1D - prerun'
171    ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 )  THEN
172       run_classification = '3D - run initialized by user'
173    ELSE
174       message_string = ' unknown action(s): ' // TRIM( initializing_actions )
175       CALL message( 'header', 'PA0191', 0, 0, 0, 6, 0 )
176    ENDIF
177    IF ( ocean )  THEN
178       run_classification = 'ocean - ' // run_classification
179    ELSE
180       run_classification = 'atmosphere - ' // run_classification
181    ENDIF
182
183!
184!-- Run-identification, date, time, host
185    host_chr = host(1:10)
186    ver_rev = TRIM( version ) // '  ' // TRIM( revision )
187    WRITE ( io, 100 )  ver_rev, TRIM( run_classification )
188    IF ( TRIM( coupling_mode ) /= 'uncoupled' )  THEN
189#if defined( __mpi2 )
190       mpi_type = 2
191#else
192       mpi_type = 1
193#endif
194       WRITE ( io, 101 )  mpi_type, coupling_mode
195    ENDIF
196    WRITE ( io, 102 )  run_date, run_identifier, run_time, runnr, &
197                       ADJUSTR( host_chr )
198#if defined( __parallel )
199    IF ( npex == -1  .AND.  pdims(2) /= 1 )  THEN
200       char1 = 'calculated'
201    ELSEIF ( ( host(1:3) == 'ibm'  .OR.  host(1:3) == 'nec'  .OR.  &
202               host(1:2) == 'lc' )  .AND.                          &
203             npex == -1  .AND.  pdims(2) == 1 )  THEN
204       char1 = 'forced'
205    ELSE
206       char1 = 'predefined'
207    ENDIF
208    IF ( threads_per_task == 1 )  THEN
209       WRITE ( io, 103 )  numprocs, pdims(1), pdims(2), TRIM( char1 )
210    ELSE
211       WRITE ( io, 104 )  numprocs*threads_per_task, numprocs, &
212                          threads_per_task, pdims(1), pdims(2), TRIM( char1 )
213    ENDIF
214    IF ( ( host(1:3) == 'ibm'  .OR.  host(1:3) == 'nec'  .OR.    &
215           host(1:2) == 'lc'   .OR.  host(1:3) == 'dec' )  .AND. &
216         npex == -1  .AND.  pdims(2) == 1 )                      &
217    THEN
218       WRITE ( io, 106 )
219    ELSEIF ( pdims(2) == 1 )  THEN
220       WRITE ( io, 107 )  'x'
221    ELSEIF ( pdims(1) == 1 )  THEN
222       WRITE ( io, 107 )  'y'
223    ENDIF
224    IF ( use_seperate_pe_for_dvrp_output )  WRITE ( io, 105 )
225    IF ( numprocs /= maximum_parallel_io_streams )  THEN
226       WRITE ( io, 108 )  maximum_parallel_io_streams
227    ENDIF
228#endif
229    WRITE ( io, 99 )
230
231!
232!-- Numerical schemes
233    WRITE ( io, 110 )
234    IF ( psolver(1:7) == 'poisfft' )  THEN
235       WRITE ( io, 111 )  TRIM( fft_method )
236       IF ( psolver == 'poisfft_hybrid' )  WRITE ( io, 138 )
237    ELSEIF ( psolver == 'sor' )  THEN
238       WRITE ( io, 112 )  nsor_ini, nsor, omega_sor
239    ELSEIF ( psolver == 'multigrid' )  THEN
240       WRITE ( io, 135 )  cycle_mg, maximum_grid_level, ngsrb
241       IF ( mg_cycles == -1 )  THEN
242          WRITE ( io, 140 )  residual_limit
243       ELSE
244          WRITE ( io, 141 )  mg_cycles
245       ENDIF
246       IF ( mg_switch_to_pe0_level == 0 )  THEN
247          WRITE ( io, 136 )  nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
248                             nzt_mg(1)
249       ELSEIF (  mg_switch_to_pe0_level /= -1 )  THEN
250          WRITE ( io, 137 )  mg_switch_to_pe0_level,            &
251                             mg_loc_ind(2,0)-mg_loc_ind(1,0)+1, &
252                             mg_loc_ind(4,0)-mg_loc_ind(3,0)+1, &
253                             nzt_mg(mg_switch_to_pe0_level),    &
254                             nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
255                             nzt_mg(1)
256       ENDIF
257    ENDIF
258    IF ( call_psolver_at_all_substeps  .AND. timestep_scheme(1:5) == 'runge' ) &
259    THEN
260       WRITE ( io, 142 )
261    ENDIF
262
263    IF ( momentum_advec == 'pw-scheme' )  THEN
264       WRITE ( io, 113 )
265    ELSEIF (momentum_advec == 'ws-scheme' ) THEN
266       WRITE ( io, 503 )
267    ELSEIF (momentum_advec == 'ups-scheme' ) THEN
268       WRITE ( io, 114 )
269       IF ( cut_spline_overshoot )  WRITE ( io, 124 )
270       IF ( overshoot_limit_u /= 0.0  .OR.  overshoot_limit_v /= 0.0  .OR. &
271            overshoot_limit_w /= 0.0 )  THEN
272          WRITE ( io, 127 )  overshoot_limit_u, overshoot_limit_v, &
273                             overshoot_limit_w
274       ENDIF
275       IF ( ups_limit_u /= 0.0  .OR.  ups_limit_v /= 0.0  .OR. &
276            ups_limit_w /= 0.0 )                               &
277       THEN
278          WRITE ( io, 125 )  ups_limit_u, ups_limit_v, ups_limit_w
279       ENDIF
280       IF ( long_filter_factor /= 0.0 )  WRITE ( io, 115 )  long_filter_factor
281    ENDIF
282    IF ( scalar_advec == 'pw-scheme' )  THEN
283       WRITE ( io, 116 )
284    ELSEIF ( scalar_advec == 'ws-scheme' )  THEN
285       WRITE ( io, 504 )
286    ELSEIF ( scalar_advec == 'ups-scheme' )  THEN
287       WRITE ( io, 117 )
288       IF ( cut_spline_overshoot )  WRITE ( io, 124 )
289       IF ( overshoot_limit_e /= 0.0  .OR.  overshoot_limit_pt /= 0.0 )  THEN
290          WRITE ( io, 128 )  overshoot_limit_e, overshoot_limit_pt
291       ENDIF
292       IF ( ups_limit_e /= 0.0  .OR.  ups_limit_pt /= 0.0 )  THEN
293          WRITE ( io, 126 )  ups_limit_e, ups_limit_pt
294       ENDIF
295    ELSE
296       WRITE ( io, 118 )
297    ENDIF
298
299    WRITE ( io, 139 )  TRIM( loop_optimization )
300
301    IF ( galilei_transformation )  THEN
302       IF ( use_ug_for_galilei_tr )  THEN
303          char1 = 'geostrophic wind'
304       ELSE
305          char1 = 'mean wind in model domain'
306       ENDIF
307       IF ( simulated_time_at_begin == simulated_time )  THEN
308          char2 = 'at the start of the run'
309       ELSE
310          char2 = 'at the end of the run'
311       ENDIF
312       WRITE ( io, 119 )  TRIM( char1 ), TRIM( char2 ), &
313                          advected_distance_x/1000.0, advected_distance_y/1000.0
314    ENDIF
315    IF ( timestep_scheme == 'leapfrog' )  THEN
316       WRITE ( io, 120 )
317    ELSEIF ( timestep_scheme == 'leapfrog+euler' )  THEN
318       WRITE ( io, 121 )
319    ELSE
320       WRITE ( io, 122 )  timestep_scheme
321    ENDIF
322    IF ( use_upstream_for_tke )  WRITE ( io, 143 )
323    IF ( rayleigh_damping_factor /= 0.0 )  THEN
324       IF ( .NOT. ocean )  THEN
325          WRITE ( io, 123 )  'above', rayleigh_damping_height, &
326               rayleigh_damping_factor
327       ELSE
328          WRITE ( io, 123 )  'below', rayleigh_damping_height, &
329               rayleigh_damping_factor
330       ENDIF
331    ENDIF
332    IF ( humidity )  THEN
333       IF ( .NOT. cloud_physics )  THEN
334          WRITE ( io, 129 )
335       ELSE
336          WRITE ( io, 130 )
337          WRITE ( io, 131 )
338          IF ( radiation )      WRITE ( io, 132 )
339          IF ( precipitation )  WRITE ( io, 133 )
340       ENDIF
341    ENDIF
342    IF ( passive_scalar )  WRITE ( io, 134 )
343    IF ( conserve_volume_flow )  THEN
344       WRITE ( io, 150 )  conserve_volume_flow_mode
345       IF ( TRIM( conserve_volume_flow_mode ) == 'bulk_velocity' )  THEN
346          WRITE ( io, 151 )  u_bulk, v_bulk
347       ENDIF
348    ELSEIF ( dp_external )  THEN
349       IF ( dp_smooth )  THEN
350          WRITE ( io, 152 )  dpdxy, dp_level_b, ', vertically smoothed.'
351       ELSE
352          WRITE ( io, 152 )  dpdxy, dp_level_b, '.'
353       ENDIF
354    ENDIF
355    IF ( large_scale_subsidence )  THEN
356        WRITE ( io, 153 )
357        WRITE ( io, 154 )
358    ENDIF
359    WRITE ( io, 99 )
360
361!
362!-- Runtime and timestep informations
363    WRITE ( io, 200 )
364    IF ( .NOT. dt_fixed )  THEN
365       WRITE ( io, 201 )  dt_max, cfl_factor
366    ELSE
367       WRITE ( io, 202 )  dt
368    ENDIF
369    WRITE ( io, 203 )  simulated_time_at_begin, end_time
370
371    IF ( time_restart /= 9999999.9  .AND. &
372         simulated_time_at_begin == simulated_time )  THEN
373       IF ( dt_restart == 9999999.9 )  THEN
374          WRITE ( io, 204 )  ' Restart at:       ',time_restart
375       ELSE
376          WRITE ( io, 205 )  ' Restart at:       ',time_restart, dt_restart
377       ENDIF
378    ENDIF
379
380    IF ( simulated_time_at_begin /= simulated_time )  THEN
381       i = MAX ( log_point_s(10)%counts, 1 )
382       IF ( ( simulated_time - simulated_time_at_begin ) == 0.0 )  THEN
383          cpuseconds_per_simulated_second = 0.0
384       ELSE
385          cpuseconds_per_simulated_second = log_point_s(10)%sum / &
386                                            ( simulated_time -    &
387                                              simulated_time_at_begin )
388       ENDIF
389       WRITE ( io, 206 )  simulated_time, log_point_s(10)%sum, &
390                          log_point_s(10)%sum / REAL( i ),     &
391                          cpuseconds_per_simulated_second
392       IF ( time_restart /= 9999999.9  .AND.  time_restart < end_time )  THEN
393          IF ( dt_restart == 9999999.9 )  THEN
394             WRITE ( io, 204 )  ' Next restart at:  ',time_restart
395          ELSE
396             WRITE ( io, 205 )  ' Next restart at:  ',time_restart, dt_restart
397          ENDIF
398       ENDIF
399    ENDIF
400
401!
402!-- Start time for coupled runs, if independent precursor runs for atmosphere
403!-- and ocean are used. In this case, coupling_start_time defines the time
404!-- when the coupling is switched on.
405    IF ( coupling_start_time /= 0.0 )  THEN
406       IF ( coupling_start_time >= simulated_time_at_begin )  THEN
407          char1 = 'Precursor run for a coupled atmosphere-ocean run'
408       ELSE
409          char1 = 'Coupled atmosphere-ocean run following independent ' // &
410                  'precursor runs'
411       ENDIF
412       WRITE ( io, 207 )  char1, coupling_start_time
413    ENDIF
414
415!
416!-- Computational grid
417    IF ( .NOT. ocean )  THEN
418       WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(nzt+1)
419       IF ( dz_stretch_level_index < nzt+1 )  THEN
420          WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
421                             dz_stretch_factor, dz_max
422       ENDIF
423    ELSE
424       WRITE ( io, 250 )  dx, dy, dz, (nx+1)*dx, (ny+1)*dy, zu(0)
425       IF ( dz_stretch_level_index > 0 )  THEN
426          WRITE ( io, 252 )  dz_stretch_level, dz_stretch_level_index, &
427                             dz_stretch_factor, dz_max
428       ENDIF
429    ENDIF
430    WRITE ( io, 254 )  nx, ny, nzt+1, MIN( nnx, nx+1 ), MIN( nny, ny+1 ), &
431                       MIN( nnz+2, nzt+2 )
432    IF ( numprocs > 1 )  THEN
433       IF ( nxa == nx  .AND.  nya == ny  .AND.  nza == nz )  THEN
434          WRITE ( io, 255 )
435       ELSE
436          WRITE ( io, 256 )  nnx-(nxa-nx), nny-(nya-ny), nzt+2
437       ENDIF
438    ENDIF
439    IF ( sloping_surface )  WRITE ( io, 260 )  alpha_surface
440
441!
442!-- Topography
443    WRITE ( io, 270 )  topography
444    SELECT CASE ( TRIM( topography ) )
445
446       CASE ( 'flat' )
447          ! no actions necessary
448
449       CASE ( 'single_building' )
450          blx = INT( building_length_x / dx )
451          bly = INT( building_length_y / dy )
452          bh  = INT( building_height / dz )
453
454          IF ( building_wall_left == 9999999.9 )  THEN
455             building_wall_left = ( nx + 1 - blx ) / 2 * dx
456          ENDIF
457          bxl = INT ( building_wall_left / dx + 0.5 )
458          bxr = bxl + blx
459
460          IF ( building_wall_south == 9999999.9 )  THEN
461             building_wall_south = ( ny + 1 - bly ) / 2 * dy
462          ENDIF
463          bys = INT ( building_wall_south / dy + 0.5 )
464          byn = bys + bly
465
466          WRITE ( io, 271 )  building_length_x, building_length_y, &
467                             building_height, bxl, bxr, bys, byn
468
469       CASE ( 'single_street_canyon' )
470          ch  = NINT( canyon_height / dz )
471          IF ( canyon_width_x /= 9999999.9 )  THEN
472!
473!--          Street canyon in y direction
474             cwx = NINT( canyon_width_x / dx )
475             IF ( canyon_wall_left == 9999999.9 )  THEN
476                canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
477             ENDIF
478             cxl = NINT( canyon_wall_left / dx )
479             cxr = cxl + cwx
480             WRITE ( io, 272 )  'y', canyon_height, ch, 'u', cxl, cxr
481
482          ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
483!
484!--          Street canyon in x direction
485             cwy = NINT( canyon_width_y / dy )
486             IF ( canyon_wall_south == 9999999.9 )  THEN
487                canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
488             ENDIF
489             cys = NINT( canyon_wall_south / dy )
490             cyn = cys + cwy
491             WRITE ( io, 272 )  'x', canyon_height, ch, 'v', cys, cyn
492          ENDIF
493
494    END SELECT
495
496    IF ( TRIM( topography ) /= 'flat' )  THEN
497       IF ( TRIM( topography_grid_convention ) == ' ' )  THEN
498          IF ( TRIM( topography ) == 'single_building' .OR.  &
499               TRIM( topography ) == 'single_street_canyon' )  THEN
500             WRITE ( io, 278 )
501          ELSEIF ( TRIM( topography ) == 'read_from_file' )  THEN
502             WRITE ( io, 279 )
503          ENDIF
504       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_edge' )  THEN
505          WRITE ( io, 278 )
506       ELSEIF ( TRIM( topography_grid_convention ) == 'cell_center' )  THEN
507          WRITE ( io, 279 )
508       ENDIF
509    ENDIF
510
511    IF ( plant_canopy ) THEN
512
513       WRITE ( io, 280 ) canopy_mode, pch_index, drag_coefficient
514       IF ( passive_scalar ) THEN
515          WRITE ( io, 281 ) scalar_exchange_coefficient,   &
516                            leaf_surface_concentration
517       ENDIF
518
519!
520!--    Heat flux at the top of vegetation
521       WRITE ( io, 282 ) cthf
522
523!
524!--    Leaf area density profile
525!--    Building output strings, starting with surface value
526       WRITE ( learde, '(F6.2)' )  lad_surface
527       gradients = '------'
528       slices = '     0'
529       coordinates = '   0.0'
530       i = 1
531       DO  WHILE ( lad_vertical_gradient_level_ind(i) /= -9999 )
532
533          WRITE (coor_chr,'(F7.2)')  lad(lad_vertical_gradient_level_ind(i))
534          learde = TRIM( learde ) // ' ' // TRIM( coor_chr )
535
536          WRITE (coor_chr,'(F7.2)')  lad_vertical_gradient(i)
537          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
538
539          WRITE (coor_chr,'(I7)')  lad_vertical_gradient_level_ind(i)
540          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
541
542          WRITE (coor_chr,'(F7.1)')  lad_vertical_gradient_level(i)
543          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
544
545          i = i + 1
546       ENDDO
547
548       WRITE ( io, 283 )  TRIM( coordinates ), TRIM( learde ), &
549                          TRIM( gradients ), TRIM( slices )
550
551    ENDIF
552
553!
554!-- Boundary conditions
555    IF ( ibc_p_b == 0 )  THEN
556       runten = 'p(0)     = 0      |'
557    ELSEIF ( ibc_p_b == 1 )  THEN
558       runten = 'p(0)     = p(1)   |'
559    ELSE
560       runten = 'p(0)     = p(1) +R|'
561    ENDIF
562    IF ( ibc_p_t == 0 )  THEN
563       roben  = 'p(nzt+1) = 0      |'
564    ELSE
565       roben  = 'p(nzt+1) = p(nzt) |'
566    ENDIF
567
568    IF ( ibc_uv_b == 0 )  THEN
569       runten = TRIM( runten ) // ' uv(0)     = -uv(1)                |'
570    ELSE
571       runten = TRIM( runten ) // ' uv(0)     = uv(1)                 |'
572    ENDIF
573    IF ( TRIM( bc_uv_t ) == 'dirichlet_0' )  THEN
574       roben  = TRIM( roben  ) // ' uv(nzt+1) = 0                     |'
575    ELSEIF ( ibc_uv_t == 0 )  THEN
576       roben  = TRIM( roben  ) // ' uv(nzt+1) = ug(nzt+1), vg(nzt+1)  |'
577    ELSE
578       roben  = TRIM( roben  ) // ' uv(nzt+1) = uv(nzt)               |'
579    ENDIF
580
581    IF ( ibc_pt_b == 0 )  THEN
582       runten = TRIM( runten ) // ' pt(0)   = pt_surface'
583    ELSEIF ( ibc_pt_b == 1 )  THEN
584       runten = TRIM( runten ) // ' pt(0)   = pt(1)'
585    ELSEIF ( ibc_pt_b == 2 )  THEN
586       runten = TRIM( runten ) // ' pt(0) = from coupled model'
587    ENDIF
588    IF ( ibc_pt_t == 0 )  THEN
589       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt_top'
590    ELSEIF( ibc_pt_t == 1 )  THEN
591       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt(nzt)'
592    ELSEIF( ibc_pt_t == 2 )  THEN
593       roben  = TRIM( roben  ) // ' pt(nzt+1) = pt(nzt) + dpt/dz_ini'
594
595    ENDIF
596
597    WRITE ( io, 300 )  runten, roben
598
599    IF ( .NOT. constant_diffusion )  THEN
600       IF ( ibc_e_b == 1 )  THEN
601          runten = 'e(0)     = e(1)'
602       ELSE
603          runten = 'e(0)     = e(1) = (u*/0.1)**2'
604       ENDIF
605       roben = 'e(nzt+1) = e(nzt) = e(nzt-1)'
606
607       WRITE ( io, 301 )  'e', runten, roben       
608
609    ENDIF
610
611    IF ( ocean )  THEN
612       runten = 'sa(0)    = sa(1)'
613       IF ( ibc_sa_t == 0 )  THEN
614          roben =  'sa(nzt+1) = sa_surface'
615       ELSE
616          roben =  'sa(nzt+1) = sa(nzt)'
617       ENDIF
618       WRITE ( io, 301 ) 'sa', runten, roben
619    ENDIF
620
621    IF ( humidity )  THEN
622       IF ( ibc_q_b == 0 )  THEN
623          runten = 'q(0)     = q_surface'
624       ELSE
625          runten = 'q(0)     = q(1)'
626       ENDIF
627       IF ( ibc_q_t == 0 )  THEN
628          roben =  'q(nzt)   = q_top'
629       ELSE
630          roben =  'q(nzt)   = q(nzt-1) + dq/dz'
631       ENDIF
632       WRITE ( io, 301 ) 'q', runten, roben
633    ENDIF
634
635    IF ( passive_scalar )  THEN
636       IF ( ibc_q_b == 0 )  THEN
637          runten = 's(0)     = s_surface'
638       ELSE
639          runten = 's(0)     = s(1)'
640       ENDIF
641       IF ( ibc_q_t == 0 )  THEN
642          roben =  's(nzt)   = s_top'
643       ELSE
644          roben =  's(nzt)   = s(nzt-1) + ds/dz'
645       ENDIF
646       WRITE ( io, 301 ) 's', runten, roben
647    ENDIF
648
649    IF ( use_surface_fluxes )  THEN
650       WRITE ( io, 303 )
651       IF ( constant_heatflux )  THEN
652          WRITE ( io, 306 )  surface_heatflux
653          IF ( random_heatflux )  WRITE ( io, 307 )
654       ENDIF
655       IF ( humidity  .AND.  constant_waterflux )  THEN
656          WRITE ( io, 311 ) surface_waterflux
657       ENDIF
658       IF ( passive_scalar  .AND.  constant_waterflux )  THEN
659          WRITE ( io, 313 ) surface_waterflux
660       ENDIF
661    ENDIF
662
663    IF ( use_top_fluxes )  THEN
664       WRITE ( io, 304 )
665       IF ( coupling_mode == 'uncoupled' )  THEN
666          WRITE ( io, 320 )  top_momentumflux_u, top_momentumflux_v
667          IF ( constant_top_heatflux )  THEN
668             WRITE ( io, 306 )  top_heatflux
669          ENDIF
670       ELSEIF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
671          WRITE ( io, 316 )
672       ENDIF
673       IF ( ocean  .AND.  constant_top_salinityflux )  THEN
674          WRITE ( io, 309 )  top_salinityflux
675       ENDIF
676       IF ( humidity  .OR.  passive_scalar )  THEN
677          WRITE ( io, 315 )
678       ENDIF
679    ENDIF
680
681    IF ( prandtl_layer )  THEN
682       WRITE ( io, 305 )  (zu(1)-zu(0)), roughness_length, kappa, &
683                          rif_min, rif_max
684       IF ( .NOT. constant_heatflux )  WRITE ( io, 308 )
685       IF ( humidity  .AND.  .NOT. constant_waterflux )  THEN
686          WRITE ( io, 312 )
687       ENDIF
688       IF ( passive_scalar  .AND.  .NOT. constant_waterflux )  THEN
689          WRITE ( io, 314 )
690       ENDIF
691    ELSE
692       IF ( INDEX(initializing_actions, 'set_1d-model_profiles') /= 0 )  THEN
693          WRITE ( io, 310 )  rif_min, rif_max
694       ENDIF
695    ENDIF
696
697    WRITE ( io, 317 )  bc_lr, bc_ns
698    IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
699       WRITE ( io, 318 )  outflow_damping_width, km_damp_max
700       IF ( turbulent_inflow )  THEN
701          WRITE ( io, 319 )  recycling_width, recycling_plane, &
702                             inflow_damping_height, inflow_damping_width
703       ENDIF
704    ENDIF
705
706!
707!-- Listing of 1D-profiles
708    WRITE ( io, 325 )  dt_dopr_listing
709    IF ( averaging_interval_pr /= 0.0 )  THEN
710       WRITE ( io, 326 )  averaging_interval_pr, dt_averaging_input_pr
711    ENDIF
712
713!
714!-- DATA output
715    WRITE ( io, 330 )
716    IF ( averaging_interval_pr /= 0.0 )  THEN
717       WRITE ( io, 326 )  averaging_interval_pr, dt_averaging_input_pr
718    ENDIF
719
720!
721!-- 1D-profiles
722    dopr_chr = 'Profile:'
723    IF ( dopr_n /= 0 )  THEN
724       WRITE ( io, 331 )
725
726       output_format = ''
727       IF ( netcdf_output )  THEN
728          IF ( netcdf_data_format == 1 )  THEN
729             output_format = 'NetCDF classic'
730          ELSE
731             output_format = 'NetCDF 64bit offset'
732          ENDIF
733       ENDIF
734       IF ( profil_output )  THEN
735          IF ( netcdf_output )  THEN
736             output_format = TRIM( output_format ) // ' and profil'
737          ELSE
738             output_format = 'profil'
739          ENDIF
740       ENDIF
741       WRITE ( io, 344 )  output_format
742
743       DO  i = 1, dopr_n
744          dopr_chr = TRIM( dopr_chr ) // ' ' // TRIM( data_output_pr(i) ) // ','
745          IF ( LEN_TRIM( dopr_chr ) >= 60 )  THEN
746             WRITE ( io, 332 )  dopr_chr
747             dopr_chr = '       :'
748          ENDIF
749       ENDDO
750
751       IF ( dopr_chr /= '' )  THEN
752          WRITE ( io, 332 )  dopr_chr
753       ENDIF
754       WRITE ( io, 333 )  dt_dopr, averaging_interval_pr, dt_averaging_input_pr
755       IF ( skip_time_dopr /= 0.0 )  WRITE ( io, 339 )  skip_time_dopr
756    ENDIF
757
758!
759!-- 2D-arrays
760    DO  av = 0, 1
761
762       i = 1
763       do2d_xy = ''
764       do2d_xz = ''
765       do2d_yz = ''
766       DO  WHILE ( do2d(av,i) /= ' ' )
767
768          l = MAX( 2, LEN_TRIM( do2d(av,i) ) )
769          do2d_mode = do2d(av,i)(l-1:l)
770
771          SELECT CASE ( do2d_mode )
772             CASE ( 'xy' )
773                ll = LEN_TRIM( do2d_xy )
774                do2d_xy = do2d_xy(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
775             CASE ( 'xz' )
776                ll = LEN_TRIM( do2d_xz )
777                do2d_xz = do2d_xz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
778             CASE ( 'yz' )
779                ll = LEN_TRIM( do2d_yz )
780                do2d_yz = do2d_yz(1:ll) // ' ' // do2d(av,i)(1:l-3) // ','
781          END SELECT
782
783          i = i + 1
784
785       ENDDO
786
787       IF ( ( ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  .OR.    &
788              ( do2d_xz /= ''  .AND.  section(1,2) /= -9999 )  .OR.    &
789              ( do2d_yz /= ''  .AND.  section(1,3) /= -9999 ) )  .AND. &
790            ( netcdf_output  .OR.  iso2d_output ) )  THEN
791
792          IF (  av == 0 )  THEN
793             WRITE ( io, 334 )  ''
794          ELSE
795             WRITE ( io, 334 )  '(time-averaged)'
796          ENDIF
797
798          IF ( do2d_at_begin )  THEN
799             begin_chr = 'and at the start'
800          ELSE
801             begin_chr = ''
802          ENDIF
803
804          output_format = ''
805          IF ( netcdf_output )  THEN
806             IF ( netcdf_data_format == 1 )  THEN
807                output_format = 'NetCDF classic'
808             ELSEIF ( netcdf_data_format == 2 )  THEN
809                output_format = 'NetCDF 64bit offset'
810             ELSEIF ( netcdf_data_format == 3 )  THEN
811                output_format = 'NetCDF4/HDF5'
812             ELSEIF ( netcdf_data_format == 4 )  THEN
813                output_format = 'NetCDF4/HDF5 clasic'
814             ENDIF
815          ENDIF
816          IF ( iso2d_output )  THEN
817             IF ( netcdf_output )  THEN
818                output_format = TRIM( output_format ) // ' and iso2d'
819             ELSE
820                output_format = 'iso2d'
821             ENDIF
822          ENDIF
823          WRITE ( io, 344 )  output_format
824
825          IF ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  THEN
826             i = 1
827             slices = '/'
828             coordinates = '/'
829!
830!--          Building strings with index and coordinate informations of the
831!--          slices
832             DO  WHILE ( section(i,1) /= -9999 )
833
834                WRITE (section_chr,'(I5)')  section(i,1)
835                section_chr = ADJUSTL( section_chr )
836                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
837
838                IF ( section(i,1) == -1 )  THEN
839                   WRITE (coor_chr,'(F10.1)')  -1.0
840                ELSE
841                   WRITE (coor_chr,'(F10.1)')  zu(section(i,1))
842                ENDIF
843                coor_chr = ADJUSTL( coor_chr )
844                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
845
846                i = i + 1
847             ENDDO
848             IF ( av == 0 )  THEN
849                WRITE ( io, 335 )  'XY', do2d_xy, dt_do2d_xy, &
850                                   TRIM( begin_chr ), 'k', TRIM( slices ), &
851                                   TRIM( coordinates )
852                IF ( skip_time_do2d_xy /= 0.0 )  THEN
853                   WRITE ( io, 339 )  skip_time_do2d_xy
854                ENDIF
855             ELSE
856                WRITE ( io, 342 )  'XY', do2d_xy, dt_data_output_av, &
857                                   TRIM( begin_chr ), averaging_interval, &
858                                   dt_averaging_input, 'k', TRIM( slices ), &
859                                   TRIM( coordinates )
860                IF ( skip_time_data_output_av /= 0.0 )  THEN
861                   WRITE ( io, 339 )  skip_time_data_output_av
862                ENDIF
863             ENDIF
864
865          ENDIF
866
867          IF ( do2d_xz /= ''  .AND.  section(1,2) /= -9999 )  THEN
868             i = 1
869             slices = '/'
870             coordinates = '/'
871!
872!--          Building strings with index and coordinate informations of the
873!--          slices
874             DO  WHILE ( section(i,2) /= -9999 )
875
876                WRITE (section_chr,'(I5)')  section(i,2)
877                section_chr = ADJUSTL( section_chr )
878                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
879
880                WRITE (coor_chr,'(F10.1)')  section(i,2) * dy
881                coor_chr = ADJUSTL( coor_chr )
882                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
883
884                i = i + 1
885             ENDDO
886             IF ( av == 0 )  THEN
887                WRITE ( io, 335 )  'XZ', do2d_xz, dt_do2d_xz, &
888                                   TRIM( begin_chr ), 'j', TRIM( slices ), &
889                                   TRIM( coordinates )
890                IF ( skip_time_do2d_xz /= 0.0 )  THEN
891                   WRITE ( io, 339 )  skip_time_do2d_xz
892                ENDIF
893             ELSE
894                WRITE ( io, 342 )  'XZ', do2d_xz, dt_data_output_av, &
895                                   TRIM( begin_chr ), averaging_interval, &
896                                   dt_averaging_input, 'j', TRIM( slices ), &
897                                   TRIM( coordinates )
898                IF ( skip_time_data_output_av /= 0.0 )  THEN
899                   WRITE ( io, 339 )  skip_time_data_output_av
900                ENDIF
901             ENDIF
902          ENDIF
903
904          IF ( do2d_yz /= ''  .AND.  section(1,3) /= -9999 )  THEN
905             i = 1
906             slices = '/'
907             coordinates = '/'
908!
909!--          Building strings with index and coordinate informations of the
910!--          slices
911             DO  WHILE ( section(i,3) /= -9999 )
912
913                WRITE (section_chr,'(I5)')  section(i,3)
914                section_chr = ADJUSTL( section_chr )
915                slices = TRIM( slices ) // TRIM( section_chr ) // '/'
916
917                WRITE (coor_chr,'(F10.1)')  section(i,3) * dx
918                coor_chr = ADJUSTL( coor_chr )
919                coordinates = TRIM( coordinates ) // TRIM( coor_chr ) // '/'
920
921                i = i + 1
922             ENDDO
923             IF ( av == 0 )  THEN
924                WRITE ( io, 335 )  'YZ', do2d_yz, dt_do2d_yz, &
925                                   TRIM( begin_chr ), 'i', TRIM( slices ), &
926                                   TRIM( coordinates )
927                IF ( skip_time_do2d_yz /= 0.0 )  THEN
928                   WRITE ( io, 339 )  skip_time_do2d_yz
929                ENDIF
930             ELSE
931                WRITE ( io, 342 )  'YZ', do2d_yz, dt_data_output_av, &
932                                   TRIM( begin_chr ), averaging_interval, &
933                                   dt_averaging_input, 'i', TRIM( slices ), &
934                                   TRIM( coordinates )
935                IF ( skip_time_data_output_av /= 0.0 )  THEN
936                   WRITE ( io, 339 )  skip_time_data_output_av
937                ENDIF
938             ENDIF
939          ENDIF
940
941       ENDIF
942
943    ENDDO
944
945!
946!-- 3d-arrays
947    DO  av = 0, 1
948
949       i = 1
950       do3d_chr = ''
951       DO  WHILE ( do3d(av,i) /= ' ' )
952
953          do3d_chr = TRIM( do3d_chr ) // ' ' // TRIM( do3d(av,i) ) // ','
954          i = i + 1
955
956       ENDDO
957
958       IF ( do3d_chr /= '' )  THEN
959          IF ( av == 0 )  THEN
960             WRITE ( io, 336 )  ''
961          ELSE
962             WRITE ( io, 336 )  '(time-averaged)'
963          ENDIF
964
965          output_format = ''
966          IF ( netcdf_output )  THEN
967             IF ( netcdf_data_format == 1 )  THEN
968                output_format = 'NetCDF classic'
969             ELSEIF ( netcdf_data_format == 2 )  THEN
970                output_format = 'NetCDF 64bit offset'
971             ELSEIF ( netcdf_data_format == 3 )  THEN
972                output_format = 'NetCDF4/HDF5'
973             ELSEIF ( netcdf_data_format == 4 )  THEN
974                output_format = 'NetCDF4/HDF5 clasic'
975             ENDIF
976          ENDIF
977          IF ( avs_output )  THEN
978             IF ( netcdf_output )  THEN
979                output_format = TRIM( output_format ) // ' and avs'
980             ELSE
981                output_format = 'avs'
982             ENDIF
983          ENDIF
984          WRITE ( io, 344 )  output_format
985
986          IF ( do3d_at_begin )  THEN
987             begin_chr = 'and at the start'
988          ELSE
989             begin_chr = ''
990          ENDIF
991          IF ( av == 0 )  THEN
992             WRITE ( io, 337 )  do3d_chr, dt_do3d, TRIM( begin_chr ), &
993                                zu(nz_do3d), nz_do3d
994          ELSE
995             WRITE ( io, 343 )  do3d_chr, dt_data_output_av,           &
996                                TRIM( begin_chr ), averaging_interval, &
997                                dt_averaging_input, zu(nz_do3d), nz_do3d
998          ENDIF
999
1000          IF ( do3d_compress )  THEN
1001             do3d_chr = ''
1002             i = 1
1003             DO WHILE ( do3d(av,i) /= ' ' )
1004
1005                SELECT CASE ( do3d(av,i) )
1006                   CASE ( 'u' )
1007                      j = 1
1008                   CASE ( 'v' )
1009                      j = 2
1010                   CASE ( 'w' )
1011                      j = 3
1012                   CASE ( 'p' )
1013                      j = 4
1014                   CASE ( 'pt' )
1015                      j = 5
1016                END SELECT
1017                WRITE ( prec, '(I1)' )  plot_3d_precision(j)%precision
1018                do3d_chr = TRIM( do3d_chr ) // ' ' // TRIM( do3d(av,i) ) // &
1019                           ':' // prec // ','
1020                i = i + 1
1021
1022             ENDDO
1023             WRITE ( io, 338 )  do3d_chr
1024
1025          ENDIF
1026
1027          IF ( av == 0 )  THEN
1028             IF ( skip_time_do3d /= 0.0 )  THEN
1029                WRITE ( io, 339 )  skip_time_do3d
1030             ENDIF
1031          ELSE
1032             IF ( skip_time_data_output_av /= 0.0 )  THEN
1033                WRITE ( io, 339 )  skip_time_data_output_av
1034             ENDIF
1035          ENDIF
1036
1037       ENDIF
1038
1039    ENDDO
1040
1041!
1042!-- masked arrays
1043    IF ( masks > 0 )  WRITE ( io, 345 )  &
1044         mask_scale_x, mask_scale_y, mask_scale_z
1045    DO  mid = 1, masks
1046       DO  av = 0, 1
1047
1048          i = 1
1049          domask_chr = ''
1050          DO  WHILE ( domask(mid,av,i) /= ' ' )
1051             domask_chr = TRIM( domask_chr ) // ' ' //  &
1052                          TRIM( domask(mid,av,i) ) // ','
1053             i = i + 1
1054          ENDDO
1055
1056          IF ( domask_chr /= '' )  THEN
1057             IF ( av == 0 )  THEN
1058                WRITE ( io, 346 )  '', mid
1059             ELSE
1060                WRITE ( io, 346 )  ' (time-averaged)', mid
1061             ENDIF
1062
1063             output_format = ''
1064             IF ( netcdf_output )  THEN
1065                IF ( netcdf_data_format == 1 )  THEN
1066                   output_format = 'NetCDF classic'
1067                ELSEIF ( netcdf_data_format == 2 )  THEN
1068                   output_format = 'NetCDF 64bit offset'
1069                ELSEIF ( netcdf_data_format == 3 )  THEN
1070                   output_format = 'NetCDF4/HDF5'
1071                ELSEIF ( netcdf_data_format == 4 )  THEN
1072                   output_format = 'NetCDF4/HDF5 clasic'
1073                ENDIF
1074             ENDIF
1075             WRITE ( io, 344 )  output_format
1076
1077             IF ( av == 0 )  THEN
1078                WRITE ( io, 347 )  domask_chr, dt_domask(mid)
1079             ELSE
1080                WRITE ( io, 348 )  domask_chr, dt_data_output_av, &
1081                                   averaging_interval, dt_averaging_input
1082             ENDIF
1083
1084             IF ( av == 0 )  THEN
1085                IF ( skip_time_domask(mid) /= 0.0 )  THEN
1086                   WRITE ( io, 339 )  skip_time_domask(mid)
1087                ENDIF
1088             ELSE
1089                IF ( skip_time_data_output_av /= 0.0 )  THEN
1090                   WRITE ( io, 339 )  skip_time_data_output_av
1091                ENDIF
1092             ENDIF
1093!
1094!--          output locations
1095             DO  dim = 1, 3
1096                IF ( mask(mid,dim,1) >= 0.0 )  THEN
1097                   count = 0
1098                   DO  WHILE ( mask(mid,dim,count+1) >= 0.0 )
1099                      count = count + 1
1100                   ENDDO
1101                   WRITE ( io, 349 )  dir(dim), dir(dim), mid, dir(dim), &
1102                                      mask(mid,dim,:count)
1103                ELSEIF ( mask_loop(mid,dim,1) < 0.0 .AND.  &
1104                         mask_loop(mid,dim,2) < 0.0 .AND.  &
1105                         mask_loop(mid,dim,3) == 0.0 )  THEN
1106                   WRITE ( io, 350 )  dir(dim), dir(dim)
1107                ELSEIF ( mask_loop(mid,dim,3) == 0.0 )  THEN
1108                   WRITE ( io, 351 )  dir(dim), dir(dim), mid, dir(dim), &
1109                                      mask_loop(mid,dim,1:2)
1110                ELSE
1111                   WRITE ( io, 351 )  dir(dim), dir(dim), mid, dir(dim), &
1112                                      mask_loop(mid,dim,1:3)
1113                ENDIF
1114             ENDDO
1115          ENDIF
1116
1117       ENDDO
1118    ENDDO
1119
1120!
1121!-- Timeseries
1122    IF ( dt_dots /= 9999999.9 )  THEN
1123       WRITE ( io, 340 )
1124
1125       output_format = ''
1126       IF ( netcdf_output )  THEN
1127          IF ( netcdf_data_format == 1 )  THEN
1128             output_format = 'NetCDF classic'
1129          ELSE
1130             output_format = 'NetCDF 64bit offset'
1131          ENDIF
1132       ENDIF
1133       IF ( profil_output )  THEN
1134          IF ( netcdf_output )  THEN
1135             output_format = TRIM( output_format ) // ' and profil'
1136          ELSE
1137             output_format = 'profil'
1138          ENDIF
1139       ENDIF
1140       WRITE ( io, 344 )  output_format
1141       WRITE ( io, 341 )  dt_dots
1142    ENDIF
1143
1144#if defined( __dvrp_graphics )
1145!
1146!-- Dvrp-output
1147    IF ( dt_dvrp /= 9999999.9 )  THEN
1148       WRITE ( io, 360 )  dt_dvrp, TRIM( dvrp_output ), TRIM( dvrp_host ), &
1149                          TRIM( dvrp_username ), TRIM( dvrp_directory )
1150       i = 1
1151       l = 0
1152       m = 0
1153       DO WHILE ( mode_dvrp(i) /= ' ' )
1154          IF ( mode_dvrp(i)(1:10) == 'isosurface' )  THEN
1155             READ ( mode_dvrp(i), '(10X,I2)' )  j
1156             l = l + 1
1157             IF ( do3d(0,j) /= ' ' )  THEN
1158                WRITE ( io, 361 )  TRIM( do3d(0,j) ), threshold(l), &
1159                                   isosurface_color(:,l)
1160             ENDIF
1161          ELSEIF ( mode_dvrp(i)(1:6) == 'slicer' )  THEN
1162             READ ( mode_dvrp(i), '(6X,I2)' )  j
1163             m = m + 1
1164             IF ( do2d(0,j) /= ' ' )  THEN
1165                WRITE ( io, 362 )  TRIM( do2d(0,j) ), &
1166                                   slicer_range_limits_dvrp(:,m)
1167             ENDIF
1168          ELSEIF ( mode_dvrp(i)(1:9) == 'particles' )  THEN
1169             WRITE ( io, 363 )  dvrp_psize
1170             IF ( particle_dvrpsize /= 'none' )  THEN
1171                WRITE ( io, 364 )  'size', TRIM( particle_dvrpsize ), &
1172                                   dvrpsize_interval
1173             ENDIF
1174             IF ( particle_color /= 'none' )  THEN
1175                WRITE ( io, 364 )  'color', TRIM( particle_color ), &
1176                                   color_interval
1177             ENDIF
1178          ENDIF
1179          i = i + 1
1180       ENDDO
1181
1182       WRITE ( io, 365 )  groundplate_color, superelevation_x, &
1183                          superelevation_y, superelevation, clip_dvrp_l, &
1184                          clip_dvrp_r, clip_dvrp_s, clip_dvrp_n
1185
1186       IF ( TRIM( topography ) /= 'flat' )  THEN
1187          WRITE ( io, 366 )  topography_color
1188          IF ( cluster_size > 1 )  THEN
1189             WRITE ( io, 367 )  cluster_size
1190          ENDIF
1191       ENDIF
1192
1193    ENDIF
1194#endif
1195
1196#if defined( __spectra )
1197!
1198!-- Spectra output
1199    IF ( dt_dosp /= 9999999.9 ) THEN
1200       WRITE ( io, 370 )
1201
1202       output_format = ''
1203       IF ( netcdf_output )  THEN
1204          IF ( netcdf_data_format == 1 )  THEN
1205             output_format = 'NetCDF classic'
1206          ELSE
1207             output_format = 'NetCDF 64bit offset'
1208          ENDIF
1209       ENDIF
1210       IF ( profil_output )  THEN
1211          IF ( netcdf_output )  THEN
1212             output_format = TRIM( output_format ) // ' and profil'
1213          ELSE
1214             output_format = 'profil'
1215          ENDIF
1216       ENDIF
1217       WRITE ( io, 344 )  output_format
1218       WRITE ( io, 371 )  dt_dosp
1219       IF ( skip_time_dosp /= 0.0 )  WRITE ( io, 339 )  skip_time_dosp
1220       WRITE ( io, 372 )  ( data_output_sp(i), i = 1,10 ),     &
1221                          ( spectra_direction(i), i = 1,10 ),  &
1222                          ( comp_spectra_level(i), i = 1,100 ), &
1223                          ( plot_spectra_level(i), i = 1,100 ), &
1224                          averaging_interval_sp, dt_averaging_input_pr
1225    ENDIF
1226#endif
1227
1228    WRITE ( io, 99 )
1229
1230!
1231!-- Physical quantities
1232    WRITE ( io, 400 )
1233
1234!
1235!-- Geostrophic parameters
1236    WRITE ( io, 410 )  omega, phi, f, fs
1237
1238!
1239!-- Other quantities
1240    WRITE ( io, 411 )  g
1241    IF ( use_reference )  THEN
1242       IF ( ocean )  THEN
1243          WRITE ( io, 412 )  prho_reference
1244       ELSE
1245          WRITE ( io, 413 )  pt_reference
1246       ENDIF
1247    ENDIF
1248
1249!
1250!-- Cloud physics parameters
1251    IF ( cloud_physics ) THEN
1252       WRITE ( io, 415 )
1253       WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
1254    ENDIF
1255
1256!-- Profile of the geostrophic wind (component ug)
1257!-- Building output strings
1258    WRITE ( ugcomponent, '(F6.2)' )  ug_surface
1259    gradients = '------'
1260    slices = '     0'
1261    coordinates = '   0.0'
1262    i = 1
1263    DO  WHILE ( ug_vertical_gradient_level_ind(i) /= -9999 )
1264     
1265       WRITE (coor_chr,'(F6.2,1X)')  ug(ug_vertical_gradient_level_ind(i))
1266       ugcomponent = TRIM( ugcomponent ) // '  ' // TRIM( coor_chr )
1267
1268       WRITE (coor_chr,'(F6.2,1X)')  ug_vertical_gradient(i)
1269       gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
1270
1271       WRITE (coor_chr,'(I6,1X)')  ug_vertical_gradient_level_ind(i)
1272       slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
1273
1274       WRITE (coor_chr,'(F6.1,1X)')  ug_vertical_gradient_level(i)
1275       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
1276
1277       IF ( i == 10 )  THEN
1278          EXIT
1279       ELSE
1280          i = i + 1
1281       ENDIF
1282
1283    ENDDO
1284
1285    WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
1286                       TRIM( gradients ), TRIM( slices )
1287
1288!-- Profile of the geostrophic wind (component vg)
1289!-- Building output strings
1290    WRITE ( vgcomponent, '(F6.2)' )  vg_surface
1291    gradients = '------'
1292    slices = '     0'
1293    coordinates = '   0.0'
1294    i = 1
1295    DO  WHILE ( vg_vertical_gradient_level_ind(i) /= -9999 )
1296
1297       WRITE (coor_chr,'(F6.2,1X)')  vg(vg_vertical_gradient_level_ind(i))
1298       vgcomponent = TRIM( vgcomponent ) // '  ' // TRIM( coor_chr )
1299
1300       WRITE (coor_chr,'(F6.2,1X)')  vg_vertical_gradient(i)
1301       gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
1302
1303       WRITE (coor_chr,'(I6,1X)')  vg_vertical_gradient_level_ind(i)
1304       slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
1305
1306       WRITE (coor_chr,'(F6.1,1X)')  vg_vertical_gradient_level(i)
1307       coordinates = TRIM( coordinates ) // '  ' // TRIM( coor_chr )
1308
1309       IF ( i == 10 )  THEN
1310          EXIT
1311       ELSE
1312          i = i + 1
1313       ENDIF
1314 
1315    ENDDO
1316
1317    WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
1318                       TRIM( gradients ), TRIM( slices )
1319
1320!
1321!-- Initial wind profiles
1322    IF ( u_profile(1) /= 9999999.9 )  WRITE ( io, 427 )
1323
1324!
1325!-- Initial temperature profile
1326!-- Building output strings, starting with surface temperature
1327    WRITE ( temperatures, '(F6.2)' )  pt_surface
1328    gradients = '------'
1329    slices = '     0'
1330    coordinates = '   0.0'
1331    i = 1
1332    DO  WHILE ( pt_vertical_gradient_level_ind(i) /= -9999 )
1333
1334       WRITE (coor_chr,'(F7.2)')  pt_init(pt_vertical_gradient_level_ind(i))
1335       temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1336
1337       WRITE (coor_chr,'(F7.2)')  pt_vertical_gradient(i)
1338       gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1339
1340       WRITE (coor_chr,'(I7)')  pt_vertical_gradient_level_ind(i)
1341       slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1342
1343       WRITE (coor_chr,'(F7.1)')  pt_vertical_gradient_level(i)
1344       coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1345
1346       IF ( i == 10 )  THEN
1347          EXIT
1348       ELSE
1349          i = i + 1
1350       ENDIF
1351
1352    ENDDO
1353
1354    WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
1355                       TRIM( gradients ), TRIM( slices )
1356
1357!
1358!-- Initial humidity profile
1359!-- Building output strings, starting with surface humidity
1360    IF ( humidity  .OR.  passive_scalar )  THEN
1361       WRITE ( temperatures, '(E8.1)' )  q_surface
1362       gradients = '--------'
1363       slices = '       0'
1364       coordinates = '     0.0'
1365       i = 1
1366       DO  WHILE ( q_vertical_gradient_level_ind(i) /= -9999 )
1367         
1368          WRITE (coor_chr,'(E8.1,4X)')  q_init(q_vertical_gradient_level_ind(i))
1369          temperatures = TRIM( temperatures ) // '  ' // TRIM( coor_chr )
1370
1371          WRITE (coor_chr,'(E8.1,4X)')  q_vertical_gradient(i)
1372          gradients = TRIM( gradients ) // '  ' // TRIM( coor_chr )
1373         
1374          WRITE (coor_chr,'(I8,4X)')  q_vertical_gradient_level_ind(i)
1375          slices = TRIM( slices ) // '  ' // TRIM( coor_chr )
1376         
1377          WRITE (coor_chr,'(F8.1,4X)')  q_vertical_gradient_level(i)
1378          coordinates = TRIM( coordinates ) // '  '  // TRIM( coor_chr )
1379
1380          IF ( i == 10 )  THEN
1381             EXIT
1382          ELSE
1383             i = i + 1
1384          ENDIF
1385
1386       ENDDO
1387
1388       IF ( humidity )  THEN
1389          WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
1390                             TRIM( gradients ), TRIM( slices )
1391       ELSE
1392          WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ), &
1393                             TRIM( gradients ), TRIM( slices )
1394       ENDIF
1395    ENDIF
1396
1397!
1398!-- Initial salinity profile
1399!-- Building output strings, starting with surface salinity
1400    IF ( ocean )  THEN
1401       WRITE ( temperatures, '(F6.2)' )  sa_surface
1402       gradients = '------'
1403       slices = '     0'
1404       coordinates = '   0.0'
1405       i = 1
1406       DO  WHILE ( sa_vertical_gradient_level_ind(i) /= -9999 )
1407
1408          WRITE (coor_chr,'(F7.2)')  sa_init(sa_vertical_gradient_level_ind(i))
1409          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1410
1411          WRITE (coor_chr,'(F7.2)')  sa_vertical_gradient(i)
1412          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1413
1414          WRITE (coor_chr,'(I7)')  sa_vertical_gradient_level_ind(i)
1415          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1416
1417          WRITE (coor_chr,'(F7.1)')  sa_vertical_gradient_level(i)
1418          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1419
1420          IF ( i == 10 )  THEN
1421             EXIT
1422          ELSE
1423             i = i + 1
1424          ENDIF
1425
1426       ENDDO
1427
1428       WRITE ( io, 425 )  TRIM( coordinates ), TRIM( temperatures ), &
1429                          TRIM( gradients ), TRIM( slices )
1430    ENDIF
1431
1432!
1433!-- Profile for the large scale vertial velocity
1434!-- Building output strings, starting with surface value
1435    IF ( large_scale_subsidence )  THEN
1436       temperatures = '   0.0'
1437       gradients = '------'
1438       slices = '     0'
1439       coordinates = '   0.0'
1440       i = 1
1441       DO  WHILE ( subs_vertical_gradient_level_i(i) /= -9999 )
1442
1443          WRITE (coor_chr,'(E10.2,7X)')  &
1444                                w_subs(subs_vertical_gradient_level_i(i))
1445          temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
1446
1447          WRITE (coor_chr,'(E10.2,7X)')  subs_vertical_gradient(i)
1448          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
1449
1450          WRITE (coor_chr,'(I10,7X)')  subs_vertical_gradient_level_i(i)
1451          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
1452
1453          WRITE (coor_chr,'(F10.2,7X)')  subs_vertical_gradient_level(i)
1454          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
1455
1456          IF ( i == 10 )  THEN
1457             EXIT
1458          ELSE
1459             i = i + 1
1460          ENDIF
1461
1462       ENDDO
1463
1464       WRITE ( io, 426 )  TRIM( coordinates ), TRIM( temperatures ), &
1465                          TRIM( gradients ), TRIM( slices )
1466    ENDIF
1467
1468!
1469!-- LES / turbulence parameters
1470    WRITE ( io, 450 )
1471
1472!--
1473! ... LES-constants used must still be added here
1474!--
1475    IF ( constant_diffusion )  THEN
1476       WRITE ( io, 451 )  km_constant, km_constant/prandtl_number, &
1477                          prandtl_number
1478    ENDIF
1479    IF ( .NOT. constant_diffusion)  THEN
1480       IF ( e_init > 0.0 )  WRITE ( io, 455 )  e_init
1481       IF ( e_min > 0.0 )  WRITE ( io, 454 )  e_min
1482       IF ( wall_adjustment )  WRITE ( io, 453 )  wall_adjustment_factor
1483       IF ( adjust_mixing_length  .AND.  prandtl_layer )  WRITE ( io, 452 )
1484    ENDIF
1485
1486!
1487!-- Special actions during the run
1488    WRITE ( io, 470 )
1489    IF ( create_disturbances )  THEN
1490       WRITE ( io, 471 )  dt_disturb, disturbance_amplitude,                   &
1491                          zu(disturbance_level_ind_b), disturbance_level_ind_b,&
1492                          zu(disturbance_level_ind_t), disturbance_level_ind_t
1493       IF ( .NOT. bc_lr_cyc  .OR.  .NOT. bc_ns_cyc )  THEN
1494          WRITE ( io, 472 )  inflow_disturbance_begin, inflow_disturbance_end
1495       ELSE
1496          WRITE ( io, 473 )  disturbance_energy_limit
1497       ENDIF
1498       WRITE ( io, 474 )  TRIM( random_generator )
1499    ENDIF
1500    IF ( pt_surface_initial_change /= 0.0 )  THEN
1501       WRITE ( io, 475 )  pt_surface_initial_change
1502    ENDIF
1503    IF ( humidity  .AND.  q_surface_initial_change /= 0.0 )  THEN
1504       WRITE ( io, 476 )  q_surface_initial_change       
1505    ENDIF
1506    IF ( passive_scalar  .AND.  q_surface_initial_change /= 0.0 )  THEN
1507       WRITE ( io, 477 )  q_surface_initial_change       
1508    ENDIF
1509
1510    IF ( particle_advection )  THEN
1511!
1512!--    Particle attributes
1513       WRITE ( io, 480 )  particle_advection_start, dt_prel, bc_par_lr, &
1514                          bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
1515                          end_time_prel, dt_sort_particles
1516       IF ( use_sgs_for_particles )  WRITE ( io, 488 )  dt_min_part
1517       IF ( random_start_position )  WRITE ( io, 481 )
1518       IF ( particles_per_point > 1 )  WRITE ( io, 489 )  particles_per_point
1519       WRITE ( io, 495 )  total_number_of_particles
1520       IF ( maximum_number_of_tailpoints /= 0 )  THEN
1521          WRITE ( io, 483 )  maximum_number_of_tailpoints
1522          IF ( minimum_tailpoint_distance /= 0 )  THEN
1523             WRITE ( io, 484 )  total_number_of_tails,      &
1524                                minimum_tailpoint_distance, &
1525                                maximum_tailpoint_age
1526          ENDIF
1527       ENDIF
1528       IF ( dt_write_particle_data /= 9999999.9 )  THEN
1529          WRITE ( io, 485 )  dt_write_particle_data
1530          output_format = ''
1531          IF ( netcdf_output )  THEN
1532             IF ( netcdf_data_format > 1 )  THEN
1533                output_format = 'netcdf (64 bit offset) and binary'
1534             ELSE
1535                output_format = 'netcdf and binary'
1536             ENDIF
1537          ELSE
1538             output_format = 'binary'
1539          ENDIF
1540          WRITE ( io, 344 )  output_format
1541       ENDIF
1542       IF ( dt_dopts /= 9999999.9 )  WRITE ( io, 494 )  dt_dopts
1543       IF ( write_particle_statistics )  WRITE ( io, 486 )
1544
1545       WRITE ( io, 487 )  number_of_particle_groups
1546
1547       DO  i = 1, number_of_particle_groups
1548          IF ( i == 1  .AND.  density_ratio(i) == 9999999.9 )  THEN
1549             WRITE ( io, 490 )  i, 0.0
1550             WRITE ( io, 492 )
1551          ELSE
1552             WRITE ( io, 490 )  i, radius(i)
1553             IF ( density_ratio(i) /= 0.0 )  THEN
1554                WRITE ( io, 491 )  density_ratio(i)
1555             ELSE
1556                WRITE ( io, 492 )
1557             ENDIF
1558          ENDIF
1559          WRITE ( io, 493 )  psl(i), psr(i), pss(i), psn(i), psb(i), pst(i), &
1560                             pdx(i), pdy(i), pdz(i)
1561          IF ( .NOT. vertical_particle_advection(i) )  WRITE ( io, 482 )
1562       ENDDO
1563
1564    ENDIF
1565
1566
1567!
1568!-- Parameters of 1D-model
1569    IF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )  THEN
1570       WRITE ( io, 500 )  end_time_1d, dt_run_control_1d, dt_pr_1d, &
1571                          mixing_length_1d, dissipation_1d
1572       IF ( damp_level_ind_1d /= nzt+1 )  THEN
1573          WRITE ( io, 502 )  zu(damp_level_ind_1d), damp_level_ind_1d
1574       ENDIF
1575    ENDIF
1576
1577!
1578!-- User-defined informations
1579    CALL user_header( io )
1580
1581    WRITE ( io, 99 )
1582
1583!
1584!-- Write buffer contents to disc immediately
1585    CALL local_flush( io )
1586
1587!
1588!-- Here the FORMATs start
1589
1590 99 FORMAT (1X,78('-'))
1591100 FORMAT (/1X,'***************************',9X,42('-')/        &
1592            1X,'* ',A,' *',9X,A/                               &
1593            1X,'***************************',9X,42('-'))
1594101 FORMAT (37X,'coupled run using MPI-',I1,': ',A/ &
1595            37X,42('-'))
1596102 FORMAT (/' Date:              ',A8,9X,'Run:       ',A20/      &
1597            ' Time:              ',A8,9X,'Run-No.:   ',I2.2/     &
1598            ' Run on host:     ',A10)
1599#if defined( __parallel )
1600103 FORMAT (' Number of PEs:',8X,I5,9X,'Processor grid (x,y): (',I3,',',I3, &
1601              ')',1X,A)
1602104 FORMAT (' Number of PEs:',8X,I5,9X,'Tasks:',I4,'   threads per task:',I4/ &
1603              37X,'Processor grid (x,y): (',I3,',',I3,')',1X,A)
1604105 FORMAT (37X,'One additional PE is used to handle'/37X,'the dvrp output!')
1605106 FORMAT (37X,'A 1d-decomposition along x is forced'/ &
1606            37X,'because the job is running on an SMP-cluster')
1607107 FORMAT (37X,'A 1d-decomposition along ',A,' is used')
1608108 FORMAT (37X,'Max. # of parallel I/O streams is ',I5)
1609#endif
1610110 FORMAT (/' Numerical Schemes:'/ &
1611             ' -----------------'/)
1612111 FORMAT (' --> Solve perturbation pressure via FFT using ',A,' routines')
1613112 FORMAT (' --> Solve perturbation pressure via SOR-Red/Black-Schema'/ &
1614            '     Iterations (initial/other): ',I3,'/',I3,'  omega = ',F5.3)
1615113 FORMAT (' --> Momentum advection via Piascek-Williams-Scheme (Form C3)', &
1616                  ' or Upstream')
1617114 FORMAT (' --> Momentum advection via Upstream-Spline-Scheme')
1618115 FORMAT ('     Tendencies are smoothed via Long-Filter with factor ',F5.3) 
1619116 FORMAT (' --> Scalar advection via Piascek-Williams-Scheme (Form C3)', &
1620                  ' or Upstream')
1621117 FORMAT (' --> Scalar advection via Upstream-Spline-Scheme')
1622118 FORMAT (' --> Scalar advection via Bott-Chlond-Scheme')
1623119 FORMAT (' --> Galilei-Transform applied to horizontal advection', &
1624            '     Translation velocity = ',A/ &
1625            '     distance advected ',A,':  ',F8.3,' km(x)  ',F8.3,' km(y)')
1626120 FORMAT (' --> Time differencing scheme: leapfrog only (no euler in case', &
1627                  ' of timestep changes)')
1628121 FORMAT (' --> Time differencing scheme: leapfrog + euler in case of', &
1629                  ' timestep changes')
1630122 FORMAT (' --> Time differencing scheme: ',A)
1631123 FORMAT (' --> Rayleigh-Damping active, starts ',A,' z = ',F8.2,' m'/ &
1632            '     maximum damping coefficient: ',F5.3, ' 1/s')
1633124 FORMAT ('     Spline-overshoots are being suppressed')
1634125 FORMAT ('     Upstream-Scheme is used if Upstream-differences fall short', &
1635                  ' of'/                                                       &
1636            '     delta_u = ',F6.4,4X,'delta_v = ',F6.4,4X,'delta_w = ',F6.4)
1637126 FORMAT ('     Upstream-Scheme is used if Upstream-differences fall short', &
1638                  ' of'/                                                       &
1639            '     delta_e = ',F6.4,4X,'delta_pt = ',F6.4)
1640127 FORMAT ('     The following absolute overshoot differences are tolerated:'/&
1641            '     delta_u = ',F6.4,4X,'delta_v = ',F6.4,4X,'delta_w = ',F6.4)
1642128 FORMAT ('     The following absolute overshoot differences are tolerated:'/&
1643            '     delta_e = ',F6.4,4X,'delta_pt = ',F6.4)
1644129 FORMAT (' --> Additional prognostic equation for the specific humidity')
1645130 FORMAT (' --> Additional prognostic equation for the total water content')
1646131 FORMAT (' --> Parameterization of condensation processes via (0%-or100%)')
1647132 FORMAT (' --> Parameterization of long-wave radiation processes via'/ &
1648            '     effective emissivity scheme')
1649133 FORMAT (' --> Precipitation parameterization via Kessler-Scheme')
1650134 FORMAT (' --> Additional prognostic equation for a passive scalar')
1651135 FORMAT (' --> Solve perturbation pressure via multigrid method (', &
1652                  A,'-cycle)'/ &
1653            '     number of grid levels:                   ',I2/ &
1654            '     Gauss-Seidel red/black iterations:       ',I2)
1655136 FORMAT ('     gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1656                  I3,')')
1657137 FORMAT ('     level data gathered on PE0 at level:     ',I2/ &
1658            '     gridpoints of coarsest subdomain (x,y,z): (',I3,',',I3,',', &
1659                  I3,')'/ &
1660            '     gridpoints of coarsest domain (x,y,z):    (',I3,',',I3,',', &
1661                  I3,')')
1662138 FORMAT ('     Using hybrid version for 1d-domain-decomposition')
1663139 FORMAT (' --> Loop optimization method: ',A)
1664140 FORMAT ('     maximum residual allowed:                ',E10.3)
1665141 FORMAT ('     fixed number of multigrid cycles:        ',I4)
1666142 FORMAT ('     perturbation pressure is calculated at every Runge-Kutta ', &
1667                  'step')
1668143 FORMAT ('     Euler/upstream scheme is used for the SGS turbulent ', &
1669                  'kinetic energy')
1670150 FORMAT (' --> Volume flow at the right and north boundary will be ', &
1671                  'conserved'/ &
1672            '     using the ',A,' mode')
1673151 FORMAT ('     with u_bulk = ',F7.3,' m/s and v_bulk = ',F7.3,' m/s')
1674152 FORMAT (' --> External pressure gradient directly prescribed by the user:',&
1675           /'     ',2(1X,E12.5),'Pa/m in x/y direction', &
1676           /'     starting from dp_level_b =', F8.3, 'm', A /)
1677153 FORMAT (' --> Large-scale vertical motion is used in the ', &
1678                  'prognostic equation for')
1679154 FORMAT ('     the potential temperature')
1680200 FORMAT (//' Run time and time step information:'/ &
1681             ' ----------------------------------'/)
1682201 FORMAT ( ' Timestep:          variable     maximum value: ',F6.3,' s', &
1683             '    CFL-factor: ',F4.2)
1684202 FORMAT ( ' Timestep:       dt = ',F6.3,' s'/)
1685203 FORMAT ( ' Start time:       ',F9.3,' s'/ &
1686             ' End time:         ',F9.3,' s')
1687204 FORMAT ( A,F9.3,' s')
1688205 FORMAT ( A,F9.3,' s',5X,'restart every',17X,F9.3,' s')
1689206 FORMAT (/' Time reached:     ',F9.3,' s'/ &
1690             ' CPU-time used:    ',F9.3,' s     per timestep:               ', &
1691               '  ',F9.3,' s'/                                                 &
1692             '                                   per second of simulated tim', &
1693               'e: ',F9.3,' s')
1694207 FORMAT ( A/' Coupling start time:',F9.3,' s')
1695250 FORMAT (//' Computational grid and domain size:'/ &
1696              ' ----------------------------------'// &
1697              ' Grid length:      dx =    ',F7.3,' m    dy =    ',F7.3, &
1698              ' m    dz =    ',F7.3,' m'/ &
1699              ' Domain size:       x = ',F10.3,' m     y = ',F10.3, &
1700              ' m  z(u) = ',F10.3,' m'/)
1701252 FORMAT (' dz constant up to ',F10.3,' m (k=',I4,'), then stretched by', &
1702              ' factor: ',F5.3/ &
1703            ' maximum dz not to be exceeded is dz_max = ',F10.3,' m'/)
1704254 FORMAT (' Number of gridpoints (x,y,z):  (0:',I4,', 0:',I4,', 0:',I4,')'/ &
1705            ' Subdomain size (x,y,z):        (  ',I4,',   ',I4,',   ',I4,')'/)
1706255 FORMAT (' Subdomains have equal size')
1707256 FORMAT (' Subdomains at the upper edges of the virtual processor grid ', &
1708              'have smaller sizes'/                                          &
1709            ' Size of smallest subdomain:    (  ',I4,',   ',I4,',   ',I4,')')
1710260 FORMAT (/' The model has a slope in x-direction. Inclination angle: ',F6.2,&
1711             ' degrees')
1712270 FORMAT (//' Topography informations:'/ &
1713              ' -----------------------'// &
1714              1X,'Topography: ',A)
1715271 FORMAT (  ' Building size (x/y/z) in m: ',F5.1,' / ',F5.1,' / ',F5.1/ &
1716              ' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
1717                ' / ',I4)
1718272 FORMAT (  ' Single quasi-2D street canyon of infinite length in ',A, &
1719              ' direction' / &
1720              ' Canyon height: ', F6.2, 'm, ch = ', I4, '.'      / &
1721              ' Canyon position (',A,'-walls): cxl = ', I4,', cxr = ', I4, '.')
1722278 FORMAT (' Topography grid definition convention:'/ &
1723            ' cell edge (staggered grid points'/  &
1724            ' (u in x-direction, v in y-direction))' /)
1725279 FORMAT (' Topography grid definition convention:'/ &
1726            ' cell center (scalar grid points)' /)
1727280 FORMAT (//' Vegetation canopy (drag) model:'/ &
1728              ' ------------------------------'// &
1729              ' Canopy mode: ', A / &
1730              ' Canopy top: ',I4 / &
1731              ' Leaf drag coefficient: ',F6.2 /)
1732281 FORMAT (/ ' Scalar_exchange_coefficient: ',F6.2 / &
1733              ' Scalar concentration at leaf surfaces in kg/m**3: ',F6.2 /)
1734282 FORMAT (' Predefined constant heatflux at the top of the vegetation: ',F6.2,' K m/s')
1735283 FORMAT (/ ' Characteristic levels of the leaf area density:'// &
1736              ' Height:              ',A,'  m'/ &
1737              ' Leaf area density:   ',A,'  m**2/m**3'/ &
1738              ' Gradient:            ',A,'  m**2/m**4'/ &
1739              ' Gridpoint:           ',A)
1740               
1741300 FORMAT (//' Boundary conditions:'/ &
1742             ' -------------------'// &
1743             '                     p                    uv             ', &
1744             '                   pt'// &
1745             ' B. bound.: ',A/ &
1746             ' T. bound.: ',A)
1747301 FORMAT (/'                     ',A// &
1748             ' B. bound.: ',A/ &
1749             ' T. bound.: ',A)
1750303 FORMAT (/' Bottom surface fluxes are used in diffusion terms at k=1')
1751304 FORMAT (/' Top surface fluxes are used in diffusion terms at k=nzt')
1752305 FORMAT (//'    Prandtl-Layer between bottom surface and first ', &
1753               'computational u,v-level:'// &
1754             '       zp = ',F6.2,' m   z0 = ',F6.4,' m   kappa = ',F4.2/ &
1755             '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
1756306 FORMAT ('       Predefined constant heatflux:   ',F9.6,' K m/s')
1757307 FORMAT ('       Heatflux has a random normal distribution')
1758308 FORMAT ('       Predefined surface temperature')
1759309 FORMAT ('       Predefined constant salinityflux:   ',F9.6,' psu m/s')
1760310 FORMAT (//'    1D-Model:'// &
1761             '       Rif value range:   ',F6.2,' <= rif <=',F6.2)
1762311 FORMAT ('       Predefined constant humidity flux: ',E10.3,' m/s')
1763312 FORMAT ('       Predefined surface humidity')
1764313 FORMAT ('       Predefined constant scalar flux: ',E10.3,' kg/(m**2 s)')
1765314 FORMAT ('       Predefined scalar value at the surface')
1766315 FORMAT ('       Humidity / scalar flux at top surface is 0.0')
1767316 FORMAT ('       Sensible heatflux and momentum flux from coupled ', &
1768                    'atmosphere model')
1769317 FORMAT (//' Lateral boundaries:'/ &
1770            '       left/right:  ',A/    &
1771            '       north/south: ',A)
1772318 FORMAT (/'       outflow damping layer width: ',I3,' gridpoints with km_', &
1773                    'max =',F5.1,' m**2/s')
1774319 FORMAT ('       turbulence recycling at inflow switched on'/ &
1775            '       width of recycling domain: ',F7.1,' m   grid index: ',I4/ &
1776            '       inflow damping height: ',F6.1,' m   width: ',F6.1,' m')
1777320 FORMAT ('       Predefined constant momentumflux:  u: ',F9.6,' m**2/s**2'/ &
1778            '                                          v: ',F9.6,' m**2/s**2')
1779325 FORMAT (//' List output:'/ &
1780             ' -----------'//  &
1781            '    1D-Profiles:'/    &
1782            '       Output every             ',F8.2,' s')
1783326 FORMAT ('       Time averaged over       ',F8.2,' s'/ &
1784            '       Averaging input every    ',F8.2,' s')
1785330 FORMAT (//' Data output:'/ &
1786             ' -----------'/)
1787331 FORMAT (/'    1D-Profiles:')
1788332 FORMAT (/'       ',A)
1789333 FORMAT ('       Output every             ',F8.2,' s',/ &
1790            '       Time averaged over       ',F8.2,' s'/ &
1791            '       Averaging input every    ',F8.2,' s')
1792334 FORMAT (/'    2D-Arrays',A,':')
1793335 FORMAT (/'       ',A2,'-cross-section  Arrays: ',A/ &
1794            '       Output every             ',F8.2,' s  ',A/ &
1795            '       Cross sections at ',A1,' = ',A/ &
1796            '       scalar-coordinates:   ',A,' m'/)
1797336 FORMAT (/'    3D-Arrays',A,':')
1798337 FORMAT (/'       Arrays: ',A/ &
1799            '       Output every             ',F8.2,' s  ',A/ &
1800            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
1801338 FORMAT ('       Compressed data output'/ &
1802            '       Decimal precision: ',A/)
1803339 FORMAT ('       No output during initial ',F8.2,' s')
1804340 FORMAT (/'    Time series:')
1805341 FORMAT ('       Output every             ',F8.2,' s'/)
1806342 FORMAT (/'       ',A2,'-cross-section  Arrays: ',A/ &
1807            '       Output every             ',F8.2,' s  ',A/ &
1808            '       Time averaged over       ',F8.2,' s'/ &
1809            '       Averaging input every    ',F8.2,' s'/ &
1810            '       Cross sections at ',A1,' = ',A/ &
1811            '       scalar-coordinates:   ',A,' m'/)
1812343 FORMAT (/'       Arrays: ',A/ &
1813            '       Output every             ',F8.2,' s  ',A/ &
1814            '       Time averaged over       ',F8.2,' s'/ &
1815            '       Averaging input every    ',F8.2,' s'/ &
1816            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
1817344 FORMAT ('       Output format: ',A/)
1818345 FORMAT (/'    Scaling lengths for output locations of all subsequent mask IDs:',/ &
1819            '       mask_scale_x (in x-direction): ',F9.3, ' m',/ &
1820            '       mask_scale_y (in y-direction): ',F9.3, ' m',/ &
1821            '       mask_scale_z (in z-direction): ',F9.3, ' m' )
1822346 FORMAT (/'    Masked data output',A,' for mask ID ',I2, ':')
1823347 FORMAT ('       Variables: ',A/ &
1824            '       Output every             ',F8.2,' s')
1825348 FORMAT ('       Variables: ',A/ &
1826            '       Output every             ',F8.2,' s'/ &
1827            '       Time averaged over       ',F8.2,' s'/ &
1828            '       Averaging input every    ',F8.2,' s')
1829349 FORMAT (/'       Output locations in ',A,'-direction in multiples of ', &
1830            'mask_scale_',A,' predefined by array mask_',I2.2,'_',A,':'/ &
1831            13('       ',8(F8.2,',')/) )
1832350 FORMAT (/'       Output locations in ',A,'-direction: ', &
1833            'all gridpoints along ',A,'-direction (default).' )
1834351 FORMAT (/'       Output locations in ',A,'-direction in multiples of ', &
1835            'mask_scale_',A,' constructed from array mask_',I2.2,'_',A,'_loop:'/ &
1836            '          loop begin:',F8.2,', end:',F8.2,', stride:',F8.2 )
1837#if defined( __dvrp_graphics )
1838360 FORMAT ('    Plot-Sequence with dvrp-software:'/ &
1839            '       Output every      ',F7.1,' s'/ &
1840            '       Output mode:      ',A/ &
1841            '       Host / User:      ',A,' / ',A/ &
1842            '       Directory:        ',A// &
1843            '       The sequence contains:')
1844361 FORMAT (/'       Isosurface of "',A,'"    Threshold value: ', E12.3/ &
1845            '          Isosurface color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
1846362 FORMAT (/'       Slicer plane ',A/ &
1847            '       Slicer limits: [',F6.2,',',F6.2,']')
1848363 FORMAT (/'       Particles'/ &
1849            '          particle size:  ',F7.2,' m')
1850364 FORMAT ('          particle ',A,' controlled by "',A,'" with interval [', &
1851                       F6.2,',',F6.2,']')
1852365 FORMAT (/'       Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &
1853            '       Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &
1854                     ')'/ &
1855            '       Clipping limits: from x = ',F9.1,' m to x = ',F9.1,' m'/ &
1856            '                        from y = ',F9.1,' m to y = ',F9.1,' m')
1857366 FORMAT (/'       Topography color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)')
1858367 FORMAT ('       Polygon reduction for topography: cluster_size = ', I1)
1859#endif
1860#if defined( __spectra )
1861370 FORMAT ('    Spectra:')
1862371 FORMAT ('       Output every ',F7.1,' s'/)
1863372 FORMAT ('       Arrays:     ', 10(A5,',')/                         &
1864            '       Directions: ', 10(A5,',')/                         &
1865            '       height levels  k = ', 20(I3,',')/                  &
1866            '                          ', 20(I3,',')/                  &
1867            '                          ', 20(I3,',')/                  &
1868            '                          ', 20(I3,',')/                  &
1869            '                          ', 19(I3,','),I3,'.'/           &
1870            '       height levels selected for standard plot:'/        &
1871            '                      k = ', 20(I3,',')/                  &
1872            '                          ', 20(I3,',')/                  &
1873            '                          ', 20(I3,',')/                  &
1874            '                          ', 20(I3,',')/                  &
1875            '                          ', 19(I3,','),I3,'.'/           &
1876            '       Time averaged over ', F7.1, ' s,' /                &
1877            '       Profiles for the time averaging are taken every ', &
1878                    F6.1,' s')
1879#endif
1880400 FORMAT (//' Physical quantities:'/ &
1881              ' -------------------'/)
1882410 FORMAT ('    Angular velocity    :   omega = ',E9.3,' rad/s'/  &
1883            '    Geograph. latitude  :   phi   = ',F4.1,' degr'/   &
1884            '    Coriolis parameter  :   f     = ',F9.6,' 1/s'/    &
1885            '                            f*    = ',F9.6,' 1/s')
1886411 FORMAT (/'    Gravity             :   g     = ',F4.1,' m/s**2')
1887412 FORMAT (/'    Reference density in buoyancy terms: ',F8.3,' kg/m**3')
1888413 FORMAT (/'    Reference temperature in buoyancy terms: ',F8.4,' K')
1889415 FORMAT (/'    Cloud physics parameters:'/ &
1890             '    ------------------------'/)
1891416 FORMAT ('        Surface pressure   :   p_0   = ',F7.2,' hPa'/      &
1892            '        Gas constant       :   R     = ',F5.1,' J/(kg K)'/ &
1893            '        Density of air     :   rho_0 = ',F5.3,' kg/m**3'/  &
1894            '        Specific heat cap. :   c_p   = ',F6.1,' J/(kg K)'/ &
1895            '        Vapourization heat :   L_v   = ',E8.2,' J/kg')
1896420 FORMAT (/'    Characteristic levels of the initial temperature profile:'// &
1897            '       Height:        ',A,'  m'/ &
1898            '       Temperature:   ',A,'  K'/ &
1899            '       Gradient:      ',A,'  K/100m'/ &
1900            '       Gridpoint:     ',A)
1901421 FORMAT (/'    Characteristic levels of the initial humidity profile:'// &
1902            '       Height:      ',A,'  m'/ &
1903            '       Humidity:    ',A,'  kg/kg'/ &
1904            '       Gradient:    ',A,'  (kg/kg)/100m'/ &
1905            '       Gridpoint:   ',A)
1906422 FORMAT (/'    Characteristic levels of the initial scalar profile:'// &
1907            '       Height:                  ',A,'  m'/ &
1908            '       Scalar concentration:    ',A,'  kg/m**3'/ &
1909            '       Gradient:                ',A,'  (kg/m**3)/100m'/ &
1910            '       Gridpoint:               ',A)
1911423 FORMAT (/'    Characteristic levels of the geo. wind component ug:'// &
1912            '       Height:      ',A,'  m'/ &
1913            '       ug:          ',A,'  m/s'/ &
1914            '       Gradient:    ',A,'  1/100s'/ &
1915            '       Gridpoint:   ',A)
1916424 FORMAT (/'    Characteristic levels of the geo. wind component vg:'// &
1917            '       Height:      ',A,'  m'/ &
1918            '       vg:          ',A,'  m/s'/ &
1919            '       Gradient:    ',A,'  1/100s'/ &
1920            '       Gridpoint:   ',A)
1921425 FORMAT (/'    Characteristic levels of the initial salinity profile:'// &
1922            '       Height:     ',A,'  m'/ &
1923            '       Salinity:   ',A,'  psu'/ &
1924            '       Gradient:   ',A,'  psu/100m'/ &
1925            '       Gridpoint:  ',A)
1926426 FORMAT (/'    Characteristic levels of the subsidence/ascent profile:'// &
1927            '       Height:      ',A,'  m'/ &
1928            '       w_subs:      ',A,'  m/s'/ &
1929            '       Gradient:    ',A,'  (m/s)/100m'/ &
1930            '       Gridpoint:   ',A)
1931427 FORMAT (/'    Initial wind profiles (u,v) are interpolated from given'// &
1932                  ' profiles')
1933450 FORMAT (//' LES / Turbulence quantities:'/ &
1934              ' ---------------------------'/)
1935451 FORMAT ('   Diffusion coefficients are constant:'/ &
1936            '   Km = ',F6.2,' m**2/s   Kh = ',F6.2,' m**2/s   Pr = ',F5.2)
1937452 FORMAT ('   Mixing length is limited to the Prandtl mixing lenth.')
1938453 FORMAT ('   Mixing length is limited to ',F4.2,' * z')
1939454 FORMAT ('   TKE is not allowed to fall below ',E9.2,' (m/s)**2')
1940455 FORMAT ('   initial TKE is prescribed as ',E9.2,' (m/s)**2')
1941470 FORMAT (//' Actions during the simulation:'/ &
1942              ' -----------------------------'/)
1943471 FORMAT ('    Disturbance impulse (u,v) every :   ',F6.2,' s'/            &
1944            '    Disturbance amplitude           :     ',F4.2, ' m/s'/       &
1945            '    Lower disturbance level         : ',F8.2,' m (GP ',I4,')'/  &
1946            '    Upper disturbance level         : ',F8.2,' m (GP ',I4,')')
1947472 FORMAT ('    Disturbances continued during the run from i/j =',I4, &
1948                 ' to i/j =',I4)
1949473 FORMAT ('    Disturbances cease as soon as the disturbance energy exceeds',&
1950                 1X,F5.3, ' m**2/s**2')
1951474 FORMAT ('    Random number generator used    : ',A/)
1952475 FORMAT ('    The surface temperature is increased (or decreased, ', &
1953                 'respectively, if'/ &
1954            '    the value is negative) by ',F5.2,' K at the beginning of the',&
1955                 ' 3D-simulation'/)
1956476 FORMAT ('    The surface humidity is increased (or decreased, ',&
1957                 'respectively, if the'/ &
1958            '    value is negative) by ',E8.1,' kg/kg at the beginning of', &
1959                 ' the 3D-simulation'/)
1960477 FORMAT ('    The scalar value is increased at the surface (or decreased, ',&
1961                 'respectively, if the'/ &
1962            '    value is negative) by ',E8.1,' kg/m**3 at the beginning of', &
1963                 ' the 3D-simulation'/)
1964480 FORMAT ('    Particles:'/ &
1965            '    ---------'// &
1966            '       Particle advection is active (switched on at t = ', F7.1, &
1967                    ' s)'/ &
1968            '       Start of new particle generations every  ',F6.1,' s'/ &
1969            '       Boundary conditions: left/right: ', A, ' north/south: ', A/&
1970            '                            bottom:     ', A, ' top:         ', A/&
1971            '       Maximum particle age:                 ',F9.1,' s'/ &
1972            '       Advection stopped at t = ',F9.1,' s'/ &
1973            '       Particles are sorted every ',F9.1,' s'/)
1974481 FORMAT ('       Particles have random start positions'/)
1975482 FORMAT ('          Particles are advected only horizontally'/)
1976483 FORMAT ('       Particles have tails with a maximum of ',I3,' points')
1977484 FORMAT ('            Number of tails of the total domain: ',I10/ &
1978            '            Minimum distance between tailpoints: ',F8.2,' m'/ &
1979            '            Maximum age of the end of the tail:  ',F8.2,' s')
1980485 FORMAT ('       Particle data are written on file every ', F9.1, ' s')
1981486 FORMAT ('       Particle statistics are written on file'/)
1982487 FORMAT ('       Number of particle groups: ',I2/)
1983488 FORMAT ('       SGS velocity components are used for particle advection'/ &
1984            '          minimum timestep for advection: ', F7.5/)
1985489 FORMAT ('       Number of particles simultaneously released at each ', &
1986                    'point: ', I5/)
1987490 FORMAT ('       Particle group ',I2,':'/ &
1988            '          Particle radius: ',E10.3, 'm')
1989491 FORMAT ('          Particle inertia is activated'/ &
1990            '             density_ratio (rho_fluid/rho_particle) = ',F5.3/)
1991492 FORMAT ('          Particles are advected only passively (no inertia)'/)
1992493 FORMAT ('          Boundaries of particle source: x:',F8.1,' - ',F8.1,' m'/&
1993            '                                         y:',F8.1,' - ',F8.1,' m'/&
1994            '                                         z:',F8.1,' - ',F8.1,' m'/&
1995            '          Particle distances:  dx = ',F8.1,' m  dy = ',F8.1, &
1996                       ' m  dz = ',F8.1,' m'/)
1997494 FORMAT ('       Output of particle time series in NetCDF format every ', &
1998                    F8.2,' s'/)
1999495 FORMAT ('       Number of particles in total domain: ',I10/)
2000500 FORMAT (//' 1D-Model parameters:'/                           &
2001              ' -------------------'//                           &
2002            '    Simulation time:                   ',F8.1,' s'/ &
2003            '    Run-controll output every:         ',F8.1,' s'/ &
2004            '    Vertical profile output every:     ',F8.1,' s'/ &
2005            '    Mixing length calculation:         ',A/         &
2006            '    Dissipation calculation:           ',A/)
2007502 FORMAT ('    Damping layer starts from ',F7.1,' m (GP ',I4,')'/)
2008503 FORMAT (' --> Momentum advection via Wicker-Skamarock-Scheme 5th order')
2009504 FORMAT (' --> Scalar advection via Wicker-Skamarock-Scheme 5th order')
2010
2011
2012 END SUBROUTINE header
Note: See TracBrowser for help on using the repository browser.