Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Changeset 198 for palm


Ignore:
Timestamp:
Sep 17, 2008 8:55:28 AM (16 years ago)
Author:
raasch
Message:

file headers updated for the next release 3.5

Location:
palm/trunk
Files:
43 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/DOC/tec/technical_documentation.html

    r197 r198  
    23512351prognostic_equations</td></tr><tr><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix: summation of<span style="font-family: Courier New,Courier,monospace;"> sums_l_l </span>(<span style="font-family: Courier New,Courier,monospace;">diffusivities</span>).<br><br>Several bugfixes in the ocean part: Initial density<span style="font-family: Courier New,Courier,monospace;"> rho </span>is calculated (<span style="font-family: Courier New,Courier,monospace;">init_ocean</span>). Error in initializing<span style="font-family: Courier New,Courier,monospace;"> u_init </span>and<span style="font-family: Courier New,Courier,monospace;"> v_init </span>removed (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>). Calculation of density flux now starts from<span style="font-family: Courier New,Courier,monospace;"> nzb+1 </span>(<span style="font-family: Courier New,Courier,monospace;">production_e</span>).<br><br>Bugfix: <span style="font-family: Courier New,Courier,monospace;">pleft</span>/<span style="font-family: Courier New,Courier,monospace;">pright</span> changed to <span style="font-family: Courier New,Courier,monospace;">pnorth</span>/<span style="font-family: Courier New,Courier,monospace;">psouth</span> in sendrecv of particle tail numbers along y, small bugfixes in the SGS part (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>).<br><br>Bugfix:<span style="font-family: Courier New,Courier,monospace;"> model_string </span>needed a default value (<span style="font-family: Courier New,Courier,monospace;">combine_plot_fields</span>).<br><br>Bugfix: wavenumber calculation for even<span style="font-family: Courier New,Courier,monospace;"> nx </span>in routines<span style="font-family: Courier New,Courier,monospace;"> maketri </span>(<span style="font-family: Courier New,Courier,monospace;">poisfft</span>).<br><br>Bugfix: assignment of fluxes at walls.<br><br>Bugfix: absolute value of<span style="font-family: Courier New,Courier,monospace;"> f </span>must be used when calculating the Blackadar mixing length (<span style="font-family: Courier New,Courier,monospace;">init_1d_model</span>).</td><td style="vertical-align: top;">advec_particles,
    23522352check_parameters, combine_plot_fields, diffusion_s, diffusivities,
    2353 init_ocean, init_1d_model, poisfft, production_e</td></tr><tr><td style="vertical-align: top;">16/09/08</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.5</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;"></td><td style="vertical-align: top;"></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;"></td><td align="undefined" valign="undefined"></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;"></td><td style="vertical-align: top;"></td></tr>
     2353init_ocean, init_1d_model, poisfft, production_e</td></tr><tr><td style="vertical-align: top;">16/09/08</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.5</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">The
     2354turbulence recycling method of Kataoka and Mizuno (2002) has been
     2355implemented. A precursor run must be used to initialize the main run
     2356with a quasi-stationary turbulence state (see <span style="font-family: monospace;">initializing_actions</span> = <span style="font-style: italic;">'read_data_for_recycling'</span>).
     2357The precursor run may have a smaller domain than the main run, which
     2358helps to spare CPU time. New parameters for steering the recycling
     2359method are <span style="font-family: monospace;">turbulent_inflow</span>, <span style="font-family: monospace;">recycling_width</span>, <span style="font-family: monospace;">inflow_damping_height</span>, and <span style="font-family: monospace;">inflow_damping_width</span>.<br><br>One
     2360new feature related with the recycling method is that within a job
     2361chain, the number of processors or the virtual processor grid is now
     2362allowed to&nbsp;change. Such a change modifies the domain decomposition
     2363and requires that the PEs within a restart job have to know from which
     2364file(s) they have to read their restart data (up to now there has been
     2365a one to one relation between the restart files and the PEs, i.e.
     2366processor 4 only reads data from file "_0004"). Informations about the
     2367subdomain index limits of all restart files are now written by PE0 (to
     2368file "_0000"). Also, all variables containing global values which do
     2369not depend on the subdomain are written to file "_0000" only. In a
     2370restart run, this file is opened by all PEs to read these global values
     2371and the informations about the subdomain index limits. Depending on
     2372these index limit information it is decided which restart file(s)
     2373has/have to be opened in order to read the complete data necessary for
     2374the current subdomain.<br><br>The plant canopy model has been extended
     2375to the prognostic equation for&nbsp;the potential temperature
     2376and&nbsp;passive scalar, as proposed by Watanabe (2004) and Shaw and
     2377Schumann (1992).<br><br>Particle advection adjusted for ocean runs.<br><br>Grid
     2378coarsening (multigrid method) is limited to subdomains (i.e. no
     2379collection of data on PE0 for further coarsening beyond subdomains) if <span style="font-family: monospace;">mg_switch_to_pe0_level</span> = <span style="font-style: italic;">-1</span>.<br><br>User-defined spectra are allowed.<br><br>Output of <span style="font-family: monospace;">q*2</span> profile and <span style="font-family: monospace;">w"u"0</span>, <span style="font-family: monospace;">w"v"0</span>, <span style="font-family: monospace;">w"q"0</span> timeseries are added.<br><br>Argument<span style="font-family: monospace;"> -c </span>introduced to most of the subjob calls, which allows the user to choose his own job catalog by setting the environment variable<span style="font-family: monospace;"> job_catalog </span>in the configuration file (default is <span style="font-family: monospace;">~/job_queue</span>).<br><br>This release contains first adjustments for the SGI-ICE system at HLRN.<br><br><span style="font-weight: bold;">Attention: On the SGI-ICE,&nbsp;runs are still aborted or are hanging for unknown reasons from time to time.</span></td><td style="vertical-align: top;">advec_particles, calc_spectra, check_open, check_parameters, data_output_spectra,
     2380header, init_3d_model, init_particles, init_pegrid, modules, netcdf, parin, particle_boundary_conds, plant_canopy_model,
     2381
     2382prognostic_equations, read_3d_binary, read_var_list, time_integration,
     2383user_interface, write_3d_binary, write_var_list<br><br><span style="font-weight: bold;">new:</span><br>inflow_turbulence<br></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">C</td><td style="vertical-align: top;">Allow 100 spectra levels instead of 10 for consistency with define_netcdf_header. (<span style="font-family: monospace;">calc_spectra</span>, <span style="font-family: monospace;">data_output_spectra</span>, <span style="font-family: monospace;">header</span>, <span style="font-family: monospace;">modules</span>, <span style="font-family: monospace;">netcdf</span>)<br><br>Modification
     2384of the integrated version of the profile function for momentum for
     2385unstable stratification; more consistent flux definitions. (<span style="font-family: monospace;">wall_fluxes</span>, <span style="font-family: monospace;">production_e</span>)<br><br>Strict grid matching along z is not needed for mg-solver. (<span style="font-family: monospace;">check_parameters</span>)<br><br>Sharp
     2386lateral edges of the plant canopy introduced. Before, smaller values of
     2387the leaf area density than those applied in the prognostic equation for
     2388the vertical velocity were used in the prognostic equations for u and v
     2389at the boundaries of the plant_canopy. (<span style="font-family: monospace;">init_3d_model</span>)<br><br>-<span style="font-family: monospace;">cross_ts_*</span> (<span style="font-family: monospace;">read_var_list</span>, <span style="font-family: monospace;">write_var_list</span>), -<span style="font-family: monospace;">dopr_time_count</span>, old profil parameters (<span style="font-family: monospace;">cross_...</span>, <span style="font-family: monospace;">dopr_crossindex</span>,&nbsp;<span style="font-family: monospace;">profile_***</span>) removed. (<span style="font-family: monospace;">read_3d_binary</span>, <span style="font-family: monospace;">write_3d_binary</span>)<br><br>Reading/writing of <span style="font-family: monospace;">hom_sum</span>, <span style="font-family: monospace;">volume_flow_area</span>, <span style="font-family: monospace;">volume_flow_initial</span> is moved from <span style="font-family: monospace;">..._3d_binary</span> to <span style="font-family: monospace;">..._var_list</span>, <span style="font-family: monospace;">hom_sum</span> is now allocated in <span style="font-family: monospace;">parin</span> instead of <span style="font-family: monospace;">init_3d_model</span>.<br><br><span style="font-family: monospace;">npex</span>, <span style="font-family: monospace;">npey</span> moved from inipar- to d3par-namelist. (<span style="font-family: monospace;">parin</span>)<br><br><span style="font-family: monospace;">myid_char14</span>
     2390removed. This was used, since on IBM machines the process rank could be
     2391changed when the final communicator was defined, in order to save the
     2392preliminary processor-id for opening the binary output file for
     2393restarts (unit 14), because otherwise a mismatch occured when reading
     2394this file in a restart job. This mechanism is not needed any more since
     2395the subdomain informations are now contained in file _0000. (<span style="font-family: monospace;">parin</span>, <span style="font-family: monospace;">check_open</span>)<br><br>Transpositions for the 2D domain decomposition have been optimized by using <span style="font-family: monospace;">f_inv</span>
     2396as an automatic array instead of providing the memory by a dummy
     2397argument. This spares one copy loop per transposition. Order of indices
     2398in the 3D loops in some of the transpose routines have been rearranged
     2399for better cache utilization. Both have been suggested by Roland
     2400Richter (SGI) as part of the HLRN-II benchmark process. (<span style="font-family: monospace;">transpose</span>)<br><br>Workaround for getting information about the coupling mode, because the<span style="font-family: monospace;"> -env </span>Option (mpiexec) is not available on many systems (SGI, lcxt4 at&nbsp; Bergen Center for Computational Science). (<span style="font-family: monospace;">palm</span>, <span style="font-family: monospace;">mrun</span>)<br><br>Format adjustments in order to allow CPU# &gt; 999. (<span style="font-family: monospace;">cpu_statistics</span>)<br><br>Provisional solution for run_control_1d output: add 'CALL check_open( 15 )'. (<span style="font-family: monospace;">init_1d_model</span>)<br><br>Subdomain data are read into temporary arrays <span style="font-family: monospace;">pf_tmp</span> / <span style="font-family: monospace;">pf3d_tmp</span> in order to avoid INTEL compiler warnings about (automatic) creation of temporary arrays. (<span style="font-family: monospace;">combine_plot_fields</span>)<br><br>Time coordinate t=0 stored on netcdf-file only if an output is required for this time for at least one of the profiles.<br><br>Output of spline timeseries <span style="font-family: monospace;">splptx</span> etc. removed. (<span style="font-family: monospace;">flow_statistics</span>)</td><td style="vertical-align: top;">advec_s_ups,
     2401advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open,
     2402check_parameters, cpu_statistics, data_output_profiles,
     2403data_output_spectra, flow_statistics, header, init_1d_model,
     2404init_3d_model, modules, netcdf, palm, parin, poisfft, production_e,
     2405read_3d_binary, read_var_list, transpose, wall_fluxes, write_3d_binary, write_var_list<br></td></tr><tr><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td align="undefined" valign="undefined">&nbsp;</td><td style="vertical-align: top;">E</td><td style="vertical-align: top;">Bugfix:
     2406modification of the calculation of the vertical turbulent momentum
     2407fluxes u'w' and v'w'; the first usws that is computed corresponds to
     2408-u'w'/u* and not as priorily assumed to (-u'w')**0.5, the first vsws
     2409that is computed corresponds to -v'w'/u* and not as priorily assumed to
     2410(-v'w')**0.5. Therefore, the intermediate result for usws has to be
     2411multiplied by -u* instead by itself in order to get u'w'. Accordingly,
     2412the intermediate result for vsws has to be multiplied by -u* instead by
     2413itself in order to get v'w'. As u* is calculated for the position of a
     2414scalar an additional interpolation of u* to the position of u and v,
     2415respectively, is necessary. As u* is not determined for the ghost
     2416points on each PE, an additional exchange of information from
     2417neighbouring PEs is necessary. Tests showed that this
     2418decreases u* by some 10% and increases TKE and momentum fluxes by some
     241910% because friction was underestimated and momentum fluxes were
     2420wrongly calculated due to the bug. (<span style="font-family: monospace;">prandtl_fluxes</span>, <span style="font-family: monospace;">wall_fluxes</span>)<br><br><span style="font-weight: bold;">Attention:</span> <br><span style="font-weight: bold;">
     2421The content of the testresult-file (run control output) has changed
     2422due to the above bugfix.</span><br><br>Bugfix:
     2423change definition of us_wall from 1D to 2D.&nbsp;<br><br>Bugfix: exponent = 1.0 / ( ny + 1.0 ) in calc_spectra_x. (<span style="font-family: monospace;">calc_spectra</span>)<br><br>Bugfix: calculation of horizontal fluxes at vertical walls. (<span style="font-family: monospace;">diffusion_s</span>)<br><br>Bugfix: zero assignments to tendency arrays in case of restarts. (<span style="font-family: monospace;">init_3d_model</span>)<br><br>Bugfix: wrong grid level used in multigrid solver, routine <span style="font-family: monospace;">restrict</span>. (<span style="font-family: monospace;">poismg</span>)<br><br>Bugfix:
     2424Construction of vertical profiles when 10 gradients have been specified
     2425in the parameter list (ug, vg, pt, q, sa, lad). (<span style="font-family: monospace;">check_parameters</span>)<br><br>Bugfix: Determination of the number of vertical levels for which spectra shall be output - 100 replaced by 10. (<span style="font-family: monospace;">netcdf</span>)<br><br>Bugfix:
     2426Lengths of those characters that are needed for the output of the
     2427characteristic levels of potential temperature, the geostrophic wind,
     2428scalar concentration, humidity and leaf area density to the header file
     2429has been increased. (<span style="font-family: monospace;">header</span>)<br><br>Bugfix: add ')' in <span style="font-family: monospace;">user_read_restart_data</span>.<br><br>The leaf area density is explicitly set to its surface value at k=0. (<span style="font-family: monospace;">check_parameters</span>)<br><br>Introduce<span style="font-family: monospace;"> prefix_chr</span> to ensure unique dvrp_file path.<br><br>small bugfixes for <span style="font-family: monospace;">user_interface</span> sample code (comments):<br>- initialize ustvst with 0.0 as it is now computed only until <span style="font-family: monospace;">nxr</span> and <span style="font-family: monospace;">nyn</span><br>- two ALLOCATE statements moved from <span style="font-family: monospace;">user_read_restart_data</span> back to <span style="font-family: monospace;">user_init</span><br>- remove 'READ (13) u2_av' statement in <span style="font-family: monospace;">user_read_restart_data</span><br><br>Bugfix: remove IF statement in <span style="font-family: monospace;">plant_canopy_model_ij</span>.<br><br>Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr) (like other scalars). (<span style="font-family: monospace;">flow_statistics</span>)<br><br>Bugfix: <span style="font-family: monospace;">dopr_time_count</span>
     2430was written on the binary file, which caused that NetCDF files newly
     2431created by restart files (no append of existing files!) contained
     2432uneccessary time levels. (<span style="font-family: monospace;">read_3d_binary</span>, <span style="font-family: monospace;">write_3d_binary</span>)<br><br>Bugfix: extra '*' removed in user_statistics sample code. (<span style="font-family: monospace;">user_interface</span>)<br><br>Bugfix: a stop command was missing in some cases of the parallel branch. (<span style="font-family: monospace;">local_stop</span>)<br><br>Bugfix in volume flow control for non-cyclic boundary conditions. (<span style="font-family: monospace;">pres</span>)<br><br>Bugfix: misplaced #endif directives. (<span style="font-family: monospace;">combine_plot_fields</span>)<br><br>Data are collected from PE0 in an ordered sequence which seems to avoid hanging of processes on SGI-ICE. (<span style="font-family: monospace;">cpu_statistics</span>)</td><td style="vertical-align: top;">calc_spectra,
     2433check_parameters, cpu_statistics, diffusion_s, flow_statistics,
     2434init_3d_model, init_dvrp, local_stop, plant_canopy_model, poismg,
     2435prandtl_fluxes, pres, read_3d_binary, user_interface, wall_fluxes,
     2436write_3d_binary<br></td></tr>
    23542437</tbody>
    23552438</table>&nbsp;<b><blink>Attention:</blink></b>

  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r197 r198  
    11New:
    22---
    3 The plant canopy model has been extended to the prognostic equation for   
    4 the potential temperature and the prognostic equation for the concentration
    5 of a passive scalar, as proposed by Watanabe (2004) and Shaw and Schumann
    6 (1992). (prognostic_equations, plant_canopy_model) 
    73
    8 Within a job chain, the number of processors or the virtual processor grid
    9 is now allowed to be changed. This modifies the domain decomposition and
    10 requires that the PEs within a restart job have to know from which file(s) they
    11 have to read their restart data (up to now there is a one to one relation
    12 between the restart files and the PEs, i.e. processor 4 only reads data
    13 from file "_0004"). Informations about the subdomain index limits of all
    14 restart files are now written by PE0 (to file "_0000"). Also, all variables
    15 containing global values which do not depend on the subdomain are written
    16 to file "_0000" only. In a restart run, this file is opened by all PEs to
    17 read these global values and the informations about the subdomain index
    18 limits. Depending on these index limit information it is decided which
    19 restart file(s) has/have to be opened in order to read the complete data
    20 necessary for the current subdomain.
    21 
    22 The turbulence recycling method of Kataoka and Mizuno (2002) has been
    23 implemented. A pre-run with smaller model domain can be used to initialize
    24 the real run, in order to spare the CPU time needed to get the turbulence
    25 into a quasi-stationary state (see initializing_actions =
    26 'read_data_for_recycling').
    27 
    28 Particle advection adjusted for ocean runs.
    29 
    30 Multigrid levels are limited by subdomains if mg_switch_to_pe0_level = -1.
    31 
    32 User-defined spectra.
    33 
    34 Output of q*2 profile added.
    35 
    36 Argument -c introduced to most of the subjob calls, which allows the user to
    37 choose his own job catalog by setting job_catalog in the configuration file
    38 (default is ~/job_queue). Workaround for mpixec with -env option.
    39 Adjustments for lcxt4 (Bergen Center for Computational Science) (mrun, subjob)
    40 
    41 advec_particles, calc_spectra, check_open, check_parameters, data_output_spectra, header, init_particles, init_pegrid, init_3d_model, modules, netcdf, parin, particle_boundary_conds, plant_canopy_model, prognostic_equations, read_var_list, read_3d_binary, time_integration, user_interface, write_var_list, write_3d_binary
    42 
    43 New: inflow_turbulence
    444
    455
    466Changed:
    477-------
    48 Allow 100 spectra levels instead of 10 for consistency with
    49 define_netcdf_header (calc_spectra, data_output_spectra, header, modules,
    50 netcdf)
    518
    52 Modification of the integrated version of the profile function for momentum
    53 for unstable stratification; more consistent flux definitions. (wall_fluxes,
    54 production_e)
    55 
    56 Strict grid matching along z is not needed for mg-solver. (check_parameters)
    57 
    58 Sharp lateral edges of the plant canopy introduced. Before smaller
    59 values of the leaf area density than those applied in the prognostic
    60 equation for the vertical velocity were used in the prognostic equations
    61 for u and v at the boundaries of the plant_canopy. (init_3d_model) 
    62 
    63 -cross_ts_* (read_var_list, write_var_list)
    64 -dopr_time_count, old profil parameters (cross_..., dopr_crossindex,
    65  profile_***) removed (read_3d_binary, write_3d_binary)
    66 
    67 hom_sum moved, volume_flow_area, volume_flow_initial moved from ..._3d_binary
    68 to ..._var_list, hom_sum is now allocated in parin instead of init_3d_model
    69 
    70 npex, npey moved from inipar- to d3par-namelist (parin)
    71 
    72 myid_char14 removed. This was used since on IBM machines the process rank may
    73 be changed when the final communicator is defined, in order to save the
    74 preliminary processor-id for opening the binary output file for restarts
    75 (unit 14), because otherwise a mismatch occured when reading this file in a
    76 restart job. This mechanism is not needed any more since the subdomain
    77 informations are now contained in file _0000. (parin, check_open)
    78 
    79 Transpositions for the 2D domain decomposition have been optimized by using
    80 f_inv as an automatic array instead of providing the memory by a dummy argument.
    81 This spares one copy loop per transposition. Order of indices in the 3D loops
    82 in some of the transpose routines have been rearranged for better cache utilization.
    83 Both have been suggested by Roland Richter (SGI) as part of the
    84 HLRN-II benchmark process. (transpose)
    85 
    86 Workaround for getting information about the coupling mode. (palm)
    87 
    88 Format adjustments in order to allow CPU# > 999. (cpu_statistics)
    89 
    90 Provisional solution for run_control_1d output: add 'CALL check_open( 15 )'
    91 (init_1d_model)
    92 
    93 Subdomain data are read into temporary arrays pf_tmp/pf3d_tmp in order to
    94 avoid INTEL compiler warnings about (automatic) creation of temporary arrays
    95 (combine_plot_fields)
    96 
    97 Time coordinate t=0 stored on netcdf-file only if an output is required for
    98 this time for at least one of the profiles
    99 
    100 Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0) added
    101 (flow_statistics)
    102 
    103 advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open, check_parameters, cpu_statistics, data_output_profiles, data_output_spectra, flow_statistics, header, init_1d_model, init_3d_model, modules, netcdf, palm, parin, poisfft, production_e, read_var_list, read_3d_binary, transpose, wall_fluxes, write_var_list, write_3d_binary
    1049
    10510
    10611Errors:
    10712------
    108 Bugfix: exponent = 1.0 / ( ny + 1.0 ) in calc_spectra_x (calc_spectra)
    10913
    110 Bugfix: change definition of us_wall from 1D to 2D. Tests showed that this decreases
    111 u* by some 10% and increases TKE and momentum fluxes by some 10% because friction
    112 was underestimated and momentum fluxes were wrongly calculated due to the bug.
    113 (prandtl_fluxes, wall_fluxes)
    114 
    115 Bugfix: calculation of horizontal fluxes at vertical walls (diffusion_s)
    116 
    117 Bugfix: zero assignments to tendency arrays in case of restarts (init_3d_model)
    118 
    119 Bugfix: wrong grid level used in multigrid solver, routine restrict (poismg)
    120 
    121 Bugfix: Construction of vertical profiles when 10 gradients have been
    122 specified in the parameter list (ug, vg, pt, q, sa, lad). (check_parameters)
    123 
    124 Bugfix: Determination of the number of vertical levels for which spectra
    125 shall be output - 100 replaced by 10. (netcdf)
    126 
    127 Bugfix: Lengths of those characters that are needed for the output of the
    128 characteristic levels of potential temperature, the geostrophic wind,
    129 scalar concentration, humidity and leaf area density to the header file
    130 has been increased. (header)
    131 
    132 Bugfix: add ')' in user_read_restart_data.
    133 
    134 The leaf area density is explicitly set to its surface value at k=0. (check_parameters)
    135 
    136 Introduce prefix_chr to ensure unique dvrp_file path.
    137 
    138 small bugfixes for user_interface sample code (comments):
    139 - initialize ustvst with 0.0 as it is now computed only until nxr and nyn
    140 - two ALLOCATE statements moved from user_read_restart_data back to user_init
    141 - remove 'READ (13) u2_av' statement in user_read_restart_data
    142 
    143 Bugfix: remove IF statement in plant_canopy_model_ij (plant_canopy_model)
    144 
    145 Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr)
    146 (like other scalars) (flow_statistics)
    147 
    148 Bugfix: dopr_time_count was written on the binary file, which caused that
    149 NetCDF files newly created by restart files (no append of existing files!)
    150 contained uneccessary time levels. (read_3d_binary, write_3d_binary)
    151 
    152 Bugfix: extra '*' removed in user_statistics sample code (user_interface)
    153 
    154 Bugfix: a stop command was missing in some cases of the parallel branch (local_stop)
    155 
    156 Bugfix in volume flow control for non-cyclic boundary conditions (pres)
    157 
    158 Bugfix: misplaced #endif directives (combine_plot_fields)
    159 
    160 data are collected from PE0 in an ordered sequence which seems to avoid hanging of processes on SGI-ICE (cpu_statistics)
    161 
    162 calc_spectra, check_parameters, cpu_statistics, diffusion_s, flow_statistics, init_dvrp, init_3d_model, local_stop, plant_canopy_model, poismg, prandtl_fluxes, pres, read_3d_binary, user_interface, wall_fluxes, write_3d_binary
    163 

  • palm/trunk/SOURCE/advec_particles.f90

    r150 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Bottom boundary condition and vertical index calculations adjusted for
    7 ! ocean runs.
    86! TEST: PRINT statements on unit 9 (commented out)
    97!
    … …  
    119! -----------------
    1210! $Id$
     11!
     12! 150 2008-02-29 08:19:58Z raasch
     13! Bottom boundary condition and vertical index calculations adjusted for
     14! ocean runs.
    1315!
    1416! 119 2007-10-17 10:27:13Z raasch

  • palm/trunk/SOURCE/advec_s_ups.f90

    r164 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Arguments removed from transpose routines
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 164 2008-05-15 08:46:15Z raasch
     13! Arguments removed from transpose routines
     14!
     15! February 2007
    1116! RCS Log replace by Id keyword, revision history cleaned up
    1217!

  • palm/trunk/SOURCE/advec_u_ups.f90

    r164 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Arguments removed from transpose routines
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 164 2008-05-15 08:46:15Z raasch
     13! Arguments removed from transpose routines
     14!
     15! February 2007
    1116! RCS Log replace by Id keyword, revision history cleaned up
    1217!

  • palm/trunk/SOURCE/advec_v_ups.f90

    r164 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Arguments removed from transpose routines
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 164 2008-05-15 08:46:15Z raasch
     13! Arguments removed from transpose routines
     14!
     15! February 2007
    1116! RCS Log replace by Id keyword, revision history cleaned up
    1217!

  • palm/trunk/SOURCE/advec_w_ups.f90

    r164 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Arguments removed from transpose routines
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 164 2008-05-15 08:46:15Z raasch
     13! Arguments removed from transpose routines
     14!
     15! February 2007
    1116! RCS Log replace by Id keyword, revision history cleaned up
    1217!

  • palm/trunk/SOURCE/calc_spectra.f90

    r192 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 192 2008-08-27 16:51:49Z letzel
    613! bugfix in calc_spectra_x: exponent = 1.0 / ( ny + 1.0 )
    714! allow 100 spectra levels instead of 10 for consistency with
    … …  
    916! user-defined spectra, arguments removed from transpose routines
    1017!
    11 ! Former revisions:
    12 ! -----------------
    13 ! $Id$
     18! February 2007
    1419! RCS Log replace by Id keyword, revision history cleaned up
    1520!

  • palm/trunk/SOURCE/check_open.f90

    r146 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 146 2008-01-17 13:08:34Z raasch
    613! First opening of unit 13 openes file _0000 on all PEs (parallel version)
    714! because only this file contains the global variables,
    815! myid_char_14 removed
    9 !
    10 ! Former revisions:
    11 ! -----------------
    12 ! $Id$
    1316!
    1417! 120 2007-10-17 11:54:43Z raasch

  • palm/trunk/SOURCE/check_parameters.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
    613! Bug fix: Construction of vertical profiles when 10 gradients have been
    714! specified in the parameter list (ug, vg, pt, q, sa, lad)
    … …  
    1219! check of turbulent_inflow and calculation of recycling_plane.
    1320! q*2 profile added
    14 !
    15 ! Former revisions:
    16 ! -----------------
    17 ! $Id$
    1821!
    1922! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/cpu_statistics.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
    613! Format adjustments in order to allow CPU# > 999,
    714! data are collected from PE0 in an ordered sequence which seems to avoid
    815! hanging of processes on SGI-ICE
    9 !
    10 ! Former revisions:
    11 ! -----------------
    12 ! $Id$
    1316!
    1417! 82 2007-04-16 15:40:52Z raasch

  • palm/trunk/SOURCE/data_output_profiles.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Time coordinate t=0 stored on netcdf-file only if an output is required for
    7 ! this time for at least one of the profiles
     6!
    87!
    98! Former revisions:
    109! -----------------
    1110! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
     13! Time coordinate t=0 stored on netcdf-file only if an output is required for
     14! this time for at least one of the profiles
     15!
     16! February 2007
    1217! RCS Log replace by Id keyword, revision history cleaned up
    1318!

  • palm/trunk/SOURCE/data_output_spectra.f90

    r189 r198  
    44! Actual revisions:
    55! -----------------
    6 ! allow 100 spectra levels instead of 10 for consistency with
    7 ! define_netcdf_header
    8 ! user-defined spectra
     6!
    97!
    108! Former revisions:
    119! -----------------
    1210! $Id$
     11!
     12! 189 2008-08-13 17:09:26Z letzel
     13! allow 100 spectra levels instead of 10 for consistency with
     14! define_netcdf_header, +user-defined spectra
     15!
     16! February 2007
    1317! RCS Log replace by Id keyword, revision history cleaned up
    1418!

  • palm/trunk/SOURCE/diffusion_s.f90

    r183 r198  
    44! Actual revisions:
    55! -----------------
    6 ! bugfix: calculation of fluxes at vertical surfaces
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 183 2008-08-04 15:39:12Z letzel
     13! bugfix: calculation of fluxes at vertical surfaces
    1114!
    1215! 129 2007-10-30 12:12:24Z letzel

  • palm/trunk/SOURCE/flow_statistics.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0)
    7 ! added
    8 ! Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr)
    9 ! (like other scalars)
     6!
    107!
    118! Former revisions:
    129! -----------------
    1310! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
     13! Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0)
     14! added,
     15! bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr)
     16! (like other scalars)
    1417!
    1518! 133 2007-11-20 10:10:53Z letzel

  • palm/trunk/SOURCE/header.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
    6 ! TEST: mg_switch_to_pe0_level = -1!!!!!!!!!
    7 ! allow 100 spectra levels instead of 10 for consistency with
    8 ! define_netcdf_header
    9 ! Bugfix in the output of the characteristic levels of potential temperature,
    10 ! geostrophic wind, scalar concentration, humidity and leaf area density
    11 ! Output of turbulence recycling informations
     6!
    127!
    138! Former revisions:
    149! -----------------
    1510! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
     13! allow 100 spectra levels instead of 10 for consistency with
     14! define_netcdf_header,
     15! bugfix in the output of the characteristic levels of potential temperature,
     16! geostrophic wind, scalar concentration, humidity and leaf area density,
     17! output of turbulence recycling informations
    1618!
    1719! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/inflow_turbulence.f90

    r163 r198  
    1010! $Id$
    1111!
     12! Initial version (2008/03/07)
    1213!
    1314! Description:

  • palm/trunk/SOURCE/init_1d_model.f90

    r184 r198  
    44! Actual revisions:
    55! -----------------
    6 ! provisional solution for run_control_1d output: add 'CALL check_open( 15 )'
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 184 2008-08-04 15:53:39Z letzel
     13! provisional solution for run_control_1d output: add 'CALL check_open( 15 )'
    1114!
    1215! 135 2007-11-22 12:24:23Z raasch

  • palm/trunk/SOURCE/init_3d_model.f90

    r181 r198  
    77! Actual revisions:
    88! -----------------
    9 ! Bugfix: zero assignments to tendency arrays in case of restarts
    10 ! Further amendments to and modifications in the initialisation of the plant
    11 ! canopy model
    12 ! Allocation of hom_sum moved to parin, initialization of spectrum_x|y directly
     9!
     10!
     11! Former revisions:
     12! -----------------
     13! $Id$
     14!
     15! 181 2008-07-30 07:07:47Z raasch
     16! bugfix: zero assignments to tendency arrays in case of restarts,
     17! further extensions and modifications in the initialisation of the plant
     18! canopy model,
     19! allocation of hom_sum moved to parin, initialization of spectrum_x|y directly
    1320! after allocating theses arrays,
    1421! read data for recycling added as new initialization option,
    1522! dummy allocation for diss
    16 !
    17 ! Former revisions:
    18 ! -----------------
    19 ! $Id$
    2023!
    2124! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/init_dvrp.f90

    r155 r198  
    55! -----------------
    66!
    7 ! introduce prefix_chr to ensure unique dvrp_file path
    87! TEST: print* statements
    98! ToDo: checking of mode_dvrp for legal values is not correct
    … …  
    1211! -----------------
    1312! $Id$
     13!
     14! 155 2008-03-28 10:56:30Z letzel
     15! introduce prefix_chr to ensure unique dvrp_file path
    1416!
    1517! 130 2007-11-13 14:08:40Z letzel

  • palm/trunk/SOURCE/init_particles.f90

    r150 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Set offset_ocean_* needed for calculating vertical indices within ocean runs
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 150 2008-02-29 08:19:58Z raasch
     13! Setting offset_ocean_* needed for calculating vertical indices within ocean
     14! runs
    1115!
    1216! 117 2007-10-11 03:27:59Z raasch

  • palm/trunk/SOURCE/init_pegrid.f90

    r197 r198  
    55! -----------------
    66! ATTENTION: nnz_x undefined problem still has to be solved!!!!!!!!
     7! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
     8!
     9! Former revisions:
     10! -----------------
     11! $Id$
     12!
     13! 197 2008-09-16 15:29:03Z raasch
    714! multigrid levels are limited by subdomains if mg_switch_to_pe0_level = -1,
    815! nz is used instead nnz for calculating mg-levels
    916! Collect on PE0 horizontal index bounds from all other PEs,
    1017! broadcast the id of the inflow PE (using the respective communicator)
    11 ! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
    12 !
    13 ! Former revisions:
    14 ! -----------------
    15 ! $Id$
    1618!
    1719! 114 2007-10-10 00:03:15Z raasch

  • palm/trunk/SOURCE/local_stop.f90

    r147 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: a stop command was missing in some cases of the parallel branch
     6!
    77!
    88!
    … …  
    1010! -----------------
    1111! $Id$
     12!
     13! 147 2008-02-01 12:41:46Z raasch
     14! Bugfix: a stop command was missing in some cases of the parallel branch
    1215!
    1316! 108 2007-08-24 15:10:38Z letzel

  • palm/trunk/SOURCE/modules.f90

    r197 r198  
    55! Actual revisions:
    66! -----------------
     7!
     8!
     9! Former revisions:
     10! -----------------
     11! $Id$
     12!
     13! 197 2008-09-16 15:29:03Z raasch
    714! allow 100 spectra levels instead of 10 for consistency with
    8 ! define_netcdf_header
    9 ! +canopy_heat_flux, cthf, lai,
     15! define_netcdf_header, +canopy_heat_flux, cthf, lai,
    1016! +leaf_surface_concentration, scalar_exchange_coefficient, sec, sls
    1117! +hor_index_bounds, hor_index_bounds_previous_run, id_inflow, id_recycling,
    … …  
    1319! ny_on_file, offset_ocean_*, recycling_plane, recycling_width, turbulent_inflow
    1420! -myid_char_14
    15 !
    16 ! Former revisions:
    17 ! -----------------
    18 ! $Id$
    1921!
    2022! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/netcdf.f90

    r189 r198  
    77! Current revisions:
    88! ------------------
    9 ! consistently allow 100 spectra levels instead of 10
    10 ! bug fix in the determination of the number of output heights for spectra
    11 !
    12 ! user-defined spectra
     9!
    1310!
    1411! Former revisions:
    1512! -----------------
    1613! $Id$
     14!
     15! 189 2008-08-13 17:09:26Z letzel
     16! consistently allow 100 spectra levels instead of 10
     17! bug fix in the determination of the number of output heights for spectra,
     18! +user-defined spectra
    1719!
    1820! 97 2007-06-21 08:23:15Z raasch

  • palm/trunk/SOURCE/palm.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Workaround for getting information about the coupling mode
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
     13! Workaround for getting information about the coupling mode
    1114!
    1215! 108 2007-08-24 15:10:38Z letzel

  • palm/trunk/SOURCE/parin.f90

    r197 r198  
    44! Actual revisions:
    55! -----------------
    6 ! +cthf
    7 ! +leaf_surface_concentration, scalar_exchange_coefficient
    8 ! +inflow_damping_height, inflow_damping_width, recycling_width,
    9 ! turbulent_inflow in inipar
    10 ! -skip_time_dosp in d3par,
    11 ! Allocation of hom_sum moved from init_3d_model to here,
    12 ! npex, npey moved from inipar to d3par, setting of myid_char_14 removed,
    13 ! lad is allways allocated
     6!
    147!
    158! Former revisions:
    169! -----------------
    1710! $Id$
     11!
     12! 197 2008-09-16 15:29:03Z raasch
     13! +cthf,leaf_surface_concentration, scalar_exchange_coefficient
     14! +inflow_damping_height, inflow_damping_width, recycling_width,
     15! turbulent_inflow in inipar, -skip_time_dosp in d3par,
     16! allocation of hom_sum moved from init_3d_model to here,
     17! npex, npey moved from inipar to d3par, setting of myid_char_14 removed,
     18! lad is allways allocated
    1819!
    1920! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/particle_boundary_conds.f90

    r150 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Vertical index calculations adjusted for ocean runs.
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 150 2008-02-29 08:19:58Z raasch
     13! Vertical index calculations adjusted for ocean runs.
     14!
    1115! Initial version (2007/03/09)
    1216!

  • palm/trunk/SOURCE/plant_canopy_model.f90

    r153 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Consideration of heat sources within the forest canopy which represent the
    7 ! rate of heat input into the air from the forest leaves
    8 !
    9 ! Evaluation of sinks and sources for scalar concentration due to canopy
    10 ! elements
     6!
    117!
    128! Bugfix: remove IF statement in plant_canopy_model_ij
    … …  
    1511! -----------------
    1612! $Id$
     13!
     14! 153 2008-03-19 09:41:30Z steinfeld
     15! heat sources within the forest canopy are added, which represent the
     16! rate of heat input into the air from the forest leaves, evaluation of sinks
     17! and sources for scalar concentration due to canopy elements
    1718!
    1819! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/poisfft.f90

    r164 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Arguments removed from transpose routines
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 164 2008-05-15 08:46:15Z raasch
     13! Arguments removed from transpose routines
    1114!
    1215! 128 2007-10-26 13:11:14Z raasch

  • palm/trunk/SOURCE/poismg.f90

    r181 r198  
    88! Actual revisions:
    99! -----------------
    10 ! Bugfix: grid_level+1 has to be used in restrict for flags-array
     10!
    1111!
    1212! Former revisions:
    1313! -----------------
    1414! $Id$
     15!
     16! 181 2008-07-30 07:07:47Z raasch
     17! Bugfix: grid_level+1 has to be used in restrict for flags-array
    1518!
    1619! 114 2007-10-10 00:03:15Z raasch

  • palm/trunk/SOURCE/prandtl_fluxes.f90

    r187 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: change definition of us_wall from 1D to 2D:
    7 !         Modification of the evaluation of the vertical turbulent momentum
    8 !         fluxes u'w' and v'w'; the first usws that is computed corresponds
    9 !         to -u'w'/u* and not as priorily assumed to (-u'w')**0.5, the first
    10 !         vsws that is computed corresponds to -v'w'/u* and not as priorily
    11 !         assumed to (-v'w')**0.5. Therefore, the intermediate result for
    12 !         usws has to be multiplied by -u* instead by itself in order to
    13 !         get u'w'. Accordingly, the intermediate result for vsws has to be
    14 !         multiplied by -u* instead by itself in order to get v'w'. As u*
    15 !         is calculated for the position of a scalar an additional
    16 !         interpolation of u* to the position of u and v, respectively,
    17 !         is necessary. As u* is not determined for the ghost points on
    18 !         each PE, an additional exchange of information from neighbouring
    19 !         PEs is necessary.! 
    20 ! Change: Modification of the integrated version of the profile function for
    21 !         momentum for unstable stratification
     6!
    227!
    238! Former revisions:
    249! -----------------
    2510! $Id$
     11!
     12! 187 2008-08-06 16:25:09Z letzel
     13! Bugfix: modification of the calculation of the vertical turbulent momentum
     14! fluxes u'w' and v'w'
     15! Bugfix: change definition of us_wall from 1D to 2D
     16! Change: modification of the integrated version of the profile function for
     17! momentum for unstable stratification (does not effect results)
    2618!
    2719! 108 2007-08-24 15:10:38Z letzel

  • palm/trunk/SOURCE/pres.f90

    r151 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix in volume flow control for non-cyclic boundary conditions
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 151 2008-03-07 13:42:18Z raasch
     13! Bugfix in volume flow control for non-cyclic boundary conditions
    1114!
    1215! 106 2007-08-16 14:30:26Z raasch

  • palm/trunk/SOURCE/production_e.f90

    r187 r198  
    44! Actual revisions:
    55! -----------------
    6 ! Change: add 'minus' sign to fluxes obtained from subroutine wall_fluxes_e for
    7 !         consistency with subroutine wall_fluxes
     6!
    87!
    98! Former revisions:
    109! -----------------
    1110! $Id$
     11!
     12! 187 2008-08-06 16:25:09Z letzel
     13! Change: add 'minus' sign to fluxes obtained from subroutine wall_fluxes_e for
     14! consistency with subroutine wall_fluxes
    1215!
    1316! 124 2007-10-19 15:47:46Z raasch

  • palm/trunk/SOURCE/prognostic_equations.f90

    r153 r198  
    44! Actual revisions:
    55! -----------------
    6 ! add call of plant_canopy_model in the prognostic equation for
    7 ! the potential temperature
    8 ! add call of plant_canopy_model in the prognostic equation for
    9 ! the concentration of a passive scalar
     6!
    107!
    118! Former revisions:
    129! -----------------
    1310! $Id$
     11!
     12! 153 2008-03-19 09:41:30Z steinfeld
     13! add call of plant_canopy_model in the prognostic equation for
     14! the potential temperature and for the passive scalar
    1415!
    1516! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/read_3d_binary.f90

    r150 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 150 2008-02-29 08:19:58Z raasch
    613! Files from which restart data are to be read are determined and subsequently
    714! opened. The total domain on the restart file is allowed to be smaller than
    … …  
    1219! reading of old profil parameters (cross_..., dopr_crossindex, profile_***)
    1320! removed, initialization of spectrum_x|y removed
    14 !
    15 ! Former revisions:
    16 ! -----------------
    17 ! $Id$
    1821!
    1922! 102 2007-07-27 09:09:17Z raasch

  • palm/trunk/SOURCE/read_var_list.f90

    r173 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 173 2008-05-23 20:39:38Z raasch
    613! +cthf, leaf_surface_concentration, scalar_exchange_coefficient
    714! +numprocs_previous_run, hor_index_bounds_previous_run, inflow_damping_factor,
    … …  
    1219! read_3d_binary to here,
    1320! routines read_parts_of_var_list and skip_var_list added at the end
    14 !
    15 ! Former revisions:
    16 ! -----------------
    17 ! $Id$
    1821!
    1922! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/time_integration.f90

    r151 r198  
    44! Actual revisions:
    55! -----------------
    6 ! inflow turbulence is imposed by calling new routine inflow_turbulence
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 151 2008-03-07 13:42:18Z raasch
     13! inflow turbulence is imposed by calling new routine inflow_turbulence
    1114!
    1215! 108 2007-08-24 15:10:38Z letzel

  • palm/trunk/SOURCE/transpose.f90

    r164 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 164 2008-05-15 08:46:15Z raasch
    613! f_inv changed from subroutine argument to automatic array in order to do
    714! re-ordering from f_in to f_inv in one step, one array work is needed instead
    815! of work1 and work2
    916!
    10 ! Former revisions:
    11 ! -----------------
    12 ! $Id$
     17! February 2007
    1318! RCS Log replace by Id keyword, revision history cleaned up
    1419!

  • palm/trunk/SOURCE/user_interface.f90

    r159 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 159 2008-04-09 08:06:03Z raasch
    613! bugfix in user_read_restart_data
    714! small bugfixes for sample code (comments):
    … …  
    1017! - remove 'READ (13) u2_av' statement in user_read_restart_data
    1118! +routines user_read_restart_data, user_spectra
    12 !
    13 ! Former revisions:
    14 ! -----------------
    15 ! $Id$
    1619!
    1720! 138 2007-11-28 10:03:58Z letzel

  • palm/trunk/SOURCE/wall_fluxes.f90

    r187 r198  
    33! Actual revisions:
    44! -----------------
    5 ! Bugfix: change definition of us_wall from 1D to 2D:
    6 !         Modification of the evaluation of the vertical turbulent momentum
    7 !         fluxes u'w' and v'w'; the first usws that is computed corresponds
    8 !         to -u'w'/u* and not as priorily assumed to (-u'w')**0.5, the first
    9 !         vsws that is computed corresponds to -v'w'/u* and not as priorily
    10 !         assumed to (-v'w')**0.5. Therefore, the intermediate result for
    11 !         usws has to be multiplied by -u* instead by itself in order to
    12 !         get u'w'. Accordingly, the intermediate result for vsws has to be
    13 !         multiplied by -u* instead by itself in order to get v'w'.
    14 !         This requires the calculation of us_wall (and vel_total, u_i, v_i, ws)
    15 !         also in wall_fluxes_e.
    16 ! Bugfix: storage of rifs to rifs_wall in wall_fluxes_e removed
    17 ! Change: add 'minus' sign to fluxes produced by subroutine wall_fluxes_e for
    18 !         consistency with subroutine wall_fluxes
    19 ! Change: Modification of the integrated version of the profile function for
    20 !         momentum for unstable stratification
     5!
    216!
    227! Former revisions:
    238! -----------------
    249! $Id$
     10!
     11! 187 2008-08-06 16:25:09Z letzel
     12! Bugfix: Modification of the evaluation of the vertical turbulent momentum
     13! fluxes u'w' and v'w (see prandtl_fluxes), this requires the calculation of
     14! us_wall (and vel_total, u_i, v_i, ws) also in wall_fluxes_e.
     15! Bugfix: change definition of us_wall from 1D to 2D
     16! Bugfix: storage of rifs to rifs_wall in wall_fluxes_e removed
     17! Change: add 'minus' sign to fluxes produced by subroutine wall_fluxes_e for
     18! consistency with subroutine wall_fluxes
     19! Change: Modification of the integrated version of the profile function for
     20! momentum for unstable stratification
     21!
    2522! Initial version (2007/03/07)
    2623!

  • palm/trunk/SOURCE/write_3d_binary.f90

    r145 r198  
    44! Actual revisions:
    55! -----------------
    6 ! write_var_list is only called for pe0
    7 ! -dopr_time_count,
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 145 2008-01-09 08:17:38Z raasch
     13! write_var_list is only called for pe0, -dopr_time_count,
    814! hom_sum moved, volume_flow_area, volume_flow_initial moved to write_var_list,
    915! writing of old profil parameters (cross_..., dopr_crossindex, profile_***)
    1016! removed
    11 !
    12 ! Former revisions:
    13 ! -----------------
    14 ! $Id$
    1517!
    1618! 102 2007-07-27 09:09:17Z raasch

  • palm/trunk/SOURCE/write_var_list.f90

    r153 r198  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 153 2008-03-19 09:41:30Z steinfeld
    613! +cthf, leaf_surface_concentration, scalar_exchange_coefficient
    714! +numprocs, hor_index_bounds, inflow_damping_height, inflow_damping_width,
    … …  
    1017! hom_sum, volume_flow_area, volume_flow_initial moved from write_3d_binary
    1118! to here
    12 !
    13 ! Former revisions:
    14 ! -----------------
    15 ! $Id$
    1619!
    1720! 138 2007-11-28 10:03:58Z letzel
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |