source: palm/trunk/SOURCE/parin.f90 @ 3049

Last change on this file since 3049 was 3049, checked in by Giersch, 3 years ago

Revision history corrected

  • Property svn:keywords set to Id
File size: 44.1 KB
Line 
1!> @file parin.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: parin.f90 3049 2018-05-29 13:52:36Z Giersch $
27! Error messages revised
28!
29! 3045 2018-05-28 07:55:41Z Giersch
30! z_max_do2d removed, error messages revised
31!
32! 2995 2018-04-19 12:13:16Z Giersch
33! time_since_reference_point must be calculated/initialized before the first 
34! call of functions related to the radiation model which occur in
35! time_integration_spinup or time_integration
36!
37! 2980 2018-04-17 15:19:27Z suehring
38! Revise message call
39!
40! 2975 2018-04-16 15:22:20Z suehring
41! - Informative message when initializing_actions has been changed
42!   to set_constant_profile in child domain
43! - Change location in message call
44!
45! 2967 2018-04-13 11:22:08Z raasch
46! bugfix: missing parallel cpp-directives added
47!
48! 2941 2018-04-03 11:54:58Z kanani
49! Fix for spinup in case of restart run
50!
51! 2938 2018-03-27 15:52:42Z suehring
52! Change initialization in case child domain should be initialized with Inifor.
53!
54! 2936 2018-03-27 14:49:27Z suehring
55! inipar renamed to initialization_parameters.
56! d3par renamed to runtime_parameters.
57!
58! 2921 2018-03-22 15:05:23Z Giersch
59! Activation of spinup has been moved from lsm/usm_parin to parin itself
60!
61! 2906 2018-03-19 08:56:40Z Giersch
62! ENVIRONMENT variables read/write_svf has been added
63!
64! 2894 2018-03-15 09:17:58Z Giersch
65! read_var_list has been renamed to rrd_global, all module related _parin
66! routines are called before reading the global restart data to overwrite them
67! in case of restart runs
68!
69! 2881 2018-03-13 16:24:40Z suehring
70! Added flag for switching on/off calculation of soil moisture
71!
72! 2849 2018-03-05 10:49:33Z Giersch
73! Position of d3par namelist in parameter file is unimportant now
74!
75! 2826 2018-02-21 12:39:28Z Giersch
76! Bugfix in setting the default boundary conditions for nest domains
77!
78! 2817 2018-02-19 16:32:21Z knoop
79! Preliminary gust module interface implemented
80!
81! 2773 2018-01-30 14:12:54Z suehring
82! Nesting for chemical species implemented
83!
84! 2766 2018-01-22 17:17:47Z kanani
85! Removed preprocessor directive __chem
86!
87! 2718 2018-01-02 08:49:38Z maronga
88! Corrected "Former revisions" section
89!
90! 2696 2017-12-14 17:12:51Z kanani
91! Change in file header (GPL part)
92! Implementation of uv exposure model (FK)
93! Added rans_mode and turbulence_closure to inipar (TG)
94! Implementation of chemistry module
95! Sorting of USE list (FK)
96! Forcing implemented, and initialization with inifor (MS)
97!
98! 2600 2017-11-01 14:11:20Z raasch
99! some comments added and variables renamed concerning r2599
100!
101! 2599 2017-11-01 13:18:45Z hellstea
102! The i/o grouping is updated to work correctly also in nested runs.
103!
104! 2575 2017-10-24 09:57:58Z maronga
105! Renamed phi -> latitude, added longitude
106!
107! 2563 2017-10-19 15:36:10Z Giersch
108! Changed position where restart files are closed.
109!
110! 2550 2017-10-16 17:12:01Z boeske
111! Added complex_terrain
112!
113! 2544 2017-10-13 18:09:32Z maronga
114! Moved day_of_year_init and time_utc_init to inipar.
115!
116! 2397 2017-09-04 16:22:48Z suehring
117! Enable initialization of 3d model by user in the child domain.
118!
119! 2375 2017-08-29 14:10:28Z schwenkel
120! Added aerosol initialization for bulk microphysics
121!
122! 2372 2017-08-25 12:37:32Z sward
123! y_shift added to namelist
124!
125! 2365 2017-08-21 14:59:59Z kanani
126! Vertical grid nesting: add vnest_start_time to d3par (SadiqHuq)
127!
128! 2339 2017-08-07 13:55:26Z gronemeier
129! corrected timestamp in header
130!
131! 2338 2017-08-07 12:15:38Z gronemeier
132! Modularize 1D model
133!
134! 2310 2017-07-11 09:37:02Z gronemeier
135! Bugfix: re-arranged call for error messages for ENVPAR file
136!
137! 2304 2017-07-04 14:35:55Z suehring
138! Bugfix, enable restarts for child domain.
139!
140! 2298 2017-06-29 09:28:18Z raasch
141! -return_addres, return_username in ENVPAR, -cross_ts_uymax, cross_ts_uymin in
142! d3par
143!
144! 2296 2017-06-28 07:53:56Z maronga
145! Added parameters for model spinup
146!
147! 2292 2017-06-20 09:51:42Z schwenkel
148! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
149! includes two more prognostic equations for cloud drop concentration (nc) 
150! and cloud water content (qc).
151!
152! 2267 2017-06-09 09:33:25Z gronemeier
153! Bugfix: removed skipping of reading namelists in case of omitted d3par
154!
155! 2259 2017-06-08 09:09:11Z gronemeier
156! Implemented synthetic turbulence generator
157!
158! 2233 2017-05-30 18:08:54Z suehring
159!
160! 2232 2017-05-30 17:47:52Z suehring
161! typo corrected
162! +wall_salinityflux
163! +tunnel_height, tunnel_lenght, tunnel_width_x, tunnel_width_y,
164!  tunnel_wall_depth
165!
166! 2118 2017-01-17 16:38:49Z raasch
167! -background_communication from inipar
168!
169! 2050 2016-11-08 15:00:55Z gronemeier
170! Implement turbulent outflow condition
171!
172! 2037 2016-10-26 11:15:40Z knoop
173! Anelastic approximation implemented
174!
175! 2035 2016-10-24 15:06:17Z suehring
176! Remove check for npex and npey in nesting case
177!
178! 2011 2016-09-19 17:29:57Z kanani
179! Added flag lsf_exception to allow explicit enabling of large scale forcing
180! together with buildings on flat terrain.
181!
182! 2007 2016-08-24 15:47:17Z kanani
183! Added call to urban surface model for reading of &urban_surface_par
184!
185! 2004 2016-08-24 10:25:59Z suehring
186! Humidity and passive scalar treated separately in nesting mode
187!
188! 2000 2016-08-20 18:09:15Z knoop
189! Forced header and separation lines into 80 columns
190!
191! 1992 2016-08-12 15:14:59Z suehring
192! +top_scalarflux
193!
194! 1960 2016-07-12 16:34:24Z suehring
195! Allocation of s_init
196!
197! 1957 2016-07-07 10:43:48Z suehring
198! flight module added
199!
200! 1955 2016-07-01 12:38:59Z hellstea
201! The parameter intializating_actions is set to 'set_constant_profiles for
202! all nest domains in order to make sure that diagnostic variables are properly
203! initialized for nest domains. Prognostic variables are later initialized by
204! interpolation from the parent domain.
205!
206! 1917 2016-05-27 14:28:12Z witha
207! Initial version of purely vertical nesting introduced.
208!
209! 1914 2016-05-26 14:44:07Z witha
210! Added call to wind turbine model for reading of &wind_turbine_par
211!
212! 1849 2016-04-08 11:33:18Z hoffmann
213! Adapted for modularization of microphysics
214!
215! 1833 2016-04-07 14:23:03Z raasch
216! call of spectra_parin
217!
218! 1831 2016-04-07 13:15:51Z hoffmann
219! turbulence renamed collision_turbulence, drizzle renamed
220! cloud_water_sedimentation
221! curvature_solution_effects removed
222!
223! 1826 2016-04-07 12:01:39Z maronga
224! Added call to radiation model for reading of &radiation_par.
225! Added call to plant canopy model for reading of &canopy_par.
226!
227! 1817 2016-04-06 15:44:20Z maronga
228! Added call to land surface model for reading of &lsm_par
229!
230! 1804 2016-04-05 16:30:18Z maronga
231! Removed code for parameter file check (__check)
232!
233! 1783 2016-03-06 18:36:17Z raasch
234! +netcdf_deflate in d3par, netcdf module and variable names changed
235!
236! 1764 2016-02-28 12:45:19Z raasch
237! cpp-statements for nesting removed, explicit settings of boundary conditions
238! in nest domains,
239! bugfix: npex/npey message moved from inipar to d3par
240! bugfix: check of lateral boundary conditions from check_parameters to here,
241! because they will be already used in init_pegrid and init_grid
242!
243! 1762 2016-02-25 12:31:13Z hellstea
244! Introduction of nested domain feature
245!
246! 1691 2015-10-26 16:17:44Z maronga
247! Added parameter most_method. Renamed prandtl_layer to constant_flux_layer.
248!
249! 1682 2015-10-07 23:56:08Z knoop
250! Code annotations made doxygen readable
251!
252! 1560 2015-03-06 10:48:54Z keck
253! +recycling_yshift
254!
255! 1496 2014-12-02 17:25:50Z maronga
256! Renamed: "radiation -> "cloud_top_radiation"
257!
258! 1484 2014-10-21 10:53:05Z kanani
259! Changes due to new module structure of the plant canopy model:
260!   canopy-model related parameters moved to new package canopy_par in
261!   subroutine package_parin
262!
263! 1429 2014-07-15 12:53:45Z knoop
264! +ensemble_member_nr to prepare the random_generator for ensemble runs
265!
266! 1402 2014-05-09 14:25:13Z raasch
267! location messages modified, batch_job included in envpar-NAMELIST
268!
269! 1384 2014-05-02 14:31:06Z raasch
270! location messages added
271!
272! 1365 2014-04-22 15:03:56Z boeske
273! Usage of large scale forcing enabled:
274! +use_subsidence_tendencies
275!
276! 1361 2014-04-16 15:17:48Z hoffmann
277! +call_microphysics_at_all_substeps
278!
279! 1359 2014-04-11 17:15:14Z hoffmann
280! REAL constants provided with KIND-attribute
281!
282! 1353 2014-04-08 15:21:23Z heinze
283! REAL constants provided with KIND-attribute
284!
285! 1327 2014-03-21 11:00:16Z raasch
286! -data_output_format, do3d_compress, do3d_comp_prec
287!
288! 1320 2014-03-20 08:40:49Z raasch
289! ONLY-attribute added to USE-statements,
290! kind-parameters added to all INTEGER and REAL declaration statements,
291! kinds are defined in new module kinds,
292! old module precision_kind is removed,
293! revision history before 2012 removed,
294! comment fields (!:) to be used for variable explanations added to
295! all variable declaration statements
296!
297! 1318 2014-03-17 13:35:16Z raasch
298! +cpu_log_barrierwait in d3par
299!
300! 1301 2014-03-06 13:29:46Z heinze
301! +large_scale_subsidence
302!
303! 1241 2013-10-30 11:36:58Z heinze
304! +nudging
305! +large_scale_forcing
306!
307! 1216 2013-08-26 09:31:42Z raasch
308! +transpose_compute_overlap in inipar
309!
310! 1195 2013-07-01 12:27:57Z heinze
311! Bugfix: allocate ref_state
312!
313! 1179 2013-06-14 05:57:58Z raasch
314! +reference_state in inipar
315!
316! 1159 2013-05-21 11:58:22Z fricke
317! +use_cmax
318!
319! 1128 2013-04-12 06:19:32Z raasch
320! +background_communication in inipar
321!
322! 1115 2013-03-26 18:16:16Z hoffmann
323! unused variables removed
324!
325! 1092 2013-02-02 11:24:22Z raasch
326! unused variables removed
327!
328! 1065 2012-11-22 17:42:36Z hoffmann
329! +nc, c_sedimentation, limiter_sedimentation, turbulence
330! -mu_constant, mu_constant_value
331!
332! 1053 2012-11-13 17:11:03Z hoffmann
333! necessary expansions according to the two new prognostic equations (nr, qr)
334! of the two-moment cloud physics scheme and steering parameters:
335! +*_init, *_surface, *_surface_initial_change, *_vertical_gradient,
336! +*_vertical_gradient_level, surface_waterflux_*,
337! +cloud_scheme, drizzle, mu_constant, mu_constant_value, ventilation_effect
338!
339! 1036 2012-10-22 13:43:42Z raasch
340! code put under GPL (PALM 3.9)
341!
342! 1015 2012-09-27 09:23:24Z raasch
343! -adjust_mixing_length
344!
345! 1003 2012-09-14 14:35:53Z raasch
346! -grid_matching
347!
348! 1001 2012-09-13 14:08:46Z raasch
349! -cut_spline_overshoot, long_filter_factor, overshoot_limit_*, ups_limit_*
350!
351! 996 2012-09-07 10:41:47Z raasch
352! -use_prior_plot1d_parameters
353!
354! 978 2012-08-09 08:28:32Z fricke
355! -km_damp_max, outflow_damping_width
356! +pt_damping_factor, pt_damping_width
357! +z0h_factor
358!
359! 964 2012-07-26 09:14:24Z raasch
360! -cross_normalized_x, cross_normalized_y, cross_xtext, z_max_do1d,
361! z_max_do1d_normalized
362!
363! 940 2012-07-09 14:31:00Z raasch
364! +neutral in inipar
365!
366! 927 2012-06-06 19:15:04Z raasch
367! +masking_method in inipar
368!
369! 824 2012-02-17 09:09:57Z raasch
370! +curvature_solution_effects in inipar
371!
372! 809 2012-01-30 13:32:58Z maronga
373! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
374!
375! 807 2012-01-25 11:53:51Z maronga
376! New cpp directive "__check" implemented which is used by check_namelist_files
377!
378! Revision 1.1  1997/07/24 11:22:50  raasch
379! Initial revision
380!
381!
382! Description:
383! ------------
384!> This subroutine reads variables controling the run from the NAMELIST files
385!------------------------------------------------------------------------------!
386 SUBROUTINE parin
387 
388
389    USE arrays_3d,                                                             &
390        ONLY:  pt_init, q_init, ref_state, s_init, sa_init,                    &     
391               ug, u_init, v_init, vg
392
393    USE chemistry_model_mod,                                                   &
394        ONLY:  chem_parin
395       
396    USE chem_modules
397
398    USE control_parameters
399
400    USE cpulog,                                                                &
401        ONLY:  cpu_log_barrierwait
402
403    USE date_and_time_mod,                                                     &
404        ONLY:  day_of_year_init, time_utc_init
405
406    USE dvrp_variables,                                                        &
407        ONLY:  local_dvrserver_running
408
409    USE flight_mod,                                                            &
410        ONLY:  flight_parin
411
412    USE grid_variables,                                                        &
413        ONLY:  dx, dy
414
415    USE gust_mod,                                                              &
416        ONLY: gust_parin
417
418    USE indices,                                                               &
419        ONLY:  nx, ny, nz
420
421    USE kinds
422
423    USE land_surface_model_mod,                                                &
424        ONLY: lsm_parin
425
426    USE microphysics_mod,                                                      &
427        ONLY:  c_sedimentation, cloud_water_sedimentation,                     &
428               collision_turbulence, curvature_solution_effects_bulk,          &
429               dry_aerosol_radius, limiter_sedimentation, na_init, nc_const,   &
430               sigma_bulk, ventilation_effect
431
432    USE model_1d_mod,                                                          &
433        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
434
435    USE netcdf_interface,                                                      &
436        ONLY:  netcdf_data_format, netcdf_deflate, netcdf_precision
437
438    USE pegrid
439               
440    USE plant_canopy_model_mod,                                                &
441         ONLY: pcm_parin
442
443    USE pmc_interface,                                                         &
444        ONLY:  nested_run, nesting_mode
445
446    USE profil_parameter,                                                      &
447        ONLY:  cross_profiles, profile_columns, profile_rows
448
449    USE progress_bar,                                                          &
450        ONLY :  batch_job
451
452    USE radiation_model_mod,                                                   &
453        ONLY: radiation_parin
454
455    USE read_restart_data_mod,                                                 &
456        ONLY:  rrd_global     
457
458    USE spectra_mod,                                                           &
459        ONLY :  spectra_parin
460
461    USE statistics,                                                            &
462        ONLY:  hom, hom_sum, pr_palm, region, statistic_regions
463
464    USE synthetic_turbulence_generator_mod,                                    &
465        ONLY:  stg_parin
466
467    USE urban_surface_mod,                                                     &
468        ONLY: usm_parin
469
470    USE uv_exposure_model_mod,                                                 &
471        ONLY:  uvem_parin
472
473    USE vertical_nesting_mod,                                                  &
474        ONLY:  vnest_start_time
475
476    USE wind_turbine_model_mod,                                                &
477        ONLY:  wtm_parin
478
479
480    IMPLICIT NONE
481
482    CHARACTER (LEN=80) ::  line  !< dummy string that contains the current line of the parameter file
483
484    INTEGER(iwp) ::  global_id      !< process id with respect to MPI_COMM_WORLD
485    INTEGER(iwp) ::  global_procs   !< # of procs with respect to MPI_COMM_WORLD
486    INTEGER(iwp) ::  i              !<
487    INTEGER(iwp) ::  ioerr          !< error flag for open/read/write
488
489    NAMELIST /inipar/  aerosol_bulk, alpha_surface, approximation, bc_e_b,     &
490                       bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b, &
491             bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t,                 &
492             bottom_salinityflux, building_height, building_length_x,          &
493             building_length_y, building_wall_left, building_wall_south,       &
494             calc_soil_moisture_during_spinup,                                 &
495             call_psolver_at_all_substeps, call_microphysics_at_all_substeps,  &
496             canyon_height,                                                    &
497             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
498             canyon_wall_south, c_sedimentation, cfl_factor, cloud_droplets,   &
499             cloud_physics, cloud_scheme, cloud_top_radiation,                 &
500             cloud_water_sedimentation,                                        &
501             collective_wait, collision_turbulence, complex_terrain,           &
502             conserve_volume_flow,                                             &
503             conserve_volume_flow_mode, constant_flux_layer,                   &
504             coupling_start_time, curvature_solution_effects_bulk,             &
505             cycle_mg, damp_level_1d,                                          &
506             data_output_during_spinup,                                        &
507             day_of_year_init,                                                 &
508             dissipation_1d,                                                   &
509             dp_external, dp_level_b, dp_smooth, dpdxy, dry_aerosol_radius,    &
510             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
511             dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
512             ensemble_member_nr, e_init, e_min, fft_method,                    &
513             flux_input_mode, flux_output_mode, forcing,                       &
514             galilei_transformation, humidity,                                 &
515             inflow_damping_height, inflow_damping_width,                      &
516             inflow_disturbance_begin, inflow_disturbance_end,                 &
517             initializing_actions, km_constant,                                &
518             large_scale_forcing, large_scale_subsidence, latitude,            &
519             limiter_sedimentation, longitude,                                 &
520             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
521             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
522             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
523             nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
524             outflow_source_plane, passive_scalar,                             &
525             prandtl_number, precipitation, psolver, pt_damping_factor,        &
526             pt_damping_width, pt_reference, pt_surface,                       &
527             pt_surface_initial_change, pt_vertical_gradient,                  &
528             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
529             q_vertical_gradient, q_vertical_gradient_level,                   &
530             random_generator, random_heatflux, rans_mode,                     &
531             rayleigh_damping_factor, rayleigh_damping_height,                 &
532             recycling_width, recycling_yshift,                                &
533             reference_state, residual_limit,                                  &
534             roughness_length, sa_surface,                                     &
535             sa_vertical_gradient, sa_vertical_gradient_level, scalar_advec,   &
536             scalar_rayleigh_damping, sigma_bulk,                              &
537             spinup_time, spinup_pt_amplitude, spinup_pt_mean,                 &
538             statistic_regions, subs_vertical_gradient,                        &
539             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
540             surface_scalarflux, surface_waterflux,                            &
541             s_surface, s_surface_initial_change, s_vertical_gradient,         &
542             s_vertical_gradient_level, time_utc_init, timestep_scheme,        &
543             topography, topography_grid_convention, top_heatflux,             &
544             top_momentumflux_u, top_momentumflux_v, top_salinityflux,         &
545             top_scalarflux, transpose_compute_overlap,                        &
546             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
547             tunnel_wall_depth, turbulence_closure,                            &
548             turbulent_inflow, turbulent_outflow,                              &
549             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
550             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
551             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
552             uv_heights, u_bulk, u_profile, vg_surface, vg_vertical_gradient,  &
553             vg_vertical_gradient_level, v_bulk, v_profile, ventilation_effect,&
554             wall_adjustment, wall_heatflux, wall_humidityflux,                &
555             wall_salinityflux, wall_scalarflux, y_shift, zeta_max, zeta_min,  &
556             z0h_factor
557
558    NAMELIST /initialization_parameters/  aerosol_bulk, alpha_surface,         &
559             approximation, bc_e_b,                                            &
560             bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b,           &
561             bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t,                 &
562             bottom_salinityflux, building_height, building_length_x,          &
563             building_length_y, building_wall_left, building_wall_south,       &
564             calc_soil_moisture_during_spinup,                                 &
565             call_psolver_at_all_substeps, call_microphysics_at_all_substeps,  &
566             canyon_height,                                                    &
567             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
568             canyon_wall_south, c_sedimentation, cfl_factor, cloud_droplets,   &
569             cloud_physics, cloud_scheme, cloud_top_radiation,                 &
570             cloud_water_sedimentation,                                        &
571             collective_wait, collision_turbulence, complex_terrain,           &
572             conserve_volume_flow,                                             &
573             conserve_volume_flow_mode, constant_flux_layer,                   &
574             coupling_start_time, curvature_solution_effects_bulk,             &
575             cycle_mg, damp_level_1d,                                          &
576             data_output_during_spinup,                                        &
577             day_of_year_init,                                                 &
578             dissipation_1d,                                                   &
579             dp_external, dp_level_b, dp_smooth, dpdxy, dry_aerosol_radius,    &
580             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
581             dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
582             ensemble_member_nr, e_init, e_min, fft_method,                    &
583             flux_input_mode, flux_output_mode, forcing,                       &
584             galilei_transformation, humidity,                                 &
585             inflow_damping_height, inflow_damping_width,                      &
586             inflow_disturbance_begin, inflow_disturbance_end,                 &
587             initializing_actions, km_constant,                                &
588             large_scale_forcing, large_scale_subsidence, latitude,            &
589             limiter_sedimentation, longitude,                                 &
590             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
591             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
592             most_method, na_init, nc_const, netcdf_precision, neutral, ngsrb, &
593             nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
594             outflow_source_plane, passive_scalar,                             &
595             prandtl_number, precipitation, psolver, pt_damping_factor,        &
596             pt_damping_width, pt_reference, pt_surface,                       &
597             pt_surface_initial_change, pt_vertical_gradient,                  &
598             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
599             q_vertical_gradient, q_vertical_gradient_level,                   &
600             random_generator, random_heatflux, rans_mode,                     &
601             rayleigh_damping_factor, rayleigh_damping_height,                 &
602             recycling_width, recycling_yshift,                                &
603             reference_state, residual_limit,                                  &
604             roughness_length, sa_surface,                                     &
605             sa_vertical_gradient, sa_vertical_gradient_level, scalar_advec,   &
606             scalar_rayleigh_damping, sigma_bulk,                              &
607             spinup_time, spinup_pt_amplitude, spinup_pt_mean,                 &
608             statistic_regions, subs_vertical_gradient,                        &
609             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
610             surface_scalarflux, surface_waterflux,                            &
611             s_surface, s_surface_initial_change, s_vertical_gradient,         &
612             s_vertical_gradient_level, time_utc_init, timestep_scheme,        &
613             topography, topography_grid_convention, top_heatflux,             &
614             top_momentumflux_u, top_momentumflux_v, top_salinityflux,         &
615             top_scalarflux, transpose_compute_overlap,                        &
616             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
617             tunnel_wall_depth, turbulence_closure,                            &
618             turbulent_inflow, turbulent_outflow,                              &
619             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
620             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
621             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
622             uv_heights, u_bulk, u_profile, vg_surface, vg_vertical_gradient,  &
623             vg_vertical_gradient_level, v_bulk, v_profile, ventilation_effect,&
624             wall_adjustment, wall_heatflux, wall_humidityflux,                &
625             wall_salinityflux, wall_scalarflux, y_shift, zeta_max, zeta_min,  &
626             z0h_factor
627             
628    NAMELIST /d3par/  averaging_interval, averaging_interval_pr,               &
629             cpu_log_barrierwait, create_disturbances,                         &
630             cross_profiles, data_output, data_output_masks,                   &
631             data_output_pr, data_output_2d_on_each_pe, disturbance_amplitude, &
632             disturbance_energy_limit, disturbance_level_b,                    &
633             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
634             dt, dt_averaging_input, dt_averaging_input_pr,                    &
635             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
636             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
637             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
638             dt_run_control,end_time, force_print_header, mask_scale_x,        &
639             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
640             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
641             normalizing_region, npex, npey, nz_do3d,                          &
642             precipitation_amount_interval, profile_columns, profile_rows,     &
643             restart_time, section_xy, section_xz, section_yz,                 &
644             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
645             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
646             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
647             termination_time_needed, vnest_start_time
648
649    NAMELIST /runtime_parameters/  averaging_interval, averaging_interval_pr,  &
650             cpu_log_barrierwait, create_disturbances,                         &
651             cross_profiles, data_output, data_output_masks,                   &
652             data_output_pr, data_output_2d_on_each_pe, disturbance_amplitude, &
653             disturbance_energy_limit, disturbance_level_b,                    &
654             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
655             dt, dt_averaging_input, dt_averaging_input_pr,                    &
656             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
657             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
658             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
659             dt_run_control,end_time, force_print_header, mask_scale_x,        &
660             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
661             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
662             normalizing_region, npex, npey, nz_do3d,                          &
663             precipitation_amount_interval, profile_columns, profile_rows,     &
664             restart_time, section_xy, section_xz, section_yz,                 &
665             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
666             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
667             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
668             termination_time_needed, vnest_start_time
669
670    NAMELIST /envpar/  batch_job, host, local_dvrserver_running,               &
671                       maximum_cpu_time_allowed, maximum_parallel_io_streams,  &
672                       read_svf, revision, run_identifier, tasks_per_node,     &
673                       write_binary, write_svf
674
675!
676!-- First read values of environment variables (this NAMELIST file is
677!-- generated by palmrun)
678    CALL location_message( 'reading environment parameters from ENVPAR', .FALSE. )
679
680    OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', IOSTAT=ioerr )
681
682    IF ( ioerr /= 0 )  THEN
683       message_string = 'local file ENVPAR not found' //                       &
684                        '&some variables for steering may not be properly set'
685       CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
686    ELSE
687       READ ( 90, envpar, IOSTAT=ioerr )
688       IF ( ioerr < 0 )  THEN
689          message_string = 'no envpar-NAMELIST found in local file '  //       &
690                           'ENVPAR& or some variables for steering may '  //   &
691                           'not be properly set'
692          CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
693       ELSEIF ( ioerr > 0 )  THEN
694          message_string = 'errors in local file ENVPAR' //                    &
695                           '&some variables for steering may not be properly set'
696          CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
697       ENDIF
698       CLOSE ( 90 )
699    ENDIF
700
701    CALL location_message( 'finished', .TRUE. )
702!
703!-- Calculate the number of groups into which parallel I/O is split.
704!-- The default for files which are opened by all PEs (or where each
705!-- PE opens his own independent file) is, that all PEs are doing input/output
706!-- in parallel at the same time. This might cause performance or even more
707!-- severe problems depending on the configuration of the underlying file
708!-- system.
709!-- Calculation of the number of blocks and the I/O group must be based on all
710!-- PEs involved in this run. Since myid and numprocs are related to the
711!-- comm2d communicator, which gives only a subset of all PEs in case of
712!-- nested runs, that information must be inquired again from the global
713!-- communicator.
714!-- First, set the default:
715#if defined( __parallel )
716    CALL MPI_COMM_RANK( MPI_COMM_WORLD, global_id, ierr )
717    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, global_procs, ierr )
718#else
719    global_id    = 0
720    global_procs = 1
721#endif
722    IF ( maximum_parallel_io_streams == -1  .OR.                               &
723         maximum_parallel_io_streams > global_procs )  THEN
724       maximum_parallel_io_streams = global_procs
725    ENDIF
726!
727!-- Now calculate the number of io_blocks and the io_group to which the
728!-- respective PE belongs. I/O of the groups is done in serial, but in parallel
729!-- for all PEs belonging to the same group.
730    io_blocks = global_procs / maximum_parallel_io_streams
731    io_group  = MOD( global_id+1, io_blocks )
732   
733    CALL location_message( 'reading NAMELIST parameters from PARIN', .FALSE. )
734!
735!-- Data is read in parallel by groups of PEs
736    DO  i = 0, io_blocks-1
737       IF ( i == io_group )  THEN
738
739!
740!--       Open the NAMELIST-file which is send with this job
741          CALL check_open( 11 )
742
743!
744!--       Read the control parameters for initialization.
745!--       The namelist "inipar" must be provided in the NAMELIST-file.
746          READ ( 11, initialization_parameters, ERR=10, END=11 )
747
748          GOTO 12
749
750 10       message_string = 'errors in initialization_parameters & or no ' //  &
751                           'initialization_parameters-namelist ' //           &
752                           'found (CRAY-machines only)'
753          CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
754
755 11       REWIND ( 11 )
756          READ ( 11, inipar, ERR=13, END=14 )
757 
758          message_string = 'namelist inipar is deprecated and will be ' //    &
759                          'removed in near future. & Please use namelist ' // &
760                          'initialization_parameters instead'
761          CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 )
762 
763          GOTO 12
764 
765 13       message_string = 'errors in inipar & or no inipar-namelist ' //      &
766                           'found (CRAY-machines only)'
767          CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
768         
769 14       message_string = 'no initialization_parameters-namelist found'
770          CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
771
772!
773!--       Try to read runtime parameters given by the user for this run
774!--       (namelist "runtime_parameters"). The namelist "runtime_parmeters"   
775!--       can be omitted. In that case default values are used for the         
776!--       parameters.
777 12       line = ' '
778
779          REWIND ( 11 )
780          line = ' '
781          DO   WHILE ( INDEX( line, '&runtime_parameters' ) == 0 )
782             READ ( 11, '(A)', END=20 )  line
783          ENDDO
784          BACKSPACE ( 11 )
785
786!
787!--       Read namelist
788          READ ( 11, runtime_parameters )
789
790          GOTO 21
791         
792 20       REWIND ( 11 )
793          line = ' '
794          DO   WHILE ( INDEX( line, '&d3par' ) == 0 )
795             READ ( 11, '(A)', END=21 )  line
796          ENDDO
797          BACKSPACE ( 11 )
798 
799 !
800!--       Read namelist
801          READ ( 11, d3par )
802 
803          message_string = 'namelist d3par is deprecated and will be ' //      &
804                          'removed in near future. &Please use namelist ' //   &
805                          'runtime_parameters instead'
806          CALL message( 'parin', 'PA0487', 0, 1, 0, 6, 0 )
807         
808 21       CONTINUE
809
810!
811!--       Check if land surface model is used and read &lsm_par if required
812          CALL lsm_parin
813
814!
815!--       Check if urban surface model is used and read &urban_surface_par if required
816          CALL usm_parin
817
818!
819!--       Check if spectra shall be calculated and read spectra_par if required
820          CALL spectra_parin
821
822!
823!--       Check if radiation model is used and read &radiation_par if required
824          CALL radiation_parin
825
826!
827!--       Check if gust module is used and read &gust_par if required
828          CALL gust_parin
829 
830 
831!--       Check if plant canopy model is used and read &canopy_par if required
832          CALL pcm_parin
833 
834!
835!--       Read control parameters for optionally used model software packages
836          CALL package_parin
837
838!
839!--       Check if wind turbine model is used and read &wind_turbine_par if
840!--       required
841          CALL wtm_parin
842!
843!--       Check if virtual flights should be carried out and read &flight_par
844!--       if required
845          CALL flight_parin
846!
847!--       Check if synthetic turbulence generator is used and read stg_par if
848!--       required
849          CALL stg_parin
850!
851!--       Read chemistry variables
852          CALL chem_parin
853!
854!--       Check if uv exposure model is used and read &uvexposure_par
855          CALL uvem_parin
856!
857!--       Read user-defined variables
858          CALL user_parin
859
860!
861!--       If required, read control parameters from restart file (produced by
862!--       a prior run). All PEs are reading from file created by PE0 (see
863!--       check_open)
864          IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
865
866             CALL rrd_global
867!
868!--          Increment the run count
869             runnr = runnr + 1
870          ENDIF
871
872!
873!--       Activate spinup
874          IF ( land_surface .OR. urban_surface )  THEN
875             IF ( spinup_time > 0.0_wp )  THEN
876                coupling_start_time = spinup_time
877                time_since_reference_point = simulated_time - coupling_start_time
878                IF ( spinup_pt_mean == 9999999.9_wp )  THEN
879                   spinup_pt_mean = pt_surface
880                ENDIF
881                end_time = end_time + spinup_time
882                IF ( TRIM( initializing_actions ) /= 'read_restart_data' )     &
883                   spinup = .TRUE.
884             ENDIF
885          ENDIF
886
887!
888!--       In case of nested runs, explicitly set nesting boundary conditions.
889!--       This will overwrite the user settings and basic defaults.
890!--       bc_lr and bc_ns always need to be cyclic for vertical nesting.
891          IF ( nested_run )  THEN
892             IF ( nesting_mode == 'vertical' )  THEN
893                IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
894                   WRITE ( message_string, *) 'bc_lr and bc_ns were set to ,', &
895                        'cyclic for vertical nesting'
896                   CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
897                   bc_lr   = 'cyclic'
898                   bc_ns   = 'cyclic'
899                ENDIF
900                IF ( nest_domain )  THEN
901                   bc_uv_t  = 'nested'
902                   bc_pt_t  = 'nested'
903                   bc_q_t   = 'nested'
904                   bc_s_t   = 'nested'
905                   bc_cs_t  = 'nested'
906                   bc_p_t   = 'neumann' 
907                ENDIF
908!
909!--          For other nesting modes only set boundary conditions for
910!--          nested domains.
911             ELSE
912                IF ( nest_domain )  THEN
913                   bc_lr    = 'nested'
914                   bc_ns    = 'nested'
915                   bc_uv_t  = 'nested'
916                   bc_pt_t  = 'nested'
917                   bc_q_t   = 'nested'
918                   bc_s_t   = 'nested'
919                   bc_cs_t  = 'nested'
920                   bc_p_t   = 'neumann'
921                ENDIF
922             ENDIF
923          ENDIF
924
925          IF ( forcing )  THEN
926             bc_lr    = 'forcing'
927             bc_ns    = 'forcing'
928             bc_uv_t  = 'forcing'
929             bc_pt_t  = 'forcing'
930             bc_q_t   = 'forcing'
931             bc_s_t   = 'forcing'  ! scalar boundary condition is not clear
932             bc_cs_t  = 'forcing'  ! same for chemical species
933             bc_p_t   = 'neumann'
934          ENDIF
935
936!         
937!--       In case of nested runs, make sure that initializing_actions =
938!--       'set_constant_profiles' even though the constant-profiles
939!--       initializations for the prognostic variables will be overwritten
940!--       by pmci_child_initialize and pmci_parent_initialize. This is,
941!--       however, important e.g. to make sure that diagnostic variables
942!--       are set properly. An exception is made in case of restart runs and
943!--       if user decides to do everything by its own.
944          IF ( nest_domain  .AND.  .NOT. (                                     &
945               TRIM( initializing_actions ) == 'read_restart_data'      .OR.   &
946               TRIM( initializing_actions ) == 'set_constant_profiles'  .OR.   &
947               TRIM( initializing_actions ) == 'by_user' ) )  THEN
948             message_string = 'initializing_actions = ' //                     &
949                              TRIM( initializing_actions ) // ' has been ' //  &
950                              'changed to set_constant_profiles in child ' //  &
951                              'domain.' 
952             CALL message( 'parin', 'PA0492', 0, 0, 0, 6, 0 )
953
954             initializing_actions = 'set_constant_profiles'
955          ENDIF
956           
957!
958!--       Check validity of lateral boundary conditions. This has to be done
959!--       here because they are already used in init_pegrid and init_grid and
960!--       therefore cannot be check in check_parameters
961          IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
962               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested'  .AND. &
963               bc_lr /= 'forcing' )  THEN
964             message_string = 'unknown boundary condition: bc_lr = "' // &
965                              TRIM( bc_lr ) // '"'
966             CALL message( 'parin', 'PA0049', 1, 2, 0, 6, 0 )
967          ENDIF
968          IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
969               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested'  .AND. &
970               bc_ns /= 'forcing' )  THEN
971             message_string = 'unknown boundary condition: bc_ns = "' // &
972                              TRIM( bc_ns ) // '"'
973             CALL message( 'parin', 'PA0050', 1, 2, 0, 6, 0 )
974          ENDIF
975!
976!--       Set internal variables used for speed optimization in if clauses
977          IF ( bc_lr /= 'cyclic' )               bc_lr_cyc    = .FALSE.
978          IF ( bc_lr == 'dirichlet/radiation' )  bc_lr_dirrad = .TRUE.
979          IF ( bc_lr == 'radiation/dirichlet' )  bc_lr_raddir = .TRUE.
980          IF ( bc_ns /= 'cyclic' )               bc_ns_cyc    = .FALSE.
981          IF ( bc_ns == 'dirichlet/radiation' )  bc_ns_dirrad = .TRUE.
982          IF ( bc_ns == 'radiation/dirichlet' )  bc_ns_raddir = .TRUE.
983
984!
985!--       Definition of names of areas used for computing statistics. They must
986!--       be defined at this place, because they are allowed to be redefined by
987!--       the user in user_parin.
988          region = 'total domain'
989
990!
991!--       Check in case of initial run, if the grid point numbers are well
992!--       defined and allocate some arrays which are already needed in
993!--       init_pegrid or check_parameters. During restart jobs, these arrays
994!--       will be allocated in rrd_global. All other arrays are allocated
995!--       in init_3d_model.
996          IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
997
998             IF ( nx <= 0 )  THEN
999                WRITE( message_string, * ) 'no value or wrong value given',    &
1000                                           ' for nx: nx=', nx
1001                CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 )
1002             ENDIF
1003             IF ( ny <= 0 )  THEN
1004                WRITE( message_string, * ) 'no value or wrong value given',    &
1005                                           ' for ny: ny=', ny
1006                CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 )
1007             ENDIF
1008             IF ( nz <= 0 )  THEN
1009                WRITE( message_string, * ) 'no value or wrong value given',    &
1010                                           ' for nz: nz=', nz
1011                CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 )
1012             ENDIF
1013!
1014!--          ATTENTION: in case of changes to the following statement please
1015!--                  also check the allocate statement in routine rrd_global
1016             ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1),        &
1017                       ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1),         &
1018                       u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1),             &
1019                       hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions),  &
1020                       hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
1021
1022             hom = 0.0_wp
1023
1024          ENDIF
1025
1026!
1027!--       NAMELIST-file is not needed anymore
1028          CALL close_file( 11 )
1029
1030       ENDIF
1031#if defined( __parallel )
1032       CALL MPI_BARRIER( MPI_COMM_WORLD, ierr )
1033#endif
1034    ENDDO
1035
1036    CALL location_message( 'finished', .TRUE. )
1037
1038 END SUBROUTINE parin
Note: See TracBrowser for help on using the repository browser.