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

Last change on this file since 94 was 94, checked in by raasch, 15 years ago

preliminary uncomplete changes for ocean version

  • Property svn:keywords set to Id
File size: 10.4 KB
Line 
1 SUBROUTINE parin
2
3!------------------------------------------------------------------------------!
4! Actual revisions:
5! -----------------
6! +ocean, sa_surface, sa_vertical_gradient, sa_vertical_gradient_level,
7! top_salinity_flux in inipar,
8! sa_init is allocated
9!
10! Former revisions:
11! -----------------
12! $Id: parin.f90 94 2007-06-01 15:25:22Z raasch $
13!
14! 87 2007-05-22 15:46:47Z raasch
15! Size of hom increased by the maximum number of user-defined profiles,
16! var_hom renamed pr_palm
17!
18! 82 2007-04-16 15:40:52Z raasch
19! +return_addres, return_username in envpar
20!
21! 75 2007-03-22 09:54:05Z raasch
22! +dt_max, netcdf_64bit_3d, precipitation_amount_interval in d3par,
23! +loop_optimization, pt_reference in inipar, -data_output_ts,
24! moisture renamed humidity
25!
26! 20 2007-02-26 00:12:32Z raasch
27! +top_heatflux, use_top_fluxes in inipar
28!
29! -netcdf_64bit_3d
30!
31! 3 2007-02-13 11:30:58Z raasch
32! +netcdf_64bit_3d in d3par,
33! RCS Log replace by Id keyword, revision history cleaned up
34!
35! Revision 1.57  2007/02/11 13:11:22  raasch
36! Values of environment variables are now read from file ENVPAR instead of
37! reading them with a system call, + NAMELIST envpar
38!
39! Revision 1.1  1997/07/24 11:22:50  raasch
40! Initial revision
41!
42!
43! Description:
44! ------------
45! This subroutine reads variables controling the run from the NAMELIST files
46!------------------------------------------------------------------------------!
47
48    USE arrays_3d
49    USE averaging
50    USE control_parameters
51    USE grid_variables
52    USE indices
53    USE model_1d
54    USE pegrid
55    USE profil_parameter
56    USE statistics
57
58    IMPLICIT NONE
59
60    INTEGER ::  idum
61
62
63    NAMELIST /inipar/  adjust_mixing_length, alpha_surface, bc_e_b, bc_lr, &
64                       bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b, &
65             bc_q_t,bc_s_b, bc_s_t, bc_uv_b, bc_uv_t, building_height, &
66             building_length_x, building_length_y, building_wall_left, &
67             building_wall_south, cloud_droplets, cloud_physics, &
68             conserve_volume_flow, cut_spline_overshoot, damp_level_1d, &
69             dissipation_1d, dt, dt_pr_1d, dt_run_control_1d, dx, dy, dz, &
70             dz_max, dz_stretch_factor, dz_stretch_level, e_min, end_time_1d, &
71             fft_method, galilei_transformation, grid_matching, humidity, &
72             inflow_disturbance_begin, inflow_disturbance_end, &
73             initializing_actions, km_constant, km_damp_max, &
74             long_filter_factor, loop_optimization, mixing_length_1d, &
75             momentum_advec, netcdf_precision, npex, npey, nsor_ini, nx, ny, &
76             nz, ocean, omega, outflow_damping_width, overshoot_limit_e, &
77             overshoot_limit_pt, overshoot_limit_u, overshoot_limit_v, &
78             overshoot_limit_w, passive_scalar, phi, prandtl_layer, &
79             precipitation, pt_reference, pt_surface, &
80             pt_surface_initial_change, pt_vertical_gradient, &
81             pt_vertical_gradient_level, q_surface, q_surface_initial_change, &
82             q_vertical_gradient, q_vertical_gradient_level, radiation, &
83             random_generator, random_heatflux, rif_max, rif_min, &
84             roughness_length, sa_surface, sa_vertical_gradient, &
85             sa_vertical_gradient_level, scalar_advec, statistic_regions, &
86             surface_heatflux, surface_pressure, surface_scalarflux, &
87             surface_waterflux, s_surface, s_surface_initial_change, &
88             s_vertical_gradient, s_vertical_gradient_level, top_heatflux, &
89             top_salinity_flux, timestep_scheme, topography, ug_surface, &
90             ug_vertical_gradient, ug_vertical_gradient_level, ups_limit_e, &
91             ups_limit_pt, ups_limit_u, ups_limit_v, ups_limit_w, &
92             use_surface_fluxes, use_top_fluxes, use_ug_for_galilei_tr, &
93             use_upstream_for_tke, vg_surface, vg_vertical_gradient, &
94             vg_vertical_gradient_level, wall_adjustment, wall_heatflux
95
96
97    NAMELIST /d3par/   averaging_interval,  averaging_interval_pr, &
98                       call_psolver_at_all_substeps, cfl_factor, &
99                       create_disturbances, cross_normalized_x, &
100                       cross_normalized_y, cross_profiles, cross_ts_uymax, &
101                       cross_ts_uymin, cross_xtext, cycle_mg, data_output, &
102                       data_output_format, data_output_pr, &
103                       data_output_2d_on_each_pe, disturbance_amplitude, &
104                       disturbance_energy_limit, disturbance_level_b, &
105                       disturbance_level_t, do2d_at_begin, do3d_at_begin, &
106                       do3d_compress, do3d_comp_prec, dt, dt_averaging_input, &
107                       dt_averaging_input_pr, dt_data_output, &
108                       dt_data_output_av, dt_disturb, dt_dopr, &
109                       dt_dopr_listing, dt_dots, dt_do2d_xy, dt_do2d_xz, &
110                       dt_do2d_yz, dt_do3d, dt_max, dt_restart, dt_run_control,&
111                       end_time, force_print_header, mg_cycles, &
112                       mg_switch_to_pe0_level, netcdf_64bit, netcdf_64bit_3d, &
113                       ngsrb, normalizing_region, nsor, nz_do3d, omega_sor, &
114                       prandtl_number, precipitation_amount_interval, &
115                       profile_columns, profile_rows, psolver, &
116                       rayleigh_damping_factor, rayleigh_damping_height, &
117                       residual_limit, restart_time, section_xy, section_xz, &
118                       section_yz, skip_time_data_output, &
119                       skip_time_data_output_av, skip_time_dopr, &
120                       skip_time_dosp, skip_time_do2d_xy, skip_time_do2d_xz, &
121                       skip_time_do2d_yz, skip_time_do3d, &
122                       termination_time_needed, use_prior_plot1d_parameters, &
123                       z_max_do1d, z_max_do1d_normalized, z_max_do2d
124
125
126    NAMELIST /envpar/  host, maximum_cpu_time_allowed, revision, return_addres,&
127                       return_username, run_identifier, tasks_per_node, &
128                       write_binary
129
130
131#if defined( __parallel )
132!
133!-- Preliminary determination of processor-id which is needed here to open the
134!-- input files belonging to the corresponding processor and to produce
135!-- messages by PE0 only (myid and myid_char are later determined in
136!-- init_pegrid)
137    CALL MPI_COMM_RANK( comm_palm, myid, ierr )
138    WRITE (myid_char,'(''_'',I4.4)')  myid
139!
140!-- Since on IBM machines the process rank may be changed when the final
141!-- communicator is defined, save the preliminary processor-id for opening
142!-- the binary output file for restarts (unit 14), because otherwise
143!-- a mismatch occurs when reading this file in the next job
144    myid_char_14 = myid_char
145#endif
146
147!
148!-- Open the NAMELIST-file which is send with this job
149    CALL check_open( 11 )
150
151!
152!-- Read the control parameters for initialization.
153!-- The namelist "inipar" must be provided in the NAMELIST-file. If this is
154!-- not the case and the file contains - instead of "inipar" - any other
155!-- namelist, a read error is created on t3e and control is transferred
156!-- to the statement with label 10. Therefore, on t3e machines one can not
157!-- distinguish between errors produced by a wrong "inipar" namelist or
158!-- because this namelist is totally missing.
159    READ ( 11, inipar, ERR=10, END=11 )
160    GOTO 12
161 10 IF ( myid == 0 )  THEN
162       PRINT*, '+++ parin: errors in \$inipar'
163       PRINT*, '           or no \$inipar-namelist found (CRAY-machines only)' 
164    ENDIF
165    CALL local_stop
166 11 IF ( myid == 0 )  THEN
167       PRINT*, '+++ parin: no \$inipar-namelist found'
168    ENDIF
169    CALL local_stop
170
171!
172!-- If required, read control parameters from restart file (produced by
173!-- a prior run)
174 12 IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
175
176       CALL read_var_list
177!
178!--    Increment the run count
179       runnr = runnr + 1
180
181    ENDIF
182
183!
184!-- Definition of names of areas used for computing statistics. They must
185!-- be defined at this place, because they are allowed to be redefined by
186!-- the user in user_parin.
187    region = 'total domain'
188
189!
190!-- Read runtime parameters given by the user for this run (namelist "d3par").
191!-- The namelist "d3par" can be omitted. In that case, default values are
192!-- used for the parameters.
193    READ ( 11, d3par, END=20 )
194
195!
196!-- Read control parameters for optionally used model software packages
197 20 CALL package_parin
198
199!
200!-- Read user-defined variables
201    CALL user_parin
202
203!
204!-- Check in case of initial run, if the grid point numbers are well defined
205!-- and allocate some arrays which are already needed in init_pegrid or
206!-- check_parameters. During restart jobs, these arrays will be allocated
207!-- in read_var_list. All other arrays are allocated in init_3d_model.
208    IF ( TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
209
210       IF ( nx <= 0 )  THEN
211          IF ( myid == 0 )  THEN
212             PRINT*, '+++ parin: no value or wrong value given for nx: nx=', nx
213          ENDIF
214          CALL local_stop
215       ENDIF
216       IF ( ny <= 0 )  THEN
217          IF ( myid == 0 )  THEN
218             PRINT*, '+++ parin: no value or wrong value given for ny: ny=', ny
219          ENDIF
220          CALL local_stop
221       ENDIF
222       IF ( nz <= 0 )  THEN
223          IF ( myid == 0 )  THEN
224             PRINT*, '+++ parin: no value or wrong value given for nz: nz=', nz
225          ENDIF
226          CALL local_stop
227       ENDIF
228
229       ALLOCATE( ug(0:nz+1), vg(0:nz+1), &
230                 pt_init(0:nz+1), q_init(0:nz+1), sa_init(0:nz+1), &
231                 u_init(0:nz+1), v_init(0:nz+1),                   &
232                 hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions) )
233       hom = 0.0
234
235    ENDIF
236
237!
238!-- NAMELIST-file is not needed anymore
239    CALL close_file( 11 )
240
241!
242!-- Read values of environment variables (this NAMELIST file is generated by
243!-- mrun)
244    OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', ERR=30 )
245    READ ( 90, envpar, ERR=31, END=32 )
246    CLOSE ( 90 )
247    RETURN
248
249 30 IF ( myid == 0 )  THEN
250       PRINT*, '+++ parin: WARNING: local file ENVPAR not found'
251       PRINT*, '           some variables for steering may not be properly set'
252    ENDIF
253    RETURN
254
255 31 IF ( myid == 0 )  THEN
256       PRINT*, '+++ parin: WARNING: errors in local file ENVPAR'
257       PRINT*, '           some variables for steering may not be properly set'
258    ENDIF
259    RETURN
260
261 32 IF ( myid == 0 )  THEN
262       PRINT*, '+++ parin: WARNING: no envpar-NAMELIST found in local file ', &
263                           'ENVPAR'
264       PRINT*, '           some variables for steering may not be properly set'
265    ENDIF
266
267 END SUBROUTINE parin
Note: See TracBrowser for help on using the repository browser.