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

Last change on this file since 4186 was 4183, checked in by oliver.maas, 5 years ago

simplified steering of recycling of absolute values. Replaced initialization parameter recycle_absolute_quantities by recycling_method_for_thermodynamic_quantities.

  • Property svn:keywords set to Id
File size: 34.9 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-2019 Leibniz Universitaet Hannover
18!------------------------------------------------------------------------------!
19!
20! Current revisions:
21! -----------------
22!
23!
24! Former revisions:
25! -----------------
26! $Id: parin.f90 4183 2019-08-23 07:33:16Z suehring $
27! replaced recycle_absolute_quantities by recycling_method_for_thermodynamic_quantities
28!
29! 4182 2019-08-22 15:20:23Z scharf
30! Corrected "Former revisions" section
31!
32! 4176 2019-08-20 14:10:41Z oliver.maas
33! added recycle_absolute_quantities to initialization_parameters namelist
34!
35! 4173 2019-08-20 12:04:06Z gronemeier
36! add vdi_internal_controls
37!
38! 4131 2019-08-02 11:06:18Z monakurppa
39! Allocate hom and hom_sum to allow profile output for salsa variables.
40!
41! 4079 2019-07-09 18:04:41Z suehring
42! +monotonic_limiter_z
43!
44! 4022 2019-06-12 11:52:39Z suehring
45! Change default top boundary condition for pressure to Neumann in offline
46! nesting case
47!
48! 4017 2019-06-06 12:16:46Z schwenkel
49! Introduce alternative switch for debug output during timestepping
50!
51! 3885 2019-04-11 11:29:34Z kanani
52! Changes related to global restructuring of location messages and introduction
53! of additional debug messages
54!
55! 3806 2019-03-21 12:45:50Z raasch
56! additional check for lateral boundary conditions added
57!
58! 3747 2019-02-16 15:15:23Z gronemeier
59! removed setting of parameter region
60!
61! 3746 2019-02-16 12:41:27Z gronemeier
62! Removed most_method
63!
64! 3649 2019-01-02 16:52:21Z suehring
65! Delete debug-print statements
66!
67! Revision 1.1  1997/07/24 11:22:50  raasch
68! Initial revision
69!
70!
71! Description:
72! ------------
73!> This subroutine reads variables controling the run from the NAMELIST files
74!>
75!> @todo: Revise max_pr_cs (profiles for chemistry)
76!------------------------------------------------------------------------------!
77 SUBROUTINE parin
78 
79
80    USE arrays_3d,                                                             &
81        ONLY:  pt_init, q_init, ref_state, s_init, sa_init,                    &
82               ug, u_init, v_init, vg
83
84    USE chem_modules
85
86    USE control_parameters
87
88    USE cpulog,                                                                &
89        ONLY:  cpu_log_barrierwait
90
91    USE date_and_time_mod,                                                     &
92        ONLY:  date_init, day_of_year_init, time_utc_init
93
94    USE grid_variables,                                                        &
95        ONLY:  dx, dy
96
97    USE indices,                                                               &
98        ONLY:  nx, ny, nz
99
100    USE kinds
101
102    USE model_1d_mod,                                                          &
103        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
104
105    USE module_interface,                                                      &
106        ONLY:  module_interface_parin
107
108    USE netcdf_interface,                                                      &
109        ONLY:  netcdf_data_format, netcdf_deflate, netcdf_precision
110
111    USE pegrid
112
113    USE pmc_interface,                                                         &
114        ONLY:  nested_run, nesting_mode
115
116    USE profil_parameter,                                                      &
117        ONLY:  cross_profiles, profile_columns, profile_rows
118
119    USE progress_bar,                                                          &
120        ONLY :  progress_bar_disabled
121
122    USE read_restart_data_mod,                                                 &
123        ONLY:  rrd_global
124
125    USE statistics,                                                            &
126        ONLY:  hom, hom_sum, pr_palm, statistic_regions
127
128    USE turbulence_closure_mod,                                                &
129        ONLY:  rans_const_c, rans_const_sigma
130
131    USE vertical_nesting_mod,                                                  &
132        ONLY:  vnest_start_time
133
134
135    IMPLICIT NONE
136
137    CHARACTER (LEN=80) ::  line  !< dummy string that contains the current line of the parameter file
138
139    INTEGER(iwp) ::  global_id      !< process id with respect to MPI_COMM_WORLD
140    INTEGER(iwp) ::  global_procs   !< # of procs with respect to MPI_COMM_WORLD
141    INTEGER(iwp) ::  i              !<
142    INTEGER(iwp) ::  ioerr          !< error flag for open/read/write
143
144    NAMELIST /inipar/  alpha_surface, approximation, bc_e_b,     &
145                       bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b, &
146             bc_q_t,bc_s_b, bc_s_t, bc_uv_b, bc_uv_t,                 &
147             building_height, building_length_x,          &
148             building_length_y, building_wall_left, building_wall_south,       &
149             calc_soil_moisture_during_spinup,                                 &
150             call_psolver_at_all_substeps,  &
151             canyon_height,                                                    &
152             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
153             canyon_wall_south, cfl_factor, cloud_droplets,   &
154             collective_wait, complex_terrain,           &
155             conserve_volume_flow,                                             &
156             conserve_volume_flow_mode, constant_flux_layer,                   &
157             coupling_start_time,             &
158             cycle_mg, damp_level_1d,                                          &
159             data_output_during_spinup,                                        &
160             date_init,                                                        &
161             day_of_year_init,                                                 &
162             dissipation_1d,                                                   &
163             dp_external, dp_level_b, dp_smooth, dpdxy,    &
164             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
165             dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
166             dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
167             e_min, fft_method, flux_input_mode, flux_output_mode,             &
168             galilei_transformation, humidity,                                 &
169             inflow_damping_height, inflow_damping_width,                      &
170             inflow_disturbance_begin, inflow_disturbance_end,                 &
171             initializing_actions, km_constant,                                &
172             large_scale_forcing, large_scale_subsidence, latitude,            &
173             longitude,                                 &
174             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
175             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
176             monotonic_limiter_z,                                              &
177             netcdf_precision, neutral, ngsrb,                                 &
178             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
179             omega_sor, outflow_source_plane, passive_scalar,                  &
180             prandtl_number, psolver, pt_damping_factor,        &
181             pt_damping_width, pt_reference, pt_surface,                       &
182             pt_surface_initial_change, pt_vertical_gradient,                  &
183             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
184             q_vertical_gradient, q_vertical_gradient_level,                   &
185             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
186             rayleigh_damping_factor, rayleigh_damping_height,                 &
187             recycling_width, recycling_yshift,                                &
188             reference_state, residual_limit,                                  &
189             roughness_length,                                                 &
190             scalar_advec,   &
191             scalar_rayleigh_damping,                              &
192             spinup_time, spinup_pt_amplitude, spinup_pt_mean,                 &
193             statistic_regions, subs_vertical_gradient,                        &
194             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
195             surface_scalarflux, surface_waterflux,                            &
196             s_surface, s_surface_initial_change, s_vertical_gradient,         &
197             s_vertical_gradient_level, time_utc_init, timestep_scheme,        &
198             topography, topography_grid_convention, top_heatflux,             &
199             top_momentumflux_u, top_momentumflux_v,                           &
200             top_scalarflux, transpose_compute_overlap,                        &
201             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
202             tunnel_wall_depth, turbulence_closure,                            &
203             turbulent_inflow, turbulent_outflow,                              &
204             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
205             use_free_convection_scaling,                                      &
206             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
207             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
208             uv_heights, u_bulk, u_profile, vdi_checks, vg_surface,            &
209             vg_vertical_gradient,  &
210             vg_vertical_gradient_level, v_bulk, v_profile,&
211             wall_adjustment, wall_heatflux, wall_humidityflux,                &
212             wall_scalarflux, y_shift, zeta_max, zeta_min,  &
213             z0h_factor
214
215    NAMELIST /initialization_parameters/  alpha_surface,                       &
216             approximation, bc_e_b,                                            &
217             bc_lr, bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b,           &
218             bc_q_t,bc_s_b, bc_s_t, bc_uv_b, bc_uv_t,                          &
219             building_height, building_length_x,                               &
220             building_length_y, building_wall_left, building_wall_south,       &
221             calc_soil_moisture_during_spinup,                                 &
222             call_psolver_at_all_substeps,                                     &
223             canyon_height,                                                    &
224             canyon_width_x, canyon_width_y, canyon_wall_left,                 &
225             canyon_wall_south, cfl_factor, cloud_droplets,                    &
226             collective_wait, complex_terrain,                                 &
227             conserve_volume_flow,                                             &
228             conserve_volume_flow_mode, constant_flux_layer,                   &
229             coupling_start_time,                                              &
230             cycle_mg, damp_level_1d,                                          &
231             data_output_during_spinup,                                        &
232             date_init,                                                        &
233             day_of_year_init,                                                 &
234             dissipation_1d,                                                   &
235             dp_external, dp_level_b, dp_smooth, dpdxy,                        &
236             dt, dt_pr_1d, dt_run_control_1d, dt_spinup, dx, dy, dz, dz_max,   &
237             dz_stretch_factor, dz_stretch_level, dz_stretch_level_start,      &
238             dz_stretch_level_end, end_time_1d, ensemble_member_nr, e_init,    &
239             e_min, fft_method, flux_input_mode, flux_output_mode,             &
240             galilei_transformation, humidity,                                 &
241             inflow_damping_height, inflow_damping_width,                      &
242             inflow_disturbance_begin, inflow_disturbance_end,                 &
243             initializing_actions, km_constant,                                &
244             large_scale_forcing, large_scale_subsidence, latitude,            &
245             longitude,                                                        &
246             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
247             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
248             monotonic_limiter_z,                                              &
249             netcdf_precision, neutral, ngsrb,                                 &
250             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
251             omega_sor, outflow_source_plane, passive_scalar,                  &
252             prandtl_number, psolver, pt_damping_factor,                       &
253             pt_damping_width, pt_reference, pt_surface,                       &
254             pt_surface_initial_change, pt_vertical_gradient,                  &
255             pt_vertical_gradient_level, q_surface, q_surface_initial_change,  &
256             q_vertical_gradient, q_vertical_gradient_level,                   &
257             random_generator, random_heatflux, rans_const_c, rans_const_sigma,&
258             rayleigh_damping_factor, rayleigh_damping_height,                 &
259             recycling_width, recycling_yshift,                                &
260             recycling_method_for_thermodynamic_quantities,                    &
261             reference_state, residual_limit,                                  &
262             roughness_length, scalar_advec,                                   &
263             scalar_rayleigh_damping,                                          &
264             spinup_time, spinup_pt_amplitude, spinup_pt_mean,                 &
265             statistic_regions, subs_vertical_gradient,                        &
266             subs_vertical_gradient_level, surface_heatflux, surface_pressure, &
267             surface_scalarflux, surface_waterflux,                            &
268             s_surface, s_surface_initial_change, s_vertical_gradient,         &
269             s_vertical_gradient_level, time_utc_init, timestep_scheme,        &
270             topography, topography_grid_convention, top_heatflux,             &
271             top_momentumflux_u, top_momentumflux_v,                           &
272             top_scalarflux, transpose_compute_overlap,                        &
273             tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,     &
274             tunnel_wall_depth, turbulence_closure,                            &
275             turbulent_inflow, turbulent_outflow,                              &
276             use_subsidence_tendencies, ug_surface, ug_vertical_gradient,      &
277             ug_vertical_gradient_level, use_surface_fluxes, use_cmax,         &
278             use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke,      &
279             use_free_convection_scaling,                                      &
280             uv_heights, u_bulk, u_profile, vdi_checks,                        &
281             vg_surface, vg_vertical_gradient,  &
282             vg_vertical_gradient_level, v_bulk, v_profile,                    &
283             wall_adjustment, wall_heatflux, wall_humidityflux,                &
284             wall_scalarflux, y_shift, zeta_max, zeta_min, z0h_factor
285             
286    NAMELIST /d3par/  averaging_interval, averaging_interval_pr,               &
287             cpu_log_barrierwait, create_disturbances,                         &
288             cross_profiles, data_output, data_output_masks,                   &
289             data_output_pr, data_output_2d_on_each_pe,                        &
290             debug_output,                                                     &
291             debug_output_timestep,                                            &
292             disturbance_amplitude,                                            &
293             disturbance_energy_limit, disturbance_level_b,                    &
294             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
295             dt, dt_averaging_input, dt_averaging_input_pr,                    &
296             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
297             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
298             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
299             dt_run_control,end_time, force_print_header, mask_k_over_surface, &
300             mask_scale_x,                                                     &
301             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
302             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
303             normalizing_region, npex, npey, nz_do3d,                          &
304             profile_columns, profile_rows,     &
305             restart_time, section_xy, section_xz, section_yz,                 &
306             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
307             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
308             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
309             termination_time_needed, vnest_start_time
310
311    NAMELIST /runtime_parameters/  averaging_interval, averaging_interval_pr,  &
312             cpu_log_barrierwait, create_disturbances,                         &
313             cross_profiles, data_output, data_output_masks,                   &
314             data_output_pr, data_output_2d_on_each_pe,                        &
315             debug_output,                                                     &
316             debug_output_timestep,                                            &
317             disturbance_amplitude,                                            &
318             disturbance_energy_limit, disturbance_level_b,                    &
319             disturbance_level_t, do2d_at_begin, do3d_at_begin,                &
320             dt, dt_averaging_input, dt_averaging_input_pr,                    &
321             dt_coupling, dt_data_output, dt_data_output_av, dt_disturb,       &
322             dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy,         &
323             dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart,              &
324             dt_run_control,end_time, force_print_header, mask_k_over_surface, &
325             mask_scale_x,                                                     &
326             mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop,  &
327             mask_y_loop, mask_z_loop, netcdf_data_format, netcdf_deflate,     &
328             normalizing_region, npex, npey, nz_do3d,                          &
329             profile_columns, profile_rows,     &
330             restart_time, section_xy, section_xz, section_yz,                 &
331             skip_time_data_output, skip_time_data_output_av, skip_time_dopr,  &
332             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
333             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
334             termination_time_needed, vnest_start_time
335
336    NAMELIST /envpar/  progress_bar_disabled, host,                            &
337                       maximum_cpu_time_allowed, maximum_parallel_io_streams,  &
338                       read_svf, revision, run_identifier, tasks_per_node,     &
339                       write_binary, write_svf
340
341!
342!-- First read values of environment variables (this NAMELIST file is
343!-- generated by palmrun)
344    CALL location_message( 'reading environment parameters from ENVPAR', 'start' )
345
346    OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', IOSTAT=ioerr )
347
348    IF ( ioerr /= 0 )  THEN
349       message_string = 'local file ENVPAR not found' //                       &
350                        '&some variables for steering may not be properly set'
351       CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
352    ELSE
353       READ ( 90, envpar, IOSTAT=ioerr )
354       IF ( ioerr < 0 )  THEN
355          message_string = 'no envpar-NAMELIST found in local file '  //       &
356                           'ENVPAR& or some variables for steering may '  //   &
357                           'not be properly set'
358          CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
359       ELSEIF ( ioerr > 0 )  THEN
360          message_string = 'errors in local file ENVPAR' //                    &
361                           '&some variables for steering may not be properly set'
362          CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
363       ENDIF
364       CLOSE ( 90 )
365    ENDIF
366
367    CALL location_message( 'reading environment parameters from ENVPAR', 'finished' )
368!
369!-- Calculate the number of groups into which parallel I/O is split.
370!-- The default for files which are opened by all PEs (or where each
371!-- PE opens his own independent file) is, that all PEs are doing input/output
372!-- in parallel at the same time. This might cause performance or even more
373!-- severe problems depending on the configuration of the underlying file
374!-- system.
375!-- Calculation of the number of blocks and the I/O group must be based on all
376!-- PEs involved in this run. Since myid and numprocs are related to the
377!-- comm2d communicator, which gives only a subset of all PEs in case of
378!-- nested runs, that information must be inquired again from the global
379!-- communicator.
380!-- First, set the default:
381#if defined( __parallel )
382    CALL MPI_COMM_RANK( MPI_COMM_WORLD, global_id, ierr )
383    CALL MPI_COMM_SIZE( MPI_COMM_WORLD, global_procs, ierr )
384#else
385    global_id    = 0
386    global_procs = 1
387#endif
388    IF ( maximum_parallel_io_streams == -1  .OR.                               &
389         maximum_parallel_io_streams > global_procs )  THEN
390       maximum_parallel_io_streams = global_procs
391    ENDIF
392!
393!-- Now calculate the number of io_blocks and the io_group to which the
394!-- respective PE belongs. I/O of the groups is done in serial, but in parallel
395!-- for all PEs belonging to the same group.
396    io_blocks = global_procs / maximum_parallel_io_streams
397    io_group  = MOD( global_id+1, io_blocks )
398   
399    CALL location_message( 'reading NAMELIST parameters from PARIN', 'start' )
400!
401!-- Data is read in parallel by groups of PEs
402    DO  i = 0, io_blocks-1
403       IF ( i == io_group )  THEN
404
405!
406!--       Open the NAMELIST-file which is send with this job
407          CALL check_open( 11 )
408
409!
410!--       Read the control parameters for initialization.
411!--       The namelist "inipar" must be provided in the NAMELIST-file.
412          READ ( 11, initialization_parameters, ERR=10, END=11 )
413          GOTO 14
414         
415 10       BACKSPACE( 11 )
416          READ( 11 , '(A)') line
417          CALL parin_fail_message( 'initialization_parameters', line )
418
419 11       REWIND ( 11 )
420          READ ( 11, inipar, ERR=12, END=13 )
421 
422          message_string = 'namelist inipar is deprecated and will be ' //    &
423                          'removed in near future. & Please use namelist ' // &
424                          'initialization_parameters instead'
425          CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 )
426 
427          GOTO 14
428 
429 12       BACKSPACE( 11 )
430          READ( 11 , '(A)') line
431          CALL parin_fail_message( 'inipar', line )
432
433 13       message_string = 'no initialization_parameters-namelist found'
434          CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
435
436!
437!--       Try to read runtime parameters given by the user for this run
438!--       (namelist "runtime_parameters"). The namelist "runtime_parmeters"   
439!--       can be omitted. In that case default values are used for the         
440!--       parameters.
441 14       line = ' '
442
443          REWIND ( 11 )
444          line = ' '
445          DO WHILE ( INDEX( line, '&runtime_parameters' ) == 0 )
446             READ ( 11, '(A)', END=16 )  line
447          ENDDO
448          BACKSPACE ( 11 )
449
450!
451!--       Read namelist
452          READ ( 11, runtime_parameters, ERR = 15 )
453          GOTO 18
454
455 15       BACKSPACE( 11 )
456          READ( 11 , '(A)') line
457          CALL parin_fail_message( 'runtime_parameters', line )
458
459 16       REWIND ( 11 )
460          line = ' '
461          DO WHILE ( INDEX( line, '&d3par' ) == 0 )
462             READ ( 11, '(A)', END=18 )  line
463          ENDDO
464          BACKSPACE ( 11 )
465
466!
467!--       Read namelist
468          READ ( 11, d3par, ERR = 17, END = 18 )
469
470          message_string = 'namelist d3par is deprecated and will be ' //      &
471                          'removed in near future. &Please use namelist ' //   &
472                          'runtime_parameters instead'
473          CALL message( 'parin', 'PA0487', 0, 1, 0, 6, 0 )
474
475          GOTO 18
476
477 17       BACKSPACE( 11 )
478          READ( 11 , '(A)') line
479          CALL parin_fail_message( 'd3par', line )
480
481 18       CONTINUE
482
483!
484!--       Check for module namelists and read them
485          CALL module_interface_parin
486
487!
488!--       If required, read control parameters from restart file (produced by
489!--       a prior run). All PEs are reading from file created by PE0 (see
490!--       check_open)
491          IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
492
493             CALL rrd_global
494!
495!--          Increment the run count
496             runnr = runnr + 1
497!
498!--          In case of a restart run, the number of user-defined profiles on
499!--          the restart file (already stored in max_pr_user) has to match the
500!--          one given for the current run. max_pr_user_tmp is calculated in
501!--          user_parin and max_pr_user is read in via rrd_global.
502             IF ( max_pr_user /= max_pr_user_tmp )  THEN
503                WRITE( message_string, * ) 'the number of user-defined ',      &
504                      'profiles given in data_output_pr (', max_pr_user_tmp,   &
505                      ') does not match the one ',                             &
506                      'found in the restart file (', max_pr_user, ')'
507                CALL message( 'user_parin', 'UI0009', 1, 2, 0, 6, 0 )
508             ENDIF
509          ELSE
510             max_pr_user = max_pr_user_tmp
511          ENDIF
512!
513!--       Activate spinup
514          IF ( land_surface .OR. urban_surface )  THEN
515             IF ( spinup_time > 0.0_wp )  THEN
516                coupling_start_time = spinup_time
517                time_since_reference_point = simulated_time - coupling_start_time
518                IF ( spinup_pt_mean == 9999999.9_wp )  THEN
519                   spinup_pt_mean = pt_surface
520                ENDIF
521                end_time = end_time + spinup_time
522                IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
523                   spinup = .TRUE.
524                ELSEIF (        TRIM( initializing_actions ) == 'read_restart_data'      &
525                         .AND.  time_since_reference_point > 0.0_wp )  THEN
526                   data_output_during_spinup = .FALSE.  !< required for correct ntdim calculation
527                                                        !< in check_parameters for restart run
528                ENDIF
529             ENDIF
530          ENDIF
531
532!
533!--       In case of nested runs, explicitly set nesting boundary conditions.
534!--       This will overwrite the user settings and basic defaults.
535!--       bc_lr and bc_ns always need to be cyclic for vertical nesting.
536          IF ( nested_run )  THEN
537             IF ( nesting_mode == 'vertical' )  THEN
538                IF (bc_lr /= 'cyclic' .OR. bc_ns /= 'cyclic' )  THEN
539                   WRITE ( message_string, *) 'bc_lr and bc_ns were set to ,', &
540                        'cyclic for vertical nesting'
541                   CALL message( 'parin', 'PA0428', 0, 0, 0, 6, 0 )
542                   bc_lr   = 'cyclic'
543                   bc_ns   = 'cyclic'
544                ENDIF
545                IF ( child_domain )  THEN
546                   bc_uv_t  = 'nested'
547                   bc_pt_t  = 'nested'
548                   bc_q_t   = 'nested'
549                   bc_s_t   = 'nested'
550                   bc_cs_t  = 'nested'
551                   bc_p_t   = 'neumann' 
552                ENDIF
553!
554!--          For other nesting modes only set boundary conditions for
555!--          nested domains.
556             ELSE
557                IF ( child_domain )  THEN
558                   bc_lr    = 'nested'
559                   bc_ns    = 'nested'
560                   bc_uv_t  = 'nested'
561                   bc_pt_t  = 'nested'
562                   bc_q_t   = 'nested'
563                   bc_s_t   = 'nested'
564                   bc_cs_t  = 'nested'
565                   bc_p_t   = 'neumann'
566                ENDIF
567             ENDIF
568          ENDIF
569!
570!--       Set boundary conditions also in case the model is offline-nested in
571!--       larger-scale models.
572          IF ( nesting_offline )  THEN
573             bc_lr    = 'nesting_offline'
574             bc_ns    = 'nesting_offline'
575             bc_uv_t  = 'nesting_offline'
576             bc_pt_t  = 'nesting_offline'
577             bc_q_t   = 'nesting_offline'
578           !  bc_s_t   = 'nesting_offline'  ! scalar boundary condition is not clear yet
579           !  bc_cs_t  = 'nesting_offline'  ! same for chemical species
580             bc_p_t   = 'neumann'
581          ENDIF
582
583!         
584!--       In case of nested runs, make sure that initializing_actions =
585!--       'set_constant_profiles' even though the constant-profiles
586!--       initializations for the prognostic variables will be overwritten
587!--       by pmci_child_initialize and pmci_parent_initialize. This is,
588!--       however, important e.g. to make sure that diagnostic variables
589!--       are set properly. An exception is made in case of restart runs and
590!--       if user decides to do everything by its own.
591          IF ( child_domain  .AND.  .NOT. (                                    &
592               TRIM( initializing_actions ) == 'read_restart_data'      .OR.   &
593               TRIM( initializing_actions ) == 'set_constant_profiles'  .OR.   &
594               TRIM( initializing_actions ) == 'by_user' ) )  THEN
595             message_string = 'initializing_actions = ' //                     &
596                              TRIM( initializing_actions ) // ' has been ' //  &
597                              'changed to set_constant_profiles in child ' //  &
598                              'domain.' 
599             CALL message( 'parin', 'PA0492', 0, 0, 0, 6, 0 )
600
601             initializing_actions = 'set_constant_profiles'
602          ENDIF           
603!
604!--       Check validity of lateral boundary conditions. This has to be done
605!--       here because they are already used in init_pegrid and init_grid and
606!--       therefore cannot be check in check_parameters
607          IF ( bc_lr /= 'cyclic'  .AND.  bc_lr /= 'dirichlet/radiation'  .AND. &
608               bc_lr /= 'radiation/dirichlet'  .AND.  bc_lr /= 'nested'  .AND. &
609               bc_lr /= 'nesting_offline' )  THEN
610             message_string = 'unknown boundary condition: bc_lr = "' // &
611                              TRIM( bc_lr ) // '"'
612             CALL message( 'parin', 'PA0049', 1, 2, 0, 6, 0 )
613          ENDIF
614          IF ( bc_ns /= 'cyclic'  .AND.  bc_ns /= 'dirichlet/radiation'  .AND. &
615               bc_ns /= 'radiation/dirichlet'  .AND.  bc_ns /= 'nested'  .AND. &
616               bc_ns /= 'nesting_offline' )  THEN
617             message_string = 'unknown boundary condition: bc_ns = "' // &
618                              TRIM( bc_ns ) // '"'
619             CALL message( 'parin', 'PA0050', 1, 2, 0, 6, 0 )
620          ENDIF
621!
622!--       Set internal variables used for speed optimization in if clauses
623          IF ( bc_lr /= 'cyclic' )               bc_lr_cyc    = .FALSE.
624          IF ( bc_lr == 'dirichlet/radiation' )  bc_lr_dirrad = .TRUE.
625          IF ( bc_lr == 'radiation/dirichlet' )  bc_lr_raddir = .TRUE.
626          IF ( bc_ns /= 'cyclic' )               bc_ns_cyc    = .FALSE.
627          IF ( bc_ns == 'dirichlet/radiation' )  bc_ns_dirrad = .TRUE.
628          IF ( bc_ns == 'radiation/dirichlet' )  bc_ns_raddir = .TRUE.
629!
630!--       Radiation-Dirichlet conditions are allowed along one of the horizontal directions only.
631!--       In general, such conditions along x and y may work, but require a) some code changes
632!--       (e.g. concerning allocation of c_u, c_v ... arrays), and b) a careful model setup by the
633!--       user, in order to guarantee that there is a clearly defined outflow at two sides.
634!--       Otherwise, the radiation condition may produce wrong results.
635          IF ( ( bc_lr_dirrad .OR. bc_lr_raddir )  .AND.  ( bc_ns_dirrad .OR. bc_ns_raddir ) )  THEN
636             message_string = 'bc_lr = "' // TRIM( bc_lr ) // '" and bc_ns = "' //                 &
637                              TRIM( bc_ns ) // '" are not allowed to be set at the same time'
638             CALL message( 'parin', 'PA0589', 1, 2, 0, 6, 0 )
639          ENDIF
640!
641!--       Check in case of initial run, if the grid point numbers are well
642!--       defined and allocate some arrays which are already needed in
643!--       init_pegrid or check_parameters. During restart jobs, these arrays
644!--       will be allocated in rrd_global. All other arrays are allocated
645!--       in init_3d_model.
646          IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
647
648             IF ( nx <= 0 )  THEN
649                WRITE( message_string, * ) 'no value or wrong value given',    &
650                                           ' for nx: nx=', nx
651                CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 )
652             ENDIF
653             IF ( ny <= 0 )  THEN
654                WRITE( message_string, * ) 'no value or wrong value given',    &
655                                           ' for ny: ny=', ny
656                CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 )
657             ENDIF
658             IF ( nz <= 0 )  THEN
659                WRITE( message_string, * ) 'no value or wrong value given',    &
660                                           ' for nz: nz=', nz
661                CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 )
662             ENDIF
663
664!
665!--          As a side condition, routine flow_statistics require at least 14
666!--          vertical grid levels (they are used to store time-series data)
667!>           @todo   Remove this restriction
668             IF ( nz < 14 )  THEN
669                WRITE( message_string, * ) 'nz >= 14 is required'
670                CALL message( 'parin', 'PA0362', 1, 2, 0, 6, 0 )
671             ENDIF
672
673!
674!--          ATTENTION: in case of changes to the following statement please
675!--                  also check the allocate statement in routine rrd_global
676             ALLOCATE( pt_init(0:nz+1), q_init(0:nz+1), s_init(0:nz+1),        &
677                       ref_state(0:nz+1), sa_init(0:nz+1), ug(0:nz+1),         &
678                       u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1),             &
679                    hom(0:nz+1,2,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions),  &
680                    hom_sum(0:nz+1,pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0:statistic_regions) )
681
682             hom = 0.0_wp
683
684          ENDIF
685
686!
687!--       NAMELIST-file is not needed anymore
688          CALL close_file( 11 )
689
690       ENDIF
691#if defined( __parallel )
692       CALL MPI_BARRIER( MPI_COMM_WORLD, ierr )
693#endif
694    ENDDO
695
696    CALL location_message( 'reading NAMELIST parameters from PARIN', 'finished' )
697
698 END SUBROUTINE parin
Note: See TracBrowser for help on using the repository browser.