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

file headers updated for the next release 3.5

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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 
Note: See TracChangeset for help on using the changeset viewer.