[1] | 1 | SUBROUTINE parin |
---|
| 2 | |
---|
| 3 | !------------------------------------------------------------------------------! |
---|
[257] | 4 | ! Current revisions: |
---|
[1] | 5 | ! ----------------- |
---|
[684] | 6 | ! |
---|
[601] | 7 | ! |
---|
| 8 | ! Former revisions: |
---|
| 9 | ! ----------------- |
---|
| 10 | ! $Id: parin.f90 684 2011-02-09 14:49:31Z raasch $ |
---|
| 11 | ! |
---|
[684] | 12 | ! 683 2011-02-09 14:25:15Z raasch |
---|
| 13 | ! +synchronous_exchange in d3par |
---|
| 14 | ! |
---|
[668] | 15 | ! 667 2010-12-23 12:06:00Z suehring/gryschka |
---|
| 16 | ! Steering parameter dissipation_control added in inipar. (commented out) |
---|
| 17 | ! |
---|
[623] | 18 | ! 622 2010-12-10 08:08:13Z raasch |
---|
| 19 | ! +collective_wait in inipar |
---|
| 20 | ! |
---|
[601] | 21 | ! 600 2010-11-24 16:10:51Z raasch |
---|
[600] | 22 | ! parameters moved from d3par to inipar: call_psolver_at_all_substeps, |
---|
| 23 | ! cfl_factor, cycle_mg, mg_cycles, mg_switch_to_pe0_level, ngsrb, nsor, |
---|
| 24 | ! omega_sor, prandtl_number, psolver, rayleigh_damping_factor, |
---|
| 25 | ! rayleigh_damping_height, residual_limit |
---|
[392] | 26 | ! |
---|
[581] | 27 | ! 580 2010-10-05 13:59:11Z heinze |
---|
| 28 | ! Renaming of ws_vertical_gradient to subs_vertical_gradient and |
---|
| 29 | ! ws_vertical_gradient_level to subs_vertical_gradient_level |
---|
| 30 | ! |
---|
[554] | 31 | ! 553 2010-09-01 14:09:06Z weinreis |
---|
| 32 | ! parameters for masked output are replaced by arrays |
---|
| 33 | ! |
---|
[494] | 34 | ! 493 2010-03-01 08:30:24Z raasch |
---|
| 35 | ! +netcdf_data_format in d3par, -netcdf_64bit, -netcdf_64bit_3d |
---|
| 36 | ! |
---|
[482] | 37 | ! 449 2010-02-02 11:23:59Z raasch |
---|
| 38 | ! +wall_humidityflux, wall_scalarflux |
---|
| 39 | ! +ws_vertical_gradient, ws_vertical_gradient_level |
---|
| 40 | ! |
---|
[449] | 41 | ! 410 2009-12-04 17:05:40Z letzel |
---|
| 42 | ! masked data output: + dt_domask, mask_01~20_x|y|z, mask_01~20_x|y|z_loop, |
---|
[493] | 43 | ! mask_scale_x|y|z, masks, skip_time_domask |
---|
[449] | 44 | ! |
---|
[392] | 45 | ! 291 2009-04-16 12:07:26Z raasch |
---|
[260] | 46 | ! +local_dvrserver_running in envpar |
---|
[257] | 47 | ! Output of messages replaced by message handling routine. |
---|
[240] | 48 | ! +canyon_height, canyon_width_x, canyon_width_y, canyon_wall_left, |
---|
[291] | 49 | ! canyon_wall_south, conserve_volume_flow_mode, coupling_start_time, |
---|
| 50 | ! dp_external, dp_level_b, dp_smooth, dpdxy, u_bulk, v_bulk in inipar |
---|
[256] | 51 | ! topography_grid_convention moved from userpar |
---|
[1] | 52 | ! |
---|
[198] | 53 | ! 197 2008-09-16 15:29:03Z raasch |
---|
| 54 | ! +cthf,leaf_surface_concentration, scalar_exchange_coefficient |
---|
| 55 | ! +inflow_damping_height, inflow_damping_width, recycling_width, |
---|
| 56 | ! turbulent_inflow in inipar, -skip_time_dosp in d3par, |
---|
| 57 | ! allocation of hom_sum moved from init_3d_model to here, |
---|
| 58 | ! npex, npey moved from inipar to d3par, setting of myid_char_14 removed, |
---|
| 59 | ! lad is allways allocated |
---|
| 60 | ! |
---|
[139] | 61 | ! 138 2007-11-28 10:03:58Z letzel |
---|
| 62 | ! +canopy_mode, drag_coefficient, lad_surface, lad_vertical_gradient, |
---|
| 63 | ! lad_vertical_gradient_level, pch_index, plant_canopy, |
---|
| 64 | ! +allocation of leaf area density field |
---|
| 65 | ! |
---|
[110] | 66 | ! 108 2007-08-24 15:10:38Z letzel |
---|
| 67 | ! +e_init, top_momentumflux_u|v in inipar, +dt_coupling in d3par |
---|
| 68 | ! |
---|
[98] | 69 | ! 95 2007-06-02 16:48:38Z raasch |
---|
| 70 | ! +bc_sa_t, bottom_salinityflux, ocean, sa_surface, sa_vertical_gradient, |
---|
| 71 | ! sa_vertical_gradient_level, top_salinityflux in inipar, |
---|
| 72 | ! sa_init is allocated |
---|
| 73 | ! |
---|
[90] | 74 | ! 87 2007-05-22 15:46:47Z raasch |
---|
| 75 | ! Size of hom increased by the maximum number of user-defined profiles, |
---|
| 76 | ! var_hom renamed pr_palm |
---|
| 77 | ! |
---|
[83] | 78 | ! 82 2007-04-16 15:40:52Z raasch |
---|
| 79 | ! +return_addres, return_username in envpar |
---|
| 80 | ! |
---|
[77] | 81 | ! 75 2007-03-22 09:54:05Z raasch |
---|
| 82 | ! +dt_max, netcdf_64bit_3d, precipitation_amount_interval in d3par, |
---|
| 83 | ! +loop_optimization, pt_reference in inipar, -data_output_ts, |
---|
| 84 | ! moisture renamed humidity |
---|
| 85 | ! |
---|
[39] | 86 | ! 20 2007-02-26 00:12:32Z raasch |
---|
| 87 | ! +top_heatflux, use_top_fluxes in inipar |
---|
| 88 | ! |
---|
[7] | 89 | ! 3 2007-02-13 11:30:58Z raasch |
---|
[3] | 90 | ! +netcdf_64bit_3d in d3par, |
---|
| 91 | ! RCS Log replace by Id keyword, revision history cleaned up |
---|
[2] | 92 | ! |
---|
[1] | 93 | ! Revision 1.57 2007/02/11 13:11:22 raasch |
---|
| 94 | ! Values of environment variables are now read from file ENVPAR instead of |
---|
| 95 | ! reading them with a system call, + NAMELIST envpar |
---|
| 96 | ! |
---|
| 97 | ! Revision 1.1 1997/07/24 11:22:50 raasch |
---|
| 98 | ! Initial revision |
---|
| 99 | ! |
---|
| 100 | ! |
---|
| 101 | ! Description: |
---|
| 102 | ! ------------ |
---|
| 103 | ! This subroutine reads variables controling the run from the NAMELIST files |
---|
| 104 | !------------------------------------------------------------------------------! |
---|
| 105 | |
---|
| 106 | USE arrays_3d |
---|
| 107 | USE averaging |
---|
| 108 | USE control_parameters |
---|
[260] | 109 | USE dvrp_variables |
---|
[1] | 110 | USE grid_variables |
---|
| 111 | USE indices |
---|
| 112 | USE model_1d |
---|
| 113 | USE pegrid |
---|
| 114 | USE profil_parameter |
---|
| 115 | USE statistics |
---|
| 116 | |
---|
| 117 | IMPLICIT NONE |
---|
| 118 | |
---|
[559] | 119 | INTEGER :: idum |
---|
[1] | 120 | |
---|
| 121 | |
---|
| 122 | NAMELIST /inipar/ adjust_mixing_length, alpha_surface, bc_e_b, bc_lr, & |
---|
| 123 | bc_ns, bc_p_b, bc_p_t, bc_pt_b, bc_pt_t, bc_q_b, & |
---|
[95] | 124 | bc_q_t,bc_s_b, bc_s_t, bc_sa_t, bc_uv_b, bc_uv_t, & |
---|
[138] | 125 | bottom_salinityflux, building_height, building_length_x, & |
---|
[240] | 126 | building_length_y, building_wall_left, building_wall_south, & |
---|
[600] | 127 | call_psolver_at_all_substeps, canopy_mode, canyon_height, & |
---|
| 128 | canyon_width_x, canyon_width_y, canyon_wall_left, & |
---|
| 129 | canyon_wall_south, cfl_factor, cloud_droplets, cloud_physics, & |
---|
[622] | 130 | collective_wait, conserve_volume_flow, conserve_volume_flow_mode, & |
---|
[291] | 131 | coupling_start_time, cthf, cut_spline_overshoot, & |
---|
[667] | 132 | cycle_mg, damp_level_1d, dissipation_1d, & !dissipation_control, & |
---|
| 133 | dp_external, dp_level_b, dp_smooth, dpdxy, drag_coefficient, & |
---|
| 134 | dt, dt_pr_1d, dt_run_control_1d, dx, dy, dz, dz_max, & |
---|
| 135 | dz_stretch_factor, dz_stretch_level, e_init, e_min, end_time_1d, & |
---|
| 136 | fft_method, galilei_transformation, grid_matching, humidity, & |
---|
[151] | 137 | inflow_damping_height, inflow_damping_width, & |
---|
[94] | 138 | inflow_disturbance_begin, inflow_disturbance_end, & |
---|
[138] | 139 | initializing_actions, km_constant, km_damp_max, lad_surface, & |
---|
| 140 | lad_vertical_gradient, lad_vertical_gradient_level, & |
---|
[153] | 141 | leaf_surface_concentration, long_filter_factor, & |
---|
[600] | 142 | loop_optimization, mg_cycles, mg_switch_to_pe0_level, & |
---|
| 143 | mixing_length_1d, momentum_advec, netcdf_precision, ngsrb, nsor, & |
---|
| 144 | nsor_ini, nx, ny, nz, ocean, omega, omega_sor, & |
---|
| 145 | outflow_damping_width, overshoot_limit_e, overshoot_limit_pt, & |
---|
| 146 | overshoot_limit_u, overshoot_limit_v, overshoot_limit_w, & |
---|
| 147 | passive_scalar, pch_index, phi, plant_canopy, prandtl_layer, & |
---|
| 148 | prandtl_number, precipitation, psolver, pt_reference, pt_surface, & |
---|
| 149 | pt_surface_initial_change, pt_vertical_gradient, & |
---|
[94] | 150 | pt_vertical_gradient_level, q_surface, q_surface_initial_change, & |
---|
| 151 | q_vertical_gradient, q_vertical_gradient_level, radiation, & |
---|
[600] | 152 | random_generator, random_heatflux, rayleigh_damping_factor, & |
---|
| 153 | rayleigh_damping_height, recycling_width, residual_limit, & |
---|
| 154 | rif_max, rif_min, roughness_length, sa_surface, & |
---|
| 155 | sa_vertical_gradient, sa_vertical_gradient_level, scalar_advec, & |
---|
[153] | 156 | scalar_exchange_coefficient, statistic_regions, & |
---|
[600] | 157 | subs_vertical_gradient, subs_vertical_gradient_level, & |
---|
[94] | 158 | surface_heatflux, surface_pressure, surface_scalarflux, & |
---|
| 159 | surface_waterflux, s_surface, s_surface_initial_change, & |
---|
[199] | 160 | s_vertical_gradient, s_vertical_gradient_level, timestep_scheme, & |
---|
[256] | 161 | topography, topography_grid_convention, top_heatflux, & |
---|
[600] | 162 | top_momentumflux_u, top_momentumflux_v, top_salinityflux, & |
---|
| 163 | turbulent_inflow, ug_surface, ug_vertical_gradient, & |
---|
| 164 | ug_vertical_gradient_level, u_bulk, ups_limit_e, ups_limit_pt, & |
---|
| 165 | ups_limit_u, ups_limit_v, ups_limit_w, use_surface_fluxes, & |
---|
| 166 | use_top_fluxes, use_ug_for_galilei_tr, use_upstream_for_tke, & |
---|
| 167 | vg_surface, vg_vertical_gradient, vg_vertical_gradient_level, & |
---|
| 168 | v_bulk, wall_adjustment, wall_heatflux, wall_humidityflux, & |
---|
| 169 | wall_scalarflux |
---|
[1] | 170 | |
---|
| 171 | |
---|
[600] | 172 | NAMELIST /d3par/ averaging_interval, averaging_interval_pr, & |
---|
| 173 | create_disturbances, cross_normalized_x, cross_normalized_y, & |
---|
| 174 | cross_profiles, cross_ts_uymax, cross_ts_uymin, cross_xtext, & |
---|
| 175 | data_output, data_output_format, data_output_masks, & |
---|
| 176 | data_output_pr, data_output_2d_on_each_pe, disturbance_amplitude, & |
---|
| 177 | disturbance_energy_limit, disturbance_level_b, & |
---|
| 178 | disturbance_level_t, do2d_at_begin, do3d_at_begin, do3d_compress, & |
---|
| 179 | do3d_comp_prec, dt, dt_averaging_input, dt_averaging_input_pr, & |
---|
| 180 | dt_coupling, dt_data_output, dt_data_output_av, dt_disturb, & |
---|
| 181 | dt_domask, dt_dopr, dt_dopr_listing, dt_dots, dt_do2d_xy, & |
---|
| 182 | dt_do2d_xz, dt_do2d_yz, dt_do3d, dt_max, dt_restart, & |
---|
| 183 | dt_run_control,end_time, force_print_header, mask_scale_x, & |
---|
| 184 | mask_scale_y, mask_scale_z, mask_x, mask_y, mask_z, mask_x_loop, & |
---|
| 185 | mask_y_loop, mask_z_loop, netcdf_data_format, & |
---|
| 186 | normalizing_region, npex, npey, nz_do3d, & |
---|
| 187 | precipitation_amount_interval, profile_columns, profile_rows, & |
---|
| 188 | restart_time, section_xy, section_xz, section_yz, & |
---|
| 189 | skip_time_data_output, skip_time_data_output_av, skip_time_dopr, & |
---|
| 190 | skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz, & |
---|
[683] | 191 | skip_time_do3d, skip_time_domask, synchronous_exchange, & |
---|
| 192 | termination_time_needed, use_prior_plot1d_parameters, z_max_do1d, & |
---|
| 193 | z_max_do1d_normalized, z_max_do2d |
---|
[1] | 194 | |
---|
| 195 | |
---|
[260] | 196 | NAMELIST /envpar/ host, local_dvrserver_running, maximum_cpu_time_allowed, & |
---|
| 197 | revision, return_addres, return_username, run_identifier, & |
---|
[667] | 198 | tasks_per_node, write_binary |
---|
[1] | 199 | |
---|
| 200 | ! |
---|
| 201 | !-- Open the NAMELIST-file which is send with this job |
---|
| 202 | CALL check_open( 11 ) |
---|
| 203 | |
---|
| 204 | ! |
---|
| 205 | !-- Read the control parameters for initialization. |
---|
| 206 | !-- The namelist "inipar" must be provided in the NAMELIST-file. If this is |
---|
| 207 | !-- not the case and the file contains - instead of "inipar" - any other |
---|
| 208 | !-- namelist, a read error is created on t3e and control is transferred |
---|
| 209 | !-- to the statement with label 10. Therefore, on t3e machines one can not |
---|
| 210 | !-- distinguish between errors produced by a wrong "inipar" namelist or |
---|
| 211 | !-- because this namelist is totally missing. |
---|
| 212 | READ ( 11, inipar, ERR=10, END=11 ) |
---|
| 213 | GOTO 12 |
---|
[274] | 214 | 10 message_string = 'errors in \$inipar &or no \$inipar-namelist ' // & |
---|
| 215 | 'found (CRAY-machines only)' |
---|
[257] | 216 | CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 ) |
---|
[1] | 217 | |
---|
[257] | 218 | 11 message_string = 'no \$inipar-namelist found' |
---|
| 219 | CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 ) |
---|
[559] | 220 | |
---|
[1] | 221 | ! |
---|
| 222 | !-- If required, read control parameters from restart file (produced by |
---|
[146] | 223 | !-- a prior run). All PEs are reading from file created by PE0 (see check_open) |
---|
[559] | 224 | 12 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN |
---|
| 225 | |
---|
| 226 | CALL read_var_list |
---|
[1] | 227 | ! |
---|
[146] | 228 | !-- The restart file will be reopened when reading the subdomain data |
---|
| 229 | CALL close_file( 13 ) |
---|
| 230 | |
---|
| 231 | ! |
---|
[1] | 232 | !-- Increment the run count |
---|
| 233 | runnr = runnr + 1 |
---|
| 234 | |
---|
[87] | 235 | ENDIF |
---|
| 236 | |
---|
[1] | 237 | ! |
---|
[87] | 238 | !-- Definition of names of areas used for computing statistics. They must |
---|
| 239 | !-- be defined at this place, because they are allowed to be redefined by |
---|
| 240 | !-- the user in user_parin. |
---|
| 241 | region = 'total domain' |
---|
| 242 | |
---|
| 243 | ! |
---|
| 244 | !-- Read runtime parameters given by the user for this run (namelist "d3par"). |
---|
| 245 | !-- The namelist "d3par" can be omitted. In that case, default values are |
---|
| 246 | !-- used for the parameters. |
---|
| 247 | READ ( 11, d3par, END=20 ) |
---|
| 248 | |
---|
| 249 | ! |
---|
| 250 | !-- Read control parameters for optionally used model software packages |
---|
| 251 | 20 CALL package_parin |
---|
| 252 | |
---|
| 253 | ! |
---|
| 254 | !-- Read user-defined variables |
---|
| 255 | CALL user_parin |
---|
| 256 | |
---|
| 257 | ! |
---|
| 258 | !-- Check in case of initial run, if the grid point numbers are well defined |
---|
| 259 | !-- and allocate some arrays which are already needed in init_pegrid or |
---|
| 260 | !-- check_parameters. During restart jobs, these arrays will be allocated |
---|
| 261 | !-- in read_var_list. All other arrays are allocated in init_3d_model. |
---|
| 262 | IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN |
---|
| 263 | |
---|
[1] | 264 | IF ( nx <= 0 ) THEN |
---|
[274] | 265 | WRITE( message_string, * ) 'no value or wrong value given for nx: ',& |
---|
| 266 | 'nx=', nx |
---|
[257] | 267 | CALL message( 'parin', 'PA0273', 1, 2, 0, 6, 0 ) |
---|
[1] | 268 | ENDIF |
---|
| 269 | IF ( ny <= 0 ) THEN |
---|
[274] | 270 | WRITE( message_string, * ) 'no value or wrong value given for ny: ',& |
---|
| 271 | 'ny=', ny |
---|
[257] | 272 | CALL message( 'parin', 'PA0274', 1, 2, 0, 6, 0 ) |
---|
[1] | 273 | ENDIF |
---|
| 274 | IF ( nz <= 0 ) THEN |
---|
[274] | 275 | WRITE( message_string, * ) 'no value or wrong value given for nz: ',& |
---|
| 276 | 'nz=', nz |
---|
[257] | 277 | CALL message( 'parin', 'PA0275', 1, 2, 0, 6, 0 ) |
---|
[1] | 278 | ENDIF |
---|
[147] | 279 | ! |
---|
| 280 | !-- ATTENTION: in case of changes to the following statement please also |
---|
| 281 | !-- check the allocate statement in routine read_var_list |
---|
[146] | 282 | ALLOCATE( lad(0:nz+1),pt_init(0:nz+1), q_init(0:nz+1), sa_init(0:nz+1), & |
---|
| 283 | ug(0:nz+1), u_init(0:nz+1), v_init(0:nz+1), vg(0:nz+1), & |
---|
| 284 | hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions), & |
---|
[145] | 285 | hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) ) |
---|
[667] | 286 | |
---|
[1] | 287 | hom = 0.0 |
---|
| 288 | |
---|
| 289 | ENDIF |
---|
| 290 | |
---|
| 291 | ! |
---|
| 292 | !-- NAMELIST-file is not needed anymore |
---|
| 293 | CALL close_file( 11 ) |
---|
| 294 | |
---|
| 295 | ! |
---|
| 296 | !-- Read values of environment variables (this NAMELIST file is generated by |
---|
| 297 | !-- mrun) |
---|
| 298 | OPEN ( 90, FILE='ENVPAR', STATUS='OLD', FORM='FORMATTED', ERR=30 ) |
---|
| 299 | READ ( 90, envpar, ERR=31, END=32 ) |
---|
| 300 | CLOSE ( 90 ) |
---|
| 301 | RETURN |
---|
| 302 | |
---|
[257] | 303 | 30 message_string = 'local file ENVPAR not found' // & |
---|
| 304 | '&some variables for steering may not be properly set' |
---|
| 305 | CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 ) |
---|
[1] | 306 | RETURN |
---|
| 307 | |
---|
[257] | 308 | 31 message_string = 'errors in local file ENVPAR' // & |
---|
| 309 | '&some variables for steering may not be properly set' |
---|
| 310 | CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 ) |
---|
[1] | 311 | RETURN |
---|
| 312 | |
---|
[257] | 313 | 32 message_string = 'no envpar-NAMELIST found in local file ENVPAR' // & |
---|
| 314 | '&some variables for steering may not be properly set' |
---|
| 315 | CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 ) |
---|
[1] | 316 | |
---|
| 317 | END SUBROUTINE parin |
---|