- Timestamp:
- Oct 5, 2007 5:13:14 AM (17 years ago)
- Location:
- palm/trunk
- Files:
-
- 38 edited
-
DOC/tec/technical_documentation.html (modified) (1 diff)
-
SOURCE/CURRENT_MODIFICATIONS (modified) (1 diff)
-
SOURCE/advec_particles.f90 (modified) (2 diffs)
-
SOURCE/advec_u_pw.f90 (modified) (1 diff)
-
SOURCE/advec_u_up.f90 (modified) (1 diff)
-
SOURCE/advec_v_pw.f90 (modified) (1 diff)
-
SOURCE/advec_v_up.f90 (modified) (1 diff)
-
SOURCE/boundary_conds.f90 (modified) (1 diff)
-
SOURCE/buoyancy.f90 (modified) (2 diffs)
-
SOURCE/check_for_restart.f90 (modified) (7 diffs)
-
SOURCE/check_open.f90 (modified) (1 diff)
-
SOURCE/check_parameters.f90 (modified) (2 diffs)
-
SOURCE/coriolis.f90 (modified) (1 diff)
-
SOURCE/diffusion_u.f90 (modified) (1 diff)
-
SOURCE/diffusion_v.f90 (modified) (1 diff)
-
SOURCE/flow_statistics.f90 (modified) (1 diff)
-
SOURCE/header.f90 (modified) (1 diff)
-
SOURCE/init_3d_model.f90 (modified) (1 diff)
-
SOURCE/init_particles.f90 (modified) (1 diff)
-
SOURCE/init_pegrid.f90 (modified) (2 diffs)
-
SOURCE/init_rankine.f90 (modified) (1 diff)
-
SOURCE/local_stop.f90 (modified) (2 diffs)
-
SOURCE/modules.f90 (modified) (2 diffs)
-
SOURCE/palm.f90 (modified) (2 diffs)
-
SOURCE/parin.f90 (modified) (1 diff)
-
SOURCE/prandtl_fluxes.f90 (modified) (1 diff)
-
SOURCE/pres.f90 (modified) (1 diff)
-
SOURCE/production_e.f90 (modified) (2 diffs)
-
SOURCE/prognostic_equations.f90 (modified) (1 diff)
-
SOURCE/read_3d_binary.f90 (modified) (1 diff)
-
SOURCE/read_var_list.f90 (modified) (1 diff)
-
SOURCE/surface_coupler.f90 (modified) (2 diffs)
-
SOURCE/swap_timelevel.f90 (modified) (1 diff)
-
SOURCE/time_integration.f90 (modified) (2 diffs)
-
SOURCE/timestep.f90 (modified) (2 diffs)
-
SOURCE/user_interface.f90 (modified) (1 diff)
-
SOURCE/write_3d_binary.f90 (modified) (1 diff)
-
SOURCE/write_var_list.f90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/tec/technical_documentation.html
r98 r110 2314 2314 density is calculated from the equation of state for seawater after 2315 2315 each timestep, using the algorithm proposed by Jackett et al. (2006, J. 2316 Atmos. Oceanic Technol., <span style="font-weight: bold;">23</span>, 1709-1728).<br>So far, only the initial hydrostatic pressure is entered into this equation.</li><li>z=0 (sea surface) is assumed at the model top (vertical grid index <span style="font-family: Courier New,Courier,monospace;">k=nzt</span> on the w-grid), with negative values of z indicating the depth.</li><li>Initial profiles are constructed (e.g. from <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient">pt_vertical_gradient</a> / <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient_level">pt_vertical_gradient_level</a>) starting from the sea surface, using surface values given by <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_surface">pt_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#sa_surface">sa_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#ug_surface">ug_surface</a>, and <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#vg_surface">vg_surface</a>.</li><li>Zero salinity flux is used as default boundary condition at the bottom of the sea.</li><li>If switched on, random perturbations are by default imposed to the upper model domain from zu(nzt*2/3) to zu(nzt-3).</li></ul>Relevant new inipar-parameters to be exclusively used for steering ocean runs are <span style="font-family: Courier New,Courier,monospace;">bc_sa_t</span>, <span style="font-family: Courier New,Courier,monospace;">bottom_salinityflux</span>, <span style="font-family: Courier New,Courier,monospace;">sa_surface</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient_level</span>, and <span style="font-family: Courier New,Courier,monospace;">top_salinityflux</span>.<br><br>Salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>) and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) are included as new 2d/3d output quantities. Vertical profiles of salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>), salinity fluxes (<span style="font-family: Courier New,Courier,monospace;">w"sa"</span>, <span style="font-family: Courier New,Courier,monospace;">w*sa*</span>, <span style="font-family: Courier New,Courier,monospace;">wsa</span>), and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) can also be output.<span style="font-family: Courier New,Courier,monospace;"></span></td><td style="vertical-align: top;">advec_s_bc, average_3d_data,<br>boundary_conds, buoyancy,<br>check_parameters, data_output_2d,<br>data_output_3d, diffusion_e,<br>diffusivities, flow_statistics, header,<br>init_grid, init_3d_model, modules,<br>netcdf, parin, production_e,<br>prognostic_equations,<br>read_3d_binary, read_var_list,<br>sum_up_3d_data, swap_timelevel,<br>time_integration, user_interface,<br>write_var_list, write_3d_binary<br><br><span style="font-weight: bold;">new:</span><br>eqn_state_seawater, init_ocean</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Inipar-parameter <span style="font-family: Courier New,Courier,monospace;">use_pt_reference</span> renamed <span style="font-family: Courier New,Courier,monospace;">use_reference.</span><br>Internal variable <span style="font-family: Courier New,Courier,monospace;">hydro_press</span> renamed <span style="font-family: Courier New,Courier,monospace;">hyp</span>, routine <span style="font-family: Courier New,Courier,monospace;">calc_mean_pt_profile</span> renamed <span style="font-family: Courier New,Courier,monospace;">calc_mean_profile</span>.<br><br>The format of the <span style="font-family: Courier New,Courier,monospace;">RUN_CONTROL</span> file has been adjusted for ocean runs.<br><br>Defaults of <span style="font-family: Courier New,Courier,monospace;">.._vertical_gradient_levels</span> changed from <span style="font-style: italic;">-1.0</span> to <span style="font-style: italic;">-9999999.9</span>.</td><td style="vertical-align: top;">advec_particles, buoyancy, <br>calc_liquid_water_content, <br>check_parameters, diffusion_e,<br>diffusivities, header,<br>init_cloud_physics, modules,<br>production_e, prognostic_equations,<br>run_control</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: height above topography instead of height above level k=0 is used for calculating the mixing length (<span style="font-family: Courier New,Courier,monospace;">diffusion_e</span> and <span style="font-family: Courier New,Courier,monospace;">diffusivities</span>).<br><br>Bugfix: error in boundary condition for TKE removed (<span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>).</td><td style="vertical-align: top;">advec_s_bc, diffusion_e, diffusivities</td></tr> 2316 Atmos. Oceanic Technol., <span style="font-weight: bold;">23</span>, 1709-1728).<br>So far, only the initial hydrostatic pressure is entered into this equation.</li><li>z=0 (sea surface) is assumed at the model top (vertical grid index <span style="font-family: Courier New,Courier,monospace;">k=nzt</span> on the w-grid), with negative values of z indicating the depth.</li><li>Initial profiles are constructed (e.g. from <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient">pt_vertical_gradient</a> / <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_vertical_gradient_level">pt_vertical_gradient_level</a>) starting from the sea surface, using surface values given by <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#pt_surface">pt_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#sa_surface">sa_surface</a>, <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#ug_surface">ug_surface</a>, and <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#vg_surface">vg_surface</a>.</li><li>Zero salinity flux is used as default boundary condition at the bottom of the sea.</li><li>If switched on, random perturbations are by default imposed to the upper model domain from zu(nzt*2/3) to zu(nzt-3).</li></ul>Relevant new inipar-parameters to be exclusively used for steering ocean runs are <span style="font-family: Courier New,Courier,monospace;">bc_sa_t</span>, <span style="font-family: Courier New,Courier,monospace;">bottom_salinityflux</span>, <span style="font-family: Courier New,Courier,monospace;">sa_surface</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient</span>, <span style="font-family: Courier New,Courier,monospace;">sa_vertical_gradient_level</span>, and <span style="font-family: Courier New,Courier,monospace;">top_salinityflux</span>.<br><br>Salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>) and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) are included as new 2d/3d output quantities. Vertical profiles of salinity (<span style="font-family: Courier New,Courier,monospace;">sa</span>), salinity fluxes (<span style="font-family: Courier New,Courier,monospace;">w"sa"</span>, <span style="font-family: Courier New,Courier,monospace;">w*sa*</span>, <span style="font-family: Courier New,Courier,monospace;">wsa</span>), and potential density (<span style="font-family: Courier New,Courier,monospace;">rho</span>) can also be output.<span style="font-family: Courier New,Courier,monospace;"></span></td><td style="vertical-align: top;">advec_s_bc, average_3d_data,<br>boundary_conds, buoyancy,<br>check_parameters, data_output_2d,<br>data_output_3d, diffusion_e,<br>diffusivities, flow_statistics, header,<br>init_grid, init_3d_model, modules,<br>netcdf, parin, production_e,<br>prognostic_equations,<br>read_3d_binary, read_var_list,<br>sum_up_3d_data, swap_timelevel,<br>time_integration, user_interface,<br>write_var_list, write_3d_binary<br><br><span style="font-weight: bold;">new:</span><br>eqn_state_seawater, init_ocean</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Inipar-parameter <span style="font-family: Courier New,Courier,monospace;">use_pt_reference</span> renamed <span style="font-family: Courier New,Courier,monospace;">use_reference.</span><br>Internal variable <span style="font-family: Courier New,Courier,monospace;">hydro_press</span> renamed <span style="font-family: Courier New,Courier,monospace;">hyp</span>, routine <span style="font-family: Courier New,Courier,monospace;">calc_mean_pt_profile</span> renamed <span style="font-family: Courier New,Courier,monospace;">calc_mean_profile</span>.<br><br>The format of the <span style="font-family: Courier New,Courier,monospace;">RUN_CONTROL</span> file has been adjusted for ocean runs.<br><br>Defaults of <span style="font-family: Courier New,Courier,monospace;">.._vertical_gradient_levels</span> changed from <span style="font-style: italic;">-1.0</span> to <span style="font-style: italic;">-9999999.9</span>.</td><td style="vertical-align: top;">advec_particles, buoyancy, <br>calc_liquid_water_content, <br>check_parameters, diffusion_e,<br>diffusivities, header,<br>init_cloud_physics, modules,<br>production_e, prognostic_equations,<br>run_control</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: height above topography instead of height above level k=0 is used for calculating the mixing length (<span style="font-family: Courier New,Courier,monospace;">diffusion_e</span> and <span style="font-family: Courier New,Courier,monospace;">diffusivities</span>).<br><br>Bugfix: error in boundary condition for TKE removed (<span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>).</td><td style="vertical-align: top;">advec_s_bc, diffusion_e, diffusivities</td></tr><tr><td style="vertical-align: top;">05/10/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.4</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">This version allows runs for a coupled atmosphere-ocean LES.<br><br>Coupling frequency is controlled by new d3par-parameter <span style="font-family: Courier New,Courier,monospace;">dt_coupling</span>, the coupling mode ('<span style="font-style: italic;">atmosphere_to_ocean</span>' or '<span style="font-style: italic;">ocean_to_atmosphere</span>') for the respective processes is read from environment variable <span style="font-family: Courier New,Courier,monospace;">coupling_mode</span>, which is set by the mpiexec-command in <span style="font-weight: bold;">mrun</span>. Communication between the two models is handled by the intercommunicator <span style="font-family: Courier New,Courier,monospace;">comm_inter</span>. Local files opened by the ocean model get the additional suffic "<span style="font-family: Courier New,Courier,monospace;">_O</span>". 2317 Saturation at k=nzb is assumed for the atmosphere coupled to ocean. 2318 Usage of the coupled model is described in the new section 3.8.<br><br>A momentum flux can be set as top boundary condition using the new inipar parameters <span style="font-family: Courier New,Courier,monospace;">top_momentumflux_u|v</span>.<br><br>Non-cyclic boundary conditions can be used along all horizontal directions.<br><br>Quantities w*p* and w"e can be output as vertical profiles.<br><br>Initial profiles are reset to constant profiles in case that <span style="font-family: Courier New,Courier,monospace;">initializing_actions</span> /= '<span style="font-style: italic;">set_constant_profiles</span>'. (init_rankine)<br><br>Optionally calculate km and kh from initial TKE e_init.<br><br><span style="font-weight: bold;">Attention:</span> <br> 2319 The format of the testresult-file (run control output) has changed 2320 due to this modification.</td><td style="vertical-align: top;">boundary_conds, 2321 check_open, check_parameters, diffusion_u, diffusion_v, 2322 flow_statistics, header, init_pegrid, init_rankine, init_3d_model, 2323 local_stop, modules, palm, parin, prandtl_fluxes, pres, 2324 prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, 2325 timestep, time_integration, write_var_list, write_3d_binary<br><br><span style="font-weight: bold;">new:</span><br>surface_coupler</td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Remaining variables <span style="font-family: Courier New,Courier,monospace;">iran</span> changed to <span style="font-family: Courier New,Courier,monospace;">iran_part</span> (advec_particles, init_particles).<br><br>In case that the presure solver is not called for every Runge-Kutta substep (<span style="font-family: Courier New,Courier,monospace;">call_psolver_at_all_substeps</span> = <span style="font-style: italic;">.F.</span>), 2326 it is called after the first substep instead of the last. In that case, 2327 random perturbations are also added to the velocity field after the 2328 first substep.<br><br>Initialization of <span style="font-family: Courier New,Courier,monospace;">km</span>, <span style="font-family: Courier New,Courier,monospace;">kh</span> = 0.00001 for <span style="font-family: Courier New,Courier,monospace;">ocean</span> = <span style="font-style: italic;">.T.</span> (for <span style="font-family: Courier New,Courier,monospace;">ocean</span> = <span style="font-style: italic;">.F.</span> it remains 0.01).<br><br>Allow <span style="font-family: Courier New,Courier,monospace;">data_output_pr</span> = '<span style="font-style: italic;">q</span>', '<span style="font-style: italic;">wq</span>', '<span style="font-style: italic;">w"q"</span>', '<span style="font-style: italic;">w*q*</span>' for <span style="font-family: Courier New,Courier,monospace;">humidity</span> = <span style="font-style: italic;">.T.</span> (instead of <span style="font-family: Courier New,Courier,monospace;">cloud_physics</span> = <span style="font-style: italic;">.T.</span>).</td><td style="vertical-align: top;">advec_particles, check_for_restart, check_parameters, init_particles, init_3d_model, time_integration<br></td></tr><tr><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;"> </td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugs from code parts for non-cyclic boundary conditions are removed:<br>loops for u and v are starting from index <span style="font-family: Courier New,Courier,monospace;">nxlu</span>, <span style="font-family: Courier New,Courier,monospace;">nysv</span>, 2329 respectively. The radiation boundary condition is used for every 2330 Runge-Kutta substep. Velocity phase speeds for the radiation boundary 2331 conditions are calculated for the first Runge-Kutta substep only and 2332 reused for the further substeps. New arrays <span style="font-family: Courier New,Courier,monospace;">c_u</span>, <span style="font-family: Courier New,Courier,monospace;">c_v</span>, and <span style="font-family: Courier New,Courier,monospace;">c_w</span> 2333 are defined for this purpose. Several index errors are removed from the 2334 radiation boundary condition code parts. Upper bounds for calculating <span style="font-family: Courier New,Courier,monospace;">u_0</span> and <span style="font-family: Courier New,Courier,monospace;">v_0</span> (in production_e) are <span style="font-style: italic;">nxr+1</span> and <span style="font-style: italic;">nyn+1</span> because otherwise these values are not available in case of non-cyclic boundary conditions.<br><br>+<span style="font-family: Courier New,Courier,monospace;">dots_num_palm</span> in module user, +module netcdf_control in user_init (both in user_interface)<br><br>Bugfix: wrong sign removed from the buoyancy production term in the case <span style="font-family: Courier New,Courier,monospace;">use_reference</span> = <span style="font-style: italic;">.T.</span> (production_e)<br><br>Bugfix: Error message concerning output of particle concentration (pc) modified (check_parameters).<br><br>Bugfix: Rayleigh damping for ocean fixed.</td><td style="vertical-align: top;">advec_u_pw, 2335 advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy, 2336 check_parameters, coriolis, diffusion_u, diffusion_v, header, 2337 init_pegrid, init_3d_model, modules, production_e, 2338 prognostic_equations, user_interface<br></td></tr> 2317 2339 </tbody> 2318 2340 </table> <b><blink>Attention:</blink></b> -
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r108 r110 1 1 New: 2 2 --- 3 Allows runs for a coupled atmosphere-ocean LES,4 coupling frequency is controlled by new d3par-parameter dt_coupling,5 the coupling mode (atmosphere_to_ocean or ocean_to_atmosphere) for the6 respective processes is read from environment variable coupling_mode,7 which is set by the mpiexec-command,8 communication between the two models is done using the intercommunicator9 comm_inter,10 local files opened by the ocean model get the additional suffic "_O".11 Assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean.12 13 A momentum flux can be set as top boundary condition using the new14 inipar parameter top_momentumflux_u|v.15 16 Non-cyclic boundary conditions can be used along all horizontal directions.17 18 Quantities w*p* and w"e can be output as vertical profiles.19 20 Initial profiles are reset to constant profiles in case that initializing_actions /= 'set_constant_profiles'. (init_rankine)21 22 Optionally calculate km and kh from initial TKE e_init.23 24 boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_rankine, init_3d_model, modules, palm, parin, pres, prandtl_fluxes, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary25 26 New:27 surface_coupler28 3 29 4 30 5 Changed: 31 6 ------- 32 Remaining variables iran changed to iran_part (advec_particles, init_particles).33 34 In case that the presure solver is not called for every Runge-Kutta substep35 (call_psolver_at_all_substeps = .F.), it is called after the first substep36 instead of the last. In that case, random perturbations are also added to the37 velocity field after the first substep.38 39 Initialization of km,kh = 0.00001 for ocean = .T. (for ocean = .F. it remains 0.01).40 41 Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T. (instead of cloud_physics = .T.).42 43 Modifications to terminate coupled runs.44 45 advec_particles, check_for_restart, check_parameters, init_particles, init_3d_model, local_stop, timestep, time_integration46 7 47 8 48 9 Errors: 49 10 ------ 50 Bugs from code parts for non-cyclic boundary conditions are removed: loops for51 u and v are starting from index nxlu, nysv, respectively. The radiation boundary52 condition is used for every Runge-Kutta substep. Velocity phase speeds for53 the radiation boundary conditions are calculated for the first Runge-Kutta54 substep only and reused for the further substeps. New arrays c_u, c_v, and c_w55 are defined for this purpose. Several index errors are removed from the56 radiation boundary condition code parts. Upper bounds for calculating57 u_0 and v_0 (in production_e) are nxr+1 and nyn+1 because otherwise these58 values are not available in case of non-cyclic boundary conditions.59 60 +dots_num_palm in module user, +module netcdf_control in user_init (both in user_interface)61 62 Bugfix: wrong sign removed from the buoyancy production term in the case use_reference = .T. (production_e)63 64 Bugfix: Error message concerning output of particle concentration (pc) modified (check_parameters).65 66 Bugfix: Rayleigh damping for ocean fixed.67 68 advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, boundary_conds, buoyancy, check_parameters, coriolis, diffusion_u, diffusion_v, header, init_pegrid, init_3d_model, modules, production_e, prognostic_equations, user_interface -
palm/trunk/SOURCE/advec_particles.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! remaining variables iran changed to iran_part7 6 ! TEST: PRINT statements on unit 9 (commented out) 8 7 ! … … 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! remaining variables iran changed to iran_part 12 14 ! 13 15 ! 95 2007-06-02 16:48:38Z raasch -
palm/trunk/SOURCE/advec_u_pw.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! i loop is starting from nxlu (needed for non-cyclic boundary conditions)6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! i loop is starting from nxlu (needed for non-cyclic boundary conditions) 11 14 ! 12 15 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/advec_u_up.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! i loop is starting from nxlu (needed for non-cyclic boundary conditions)6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! i loop is starting from nxlu (needed for non-cyclic boundary conditions) 11 14 ! 12 15 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/advec_v_pw.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! j loop is starting from nysv (needed for non-cyclic boundary conditions)6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! j loop is starting from nysv (needed for non-cyclic boundary conditions) 11 14 ! 12 15 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/advec_v_up.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! j loop is starting from nysv (needed for non-cyclic boundary conditions)6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! j loop is starting from nysv (needed for non-cyclic boundary conditions) 11 14 ! 12 15 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/boundary_conds.f90
r107 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 7 ! 8 ! Former revisions: 9 ! ----------------- 10 ! $Id$ 11 ! 12 ! 107 2007-08-17 13:54:45Z raasch 6 13 ! Boundary conditions for temperature adjusted for coupled runs, 7 14 ! bugfixes for the radiation boundary conditions at the outflow: radiation 8 15 ! conditions are used for every substep, phase speeds are calculated for the 9 16 ! first Runge-Kutta substep only and then reused, several index values changed 10 !11 !12 ! Former revisions:13 ! -----------------14 ! $Id$15 17 ! 16 18 ! 95 2007-06-02 16:48:38Z raasch -
palm/trunk/SOURCE/buoyancy.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! i loop for u-component (sloping surface) is starting from nxlu (needed for 7 ! non-cyclic boundary conditions) 6 ! 8 7 ! 9 8 ! Former revisions: … … 11 10 ! $Id$ 12 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! i loop for u-component (sloping surface) is starting from nxlu (needed for 14 ! non-cyclic boundary conditions) 15 ! 13 16 ! 97 2007-06-21 08:23:15Z raasch 14 17 ! Routine reneralized to be used with temperature AND density: -
palm/trunk/SOURCE/check_for_restart.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! modifications to terminate coupled runs7 6 ! 8 7 ! … … 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 13 ! modifications to terminate coupled runs 14 ! 12 15 ! RCS Log replace by Id keyword, revision history cleaned up 13 16 ! … … 61 64 !-- Output that job will be terminated 62 65 IF ( terminate_run .AND. myid == 0 ) THEN 63 PRINT*, '*** WARNING: run will be terminated because it is running out of job cpu', &64 ' limit'66 PRINT*, '*** WARNING: run will be terminated because it is running out', & 67 ' of job cpu limit' 65 68 PRINT*, ' remaining time: ', remaining_time, ' s' 66 PRINT*, ' termination time needed:', termination_time_needed, ' s' 69 PRINT*, ' termination time needed:', termination_time_needed,& 70 ' s' 67 71 ENDIF 68 72 … … 72 76 !-- informed of another termination reason (terminate_coupled > 0) before, 73 77 !-- or vice versa (terminate_coupled_remote > 0). 74 IF ( terminate_run .AND. TRIM( coupling_mode ) /= 'uncoupled' &75 .AND. terminate_coupled == 0 &76 .AND. terminate_coupled_remote == 0 ) THEN 78 IF ( terminate_run .AND. TRIM( coupling_mode ) /= 'uncoupled' .AND. & 79 terminate_coupled == 0 .AND. terminate_coupled_remote == 0 ) THEN 80 77 81 terminate_coupled = 3 78 CALL MPI_SENDRECV( & 79 terminate_coupled, 1, MPI_INTEGER, myid, 0, & 80 terminate_coupled_remote, 1, MPI_INTEGER, myid, 0, & 81 comm_inter, status, ierr ) 82 CALL MPI_SENDRECV( terminate_coupled, 1, MPI_INTEGER, myid, 0, & 83 terminate_coupled_remote, 1, MPI_INTEGER, myid, 0, & 84 comm_inter, status, ierr ) 82 85 ENDIF 83 86 … … 101 104 102 105 IF ( myid == 0 ) THEN 103 PRINT*, '*** INFORMATIVE: run will be terminated due to user settings of' 106 PRINT*, '*** INFORMATIVE: run will be terminated due to user ', & 107 'settings of' 104 108 PRINT*, ' restart_time / dt_restart' 105 109 PRINT*, ' new restart time is: ', time_restart, ' s' … … 110 114 !-- informed of another termination reason (terminate_coupled > 0) before, 111 115 !-- or vice versa (terminate_coupled_remote > 0). 112 IF ( coupling_mode /= 'uncoupled' .AND. terminate_coupled == 0 & 113 .AND. terminate_coupled_remote == 0) THEN 116 IF ( coupling_mode /= 'uncoupled' .AND. terminate_coupled == 0 .AND. & 117 terminate_coupled_remote == 0) THEN 118 114 119 IF ( dt_restart /= 9999999.9 ) THEN 115 120 terminate_coupled = 4 … … 117 122 terminate_coupled = 5 118 123 ENDIF 119 CALL MPI_SENDRECV( &120 terminate_coupled, 1, MPI_INTEGER, myid, 0, &121 terminate_coupled_remote, 1, MPI_INTEGER, myid, 0, &122 comm_inter, status, ierr )124 CALL MPI_SENDRECV( & 125 terminate_coupled, 1, MPI_INTEGER, myid, 0, & 126 terminate_coupled_remote, 1, MPI_INTEGER, myid, 0, & 127 comm_inter, status, ierr ) 123 128 ENDIF 124 129 ELSE -
palm/trunk/SOURCE/check_open.f90
r105 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Different filenames are used in case of a coupled simulation, 7 ! coupling_char added to all relevant filenames 6 ! 8 7 ! 9 8 ! Former revisions: 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 105 2007-08-08 07:12:55Z raasch 13 ! Different filenames are used in case of a coupled simulation, 14 ! coupling_char added to all relevant filenames 12 15 ! 13 16 ! 82 2007-04-16 15:40:52Z raasch -
palm/trunk/SOURCE/check_parameters.f90
r109 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 7 ! 8 ! Former revisions: 9 ! ----------------- 10 ! $Id$ 11 ! 12 ! 109 2007-08-28 15:26:47Z letzel 6 13 ! Check coupling_mode and set default (obligatory) values (like boundary 7 14 ! conditions for temperature and fluxes) in case of coupled runs. … … 14 21 ! Rayleigh damping for ocean fixed. 15 22 ! Check and, if necessary, set default value for dt_coupling 16 !17 ! Former revisions:18 ! -----------------19 ! $Id$20 23 ! 21 24 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/coriolis.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! loops for u and v are starting from index nxlu, nysv, respectively (needed 7 ! for non-cyclic boundary conditions) 6 ! 8 7 ! 9 8 ! Former revisions: 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! loops for u and v are starting from index nxlu, nysv, respectively (needed 14 ! for non-cyclic boundary conditions) 12 15 ! 13 16 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/diffusion_u.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Momentumflux at top (uswst) included as boundary condition, 7 ! i loop is starting from nxlu (needed for non-cyclic boundary conditions) 6 ! 8 7 ! 9 8 ! Former revisions: 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! Momentumflux at top (uswst) included as boundary condition, 14 ! i loop is starting from nxlu (needed for non-cyclic boundary conditions) 12 15 ! 13 16 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/diffusion_v.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Momentumflux at top (vswst) included as boundary condition, 7 ! j loop is starting from nysv (needed for non-cyclic boundary conditions) 6 ! 8 7 ! 9 8 ! Former revisions: 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! Momentumflux at top (vswst) included as boundary condition, 14 ! j loop is starting from nysv (needed for non-cyclic boundary conditions) 12 15 ! 13 16 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/flow_statistics.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Prescribed momentum fluxes at the top surface are used, 7 ! profiles for w*p* and w"e are calculated 6 ! 8 7 ! 9 8 ! Former revisions: 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! Prescribed momentum fluxes at the top surface are used, 14 ! profiles for w*p* and w"e are calculated 12 15 ! 13 16 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/header.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 7 ! 8 ! Former revisions: 9 ! ----------------- 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 6 13 ! Output of informations for coupled model runs (boundary conditions etc.) 7 14 ! + output of momentumfluxes at the top boundary 8 15 ! Rayleigh damping for ocean, e_init 9 !10 ! Former revisions:11 ! -----------------12 ! $Id$13 16 ! 14 17 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/init_3d_model.f90
r108 r110 7 7 ! Actual revisions: 8 8 ! ----------------- 9 ! 10 ! 11 ! Former revisions: 12 ! ----------------- 13 ! $Id$ 14 ! 15 ! 108 2007-08-24 15:10:38Z letzel 9 16 ! Flux initialization in case of coupled runs, +momentum fluxes at top boundary, 10 17 ! +arrays for phase speed c_u, c_v, c_w, indices for u|v|w_m_l|r changed 11 18 ! +qswst_remote in case of atmosphere model with humidity coupled to ocean 12 ! Rayleigh damping for ocean 13 ! optionally calculate km and kh from initial TKE e_init 14 ! 15 ! Former revisions: 16 ! ----------------- 17 ! $Id$ 19 ! Rayleigh damping for ocean, optionally calculate km and kh from initial 20 ! TKE e_init 18 21 ! 19 22 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/init_particles.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! variable iran replaced by iran_part6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! variable iran replaced by iran_part 11 14 ! 12 15 ! 82 2007-04-16 15:40:52Z raasch -
palm/trunk/SOURCE/init_pegrid.f90
r108 r110 5 5 ! ----------------- 6 6 ! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values 7 ! Intercommunicator (comm_inter) and derived data type (type_xy) for8 ! coupled model runs created, assign coupling_mode_remote9 ! indices nxlu and nysv are calculated (needed for non-cyclic boundary10 ! conditions)11 7 ! 12 8 ! Former revisions: 13 9 ! ----------------- 14 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 13 ! Intercommunicator (comm_inter) and derived data type (type_xy) for 14 ! coupled model runs created, assign coupling_mode_remote, 15 ! indices nxlu and nysv are calculated (needed for non-cyclic boundary 16 ! conditions) 15 17 ! 16 18 ! 82 2007-04-16 15:40:52Z raasch … … 894 896 #endif 895 897 ! 896 !-- At the outflow, u or v, respectively, have to be calculated for one grid897 !-- point less.898 !-- At the outflow, u or v, respectively, have to be calculated for one more 899 !-- grid point. 898 900 IF ( outflow_l ) THEN 899 901 nxlu = nxl + 1 -
palm/trunk/SOURCE/init_rankine.f90
r107 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Initial profiles are reset to constant profiles6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 107 2007-08-17 13:54:45Z raasch 13 ! Initial profiles are reset to constant profiles 11 14 ! 12 15 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/local_stop.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! modifications to terminate coupled runs6 ! 7 7 ! 8 8 ! … … 10 10 ! ----------------- 11 11 ! $Id$ 12 ! 13 ! 108 2007-08-24 15:10:38Z letzel 14 ! modifications to terminate coupled runs 15 ! 12 16 ! RCS Log replace by Id keyword, revision history cleaned up 13 17 ! -
palm/trunk/SOURCE/modules.f90
r108 r110 5 5 ! Actual revisions: 6 6 ! ----------------- 7 ! 8 ! 9 ! Former revisions: 10 ! ----------------- 11 ! $Id$ 12 ! 13 ! 108 2007-08-24 15:10:38Z letzel 7 14 ! +comm_inter, constant_top_momentumflux, coupling_char, coupling_mode, 8 15 ! coupling_mode_remote, c_u, c_v, c_w, dt_coupling, e_init, humidity_remote, … … 10 17 ! terminate_coupled_remote, time_coupling, top_momentumflux_u|v, type_xy, 11 18 ! uswst*, vswst* 12 !13 ! Former revisions:14 ! -----------------15 ! $Id$16 19 ! 17 20 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/palm.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Get coupling mode from environment variable 7 ! Change localtion of debug output 6 ! 8 7 ! 9 8 ! Former revisions: 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 13 ! Get coupling mode from environment variable, change location of debug output 12 14 ! 13 15 ! 75 2007-03-22 09:54:05Z raasch … … 64 66 INTEGER :: i, run_description_header_i(80) 65 67 66 version = 'PALM 3. 3a'68 version = 'PALM 3.4' 67 69 68 70 #if defined( __parallel ) -
palm/trunk/SOURCE/parin.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +e_init, top_momentumflux_u|v in inipar, +dt_coupling in d3par6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 13 ! +e_init, top_momentumflux_u|v in inipar, +dt_coupling in d3par 11 14 ! 12 15 ! 95 2007-06-02 16:48:38Z raasch -
palm/trunk/SOURCE/prandtl_fluxes.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 13 ! assume saturation at k=nzb_s_inner(j,i) for atmosphere coupled to ocean 11 14 ! 12 15 ! 75 2007-03-22 09:54:05Z raasch -
palm/trunk/SOURCE/pres.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Volume flow conservation added for the remaining three outflow boundaries6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 13 ! Volume flow conservation added for the remaining three outflow boundaries 11 14 ! 12 15 ! 85 2007-05-11 09:35:14Z raasch -
palm/trunk/SOURCE/production_e.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 7 ! 8 ! Former revisions: 9 ! ----------------- 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 6 13 ! Bugfix: wrong sign removed from the buoyancy production term in the case 7 14 ! use_reference = .T., … … 9 16 ! not available in case of non-cyclic boundary conditions) 10 17 ! Bugfix for ocean density flux at bottom 11 !12 ! Former revisions:13 ! -----------------14 ! $Id$15 18 ! 16 19 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/prognostic_equations.f90
r106 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 7 ! 8 ! Former revisions: 9 ! ----------------- 10 ! $Id$ 11 ! 12 ! 106 2007-08-16 14:30:26Z raasch 6 13 ! +uswst, vswst as arguments in calls of diffusion_u|v, 7 14 ! loops for u and v are starting from index nxlu, nysv, respectively (needed 8 15 ! for non-cyclic boundary conditions) 9 !10 ! Former revisions:11 ! -----------------12 ! $Id$13 16 ! 14 17 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/read_3d_binary.f90
r102 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +uswst, uswst_m, vswst, vswst_m6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 102 2007-07-27 09:09:17Z raasch 13 ! +uswst, uswst_m, vswst, vswst_m 11 14 ! 12 15 ! 96 2007-06-04 08:07:41Z raasch -
palm/trunk/SOURCE/read_var_list.f90
r102 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +time_coupling, top_momentumflux_u|v6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 102 2007-07-27 09:09:17Z raasch 13 ! +time_coupling, top_momentumflux_u|v 11 14 ! 12 15 ! 95 2007-06-02 16:48:38Z raasch -
palm/trunk/SOURCE/surface_coupler.f90
r109 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! include latent heatflux and salinity flux for atmosphere runs with moisture7 ! modifications to terminate coupled runs8 6 ! 9 7 ! … … 12 10 ! $Id$ 13 11 ! 12 ! 109 2007-08-28 15:26:47Z letzel 14 13 ! Initial revision 15 14 ! -
palm/trunk/SOURCE/swap_timelevel.f90
r102 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! swaping of uswst, vswst included6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 102 2007-07-27 09:09:17Z raasch 13 ! swaping of uswst, vswst included 11 14 ! 12 15 ! 95 2007-06-02 16:48:38Z raasch -
palm/trunk/SOURCE/time_integration.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 7 ! 8 ! Former revisions: 9 ! ----------------- 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 6 13 ! Call of new routine surface_coupler, 7 14 ! presure solver is called after the first Runge-Kutta substep instead of the … … 9 16 ! random perturbation has to be added to the velocity fields also after the 10 17 ! first substep 11 ! modifications to terminate coupled runs12 !13 !14 ! Former revisions:15 ! -----------------16 ! $Id$17 18 ! 18 19 ! 97 2007-06-21 08:23:15Z raasch -
palm/trunk/SOURCE/timestep.f90
r108 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! modifications to terminate coupled runs7 6 ! 8 7 ! … … 10 9 ! ----------------- 11 10 ! $Id$ 11 ! 12 ! 108 2007-08-24 15:10:38Z letzel 13 ! modifications to terminate coupled runs 14 ! 12 15 ! RCS Log replace by Id keyword, revision history cleaned up 13 16 ! -
palm/trunk/SOURCE/user_interface.f90
r105 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +dots_num_palm in module user, +module netcdf_control in user_init6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 105 2007-08-08 07:12:55Z raasch 13 ! +dots_num_palm in module user, +module netcdf_control in user_init 11 14 ! 12 15 ! 95 2007-06-02 16:48:38Z raasch -
palm/trunk/SOURCE/write_3d_binary.f90
r102 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +uswst, uswst_m, vswst, vswst_m6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 102 2007-07-27 09:09:17Z raasch 13 ! +uswst, uswst_m, vswst, vswst_m 11 14 ! 12 15 ! 96 2007-06-04 08:07:41Z raasch -
palm/trunk/SOURCE/write_var_list.f90
r102 r110 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! +top_momentumflux_u|v, time_coupling6 ! 7 7 ! 8 8 ! Former revisions: 9 9 ! ----------------- 10 10 ! $Id$ 11 ! 12 ! 102 2007-07-27 09:09:17Z raasch 13 ! +top_momentumflux_u|v, time_coupling 11 14 ! 12 15 ! 95 2007-06-02 16:48:38Z raasch
Note: See TracChangeset
for help on using the changeset viewer.
