source: palm/trunk/SOURCE/write_var_list.f90 @ 2846

Last change on this file since 2846 was 2818, checked in by maronga, 7 years ago

improvements in spinup mechanism

  • Property svn:keywords set to Id
File size: 30.0 KB
Line 
1!> @file write_var_list.f90
2!------------------------------------------------------------------------------!
3! This file is part of the PALM model system.
4!
5! PALM is free software: you can redistribute it and/or modify it under the
6! terms of the GNU General Public License as published by the Free Software
7! Foundation, either version 3 of the License, or (at your option) any later
8! version.
9!
10! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
11! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
12! A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
13!
14! You should have received a copy of the GNU General Public License along with
15! PALM. If not, see <http://www.gnu.org/licenses/>.
16!
17! Copyright 1997-2018 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: write_var_list.f90 2818 2018-02-19 16:42:36Z raasch $
27! Added spinup_time to allow for restarts with spinup in precursor run.
28!
29! 2776 2018-01-31 10:44:42Z Giersch
30! Variable synthetic_turbulence_generator has been abbreviated
31!
32! 2718 2018-01-02 08:49:38Z maronga
33! Corrected "Former revisions" section
34!
35! 2696 2017-12-14 17:12:51Z kanani
36! Change in file header (GPL part)
37! + rans_mode, turbulence_closure (TG)
38! Implementation of chemistry module (FK)
39! Removal of chemistry module part. It is not needed and causes bugs (FK)
40!
41! 2576 2017-10-24 13:49:46Z Giersch
42! The end of the standard parameter list is again marked with *** end *** like
43! it was before revision 2563
44!
45! 2575 2017-10-24 09:57:58Z maronga
46! Renamed phi -> latitude, added longitude
47!
48! 2574 2017-10-23 09:02:53Z Giersch
49! Bugfix: Increased length of the character string *** end ***
50!
51! 2563 2017-10-19 15:36:10Z Giersch
52! Function call wtm_write_restart_data was added and the end of the standard
53! parameter list is now marked with *** end default ***. The end of the whole
54! parameter list including module parameter is marked with *** end ***
55!
56! 2372 2017-08-25 12:37:32Z sward
57! y_shift added to vars, version no. increased
58!
59! 2365 2017-08-21 14:59:59Z kanani
60! Vertical nesting implemented (SadiqHuq)
61!
62! 2339 2017-08-07 13:55:26Z gronemeier
63! corrected timestamp in header
64!
65! 2338 2017-08-07 12:15:38Z gronemeier
66! Modularize 1D model
67!
68! 2265 2017-06-08 16:58:28Z schwenkel
69! Unused variables removed.
70!
71! 2259 2017-06-08 09:09:11Z gronemeier
72! Implemented synthetic turbulence generator
73!
74! 2233 2017-05-30 18:08:54Z suehring
75!
76! 2232 2017-05-30 17:47:52Z suehring
77! Replace wall_qflux, wall_sflux by wall_humidityflux and wall_scalarflux; add
78! wall_salinityflux
79! +tunnel_height, tunnel_lenth, tunnel_width_x, tunnel_width_y,
80! tunnel_wall_depth
81!
82! 2042 2016-11-02 13:47:31Z suehring
83! Bugfix, write restart data for wall_heatflux, wall_qflux and wall_sflux
84!
85! 2000 2016-08-20 18:09:15Z knoop
86! Forced header and separation lines into 80 columns
87!
88! 1992 2016-08-12 15:14:59Z suehring
89! top scalarflux added
90!
91! 1957 2016-07-07 10:43:48Z suehring
92! flight module added
93!
94! 1849 2016-04-08 11:33:18Z hoffmann
95! Adapted for modularization of microphysics
96!
97! 1833 2016-04-07 14:23:03Z raasch
98! spectra_mod added
99!
100! 1831 2016-04-07 13:15:51Z hoffmann
101! turbulence renamed collision_turbulence, drizzle renamed
102! cloud_water_sedimentation
103!
104! 1783 2016-03-06 18:36:17Z raasch
105! netcdf module name changed + related changes
106!
107! 1705 2015-11-02 14:28:56Z maronga
108! Bugfix: two lines required swapping
109!
110! 1691 2015-10-26 16:17:44Z maronga
111! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed
112! output of prandtl_layer and rif_min, rif_max.
113!
114! 1682 2015-10-07 23:56:08Z knoop
115! Code annotations made doxygen readable
116!
117! 1585 2015-04-30 07:05:52Z maronga
118! Adapted for RRTMG
119!
120! 1551 2015-03-03 14:18:16Z maronga
121! Typo removed
122!
123! 1502 2014-12-03 18:22:31Z kanani
124! Canopy module and parameters removed (parameters are always read from
125! canopy_par NAMELIST for initial and restart runs),
126! Bugfix: added blanks in "cloud_top_radiation"-string to a total of 30
127! characters
128!
129! 1496 2014-12-02 17:25:50Z maronga
130! Renamed "radiation" -> "cloud_top_radiation"
131!
132! 1484 2014-10-21 10:53:05Z kanani
133! Changes in the course of the canopy-model modularization:
134!   parameters alpha_lad, beta_lad, lai_beta added,
135!   module plant_canopy_model_mod added,
136!   drag_coefficient, leaf_surface_concentration and scalar_exchange_coefficient
137!   renamed to canopy_drag_coeff, leaf_surface_conc and leaf_scalar_exch_coeff
138!
139! 1324 2014-03-21 09:13:16Z suehring
140! Bugfix: ONLY statement for module netcdf_control removed
141!
142! 1320 2014-03-20 08:40:49Z raasch
143! revision history before 2012 removed,
144! comment fields (!:) to be used for variable explanations added to
145! all variable declaration statements
146!
147! 1308 2014-03-13 14:58:42Z fricke
148! +do2d_xy_time_count, do2d_xz_time_count, do2d_yz_time_count,
149! +do3d_time_count
150!
151! 1241 2013-10-30 11:36:58Z heinze
152! +nudging
153! +large_scale_forcing
154!
155! 1179 2013-06-14 05:57:58Z raasch
156! +reference_state, ref_state
157!
158! 1115 2013-03-26 18:16:16Z hoffmann
159! unused variables removed
160!
161! 1065 2012-11-22 17:42:36Z hoffmann
162! +nc, c_sedimentation, turbulence, limiter_sedimentation
163! -mu_constant, mu_constant_value
164!
165! 1053 2012-11-13 17:11:03Z hoffmann
166! necessary expansions according to the two new prognostic equations (nr, qr)
167! of the two-moment cloud physics scheme:
168! +bc_*_b, bc_*_t, bc_*_t_val, *_init, *_surface, *_surface_initial_change,
169! +*_vertical_gradient, *_vertical_gradient_level, *_vertical_gradient_level_ind,
170! +surface_waterflux_*
171!
172! in addition, steering parameters parameters of the two-moment cloud physics
173! scheme:
174! +cloud_scheme, +drizzle, +mu_constant, +mu_constant_value, +ventilation_effect
175!
176! 1036 2012-10-22 13:43:42Z raasch
177! code put under GPL (PALM 3.9)
178!
179! 1015 2012-09-27 09:23:24Z raasch
180! -adjust_mixing_length
181!
182! 1003 2012-09-14 14:35:53Z raasch
183! -grid_matching
184!
185! 1001 2012-09-13 14:08:46Z raasch
186! -cut_spline_overshoot, last_dt_change, long_filter_factor, overshoot_limit_*,
187! ups_limit_*
188!
189! 978 2012-08-09 08:28:32Z fricke
190! -km_damp_max, outflow_damping_width
191! +pt_damping_factor, pt_damping_width
192! +z0h_factor
193!
194! 940 2012-07-09 14:31:00Z raasch
195! +neutral
196!
197! 927 2012-06-06 19:15:04Z raasch
198! +masking_method
199!
200! 849 2012-03-15 10:35:09Z raasch
201! first_call_advec_particles renamed first_call_lpm
202!
203! 824 2012-02-17 09:09:57Z raasch
204! +curvature_solution_effects
205!
206! Revision 1.1  1998/03/18 20:20:38  raasch
207! Initial revision
208!
209!
210! Description:
211! ------------
212!> Writing values of control variables to restart-file (binary format).
213!> This information are only written to the file opened by PE0.
214!------------------------------------------------------------------------------!
215 SUBROUTINE write_var_list
216
217
218    USE arrays_3d,                                                             &
219        ONLY:  inflow_damping_factor, mean_inflow_profiles, pt_init,           &
220               q_init, ref_state, s_init, sa_init, u_init, ug, v_init, vg
221
222    USE control_parameters
223
224    USE flight_mod,                                                            &
225        ONLY:  flight_write_restart_data
226
227    USE grid_variables,                                                        &
228        ONLY:  dx, dy
229
230    USE indices,                                                               &
231        ONLY:  nz, nx, ny
232
233    USE microphysics_mod,                                                      &
234        ONLY:  c_sedimentation, cloud_water_sedimentation,                     &
235               collision_turbulence, limiter_sedimentation, nc_const,          &
236               ventilation_effect
237
238    USE model_1d_mod,                                                          &
239        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
240
241    USE netcdf_interface,                                                      &
242        ONLY:  netcdf_precision, output_for_t0
243
244    USE particle_attributes,                                                   &
245        ONLY:  curvature_solution_effects
246
247    USE pegrid
248
249    USE radiation_model_mod,                                                   &
250        ONLY:  time_radiation
251
252    USE spectra_mod,                                                           &
253        ONLY:  average_count_sp
254
255    USE synthetic_turbulence_generator_mod,                                    &
256        ONLY:  stg_write_restart_data
257
258    USE statistics,                                                            &
259        ONLY:  statistic_regions, hom, hom_sum, u_max, u_max_ijk, v_max,       &
260               v_max_ijk, w_max, w_max_ijk
261
262    USE vertical_nesting_mod,                                                  &
263        ONLY:  vnest_init
264
265    USE wind_turbine_model_mod,                                                &
266        ONLY:  wtm_write_restart_data 
267
268    IMPLICIT NONE
269
270    CHARACTER (LEN=10) ::  binary_version   !<
271
272
273    binary_version = '4.6'
274
275    WRITE ( 14 )  binary_version
276
277    WRITE ( 14 )  'numprocs                      '
278    WRITE ( 14 )  numprocs
279    WRITE ( 14 )  'hor_index_bounds              '
280    WRITE ( 14 )  hor_index_bounds
281    WRITE ( 14 )  'nz                            '
282    WRITE ( 14 )  nz
283    WRITE ( 14 )  'max_pr_user                   '
284    WRITE ( 14 )  max_pr_user
285    WRITE ( 14 )  'statistic_regions             '
286    WRITE ( 14 )  statistic_regions
287
288!
289!-- Caution: After changes in the following parameter-list, the
290!-- -------  version number stored in the variable binary_version has to be
291!--          increased. The same changes must also be done in the parameter-
292!--          list in read_var_list.
293
294    WRITE ( 14 )  'advected_distance_x           '
295    WRITE ( 14 )  advected_distance_x
296    WRITE ( 14 )  'advected_distance_y           '
297    WRITE ( 14 )  advected_distance_y
298    WRITE ( 14 )  'alpha_surface                 '
299    WRITE ( 14 )  alpha_surface
300    WRITE ( 14 )  'average_count_pr              '
301    WRITE ( 14 )  average_count_pr
302    WRITE ( 14 )  'average_count_sp              '
303    WRITE ( 14 )  average_count_sp
304    WRITE ( 14 )  'average_count_3d              '
305    WRITE ( 14 )  average_count_3d
306    WRITE ( 14 )  'bc_e_b                        '
307    WRITE ( 14 )  bc_e_b
308    WRITE ( 14 )  'bc_lr                         '
309    WRITE ( 14 )  bc_lr
310    WRITE ( 14 )  'bc_ns                         '
311    WRITE ( 14 )  bc_ns
312    WRITE ( 14 )  'bc_p_b                        '
313    WRITE ( 14 )  bc_p_b
314    WRITE ( 14 )  'bc_p_t                        '
315    WRITE ( 14 )  bc_p_t
316    WRITE ( 14 )  'bc_pt_b                       '
317    WRITE ( 14 )  bc_pt_b
318    WRITE ( 14 )  'bc_pt_t                       '
319    WRITE ( 14 )  bc_pt_t
320    WRITE ( 14 )  'bc_pt_t_val                   '
321    WRITE ( 14 )  bc_pt_t_val
322    WRITE ( 14 )  'bc_q_b                        '
323    WRITE ( 14 )  bc_q_b
324    WRITE ( 14 )  'bc_q_t                        '
325    WRITE ( 14 )  bc_q_t
326    WRITE ( 14 )  'bc_q_t_val                    '
327    WRITE ( 14 )  bc_q_t_val
328    WRITE ( 14 )  'bc_s_b                        '
329    WRITE ( 14 )  bc_s_b
330    WRITE ( 14 )  'bc_s_t                        '
331    WRITE ( 14 )  bc_s_t
332    WRITE ( 14 )  'bc_sa_t                       '
333    WRITE ( 14 )  bc_sa_t
334    WRITE ( 14 )  'bc_uv_b                       '
335    WRITE ( 14 )  bc_uv_b
336    WRITE ( 14 )  'bc_uv_t                       '
337    WRITE ( 14 )  bc_uv_t
338    WRITE ( 14 )  'bottom_salinityflux           '
339    WRITE ( 14 )  bottom_salinityflux
340    WRITE ( 14 )  'building_height               '
341    WRITE ( 14 )  building_height
342    WRITE ( 14 )  'building_length_x             '
343    WRITE ( 14 )  building_length_x
344    WRITE ( 14 )  'building_length_y             '
345    WRITE ( 14 )  building_length_y
346    WRITE ( 14 )  'building_wall_left            '
347    WRITE ( 14 )  building_wall_left
348    WRITE ( 14 )  'building_wall_south           '
349    WRITE ( 14 )  building_wall_south
350    WRITE ( 14 )  'call_psolver_at_all_substeps  '
351    WRITE ( 14 )  call_psolver_at_all_substeps
352    WRITE ( 14 )  'canyon_height                 '
353    WRITE ( 14 )  canyon_height
354    WRITE ( 14 )  'canyon_width_x                '
355    WRITE ( 14 )  canyon_width_x
356    WRITE ( 14 )  'canyon_width_y                '
357    WRITE ( 14 )  canyon_width_y
358    WRITE ( 14 )  'canyon_wall_left              '
359    WRITE ( 14 )  canyon_wall_left
360    WRITE ( 14 )  'canyon_wall_south             '
361    WRITE ( 14 )  canyon_wall_south
362    WRITE ( 14 )  'c_sedimentation               '
363    WRITE ( 14 )  c_sedimentation
364    WRITE ( 14 )  'cfl_factor                    '
365    WRITE ( 14 )  cfl_factor
366    WRITE ( 14 )  'cloud_droplets                '
367    WRITE ( 14 )  cloud_droplets
368    WRITE ( 14 )  'cloud_physics                 '
369    WRITE ( 14 )  cloud_physics
370    WRITE ( 14 )  'cloud_scheme                  '
371    WRITE ( 14 )  cloud_scheme
372    WRITE ( 14 )  'collective_wait               '
373    WRITE ( 14 )  collective_wait
374    WRITE ( 14 )  'conserve_volume_flow          '
375    WRITE ( 14 )  conserve_volume_flow
376    WRITE ( 14 )  'conserve_volume_flow_mode     '
377    WRITE ( 14 )  conserve_volume_flow_mode
378    WRITE ( 14 )  'coupling_start_time           '
379    WRITE ( 14 )  coupling_start_time
380    WRITE ( 14 )  'constant_flux_layer           '
381    WRITE ( 14 )  constant_flux_layer
382    WRITE ( 14 )  'current_timestep_number       '
383    WRITE ( 14 )  current_timestep_number
384    WRITE ( 14 )  'curvature_solution_effects    '
385    WRITE ( 14 )  curvature_solution_effects
386    WRITE ( 14 )  'cycle_mg                      '
387    WRITE ( 14 )  cycle_mg
388    WRITE ( 14 )  'damp_level_1d                 '
389    WRITE ( 14 )  damp_level_1d
390    WRITE ( 14 )  'dissipation_1d                '
391    WRITE ( 14 )  dissipation_1d
392    WRITE ( 14 )  'do2d_xy_time_count            '
393    WRITE ( 14 )  do2d_xy_time_count
394    WRITE ( 14 )  'do2d_xz_time_count            '
395    WRITE ( 14 )  do2d_xz_time_count
396    WRITE ( 14 )  'do2d_yz_time_count            '
397    WRITE ( 14 )  do2d_yz_time_count
398    WRITE ( 14 )  'do3d_time_count               '
399    WRITE ( 14 )  do3d_time_count
400    WRITE ( 14 )  'dp_external                   '
401    WRITE ( 14 )  dp_external
402    WRITE ( 14 )  'dp_level_b                    '
403    WRITE ( 14 )  dp_level_b
404    WRITE ( 14 )  'dp_smooth                     '
405    WRITE ( 14 )  dp_smooth
406    WRITE ( 14 )  'dpdxy                         '
407    WRITE ( 14 )  dpdxy
408    WRITE ( 14 )  'cloud_water_sedimentation     '
409    WRITE ( 14 )  cloud_water_sedimentation
410    WRITE ( 14 )  'dt_pr_1d                      '
411    WRITE ( 14 )  dt_pr_1d
412    WRITE ( 14 )  'dt_run_control_1d             '
413    WRITE ( 14 )  dt_run_control_1d
414    WRITE ( 14 )  'dt_3d                         '
415    WRITE ( 14 )  dt_3d
416    WRITE ( 14 )  'dvrp_filecount                '
417    WRITE ( 14 )  dvrp_filecount
418    WRITE ( 14 )  'dx                            '
419    WRITE ( 14 )  dx
420    WRITE ( 14 )  'dy                            '
421    WRITE ( 14 )  dy
422    WRITE ( 14 )  'dz                            '
423    WRITE ( 14 )  dz
424    WRITE ( 14 )  'dz_max                        '
425    WRITE ( 14 )  dz_max
426    WRITE ( 14 )  'dz_stretch_factor             '
427    WRITE ( 14 )  dz_stretch_factor
428    WRITE ( 14 )  'dz_stretch_level              '
429    WRITE ( 14 )  dz_stretch_level
430    WRITE ( 14 )  'e_min                         '
431    WRITE ( 14 )  e_min
432    WRITE ( 14 )  'end_time_1d                   '
433    WRITE ( 14 )  end_time_1d
434    WRITE ( 14 )  'fft_method                    '
435    WRITE ( 14 )  fft_method
436    WRITE ( 14 )  'first_call_lpm                '
437    WRITE ( 14 )  first_call_lpm
438    WRITE ( 14 )  'galilei_transformation        '
439    WRITE ( 14 )  galilei_transformation
440    WRITE ( 14 )  'hom                           '
441    WRITE ( 14 )  hom
442    WRITE ( 14 )  'hom_sum                       '
443    WRITE ( 14 )  hom_sum
444    WRITE ( 14 )  'humidity                      '
445    WRITE ( 14 )  humidity
446    IF ( ALLOCATED( inflow_damping_factor ) )  THEN
447       WRITE ( 14 )  'inflow_damping_factor         '
448       WRITE ( 14 )  inflow_damping_factor
449    ENDIF
450    WRITE ( 14 )  'inflow_damping_height         '
451    WRITE ( 14 )  inflow_damping_height
452    WRITE ( 14 )  'inflow_damping_width          '
453    WRITE ( 14 )  inflow_damping_width
454    WRITE ( 14 )  'inflow_disturbance_begin      '
455    WRITE ( 14 )  inflow_disturbance_begin
456    WRITE ( 14 )  'inflow_disturbance_end        '
457    WRITE ( 14 )  inflow_disturbance_end
458    WRITE ( 14 )  'km_constant                   '
459    WRITE ( 14 )  km_constant
460    WRITE ( 14 )  'large_scale_forcing           '
461    WRITE ( 14 )  large_scale_forcing
462    WRITE ( 14 )  'large_scale_subsidence        '
463    WRITE ( 14 )  large_scale_subsidence
464    WRITE ( 14 )  'limiter_sedimentation         '
465    WRITE ( 14 )  limiter_sedimentation
466    WRITE ( 14 )  'loop_optimization             '
467    WRITE ( 14 )  loop_optimization
468    WRITE ( 14 )  'masking_method                '
469    WRITE ( 14 )  masking_method
470    IF ( ALLOCATED( mean_inflow_profiles ) )  THEN
471       WRITE ( 14 )  'mean_inflow_profiles          '
472       WRITE ( 14 )  mean_inflow_profiles
473    ENDIF
474    WRITE ( 14 )  'mg_cycles                     '
475    WRITE ( 14 )  mg_cycles
476    WRITE ( 14 )  'mg_switch_to_pe0_level        '
477    WRITE ( 14 )  mg_switch_to_pe0_level
478    WRITE ( 14 )  'mixing_length_1d              '
479    WRITE ( 14 )  mixing_length_1d
480    WRITE ( 14 )  'momentum_advec                '
481    WRITE ( 14 )  momentum_advec
482    WRITE ( 14 )  'most_method                   '
483    WRITE ( 14 )  most_method
484    WRITE ( 14 )  'nc_const                      '
485    WRITE ( 14 )  nc_const
486    WRITE ( 14 )  'netcdf_precision              '
487    WRITE ( 14 )  netcdf_precision
488    WRITE ( 14 )  'neutral                       '
489    WRITE ( 14 )  neutral
490    WRITE ( 14 )  'ngsrb                         '
491    WRITE ( 14 )  ngsrb
492    WRITE ( 14 )  'nsor                          '
493    WRITE ( 14 )  nsor
494    WRITE ( 14 )  'nsor_ini                      '
495    WRITE ( 14 )  nsor_ini
496    WRITE ( 14 )  'nudging                       '
497    WRITE ( 14 )  nudging
498    WRITE ( 14 )  'num_leg                       '
499    WRITE ( 14 )  num_leg
500    WRITE ( 14 )  'nx                            '
501    WRITE ( 14 )  nx
502    WRITE ( 14 )  'ny                            '
503    WRITE ( 14 )  ny
504    WRITE ( 14 )  'ocean                         '
505    WRITE ( 14 )  ocean
506    WRITE ( 14 )  'old_dt                        '
507    WRITE ( 14 )  old_dt
508    WRITE ( 14 )  'omega                         '
509    WRITE ( 14 )  omega
510    WRITE ( 14 )  'omega_sor                     '
511    WRITE ( 14 )  omega_sor
512    WRITE ( 14 )  'output_for_t0                 '
513    WRITE ( 14 )  output_for_t0
514    WRITE ( 14 )  'passive_scalar                '
515    WRITE ( 14 )  passive_scalar
516    WRITE ( 14 )  'latitude                      '
517    WRITE ( 14 )  latitude
518    WRITE ( 14 )  'longitude                     '
519    WRITE ( 14 )  longitude
520    WRITE ( 14 )  'prandtl_number                '
521    WRITE ( 14 )  prandtl_number
522    WRITE ( 14 )  'precipitation                 '
523    WRITE ( 14 )  precipitation
524    WRITE ( 14 )  'psolver                       '
525    WRITE ( 14 )  psolver
526    WRITE ( 14 )  'pt_damping_factor             '
527    WRITE ( 14 )  pt_damping_factor
528    WRITE ( 14 )  'pt_damping_width              '
529    WRITE ( 14 )  pt_damping_width
530    WRITE ( 14 )  'pt_init                       '
531    WRITE ( 14 )  pt_init
532    WRITE ( 14 )  'pt_reference                  '
533    WRITE ( 14 )  pt_reference
534    WRITE ( 14 )  'pt_surface                    '
535    WRITE ( 14 )  pt_surface
536    WRITE ( 14 )  'pt_surface_initial_change     '
537    WRITE ( 14 )  pt_surface_initial_change
538    WRITE ( 14 )  'pt_vertical_gradient          '
539    WRITE ( 14 )  pt_vertical_gradient
540    WRITE ( 14 )  'pt_vertical_gradient_level    '
541    WRITE ( 14 )  pt_vertical_gradient_level
542    WRITE ( 14 )  'pt_vertical_gradient_level_ind'
543    WRITE ( 14 )  pt_vertical_gradient_level_ind
544    WRITE ( 14 )  'q_init                        '
545    WRITE ( 14 )  q_init
546    WRITE ( 14 )  'q_surface                     '
547    WRITE ( 14 )  q_surface
548    WRITE ( 14 )  'q_surface_initial_change      '
549    WRITE ( 14 )  q_surface_initial_change
550    WRITE ( 14 )  'q_vertical_gradient           '
551    WRITE ( 14 )  q_vertical_gradient
552    WRITE ( 14 )  'q_vertical_gradient_level     '
553    WRITE ( 14 )  q_vertical_gradient_level
554    WRITE ( 14 )  'q_vertical_gradient_level_ind '
555    WRITE ( 14 )  q_vertical_gradient_level_ind
556    WRITE ( 14 )  'cloud_top_radiation           '
557    WRITE ( 14 )  cloud_top_radiation
558    WRITE ( 14 )  'random_generator              '
559    WRITE ( 14 )  random_generator
560    WRITE ( 14 )  'random_heatflux               '
561    WRITE ( 14 )  random_heatflux
562    WRITE ( 14 )  'rans_mode                     '
563    WRITE ( 14 )  rans_mode
564    WRITE ( 14 )  'rayleigh_damping_factor       '
565    WRITE ( 14 )  rayleigh_damping_factor
566    WRITE ( 14 )  'rayleigh_damping_height       '
567    WRITE ( 14 )  rayleigh_damping_height
568    WRITE ( 14 )  'recycling_width               '
569    WRITE ( 14 )  recycling_width
570    WRITE ( 14 )  'reference_state               '
571    WRITE ( 14 )  reference_state
572    WRITE ( 14 )  'ref_state                     '
573    WRITE ( 14 )  ref_state
574    WRITE ( 14 )  'residual_limit                '
575    WRITE ( 14 )  residual_limit
576    WRITE ( 14 )  'roughness_length              '
577    WRITE ( 14 )  roughness_length
578    WRITE ( 14 )  'runnr                         '
579    WRITE ( 14 )  runnr
580    WRITE ( 14 )  'run_coupled                   '
581    WRITE ( 14 )  run_coupled
582    WRITE ( 14 )  's_init                        '
583    WRITE ( 14 )  s_init
584    WRITE ( 14 )  's_surface                     '
585    WRITE ( 14 )  s_surface
586    WRITE ( 14 )  's_surface_initial_change      '
587    WRITE ( 14 )  s_surface_initial_change
588    WRITE ( 14 )  's_vertical_gradient           '
589    WRITE ( 14 )  s_vertical_gradient
590    WRITE ( 14 )  's_vertical_gradient_level     '
591    WRITE ( 14 )  s_vertical_gradient_level
592    WRITE ( 14 )  's_vertical_gradient_level_ind '
593    WRITE ( 14 )  s_vertical_gradient_level_ind
594    WRITE ( 14 )  'sa_init                       '
595    WRITE ( 14 )  sa_init
596    WRITE ( 14 )  'sa_surface                    '
597    WRITE ( 14 )  sa_surface
598    WRITE ( 14 )  'sa_vertical_gradient          '
599    WRITE ( 14 )  sa_vertical_gradient
600    WRITE ( 14 )  'sa_vertical_gradient_level    '
601    WRITE ( 14 )  sa_vertical_gradient_level
602    WRITE ( 14 )  'scalar_advec                  '
603    WRITE ( 14 )  scalar_advec
604    WRITE ( 14 )  'simulated_time                '
605    WRITE ( 14 )  simulated_time
606    WRITE ( 14 )  'spinup_time                   '
607    WRITE ( 14 )  spinup_time
608    WRITE ( 14 )  'surface_heatflux              '
609    WRITE ( 14 )  surface_heatflux
610    WRITE ( 14 )  'surface_pressure              '
611    WRITE ( 14 )  surface_pressure
612    WRITE ( 14 )  'surface_scalarflux            '
613    WRITE ( 14 )  surface_scalarflux
614    WRITE ( 14 )  'surface_waterflux             '
615    WRITE ( 14 )  surface_waterflux
616    WRITE ( 14 )  's_surface                     '
617    WRITE ( 14 )  s_surface
618    WRITE ( 14 )  's_surface_initial_change      '
619    WRITE ( 14 )  s_surface_initial_change
620    WRITE ( 14 )  's_vertical_gradient           '
621    WRITE ( 14 )  s_vertical_gradient
622    WRITE ( 14 )  's_vertical_gradient_level     '
623    WRITE ( 14 )  s_vertical_gradient_level
624    WRITE ( 14 )  'time_coupling                 '
625    WRITE ( 14 )  time_coupling
626    WRITE ( 14 )  'time_disturb                  '
627    WRITE ( 14 )  time_disturb
628    WRITE ( 14 )  'time_domask                   '
629    WRITE ( 14 )  time_domask
630    WRITE ( 14 )  'time_dopr                     '
631    WRITE ( 14 )  time_dopr
632    WRITE ( 14 )  'time_dopr_av                  '
633    WRITE ( 14 )  time_dopr_av
634    WRITE ( 14 )  'time_dopr_listing             '
635    WRITE ( 14 )  time_dopr_listing
636    WRITE ( 14 )  'time_dopts                    '
637    WRITE ( 14 )  time_dopts
638    WRITE ( 14 )  'time_dosp                     '
639    WRITE ( 14 )  time_dosp
640    WRITE ( 14 )  'time_dots                     '
641    WRITE ( 14 )  time_dots
642    WRITE ( 14 )  'time_do2d_xy                  '
643    WRITE ( 14 )  time_do2d_xy
644    WRITE ( 14 )  'time_do2d_xz                  '
645    WRITE ( 14 )  time_do2d_xz
646    WRITE ( 14 )  'time_do2d_yz                  '
647    WRITE ( 14 )  time_do2d_yz
648    WRITE ( 14 )  'time_do3d                     '
649    WRITE ( 14 )  time_do3d
650    WRITE ( 14 )  'time_do_av                    '
651    WRITE ( 14 )  time_do_av
652    WRITE ( 14 )  'time_do_sla                   '
653    WRITE ( 14 )  time_do_sla
654    WRITE ( 14 )  'time_dvrp                     '
655    WRITE ( 14 )  time_dvrp
656    WRITE ( 14 )  'time_radiation                '
657    WRITE ( 14 )  time_radiation
658    WRITE ( 14 )  'time_restart                  '
659    WRITE ( 14 )  time_restart
660    WRITE ( 14 )  'time_run_control              '
661    WRITE ( 14 )  time_run_control
662    WRITE ( 14 )  'time_since_reference_point    '
663    WRITE ( 14 )  time_since_reference_point
664    WRITE ( 14 )  'timestep_scheme               '
665    WRITE ( 14 )  timestep_scheme
666    WRITE ( 14 )  'topography                    '
667    WRITE ( 14 )  topography
668    WRITE ( 14 )  'topography_grid_convention    '
669    WRITE ( 14 )  topography_grid_convention
670    WRITE ( 14 )  'top_heatflux                  '
671    WRITE ( 14 )  top_heatflux
672    WRITE ( 14 )  'top_momentumflux_u            '
673    WRITE ( 14 )  top_momentumflux_u
674    WRITE ( 14 )  'top_momentumflux_v            '
675    WRITE ( 14 )  top_momentumflux_v
676    WRITE ( 14 )  'top_salinityflux              '
677    WRITE ( 14 )  top_salinityflux
678    WRITE ( 14 )  'top_scalarflux                '
679    WRITE ( 14 )  top_scalarflux
680    WRITE ( 14 )  'tsc                           '
681    WRITE ( 14 )  tsc
682    WRITE ( 14 )  'tunnel_height                 '
683    WRITE ( 14 )  tunnel_height
684    WRITE ( 14 )  'tunnel_length                 '
685    WRITE ( 14 )  tunnel_length
686    WRITE ( 14 )  'tunnel_width_x                '
687    WRITE ( 14 )  tunnel_width_x
688    WRITE ( 14 )  'tunnel_width_y                '
689    WRITE ( 14 )  tunnel_width_y
690    WRITE ( 14 )  'tunnel_wall_depth             '
691    WRITE ( 14 )  tunnel_wall_depth
692    WRITE ( 14 )  'turbulence_closure            '
693    WRITE ( 14 )  turbulence_closure
694    WRITE ( 14 )  'collision_turbulence          '
695    WRITE ( 14 )  collision_turbulence
696    WRITE ( 14 )  'turbulent_inflow              '
697    WRITE ( 14 )  turbulent_inflow
698    WRITE ( 14 )  'u_bulk                        '
699    WRITE ( 14 )  u_bulk
700    WRITE ( 14 )  'u_init                        '
701    WRITE ( 14 )  u_init
702    WRITE ( 14 )  'u_max                         '
703    WRITE ( 14 )  u_max
704    WRITE ( 14 )  'u_max_ijk                     '
705    WRITE ( 14 )  u_max_ijk
706    WRITE ( 14 )  'ug                            '
707    WRITE ( 14 )  ug
708    WRITE ( 14 )  'ug_surface                    '
709    WRITE ( 14 )  ug_surface
710    WRITE ( 14 )  'ug_vertical_gradient          '
711    WRITE ( 14 )  ug_vertical_gradient
712    WRITE ( 14 )  'ug_vertical_gradient_level    '
713    WRITE ( 14 )  ug_vertical_gradient_level
714    WRITE ( 14 )  'ug_vertical_gradient_level_ind'
715    WRITE ( 14 )  ug_vertical_gradient_level_ind
716    WRITE ( 14 )  'use_surface_fluxes            '
717    WRITE ( 14 )  use_surface_fluxes
718    WRITE ( 14 )  'use_top_fluxes                '
719    WRITE ( 14 )  use_top_fluxes
720    WRITE ( 14 )  'use_ug_for_galilei_tr         '
721    WRITE ( 14 )  use_ug_for_galilei_tr
722    WRITE ( 14 )  'use_upstream_for_tke          '
723    WRITE ( 14 )  use_upstream_for_tke
724    WRITE ( 14 )  'v_bulk                        '
725    WRITE ( 14 )  v_bulk
726    WRITE ( 14 )  'v_init                        '
727    WRITE ( 14 )  v_init
728    WRITE ( 14 )  'v_max                         '
729    WRITE ( 14 )  v_max
730    WRITE ( 14 )  'v_max_ijk                     '
731    WRITE ( 14 )  v_max_ijk
732    WRITE ( 14 )  'ventilation_effect            '
733    WRITE ( 14 )  ventilation_effect
734    WRITE ( 14 )  'vg                            '
735    WRITE ( 14 )  vg
736    WRITE ( 14 )  'vg_surface                    '
737    WRITE ( 14 )  vg_surface
738    WRITE ( 14 )  'vg_vertical_gradient          '
739    WRITE ( 14 )  vg_vertical_gradient
740    WRITE ( 14 )  'vg_vertical_gradient_level    '
741    WRITE ( 14 )  vg_vertical_gradient_level
742    WRITE ( 14 )  'vg_vertical_gradient_level_ind'
743    WRITE ( 14 )  vg_vertical_gradient_level_ind
744    WRITE ( 14 )  'virtual_flight                '
745    WRITE ( 14 )  virtual_flight
746    WRITE ( 14 )  'vnest_init                    '
747    WRITE ( 14 )  vnest_init
748    WRITE ( 14 )  'volume_flow_area              '
749    WRITE ( 14 )  volume_flow_area
750    WRITE ( 14 )  'volume_flow_initial           '
751    WRITE ( 14 )  volume_flow_initial
752    WRITE ( 14 )  'wall_adjustment               '
753    WRITE ( 14 )  wall_adjustment
754    WRITE ( 14 )  'subs_vertical_gradient        '
755    WRITE ( 14 )  subs_vertical_gradient
756    WRITE ( 14 )  'subs_vertical_gradient_level  '
757    WRITE ( 14 )  subs_vertical_gradient_level
758    WRITE ( 14 )  'subs_vertical_gradient_level_i'
759    WRITE ( 14 )  subs_vertical_gradient_level_i
760    WRITE ( 14 )  'wall_heatflux                 '
761    WRITE ( 14 )  wall_heatflux
762    WRITE ( 14 )  'wall_humidityflux             '
763    WRITE ( 14 )  wall_humidityflux
764    WRITE ( 14 )  'wall_scalarflux               '
765    WRITE ( 14 )  wall_scalarflux
766    WRITE ( 14 )  'wall_salinityflux             '
767    WRITE ( 14 )  wall_salinityflux
768    WRITE ( 14 )  'w_max                         '
769    WRITE ( 14 )  w_max
770    WRITE ( 14 )  'w_max_ijk                     '
771    WRITE ( 14 )  w_max_ijk
772    WRITE ( 14 )  'y_shift                       '
773    WRITE ( 14 )  y_shift
774    WRITE ( 14 )  'zeta_max                      '
775    WRITE ( 14 )  zeta_max
776    WRITE ( 14 )  'zeta_min                      '
777    WRITE ( 14 )  zeta_min
778    WRITE ( 14 )  'z0h_factor                    '
779    WRITE ( 14 )  z0h_factor
780
781!
782!-- Set the end-of-file mark for default parameter
783    WRITE ( 14 )  '*** end ***                   '
784
785!   
786!-- ADD NEW "CALL <module>_skip_var_list" BELOW THE LAST CALL
787!-- (The order of the calls MUST be the same as in parin, and
788!--  read/skip_var_list)!
789!
790!-- If required, write restart data for wind turbine model.
791    IF ( wind_turbine )  CALL wtm_write_restart_data
792
793!
794!-- If required, write restart data for virtual measurements.
795    IF ( virtual_flight )  CALL flight_write_restart_data
796
797!
798!-- If required, write restart data for turbulence generator.
799    IF ( syn_turb_gen )  CALL stg_write_restart_data
800
801
802
803 END SUBROUTINE write_var_list
Note: See TracBrowser for help on using the repository browser.