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

Last change on this file since 2050 was 2050, checked in by gronemeier, 7 years ago

Implement turbulent outflow condition

  • Property svn:keywords set to Id
File size: 26.7 KB
Line 
1!> @file parin.f90
2!------------------------------------------------------------------------------!
3! This file is part of PALM.
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-2016 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22! Implement turbulent outflow condition
23!
24! Former revisions:
25! -----------------
26! $Id: parin.f90 2050 2016-11-08 15:00:55Z gronemeier $
27!
28! 2037 2016-10-26 11:15:40Z knoop
29! Anelastic approximation implemented
30!
31! 2035 2016-10-24 15:06:17Z suehring
32! Remove check for npex and npey in nesting case
33!
34! 2011 2016-09-19 17:29:57Z kanani
35! Added flag lsf_exception to allow explicit enabling of large scale forcing
36! together with buildings on flat terrain.
37!
38! 2007 2016-08-24 15:47:17Z kanani
39! Added call to urban surface model for reading of &urban_surface_par
40!
41! 2004 2016-08-24 10:25:59Z suehring
42! Humidity and passive scalar treated separately in nesting mode
43!
44! 2000 2016-08-20 18:09:15Z knoop
45! Forced header and separation lines into 80 columns
46!
47! 1992 2016-08-12 15:14:59Z suehring
48! +top_scalarflux
49!
50! 1960 2016-07-12 16:34:24Z suehring
51! Allocation of s_init
52!
53! 1957 2016-07-07 10:43:48Z suehring
54! flight module added
55!
56! 1955 2016-07-01 12:38:59Z hellstea
57! The parameter intializating_actions is set to 'set_constant_profiles for
58! all nest domains in order to make sure that diagnostic variables are properly
59! initialized for nest domains. Prognostic variables are later initialized by
60! interpolation from the parent domain.
61!
62! 1917 2016-05-27 14:28:12Z witha
63! Initial version of purely vertical nesting introduced.
64!
65! 1914 2016-05-26 14:44:07Z witha
66! Added call to wind turbine model for reading of &wind_turbine_par
67!
68! 1849 2016-04-08 11:33:18Z hoffmann
69! Adapted for modularization of microphysics
70!
71! 1833 2016-04-07 14:23:03Z raasch
72! call of spectra_parin
73!
74! 1831 2016-04-07 13:15:51Z hoffmann
75! turbulence renamed collision_turbulence, drizzle renamed
76! cloud_water_sedimentation
77! curvature_solution_effects removed
78!
79! 1826 2016-04-07 12:01:39Z maronga
80! Added call to radiation model for reading of &radiation_par.
81! Added call to plant canopy model for reading of &canopy_par.
82!
83! 1817 2016-04-06 15:44:20Z maronga
84! Added call to land surface model for reading of &lsm_par
85!
86! 1804 2016-04-05 16:30:18Z maronga
87! Removed code for parameter file check (__check)
88!
89! 1783 2016-03-06 18:36:17Z raasch
90! +netcdf_deflate in d3par, netcdf module and variable names changed
91!
92! 1764 2016-02-28 12:45:19Z raasch
93! cpp-statements for nesting removed, explicit settings of boundary conditions
94! in nest domains,
95! bugfix: npex/npey message moved from inipar to d3par
96! bugfix: check of lateral boundary conditions from check_parameters to here,
97! because they will be already used in init_pegrid and init_grid
98!
99! 1762 2016-02-25 12:31:13Z hellstea
100! Introduction of nested domain feature
101!
102! 1691 2015-10-26 16:17:44Z maronga
103! Added parameter most_method. Renamed prandtl_layer to constant_flux_layer.
104!
105! 1682 2015-10-07 23:56:08Z knoop
106! Code annotations made doxygen readable
107!
108! 1560 2015-03-06 10:48:54Z keck
109! +recycling_yshift
110!
111! 1496 2014-12-02 17:25:50Z maronga
112! Renamed: "radiation -> "cloud_top_radiation"
113!
114! 1484 2014-10-21 10:53:05Z kanani
115! Changes due to new module structure of the plant canopy model:
116!   canopy-model related parameters moved to new package canopy_par in
117!   subroutine package_parin
118!
119! 1429 2014-07-15 12:53:45Z knoop
120! +ensemble_member_nr to prepare the random_generator for ensemble runs
121!
122! 1402 2014-05-09 14:25:13Z raasch
123! location messages modified, batch_job included in envpar-NAMELIST
124!
125! 1384 2014-05-02 14:31:06Z raasch
126! location messages added
127!
128! 1365 2014-04-22 15:03:56Z boeske
129! Usage of large scale forcing enabled:
130! +use_subsidence_tendencies
131!
132! 1361 2014-04-16 15:17:48Z hoffmann
133! +call_microphysics_at_all_substeps
134!
135! 1359 2014-04-11 17:15:14Z hoffmann
136! REAL constants provided with KIND-attribute
137!
138! 1353 2014-04-08 15:21:23Z heinze
139! REAL constants provided with KIND-attribute
140!
141! 1327 2014-03-21 11:00:16Z raasch
142! -data_output_format, do3d_compress, do3d_comp_prec
143!
144! 1320 2014-03-20 08:40:49Z raasch
145! ONLY-attribute added to USE-statements,
146! kind-parameters added to all INTEGER and REAL declaration statements,
147! kinds are defined in new module kinds,
148! old module precision_kind is removed,
149! revision history before 2012 removed,
150! comment fields (!:) to be used for variable explanations added to
151! all variable declaration statements
152!
153! 1318 2014-03-17 13:35:16Z raasch
154! +cpu_log_barrierwait in d3par
155!
156! 1301 2014-03-06 13:29:46Z heinze
157! +large_scale_subsidence
158!
159! 1241 2013-10-30 11:36:58Z heinze
160! +nudging
161! +large_scale_forcing
162!
163! 1216 2013-08-26 09:31:42Z raasch
164! +transpose_compute_overlap in inipar
165!
166! 1195 2013-07-01 12:27:57Z heinze
167! Bugfix: allocate ref_state
168!
169! 1179 2013-06-14 05:57:58Z raasch
170! +reference_state in inipar
171!
172! 1159 2013-05-21 11:58:22Z fricke
173! +use_cmax
174!
175! 1128 2013-04-12 06:19:32Z raasch
176! +background_communication in inipar
177!
178! 1115 2013-03-26 18:16:16Z hoffmann
179! unused variables removed
180!
181! 1092 2013-02-02 11:24:22Z raasch
182! unused variables removed
183!
184! 1065 2012-11-22 17:42:36Z hoffmann
185! +nc, c_sedimentation, limiter_sedimentation, turbulence
186! -mu_constant, mu_constant_value
187!
188! 1053 2012-11-13 17:11:03Z hoffmann
189! necessary expansions according to the two new prognostic equations (nr, qr)
190! of the two-moment cloud physics scheme and steering parameters:
191! +*_init, *_surface, *_surface_initial_change, *_vertical_gradient,
192! +*_vertical_gradient_level, surface_waterflux_*,
193! +cloud_scheme, drizzle, mu_constant, mu_constant_value, ventilation_effect
194!
195! 1036 2012-10-22 13:43:42Z raasch
196! code put under GPL (PALM 3.9)
197!
198! 1015 2012-09-27 09:23:24Z raasch
199! -adjust_mixing_length
200!
201! 1003 2012-09-14 14:35:53Z raasch
202! -grid_matching
203!
204! 1001 2012-09-13 14:08:46Z raasch
205! -cut_spline_overshoot, long_filter_factor, overshoot_limit_*, ups_limit_*
206!
207! 996 2012-09-07 10:41:47Z raasch
208! -use_prior_plot1d_parameters
209!
210! 978 2012-08-09 08:28:32Z fricke
211! -km_damp_max, outflow_damping_width
212! +pt_damping_factor, pt_damping_width
213! +z0h_factor
214!
215! 964 2012-07-26 09:14:24Z raasch
216! -cross_normalized_x, cross_normalized_y, cross_xtext, z_max_do1d,
217! z_max_do1d_normalized
218!
219! 940 2012-07-09 14:31:00Z raasch
220! +neutral in inipar
221!
222! 927 2012-06-06 19:15:04Z raasch
223! +masking_method in inipar
224!
225! 824 2012-02-17 09:09:57Z raasch
226! +curvature_solution_effects in inipar
227!
228! 809 2012-01-30 13:32:58Z maronga
229! Bugfix: replaced .AND. and .NOT. with && and ! in the preprocessor directives
230!
231! 807 2012-01-25 11:53:51Z maronga
232! New cpp directive "__check" implemented which is used by check_namelist_files
233!
234! Revision 1.1  1997/07/24 11:22:50  raasch
235! Initial revision
236!
237!
238! Description:
239! ------------
240!> This subroutine reads variables controling the run from the NAMELIST files
241!------------------------------------------------------------------------------!
242 SUBROUTINE parin
243 
244
245    USE arrays_3d,                                                             &
246        ONLY:  pt_init, q_init, ref_state, s_init, sa_init, ug, u_init, v_init,&
247               vg
248
249    USE plant_canopy_model_mod,                                                &
250         ONLY: pcm_parin 
251
252    USE control_parameters
253
254    USE cpulog,                                                                &
255        ONLY:  cpu_log_barrierwait
256
257    USE dvrp_variables,                                                        &
258        ONLY:  local_dvrserver_running
259
260    USE flight_mod,                                                            &
261        ONLY:  flight_parin
262
263    USE grid_variables,                                                        &
264        ONLY:  dx, dy
265
266    USE indices,                                                               &
267        ONLY:  nx, ny, nz
268
269    USE kinds
270
271    USE land_surface_model_mod,                                                &
272        ONLY: lsm_parin
273
274    USE microphysics_mod,                                                      &
275        ONLY:  c_sedimentation, cloud_water_sedimentation,                     &
276               collision_turbulence, limiter_sedimentation, nc_const,          &
277               ventilation_effect
278
279    USE model_1d,                                                              &
280        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
281
282    USE pegrid
283
284    USE netcdf_interface,                                                      &
285        ONLY:  netcdf_data_format, netcdf_deflate, netcdf_precision
286
287    USE pmc_interface,                                                         &
288        ONLY:  nested_run, nesting_mode
289
290    USE profil_parameter,                                                      &
291        ONLY:  cross_profiles, cross_ts_uymax, cross_ts_uymin,                 &
292               profile_columns, profile_rows
293
294    USE progress_bar,                                                          &
295        ONLY :  batch_job
296
297    USE radiation_model_mod,                                                   &
298        ONLY: radiation_parin 
299
300    USE spectra_mod,                                                           &
301        ONLY :  spectra_parin
302
303    USE statistics,                                                            &
304        ONLY:  hom, hom_sum, pr_palm, region, statistic_regions
305
306    USE urban_surface_mod,                                                     &
307        ONLY: usm_parin
308
309    USE wind_turbine_model_mod,                                                &
310        ONLY:  wtm_parin
311
312
313    IMPLICIT NONE
314
315    INTEGER(iwp) ::  i   !<
316
317
318    NAMELIST /inipar/  alpha_surface, approximation,                           &
319                       background_communication, bc_e_b, bc_lr,                &
320                       bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b,        &
321             bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t,                 &
322             bottom_salinityflux, building_height, building_length_x,          &
323             building_length_y, building_wall_left, building_wall_south,       &
324             call_psolver_at_all_substeps, call_microphysics_at_all_substeps,  &
325             canyon_height,                                                    &
326             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
327             canyon_wall_south, c_sedimentation, cfl_factor, cloud_droplets,   &
328             cloud_physics, cloud_scheme, cloud_top_radiation,                 &
329             cloud_water_sedimentation,                                        &
330             collective_wait, collision_turbulence, conserve_volume_flow,      &
331             conserve_volume_flow_mode, constant_flux_layer,                   &
332             coupling_start_time,                                              &
333             cycle_mg, damp_level_1d,                                          &
334             dissipation_1d,                                                   &
335             dp_external, dp_level_b, dp_smooth, dpdxy,                        &
336             dt, dt_pr_1d, dt_run_control_1d, dx, dy, dz, dz_max,              &
337             dz_stretch_factor, dz_stretch_level, end_time_1d,                 &
338             ensemble_member_nr, e_init, e_min, fft_method,                    &
339             flux_input_mode, flux_output_mode,                                &
340             galilei_transformation, humidity,                                 &
341             inflow_damping_height, inflow_damping_width,                      &
342             inflow_disturbance_begin, inflow_disturbance_end,                 &
343             initializing_actions, km_constant,                                &
344             large_scale_forcing, large_scale_subsidence,                      &
345             limiter_sedimentation,                                            &
346             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
347             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
348             most_method, nc_const, netcdf_precision, neutral, ngsrb,          &
349             nsor, nsor_ini, nudging, nx, ny, nz, ocean, omega, omega_sor,     &
350             outflow_source_plane, passive_scalar, phi,                        &
351             prandtl_number, precipitation, psolver, pt_damping_factor,        &
352             pt_damping_width, pt_reference, pt_surface,                       &
353             pt_surface_initial_change, pt_vertical_gradient,                  &
354             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
355             q_vertical_gradient, q_vertical_gradient_level,                   &
356             random_generator, random_heatflux,                                &
357             rayleigh_damping_factor, rayleigh_damping_height,                 &
358             recycling_width, recycling_yshift,                                &
359             reference_state, residual_limit,                                  &
360             roughness_length, sa_surface,                                     &
361             sa_vertical_gradient, sa_vertical_gradient_level, scalar_advec,   &
362             scalar_rayleigh_damping,                                          &
363             statistic_regions, subs_vertical_gradient,                        &
364             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
365             surface_scalarflux, surface_waterflux,                            &
366             s_surface, s_surface_initial_change, s_vertical_gradient,         &
367             s_vertical_gradient_level, timestep_scheme,                       &
368             topography, topography_grid_convention, top_heatflux,             &
369             top_momentumflux_u, top_momentumflux_v, top_salinityflux,         &
370             top_scalarflux, transpose_compute_overlap, turbulent_inflow,      &
371             turbulent_outflow,                                                &
372             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
373             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
374             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
375             uv_heights, u_bulk, u_profile, vg_surface, vg_vertical_gradient,  &
376             vg_vertical_gradient_level, v_bulk, v_profile, ventilation_effect,&
377             wall_adjustment, wall_heatflux, wall_humidityflux,                &
378             wall_scalarflux, zeta_max, zeta_min, z0h_factor
379     
380    NAMELIST /d3par/  averaging_interval, averaging_interval_pr,               &
381             cpu_log_barrierwait, create_disturbances,                         &
382             cross_profiles, cross_ts_uymax, cross_ts_uymin,                   &
383             data_output, data_output_masks,                                   &
384             data_output_pr, data_output_2d_on_each_pe, disturbance_amplitude, &
385             disturbance_energy_limit, disturbance_level_b,                    &
386             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
387             dt, dt_averaging_input, dt_averaging_input_pr,                    &
388             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
389             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
390             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
391             dt_run_control,end_time, force_print_header, mask_scale_x,        &
392             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
393             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
394             normalizing_region, npex, npey, nz_do3d,                          &
395             precipitation_amount_interval, profile_columns, profile_rows,     &
396             restart_time, section_xy, section_xz, section_yz,                 &
397             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
398             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
399             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
400             termination_time_needed, z_max_do2d
401
402
403    NAMELIST /envpar/  batch_job, host, local_dvrserver_running,               &
404                       maximum_cpu_time_allowed, maximum_parallel_io_streams,  &
405                       revision, return_addres, return_username,               &
406                       run_identifier, tasks_per_node, write_binary
407
408!
409!-- First read values of environment variables (this NAMELIST file is
410!-- generated by mrun)
411    CALL location_message( 'reading environment parameters from ENVPAR', .FALSE. )
412    OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', ERR=30 )
413    READ ( 90, envpar, ERR=31, END=32 )
414    CLOSE ( 90 )
415    CALL location_message( 'finished', .TRUE. )
416
417!
418!-- Calculate the number of groups into which parallel I/O is split.
419!-- The default for files which are opened by all PEs (or where each
420!-- PE opens his own independent file) is, that all PEs are doing input/output
421!-- in parallel at the same time. This might cause performance or even more
422!-- severe problems depending on the configuration of the underlying file
423!-- system.
424!-- First, set the default:
425    IF ( maximum_parallel_io_streams == -1  .OR.                               &
426         maximum_parallel_io_streams > numprocs )  THEN
427       maximum_parallel_io_streams = numprocs
428    ENDIF
429!
430!-- Now calculate the number of io_blocks and the io_group to which the
431!-- respective PE belongs. I/O of the groups is done in serial, but in parallel
432!-- for all PEs belonging to the same group.
433!-- These settings are repeated in init_pegrid for the communicator comm2d,
434!-- which is not available here
435    io_blocks = numprocs / maximum_parallel_io_streams
436    io_group  = MOD( myid+1, io_blocks )
437
438    CALL location_message( 'reading NAMELIST parameters from PARIN', .FALSE. )
439!
440!-- Data is read in parallel by groups of PEs
441    DO  i = 0, io_blocks-1
442       IF ( i == io_group )  THEN
443
444!
445!--       Open the NAMELIST-file which is send with this job
446          CALL check_open( 11 )
447
448!
449!--       Read the control parameters for initialization.
450!--       The namelist "inipar" must be provided in the NAMELIST-file.
451          READ ( 11, inipar, ERR=10, END=11 )
452
453          GOTO 12
454
455 10       message_string = 'errors in \$inipar &or no \$inipar-namelist ' //   &
456                           'found (CRAY-machines only)'
457          CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
458
459 11       message_string = 'no \$inipar-namelist found'
460          CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
461
462!
463!--       If required, read control parameters from restart file (produced by
464!--       a prior run). All PEs are reading from file created by PE0 (see
465!--       check_open)
466 12       IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
467             CALL read_var_list
468!
469!--          The restart file will be reopened when reading the subdomain data
470             CALL close_file( 13 )
471
472!
473!--          Increment the run count
474             runnr = runnr + 1
475          ENDIF
476
477!
478!--       In case of nested runs, explicitly set nesting boundary conditions.
479!--       This will overwrite the user settings. bc_lr and bc_ns always need
480!--       to be cyclic for vertical nesting.
481          IF ( nesting_mode == 'vertical' )  THEN
482             IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
483                WRITE ( message_string, *)  'bc_lr and bc_ns were set to ,',   &
484                                      'cyclic for vertical nesting'
485                CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
486                bc_lr   = 'cyclic'
487                bc_ns   = 'cyclic'
488             ENDIF
489             IF ( nest_domain )  THEN
490                bc_uv_t = 'nested'
491                bc_pt_t = 'nested'
492                bc_q_t  = 'nested'
493                bc_s_t  = 'nested'
494                bc_p_t  = 'neumann'
495             ENDIF
496          ELSE
497         
498!
499!--       For other nesting modes only set boundary conditions for
500!--       nested domains.
501             IF ( nest_domain )  THEN
502                bc_lr   = 'nested'
503                bc_ns   = 'nested'
504                bc_uv_t = 'nested'
505                bc_pt_t = 'nested'
506                bc_q_t  = 'nested'
507                bc_s_t  = 'nested'
508                bc_p_t  = 'neumann'
509             ENDIF
510          ENDIF
511
512!         
513!--       In case of nested runs, make sure that initializing_actions =
514!--       'set_constant_profiles' even though the constant-profiles
515!--       initializations for the prognostic variables will be overwritten
516!--       by pmci_child_initialize and pmci_parent_initialize. This is,
517!--       however, important e.g. to make sure that diagnostic variables
518!--       are set properly.
519          IF ( nest_domain )  THEN
520             initializing_actions = 'set_constant_profiles'
521          ENDIF
522
523           
524!
525!--       Check validity of lateral boundary conditions. This has to be done
526!--       here because they are already used in init_pegrid and init_grid and
527!--       therefore cannot be check in check_parameters
528          IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
529               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested' )  THEN
530             message_string = 'unknown boundary condition: bc_lr = "' // &
531                              TRIM( bc_lr ) // '"'
532             CALL message( 'check_parameters', 'PA0049', 1, 2, 0, 6, 0 )
533          ENDIF
534          IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
535               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested' )  THEN
536             message_string = 'unknown boundary condition: bc_ns = "' // &
537                              TRIM( bc_ns ) // '"'
538             CALL message( 'check_parameters', 'PA0050', 1, 2, 0, 6, 0 )
539          ENDIF
540
541!
542!--       Set internal variables used for speed optimization in if clauses
543          IF ( bc_lr /= 'cyclic' )               bc_lr_cyc    = .FALSE.
544          IF ( bc_lr == 'dirichlet/radiation' )  bc_lr_dirrad = .TRUE.
545          IF ( bc_lr == 'radiation/dirichlet' )  bc_lr_raddir = .TRUE.
546          IF ( bc_ns /= 'cyclic' )               bc_ns_cyc    = .FALSE.
547          IF ( bc_ns == 'dirichlet/radiation' )  bc_ns_dirrad = .TRUE.
548          IF ( bc_ns == 'radiation/dirichlet' )  bc_ns_raddir = .TRUE.
549
550!
551!--       Definition of names of areas used for computing statistics. They must
552!--       be defined at this place, because they are allowed to be redefined by
553!--       the user in user_parin.
554          region = 'total domain'
555
556!
557!--       Read runtime parameters given by the user for this run (namelist
558!--       "d3par"). The namelist "d3par" can be omitted. In that case, default
559!--       values are used for the parameters.
560          READ ( 11, d3par, END=20 )
561
562!
563!--       Check if land surface model is used and read &lsm_par if required
564          CALL lsm_parin
565
566!
567!--       Check if urban surface model is used and read &urban_surface_par if required
568          CALL usm_parin
569
570!
571!--       Check if spectra shall be calculated and read spectra_par if required
572          CALL spectra_parin
573
574!
575!--       Check if radiation model is used and read &radiation_par if required
576          CALL radiation_parin
577 
578 
579!--       Check if plant canopy model is used and read &canopy_par if required
580          CALL pcm_parin
581 
582!
583!--       Read control parameters for optionally used model software packages
584 20       CALL package_parin
585
586!
587!--       Check if wind turbine model is used and read &wind_turbine_par if
588!--       required
589          CALL wtm_parin
590!
591!--       Check if virtual flights should be carried out and read &flight_part
592!--       if required
593          CALL flight_parin
594
595!
596!--       Read user-defined variables
597          CALL user_parin
598
599!
600!--       Check in case of initial run, if the grid point numbers are well
601!--       defined and allocate some arrays which are already needed in
602!--       init_pegrid or check_parameters. During restart jobs, these arrays
603!--       will be allocated in read_var_list. All other arrays are allocated
604!--       in init_3d_model.
605          IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
606
607             IF ( nx <= 0 )  THEN
608                WRITE( message_string, * ) 'no value or wrong value given',    &
609                                           ' for nx: nx=', nx
610                CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 )
611             ENDIF
612             IF ( ny <= 0 )  THEN
613                WRITE( message_string, * ) 'no value or wrong value given',    &
614                                           ' for ny: ny=', ny
615                CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 )
616             ENDIF
617             IF ( nz <= 0 )  THEN
618                WRITE( message_string, * ) 'no value or wrong value given',    &
619                                           ' for nz: nz=', nz
620                CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 )
621             ENDIF
622!
623!--          ATTENTION: in case of changes to the following statement please
624!--                  also check the allocate statement in routine read_var_list
625             ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1),        &
626                       ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1),         &
627                       u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1),             &
628                       hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions),  &
629                       hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
630
631             hom = 0.0_wp
632
633          ENDIF
634
635!
636!--       NAMELIST-file is not needed anymore
637          CALL close_file( 11 )
638
639       ENDIF
640#if defined( __parallel )
641       CALL MPI_BARRIER( MPI_COMM_WORLD, ierr )
642#endif
643    ENDDO
644
645    CALL location_message( 'finished', .TRUE. )
646
647    RETURN
648
649 30 message_string = 'local file ENVPAR not found' //                          &
650                     '&some variables for steering may not be properly set'
651    CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
652    RETURN
653
654 31 message_string = 'errors in local file ENVPAR' //                          &
655                     '&some variables for steering may not be properly set'
656    CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
657    RETURN
658
659 32 message_string = 'no envpar-NAMELIST found in local file ENVPAR'  //       &
660                     '&some variables for steering may not be properly set'
661    CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
662
663 END SUBROUTINE parin
Note: See TracBrowser for help on using the repository browser.