| [5] | 1 | <!DOCTYPE html PUBLIC "-//w3c//dtd html 4.0 transitional//en"> |
|---|
| 2 | <html><head> |
|---|
| [392] | 3 | |
|---|
| [481] | 4 | |
|---|
| 5 | |
|---|
| [482] | 6 | |
|---|
| 7 | |
|---|
| [392] | 8 | <meta content="text/html; charset=ISO-8859-1" http-equiv="Content-Type"> <meta content="Mozilla/4.73 [en] (X11; U; HP-UX B.10.20 9000/782) [Netscape]" name="GENERATOR"> <meta content="Siegfried Raasch" name="Author"> <meta content="parallel LES model" name="Keywords"><title>PALM technical/numerical documentation</title></head><body><font size="+4">PALM |
|---|
| [5] | 9 | </font>a |
|---|
| 10 | <b>pa</b>rallelized |
|---|
| 11 | <b>L</b>ES <b>m</b>odel |
|---|
| [62] | 12 | <br><hr width="100%"> |
|---|
| 13 | <br><font size="+2">technical/numerical documentation</font> |
|---|
| 14 | <br> |
|---|
| 15 | <br> |
|---|
| 16 | <br><table nosave="" cellpadding="0" cellspacing="0"> <caption> <br> </caption><tbody> |
|---|
| [482] | 17 | </tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Current fixed release:</b></td> <td><span style="font-weight: bold;">3.7a<br> |
|---|
| [481] | 18 | </span></td> </tr> |
|---|
| [62] | 19 | <tr nosave=""> <td nosave=""><b>Last |
|---|
| 20 | change of this document</b>: </td> <td nosave=""><b>$Id: technical_documentation.html 482 2010-02-05 06:37:03Z raasch $</b></td> </tr> |
|---|
| 21 | </tbody> |
|---|
| 22 | </table><hr width="100%"> |
|---|
| [5] | 23 | <p><font size="+2">Table of contents</font> |
|---|
| [62] | 24 | <br> |
|---|
| 25 | <br> |
|---|
| 26 | </p><ul> <a href="#Kapitel1.0">1.0</a> |
|---|
| 27 | History of changes <p><a href="#Kapitel2.0">2.0</a> |
|---|
| 28 | How to change the model source |
|---|
| 29 | code </p> <p><a href="#Kapitel3.0">3.0</a> |
|---|
| [481] | 30 | Description of selected model numerics<br> |
|---|
| [5] | 31 | </p> |
|---|
| [392] | 32 | </ul><b></b><br> |
|---|
| 33 | |
|---|
| [62] | 34 | <h2><a name="Kapitel1.0"></a>1.0 |
|---|
| 35 | History of changes</h2> |
|---|
| [481] | 36 | The following table documents all changes to the PALM repository since 2009. Older changes can be found <a href="tec_doc_1998-2008.html">here</a>.<br> |
|---|
| 37 | <br> |
|---|
| [482] | 38 | Entries in the fifth column specify the type of the change: |
|---|
| [5] | 39 | change |
|---|
| 40 | of existing code (C), new code (N) or bugfix (B). More detailed |
|---|
| 41 | information |
|---|
| [481] | 42 | about the changes may sometimes be found in the header comment lines of the |
|---|
| [5] | 43 | respective |
|---|
| 44 | routines. |
|---|
| [62] | 45 | <br> |
|---|
| 46 | <br> |
|---|
| [481] | 47 | <table style="text-align: left; width: 100%; vertical-align: top;" border="1" cellpadding="2" cellspacing="2"> |
|---|
| 48 | <tbody> |
|---|
| 49 | <tr> |
|---|
| [482] | 50 | <td style="vertical-align: middle; width: 27px;"><big><span style="font-weight: bold;">date</span></big><br> |
|---|
| [392] | 51 | </td> |
|---|
| [481] | 52 | <td style="width: 38px;"><big><span style="font-weight: bold;">author</span></big><br> |
|---|
| [392] | 53 | </td> |
|---|
| [481] | 54 | <td style="width: 42px;"><big><span style="font-weight: bold;">svn revision</span></big><br> |
|---|
| [392] | 55 | </td> |
|---|
| [481] | 56 | <td style="width: 40px;"><big><span style="font-weight: bold;">latest release</span></big><br> |
|---|
| [392] | 57 | </td> |
|---|
| [481] | 58 | <td style="width: 27px;"><big><span style="font-weight: bold;">type</span></big><br> |
|---|
| 59 | </td> |
|---|
| 60 | <td style="width: 497px;"><big><span style="font-weight: bold;">description</span></big><br> |
|---|
| 61 | </td> |
|---|
| 62 | <td style="width: 189px;"><big><span style="font-weight: bold;">affected routines</span></big><br> |
|---|
| 63 | </td> |
|---|
| 64 | </tr> |
|---|
| [482] | 65 | |
|---|
| 66 | |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | |
|---|
| 70 | |
|---|
| 71 | |
|---|
| 72 | |
|---|
| 73 | |
|---|
| 74 | |
|---|
| [481] | 75 | <tr> |
|---|
| [482] | 76 | <td style="vertical-align: top; width: 27px;">05/02/10<br> |
|---|
| [481] | 77 | </td> |
|---|
| [482] | 78 | <td style="vertical-align: top; width: 38px;">SR<br> |
|---|
| [481] | 79 | </td> |
|---|
| [482] | 80 | <td style="vertical-align: top; width: 42px;">483<br> |
|---|
| [481] | 81 | </td> |
|---|
| [482] | 82 | <td style="vertical-align: top; width: 40px;">3.7a<br> |
|---|
| [481] | 83 | </td> |
|---|
| [482] | 84 | <td style="vertical-align: top; width: 27px;">N<br> |
|---|
| [481] | 85 | </td> |
|---|
| [482] | 86 | <td style="vertical-align: top; width: 497px;"><span style="font-family: Courier New,Courier,monospace;">compare_palm_logs </span>is additionally compiled with <span style="font-weight: bold;">mbuild</span> <span style="font-family: Courier New,Courier,monospace;">-u</span> . (<span style="font-family: Courier New,Courier,monospace;">Makefile</span> in <span style="font-family: Courier New,Courier,monospace;">trunk/UTIL</span>)<br> |
|---|
| 87 | <br> |
|---|
| 88 | make options (<span style="font-family: Courier New,Courier,monospace;">mopts</span>) to be set by configuration file implemented. (<span style="font-weight: bold;">mrun</span>, <span style="font-weight: bold;">mbuild</span>)<br> |
|---|
| 89 | <br> |
|---|
| 90 | <span style="font-family: Courier New,Courier,monospace;">humidity</span> = <span style="font-style: italic;">.T.</span> can be used for runs with topography. <span style="font-family: Courier New,Courier,monospace;">wall_humidityflux</span> and<br><span style="font-family: Courier New,Courier,monospace;"> |
|---|
| 91 | wall_scalarflux</span> are the corresponding new inipar parameters.<br> |
|---|
| 92 | (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>)<br> |
|---|
| 93 | <br> |
|---|
| 94 | Large scale vertical motion (subsidence/ascent) can be added as |
|---|
| 95 | additional source to the prognostic equation for potential temperature. |
|---|
| 96 | See new parameters <span style="font-family: Courier New,Courier,monospace;">large_scale_subsidence</span>, <span style="font-family: Courier New,Courier,monospace;">ws_vertical_gradient</span>, and <span style="font-family: Courier New,Courier,monospace;">ws_vertical_gradient_level</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, <span style="font-family: Courier New,Courier,monospace;">header</span>, <span style="font-family: Courier New,Courier,monospace;">Makefile</span>, <span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">parin</span>, <span style="font-family: Courier New,Courier,monospace;">prognostic_equations</span>, <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, <span style="font-family: Courier New,Courier,monospace;">subsidence</span>, <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 97 | <br> |
|---|
| 98 | A simple method for installing and running palm (with limited features) has been added. (<span style="font-family: Courier New,Courier,monospace;">Makefile</span>, new scripts <span style="font-weight: bold;">palm_simple_install</span>, <span style="font-weight: bold;">palm_simple_run</span>)<br> |
|---|
| 99 | <br> |
|---|
| 100 | Masked data output has been implemented but it is not documented in the |
|---|
| 101 | user manual so far. This feature can already use NetCDF4 file format. |
|---|
| 102 | (check_open, check_parameters, close_file, modules, netcdf, parin, |
|---|
| 103 | time_integration, user_parin, new: data_output_mask, init_masks, |
|---|
| 104 | user_data_output_mask)<br> |
|---|
| [481] | 105 | </td> |
|---|
| [482] | 106 | <td style="vertical-align: top; width: 189px;">check_open, check |
|---|
| 107 | parameters, close_file, header, init_3d_model, Makefile, mbuild, modules, mrun, |
|---|
| 108 | netcdf, parin, prognostic_equations, read_var_list, time_integration, user_parin, write_var_list<br> |
|---|
| 109 | <br><span style="font-weight: bold;"> |
|---|
| 110 | new:</span> data_output_mask, init_masks, palm_simple_install, palm_simple_run, subsidence, user_data_output_mask<br> |
|---|
| 111 | <br> |
|---|
| [481] | 112 | </td> |
|---|
| 113 | </tr> |
|---|
| 114 | <tr> |
|---|
| 115 | <td style="vertical-align: top; width: 27px;"><br> |
|---|
| 116 | </td> |
|---|
| 117 | <td style="vertical-align: top; width: 38px;"><br> |
|---|
| 118 | </td> |
|---|
| 119 | <td style="vertical-align: top; width: 42px;"><br> |
|---|
| 120 | </td> |
|---|
| 121 | <td style="vertical-align: top; width: 40px;"><br> |
|---|
| 122 | </td> |
|---|
| [482] | 123 | <td style="vertical-align: top; width: 27px;">C<br> |
|---|
| [481] | 124 | </td> |
|---|
| [482] | 125 | <td style="vertical-align: top; width: 497px;">2d-decomposition is default for Cray-XT machines. (<span style="font-family: Courier New,Courier,monospace;">init_pegrid</span>)<br> |
|---|
| 126 | <br><span style="font-family: Courier New,Courier,monospace;"> |
|---|
| 127 | var_ts</span> is replaced by <span style="font-family: Courier New,Courier,monospace;">dots_max</span>. (<span style="font-family: Courier New,Courier,monospace;">modules</span>, <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> |
|---|
| 128 | <br> |
|---|
| 129 | Every cloud droplet has now an own weighting factor and can be deleted |
|---|
| 130 | due to collisions. Condensation and collision of cloud droplets are |
|---|
| 131 | adjusted accordingly. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br> |
|---|
| 132 | <br> |
|---|
| 133 | Collision efficiency for large cloud droplets has changed according to table of Rogers and Yau. (<span style="font-family: Courier New,Courier,monospace;">collision_efficiency</span>)<br> |
|---|
| [481] | 134 | </td> |
|---|
| [482] | 135 | <td style="vertical-align: top; width: 189px;">advec_particles, collision_efficiency, init_3d_model, init_pegrid, modules<br> |
|---|
| [481] | 136 | </td> |
|---|
| 137 | </tr> |
|---|
| 138 | <tr> |
|---|
| 139 | <td style="vertical-align: top; width: 27px;"><br> |
|---|
| 140 | </td> |
|---|
| 141 | <td style="vertical-align: top; width: 38px;"><br> |
|---|
| 142 | </td> |
|---|
| 143 | <td style="vertical-align: top; width: 42px;"><br> |
|---|
| 144 | </td> |
|---|
| 145 | <td style="vertical-align: top; width: 40px;"><br> |
|---|
| 146 | </td> |
|---|
| [482] | 147 | <td style="vertical-align: top; width: 27px;">B<br> |
|---|
| [481] | 148 | </td> |
|---|
| [482] | 149 | <td style="vertical-align: top; width: 497px;">Bugfix for generating serial jobs. (<span style="font-weight: bold;">subjob</span>)<br> |
|---|
| 150 | <br> |
|---|
| 151 | Bugfix: index problem concerning gradient_level indices removed. (<span style="font-family: Courier New,Courier,monospace;">header</span>)<br> |
|---|
| 152 | <br> |
|---|
| 153 | Dimension of array <span style="font-family: Courier New,Courier,monospace;">stat</span> in cascade change to prevent type problems with mpi2 libraries. (<span style="font-family: Courier New,Courier,monospace;">poisfft_hybrid</span>)<br> |
|---|
| 154 | <br> |
|---|
| 155 | Loop was split to make runs reproducible when using ifort compiler. (<span style="font-family: Courier New,Courier,monospace;">disturb_field</span>)<br> |
|---|
| 156 | <br> |
|---|
| 157 | Bugfix: exchange of ghost points for <span style="font-family: Courier New,Courier,monospace;">prho</span> included. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)<br> |
|---|
| 158 | <br> |
|---|
| 159 | Bugfix in calculation of time-averaged surface heatfluxes. (<span style="font-family: Courier New,Courier,monospace;">sum_up_3d_data</span>)<br> |
|---|
| 160 | <br> |
|---|
| 161 | Bugfix in calculation of precipitation_rate. (<span style="font-family: Courier New,Courier,monospace;">calc_precipitation</span>)<br> |
|---|
| 162 | <br> |
|---|
| 163 | Bugfix: initial data assignments to some dvrp arrays changed due to error messages from gfortran compiler. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> |
|---|
| 164 | <br> |
|---|
| 165 | Bugfix in calculation of cloud droplet velocity. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br> |
|---|
| 166 | <br> |
|---|
| 167 | Bugfix in transfer of particles at south/left edge. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br> |
|---|
| 168 | <br> |
|---|
| 169 | Bugfix in calculation of collision efficiency. (<span style="font-family: Courier New,Courier,monospace;">collision_efficiency</span>)<span style="font-family: Courier New,Courier,monospace;"></span><span style="font-family: Courier New,Courier,monospace;"></span><br> |
|---|
| [481] | 170 | </td> |
|---|
| [482] | 171 | <td style="vertical-align: top; width: 189px;">advec_particles, |
|---|
| 172 | calc_precipitation, collision_efficiency, disturb_field, header, |
|---|
| 173 | modules, poisfft_hybrid, subjob, sum_up_3d_data, |
|---|
| 174 | time_integration<br> |
|---|
| 175 | <br> |
|---|
| [481] | 176 | </td> |
|---|
| 177 | </tr> |
|---|
| 178 | <tr> |
|---|
| 179 | <td style="vertical-align: top; width: 27px;">23/09/09<br> |
|---|
| 180 | </td> |
|---|
| 181 | <td style="vertical-align: top; width: 38px;">SR<br> |
|---|
| 182 | </td> |
|---|
| 183 | <td style="vertical-align: top; width: 42px;">392<br> |
|---|
| 184 | </td> |
|---|
| 185 | <td style="vertical-align: top; width: 40px;">3.7<br> |
|---|
| 186 | </td> |
|---|
| 187 | <td style="vertical-align: top; width: 27px;">N<br> |
|---|
| 188 | </td> |
|---|
| 189 | <td style="vertical-align: top; width: 497px;"><span style="font-weight: bold;">Most |
|---|
| 190 | important changes:</span> |
|---|
| [392] | 191 | atmosphere-ocean coupling has been improved and steering of |
|---|
| 192 | dvr-software is more user-friendly now. Beside that, the release |
|---|
| 193 | contains a lot of optimizations, changes, and bugfixes which have |
|---|
| 194 | accumulated over the last months.<br> |
|---|
| [481] | 195 | |
|---|
| [392] | 196 | <br> |
|---|
| [481] | 197 | |
|---|
| [392] | 198 | The atmosphere-ocean coupling allows to use independent precursor runs |
|---|
| 199 | in order to account for different spin-up times. The time when coupling |
|---|
| [481] | 200 | has to be started is given by new inipar parameter <span style="font-family: Courier New,Courier,monospace;">coupling_start_time</span>. |
|---|
| 201 | The precursor ocean run has to be started using new <span style="font-weight: bold;">mrun</span> option "<span style="font-family: Courier New,Courier,monospace;">-y</span>" in |
|---|
| 202 | order to add the appendix "<span style="font-family: Courier New,Courier,monospace;">_O</span>" to all |
|---|
| 203 | output files. (<span style="font-family: Courier New,Courier,monospace;">check_for_restart</span>, |
|---|
| 204 | <span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 205 | <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, |
|---|
| 206 | <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>, |
|---|
| 207 | <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, |
|---|
| 208 | <span style="font-family: Courier New,Courier,monospace;">data_output_ptseries</span>, |
|---|
| 209 | <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>, |
|---|
| 210 | <span style="font-family: Courier New,Courier,monospace;">data_output_tseries</span>, |
|---|
| 211 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 212 | <span style="font-family: Courier New,Courier,monospace;">init_coupling</span>, |
|---|
| 213 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 214 | <span style="font-family: Courier New,Courier,monospace;">mrun</span>, |
|---|
| 215 | <span style="font-family: Courier New,Courier,monospace;">parin</span>, |
|---|
| 216 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 217 | <span style="font-family: Courier New,Courier,monospace;">surface_coupler</span>, |
|---|
| 218 | <span style="font-family: Courier New,Courier,monospace;">time_integration</span>, |
|---|
| 219 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| [392] | 220 | |
|---|
| 221 | <br> |
|---|
| [481] | 222 | |
|---|
| [392] | 223 | Clipping of dvrp output implemented. Default colourtable for particles |
|---|
| 224 | implemented, particle attributes (color, dvrp_size) can be set with new |
|---|
| [481] | 225 | parameters <span style="font-family: Courier New,Courier,monospace;">particle_color</span>, |
|---|
| 226 | <span style="font-family: Courier New,Courier,monospace;">particle_dvrpsize</span>, |
|---|
| 227 | <span style="font-family: Courier New,Courier,monospace;">color_interval</span>, |
|---|
| 228 | <span style="font-family: Courier New,Courier,monospace;">dvrpsize_interval</span>. |
|---|
| 229 | Slicer attributes (dvrp) are set with new routine <span style="font-family: Courier New,Courier,monospace;">set_slicer_attributes_dvrp</span> |
|---|
| 230 | and are controlled with existing parameters <span style="font-family: Courier New,Courier,monospace;">slicer_range_limits</span>. |
|---|
| 231 | (<span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>, |
|---|
| 232 | <span style="font-family: Courier New,Courier,monospace;">data_output_dvrp</span>, |
|---|
| 233 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 234 | <span style="font-family: Courier New,Courier,monospace;">set_slicer_attributes_dvrp</span>, |
|---|
| 235 | <span style="font-family: Courier New,Courier,monospace;">user_data_output_dvrp</span>)<br> |
|---|
| [392] | 236 | |
|---|
| 237 | <br> |
|---|
| [481] | 238 | |
|---|
| [392] | 239 | Polygon reduction for topography and ground plate isosurface. Reduction |
|---|
| 240 | level for buildings can be chosen with parameter cluster_size. (<span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>)<br> |
|---|
| 241 | |
|---|
| 242 | <br> |
|---|
| 243 | |
|---|
| [481] | 244 | Variables <span style="font-family: Courier New,Courier,monospace;">bc_lr</span> |
|---|
| 245 | / <span style="font-family: Courier New,Courier,monospace;">bc_ns</span> |
|---|
| 246 | in most subroutines replaced by LOGICAL variables <span style="font-family: Courier New,Courier,monospace;">bc_lr_cyc</span>,<br> |
|---|
| 247 | |
|---|
| 248 | <span style="font-family: Courier New,Courier,monospace;">bc_ns_cyc</span> |
|---|
| 249 | for speed optimization. This gives a significant speedup of more than |
|---|
| 250 | 10%. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 251 | <span style="font-family: Courier New,Courier,monospace;">diffusion_u</span>, |
|---|
| 252 | <span style="font-family: Courier New,Courier,monospace;">diffusion_v</span>, |
|---|
| 253 | <span style="font-family: Courier New,Courier,monospace;">diffusion_w</span>, |
|---|
| 254 | <span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> |
|---|
| 255 | |
|---|
| [392] | 256 | <br> |
|---|
| [481] | 257 | |
|---|
| [392] | 258 | Scripts have been adapted for machine lck (Yonsei Univ.). (<span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">mbuild</span>, <span style="font-family: Courier New,Courier,monospace;">subjob</span>)<br> |
|---|
| [481] | 259 | |
|---|
| [392] | 260 | <br> |
|---|
| [481] | 261 | |
|---|
| 262 | Additional timestep criterion in case of simulations with plant canopy. |
|---|
| 263 | (<span style="font-family: Courier New,Courier,monospace;">timestep</span>)<br> |
|---|
| 264 | |
|---|
| [392] | 265 | <br> |
|---|
| [481] | 266 | |
|---|
| 267 | Check for illegal entries in <span style="font-family: Courier New,Courier,monospace;">section_xy</span>|<span style="font-family: Courier New,Courier,monospace;">xz</span>|<span style="font-family: Courier New,Courier,monospace;">yz</span> that |
|---|
| 268 | exceed <span style="font-family: Courier New,Courier,monospace;">nz+1</span>|<span style="font-family: Courier New,Courier,monospace;">ny+1</span>|<span style="font-family: Courier New,Courier,monospace;">nx+1</span>. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br> |
|---|
| 269 | |
|---|
| [392] | 270 | <br> |
|---|
| [481] | 271 | |
|---|
| 272 | External pressure gradient can be used as driving force. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 273 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 274 | <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, |
|---|
| 275 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 276 | <span style="font-family: Courier New,Courier,monospace;">parin</span>, |
|---|
| 277 | <span style="font-family: Courier New,Courier,monospace;">prognostic_equations</span>, |
|---|
| 278 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 279 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 280 | |
|---|
| [392] | 281 | <br> |
|---|
| [481] | 282 | |
|---|
| 283 | New topography case '<span style="font-style: italic;">single_street_canyon</span>'. |
|---|
| 284 | (<span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 285 | <span style="font-family: Courier New,Courier,monospace;">init_grid</span>, |
|---|
| 286 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 287 | <span style="font-family: Courier New,Courier,monospace;">parin</span>, |
|---|
| 288 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 289 | <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span>, |
|---|
| 290 | <span style="font-family: Courier New,Courier,monospace;">user_header</span>, |
|---|
| 291 | <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>, |
|---|
| 292 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 293 | |
|---|
| [392] | 294 | <br> |
|---|
| [481] | 295 | |
|---|
| 296 | Option to predefine a target bulk velocity for <span style="font-family: Courier New,Courier,monospace;">conserve_volume_flow</span>. |
|---|
| 297 | (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 298 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 299 | <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, |
|---|
| 300 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 301 | <span style="font-family: Courier New,Courier,monospace;">parin</span>, |
|---|
| 302 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 303 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 304 | |
|---|
| [392] | 305 | <br> |
|---|
| [481] | 306 | |
|---|
| 307 | Option for user defined 2D data output in xy cross sections at z=<span style="font-family: Courier New,Courier,monospace;">nzb</span>+1. (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, |
|---|
| 308 | <span style="font-family: Courier New,Courier,monospace;">user_data_output_2d</span>)<br> |
|---|
| 309 | |
|---|
| [392] | 310 | <br> |
|---|
| [481] | 311 | |
|---|
| 312 | xy cross section output of surface heatfluxes (sensible and latent). (<span style="font-family: Courier New,Courier,monospace;">average_3d_data</span>, |
|---|
| 313 | <span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 314 | <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, |
|---|
| 315 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 316 | <span style="font-family: Courier New,Courier,monospace;">read_3d_binary</span>, |
|---|
| 317 | <span style="font-family: Courier New,Courier,monospace;">sum_up_3d_data</span>, |
|---|
| 318 | <span style="font-family: Courier New,Courier,monospace;">write_3d_binary</span>)</td> |
|---|
| 319 | <td style="vertical-align: top; width: 189px;">average_3d_data, |
|---|
| [392] | 320 | check_for_restart, check_parameters, data_output_2d, data_output_3d, |
|---|
| 321 | data_output_dvrp, data_output_profiles, data_output_ptseries, |
|---|
| 322 | data_output_spectra, data_output_tseries, diffusion_u, diffusion_v, |
|---|
| 323 | diffusion_w, init_coupling, init_dvrp, init_grid, init_3d_model, |
|---|
| 324 | header, mbuild, modules, mrun, package_parin, parin, |
|---|
| 325 | prognostic_equations, read_3d_binary, read_var_list, subjob, |
|---|
| 326 | sum_up_3d_data, surface_coupler, timestep, time_integration, |
|---|
| 327 | user_check_parameters, user_data_output_2d, user_data_output_dvrp, |
|---|
| 328 | user_header, user_init_grid, write_3d_binary, write_var_list<br> |
|---|
| [481] | 329 | |
|---|
| [392] | 330 | <br> |
|---|
| [481] | 331 | |
|---|
| 332 | <span style="font-weight: bold;">new:</span> |
|---|
| 333 | set_particle_attributes, set_slicer_attributes_dvrp</td> |
|---|
| [392] | 334 | </tr> |
|---|
| 335 | <tr> |
|---|
| [481] | 336 | <td style="vertical-align: top; width: 27px;"><br> |
|---|
| [392] | 337 | </td> |
|---|
| [481] | 338 | <td style="vertical-align: top; width: 38px;"><br> |
|---|
| [392] | 339 | </td> |
|---|
| [481] | 340 | <td style="vertical-align: top; width: 42px;"><br> |
|---|
| [392] | 341 | </td> |
|---|
| [481] | 342 | <td style="vertical-align: top; width: 40px;"><br> |
|---|
| [392] | 343 | </td> |
|---|
| [481] | 344 | <td style="vertical-align: top; width: 27px;">C<br> |
|---|
| 345 | </td> |
|---|
| 346 | <td style="vertical-align: top; width: 497px;">Output of messages replaced by |
|---|
| 347 | message handling routine. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>, |
|---|
| 348 | <span style="font-family: Courier New,Courier,monospace;">advec_s_bc</span>, |
|---|
| 349 | <span style="font-family: Courier New,Courier,monospace;">buoyancy</span>, |
|---|
| 350 | <span style="font-family: Courier New,Courier,monospace;">calc_spectra</span>, |
|---|
| 351 | <span style="font-family: Courier New,Courier,monospace;">check_for_restart</span>, |
|---|
| 352 | <span style="font-family: Courier New,Courier,monospace;">check_open</span>, |
|---|
| 353 | <span style="font-family: Courier New,Courier,monospace;">coriolis</span>, |
|---|
| 354 | <span style="font-family: Courier New,Courier,monospace;">cpu_log</span>, |
|---|
| 355 | <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, |
|---|
| 356 | <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>, |
|---|
| 357 | <span style="font-family: Courier New,Courier,monospace;">data_output_dvrp</span>, |
|---|
| 358 | <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, |
|---|
| 359 | <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>, |
|---|
| 360 | <span style="font-family: Courier New,Courier,monospace;">fft_xy</span>, |
|---|
| 361 | <span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>, |
|---|
| 362 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 363 | <span style="font-family: Courier New,Courier,monospace;">init_1d_model</span>, |
|---|
| 364 | <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>, |
|---|
| 365 | <span style="font-family: Courier New,Courier,monospace;">init_dvrp</span>, |
|---|
| 366 | <span style="font-family: Courier New,Courier,monospace;">init_grid</span>, |
|---|
| 367 | <span style="font-family: Courier New,Courier,monospace;">init_particles</span>, |
|---|
| 368 | <span style="font-family: Courier New,Courier,monospace;">init_pegrid</span>, |
|---|
| 369 | <span style="font-family: Courier New,Courier,monospace;">netcdf</span>, |
|---|
| 370 | <span style="font-family: Courier New,Courier,monospace;">parin</span>, |
|---|
| 371 | <span style="font-family: Courier New,Courier,monospace;">plant_canopy_model</span>, |
|---|
| 372 | <span style="font-family: Courier New,Courier,monospace;">poisfft_hybrid</span>, |
|---|
| 373 | <span style="font-family: Courier New,Courier,monospace;">poismg</span>, |
|---|
| 374 | <span style="font-family: Courier New,Courier,monospace;">read_3d_binary</span>, |
|---|
| 375 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 376 | <span style="font-family: Courier New,Courier,monospace;">surface_coupler</span>, |
|---|
| 377 | <span style="font-family: Courier New,Courier,monospace;">temperton_fft</span>, |
|---|
| 378 | <span style="font-family: Courier New,Courier,monospace;">timestep</span>, |
|---|
| 379 | <span style="font-family: Courier New,Courier,monospace;">user_actions</span>, |
|---|
| 380 | <span style="font-family: Courier New,Courier,monospace;">user_data_output_dvrp</span>, |
|---|
| 381 | <span style="font-family: Courier New,Courier,monospace;">user_dvrp_coltab</span>, |
|---|
| 382 | <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>, |
|---|
| 383 | <span style="font-family: Courier New,Courier,monospace;">user_init_plant_canopy</span>, |
|---|
| 384 | <span style="font-family: Courier New,Courier,monospace;">user_parin</span>, |
|---|
| 385 | <span style="font-family: Courier New,Courier,monospace;">user_read_restart_data</span>, |
|---|
| 386 | <span style="font-family: Courier New,Courier,monospace;">user_spectra</span>)<br> |
|---|
| [392] | 387 | |
|---|
| 388 | <br> |
|---|
| 389 | |
|---|
| [481] | 390 | Output of NetCDF messages with aid of message handling routine. (<span style="font-family: Courier New,Courier,monospace;">check_open</span>, |
|---|
| 391 | <span style="font-family: Courier New,Courier,monospace;">close_file</span>, |
|---|
| 392 | <span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, |
|---|
| 393 | <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>, |
|---|
| 394 | <span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, |
|---|
| 395 | <span style="font-family: Courier New,Courier,monospace;">data_output_ptseries</span>, |
|---|
| 396 | <span style="font-family: Courier New,Courier,monospace;">data_output_spectra</span>, |
|---|
| 397 | <span style="font-family: Courier New,Courier,monospace;">data_output_tseries</span>, |
|---|
| 398 | <span style="font-family: Courier New,Courier,monospace;">netcdf</span><span style="font-family: Courier New,Courier,monospace;"></span>)<br> |
|---|
| 399 | |
|---|
| [392] | 400 | <br> |
|---|
| 401 | |
|---|
| [481] | 402 | __lcmuk changed to __lc to avoid problems with Intel compiler on |
|---|
| 403 | sgi-ice. (<span style="font-family: Courier New,Courier,monospace;">poisfft</span>)<br> |
|---|
| 404 | |
|---|
| [392] | 405 | <br> |
|---|
| [481] | 406 | |
|---|
| 407 | For extended NetCDF files, the updated title attribute includes an |
|---|
| 408 | update of time_average_text where appropriate. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> |
|---|
| 409 | |
|---|
| [392] | 410 | <br> |
|---|
| [481] | 411 | |
|---|
| [392] | 412 | In case of restart runs without extension, initial profiles are not |
|---|
| 413 | written to NetCDF-file anymore. So far, time axis always started at |
|---|
| [481] | 414 | t=0. (<span style="font-family: Courier New,Courier,monospace;">data_output_profiles</span>, |
|---|
| 415 | <span style="font-family: Courier New,Courier,monospace;">modules</span>, |
|---|
| 416 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 417 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 418 | |
|---|
| [392] | 419 | <br> |
|---|
| [481] | 420 | |
|---|
| 421 | <span style="font-family: Courier New,Courier,monospace;">initializing_actions</span>='<span style="font-style: italic;">read_data_for_recycling</span>' renamed to |
|---|
| 422 | '<span style="font-style: italic;">cyclic_fill</span>'. A cyclic fill |
|---|
| 423 | can now <br> |
|---|
| 424 | |
|---|
| 425 | be used independent of <span style="font-family: Courier New,Courier,monospace;">turbulent_inflow</span>. |
|---|
| 426 | (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 427 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 428 | <span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> |
|---|
| 429 | |
|---|
| [392] | 430 | <br> |
|---|
| [481] | 431 | |
|---|
| [392] | 432 | 2 NetCDF error numbers changed. (<span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>)<br> |
|---|
| [481] | 433 | |
|---|
| [392] | 434 | <br> |
|---|
| [481] | 435 | |
|---|
| [392] | 436 | A link to the webspage appendix_a.html is printed for further |
|---|
| 437 | information about possible reasons of the error that appeared. (message)<br> |
|---|
| [481] | 438 | |
|---|
| [392] | 439 | <br> |
|---|
| [481] | 440 | |
|---|
| [392] | 441 | Temperature gradient criterion for estimating the boundary layer height |
|---|
| 442 | replaced by the gradient criterion of Sullivan et al. (1998). (<span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>)<br> |
|---|
| [481] | 443 | |
|---|
| [392] | 444 | <br> |
|---|
| [481] | 445 | |
|---|
| 446 | NetCDF unit attribute in timeseries output in case of statistic regions |
|---|
| 447 | added. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> |
|---|
| 448 | |
|---|
| [392] | 449 | <br> |
|---|
| [481] | 450 | |
|---|
| 451 | Maximum number of tails is calculated from maximum number of particles |
|---|
| 452 | and <span style="font-family: Courier New,Courier,monospace;">skip_particles_for_tail</span>. |
|---|
| 453 | (<span style="font-family: Courier New,Courier,monospace;">init_particles</span>)<br> |
|---|
| 454 | |
|---|
| [392] | 455 | <br> |
|---|
| [481] | 456 | |
|---|
| 457 | Value of <span style="font-family: Courier New,Courier,monospace;">vertical_particle_advection</span> |
|---|
| 458 | is allowed to differ for each particle group. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>, |
|---|
| 459 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 460 | <span style="font-family: Courier New,Courier,monospace;">modules</span>) |
|---|
| [392] | 461 | <br> |
|---|
| [481] | 462 | |
|---|
| [392] | 463 | <br> |
|---|
| [481] | 464 | |
|---|
| 465 | First constant in array <span style="font-family: Courier New,Courier,monospace;">den</span> also |
|---|
| 466 | defined as type double. (<span style="font-family: Courier New,Courier,monospace;">eqn_state_seawater</span>)<br> |
|---|
| 467 | |
|---|
| [392] | 468 | <br> |
|---|
| [481] | 469 | |
|---|
| 470 | Parameter <span style="font-family: Courier New,Courier,monospace;">dvrp_psize</span> |
|---|
| 471 | moved from <span style="font-family: Courier New,Courier,monospace;">particles_par</span> |
|---|
| 472 | to <span style="font-family: Courier New,Courier,monospace;">dvrp_graphics_par</span>. |
|---|
| 473 | (<span style="font-family: Courier New,Courier,monospace;">package_parin</span>)<br> |
|---|
| 474 | |
|---|
| [392] | 475 | <br> |
|---|
| [481] | 476 | |
|---|
| 477 | <span style="font-family: Courier New,Courier,monospace;">topography_grid_convention</span> |
|---|
| 478 | moved from <span style="font-family: Courier New,Courier,monospace;">userpar</span> |
|---|
| 479 | to <span style="font-family: Courier New,Courier,monospace;">inipar</span> |
|---|
| 480 | (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>, |
|---|
| 481 | <span style="font-family: Courier New,Courier,monospace;">header</span>, |
|---|
| 482 | <span style="font-family: Courier New,Courier,monospace;">parin</span>, |
|---|
| 483 | <span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 484 | <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span>, |
|---|
| 485 | <span style="font-family: Courier New,Courier,monospace;">user_header</span>, |
|---|
| 486 | <span style="font-family: Courier New,Courier,monospace;">user_init_grid</span>, |
|---|
| 487 | <span style="font-family: Courier New,Courier,monospace;">user_parin</span>, |
|---|
| 488 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 489 | |
|---|
| [392] | 490 | <br> |
|---|
| [481] | 491 | |
|---|
| 492 | Default value of <span style="font-family: Courier New,Courier,monospace;">grid_matching</span> |
|---|
| 493 | changed to '<span style="font-style: italic;">strict</span>'. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> |
|---|
| 494 | |
|---|
| 495 | <br> |
|---|
| 496 | |
|---|
| 497 | Adjustments for runs on lcxt4 (necessary due to a software update on |
|---|
| 498 | CRAY) and for coupled runs on ibmy. (<span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">subjob</span>)</td> |
|---|
| 499 | <td style="vertical-align: top; width: 189px;">advec_particles, |
|---|
| 500 | advec_s_bc, |
|---|
| [392] | 501 | buoyancy, calc_spectra, check_for_restart, check_open, |
|---|
| 502 | check_parameters, close_file, coriolis, cpu_log, data_output_2d, |
|---|
| 503 | data_output_3d, data_output_dvrp, data_output_profiles, |
|---|
| 504 | data_output_ptseries, data_output_spectra, data_output_tseries, |
|---|
| 505 | eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, |
|---|
| 506 | init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, |
|---|
| 507 | message, mrun, netcdf, package_parin, parin, plant_canopy_model, |
|---|
| 508 | poisfft, poisfft_hybrid, poismg, read_3d_binary, read_var_list, subjob, |
|---|
| [481] | 509 | surface_coupler, temperton_fft, timestep, |
|---|
| 510 | user_actions,user_check_parameters, user_data_output_dvrp, |
|---|
| 511 | user_dvrp_coltab, user_header, user_init_grid, user_init_plant_canopy, |
|---|
| 512 | user_parin, user_read_restart_data, user_spectra, write_var_list </td> |
|---|
| [392] | 513 | </tr> |
|---|
| 514 | <tr> |
|---|
| [481] | 515 | <td style="vertical-align: top; width: 27px;"><br> |
|---|
| [392] | 516 | </td> |
|---|
| [481] | 517 | <td style="vertical-align: top; width: 38px;"><br> |
|---|
| [392] | 518 | </td> |
|---|
| [481] | 519 | <td style="vertical-align: top; width: 42px;"><br> |
|---|
| [392] | 520 | </td> |
|---|
| [481] | 521 | <td style="vertical-align: top; width: 40px;"><br> |
|---|
| [392] | 522 | </td> |
|---|
| [482] | 523 | <td style="vertical-align: top; width: 27px;">B<br> |
|---|
| [481] | 524 | </td> |
|---|
| 525 | <td style="vertical-align: top; width: 497px;">Bugfix: Initial hydrostatic |
|---|
| 526 | pressure profile in case of ocean runs is now calculated in 5 iteration |
|---|
| 527 | steps. (<span style="font-family: Courier New,Courier,monospace;">init_ocean</span>)<br> |
|---|
| 528 | |
|---|
| [392] | 529 | <br> |
|---|
| [481] | 530 | |
|---|
| 531 | Bugfix: wrong sign in buoyancy production of ocean part in case of not |
|---|
| 532 | using the reference density (only in 3D routine <span style="font-family: Courier New,Courier,monospace;">production_e</span>). |
|---|
| 533 | (<span style="font-family: Courier New,Courier,monospace;">production_e</span>)<br> |
|---|
| 534 | |
|---|
| [392] | 535 | <br> |
|---|
| [481] | 536 | |
|---|
| [392] | 537 | Bugfix: output of averaged 2d/3d quantities requires that an avaraging |
|---|
| 538 | interval has been set, respective error message is included. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br> |
|---|
| [481] | 539 | |
|---|
| [392] | 540 | <br> |
|---|
| [481] | 541 | |
|---|
| 542 | Bugfix: Output on unit 14 only if requested by <span style="font-family: Courier New,Courier,monospace;">write_binary</span>. |
|---|
| 543 | (<span style="font-family: Courier New,Courier,monospace;">user_last_actions</span>)<br> |
|---|
| 544 | |
|---|
| [392] | 545 | <br> |
|---|
| [481] | 546 | |
|---|
| 547 | Bugfix to avoid zero division by <span style="font-family: Courier New,Courier,monospace;">km_neutral</span>. |
|---|
| 548 | (<span style="font-family: Courier New,Courier,monospace;">production_e</span>)<br> |
|---|
| 549 | |
|---|
| [392] | 550 | <br> |
|---|
| [481] | 551 | |
|---|
| [392] | 552 | Bugfix for extended NetCDF files: In order to avoid 'data mode' errors |
|---|
| [481] | 553 | if updated attributes are larger than their original size, <span style="font-family: Courier New,Courier,monospace;">NF90_PUT_ATT</span> |
|---|
| 554 | is called in 'define mode' enclosed by <span style="font-family: Courier New,Courier,monospace;">NF90_REDEF</span> |
|---|
| 555 | and <span style="font-family: Courier New,Courier,monospace;">NF90_ENDDEF</span> |
|---|
| [392] | 556 | calls. This implies a possible performance loss; an alternative |
|---|
| 557 | strategy would be to ensure equal attribute size in a job chain. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> |
|---|
| [481] | 558 | |
|---|
| [392] | 559 | <br> |
|---|
| [481] | 560 | |
|---|
| [392] | 561 | Bugfix: correction of initial volume flow for non-flat topography. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> |
|---|
| [481] | 562 | |
|---|
| [392] | 563 | <br> |
|---|
| [481] | 564 | |
|---|
| [392] | 565 | Bugfix: zero initialization of arrays within buildings for '<span style="font-style: italic;">cyclic_fill</span>'. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> |
|---|
| [481] | 566 | |
|---|
| [392] | 567 | <br> |
|---|
| [481] | 568 | |
|---|
| 569 | Bugfix: <span style="font-family: Courier New,Courier,monospace;">to_be_resorted</span> |
|---|
| 570 | => <span style="font-family: Courier New,Courier,monospace;">s_av</span> |
|---|
| 571 | for time-averaged scalars. (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>, |
|---|
| 572 | <span style="font-family: Courier New,Courier,monospace;">data_output_3d</span>)<br> |
|---|
| 573 | |
|---|
| [392] | 574 | <br> |
|---|
| [481] | 575 | |
|---|
| 576 | Bugfix: avoid that <span style="font-family: Courier New,Courier,monospace;">ngp_2dh_s_inner</span> |
|---|
| 577 | becomes zero. (<span style="font-family: Courier New,Courier,monospace;">init_3_model</span>)<br> |
|---|
| 578 | |
|---|
| [392] | 579 | <br> |
|---|
| [481] | 580 | |
|---|
| [392] | 581 | Typographical error: unit of wpt in dots_unit. (<span style="font-family: Courier New,Courier,monospace;">modules</span>)<br> |
|---|
| [481] | 582 | |
|---|
| [392] | 583 | <br> |
|---|
| [481] | 584 | |
|---|
| [392] | 585 | Bugfix: error in check, if particles moved further than one subdomain |
|---|
| 586 | length. This check must not be applied for newly released particles. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>)<br> |
|---|
| [481] | 587 | |
|---|
| [392] | 588 | <br> |
|---|
| [481] | 589 | |
|---|
| [392] | 590 | Bugfix: several tail counters are initialized, <span style="font-family: Courier New,Courier,monospace;">particle_tail_coordinates</span> |
|---|
| 591 | is only written to file if its third index is > 0, arrays for tails |
|---|
| 592 | are allocated with a minimum size of 10 tails if there is no tail |
|---|
| [481] | 593 | initially. (<span style="font-family: Courier New,Courier,monospace;">advec_particles</span>, |
|---|
| 594 | <span style="font-family: Courier New,Courier,monospace;">init_particles</span><span style="font-family: Courier New,Courier,monospace;"></span>)<br> |
|---|
| 595 | |
|---|
| [392] | 596 | <br> |
|---|
| [481] | 597 | |
|---|
| [392] | 598 | Bugfix: pressure included for profile output. (<span style="font-family: Courier New,Courier,monospace;">check_parameters</span>)<br> |
|---|
| [481] | 599 | |
|---|
| [392] | 600 | <br> |
|---|
| [481] | 601 | |
|---|
| 602 | Bugfix: Type of <span style="font-family: Courier New,Courier,monospace;">count</span> and <span style="font-family: Courier New,Courier,monospace;">count_rate</span> |
|---|
| 603 | changed to default <span style="font-family: Courier New,Courier,monospace;">INTEGER</span> on |
|---|
| 604 | NEC machines. (<span style="font-family: Courier New,Courier,monospace;">cpu_log</span>)<br> |
|---|
| 605 | |
|---|
| [392] | 606 | <br> |
|---|
| [481] | 607 | |
|---|
| 608 | Bugfix: output of particle time series only if particle advection is |
|---|
| 609 | switched on. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)<br> |
|---|
| 610 | |
|---|
| [392] | 611 | <br> |
|---|
| [481] | 612 | |
|---|
| 613 | Bugfix: qsws was calculated in case of <span style="font-family: Courier New,Courier,monospace;">constant heatflux</span> |
|---|
| 614 | = <span style="font-family: Courier New,Courier,monospace;">.FALSE.</span>. |
|---|
| 615 | (<span style="font-family: Courier New,Courier,monospace;">prandtl_fluxes</span>)<br> |
|---|
| 616 | |
|---|
| [392] | 617 | <br> |
|---|
| [481] | 618 | |
|---|
| 619 | Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and |
|---|
| 620 | z0). (<span style="font-family: Courier New,Courier,monospace;">data_output_2d</span>)<br> |
|---|
| 621 | |
|---|
| [392] | 622 | <br> |
|---|
| [481] | 623 | |
|---|
| [392] | 624 | Typographical errors. (<span style="font-family: Courier New,Courier,monospace;">netcdf</span>)<br> |
|---|
| [481] | 625 | |
|---|
| [392] | 626 | <br> |
|---|
| [481] | 627 | |
|---|
| 628 | If the inversion height calculated by the prerun is zero, |
|---|
| 629 | inflow_damping_height must be explicitly specified. (<span style="font-family: Courier New,Courier,monospace;">init_3d_model</span>)<br> |
|---|
| 630 | |
|---|
| [392] | 631 | <br> |
|---|
| [481] | 632 | |
|---|
| [392] | 633 | Small bugfix concerning 3d 64bit netcdf output format. (<span style="font-family: Courier New,Courier,monospace;">header</span>)<br> |
|---|
| [481] | 634 | |
|---|
| [392] | 635 | <br> |
|---|
| [481] | 636 | |
|---|
| [392] | 637 | Bugfix: <span style="font-family: Courier New,Courier,monospace;">dt_fixed</span> |
|---|
| 638 | removed from the restart file, because otherwise, no change from a |
|---|
| [481] | 639 | fixed to a variable timestep would be possible in restart runs. (<span style="font-family: Courier New,Courier,monospace;">read_var_list</span>, |
|---|
| 640 | <span style="font-family: Courier New,Courier,monospace;">write_var_list</span>)<br> |
|---|
| 641 | |
|---|
| [392] | 642 | <br> |
|---|
| [481] | 643 | |
|---|
| 644 | Bugfix: initial setting of <span style="font-family: Courier New,Courier,monospace;">time_coupling</span> |
|---|
| 645 | in coupled restart runs. (<span style="font-family: Courier New,Courier,monospace;">time_integration</span>)</td> |
|---|
| 646 | <td style="vertical-align: top; width: 189px;">advec_particles, |
|---|
| [392] | 647 | check_parameters, cpu_log, data_output_2d, data_output_3d, header, |
|---|
| 648 | init_3d_model, init_particles, init_ocean, modules, netcdf, |
|---|
| 649 | prandtl_fluxes, production_e, read_var_list, time_integration, |
|---|
| 650 | user_last_actions, write_var_list<br> |
|---|
| [481] | 651 | |
|---|
| 652 | </td> |
|---|
| 653 | </tr> |
|---|
| 654 | <tr> |
|---|
| 655 | <td style="vertical-align: top; width: 27px;">02/02/09<br> |
|---|
| 656 | </td> |
|---|
| 657 | <td style="vertical-align: top; width: 38px;">SR<br> |
|---|
| 658 | </td> |
|---|
| 659 | <td style="vertical-align: top; width: 42px;">228<br> |
|---|
| 660 | </td> |
|---|
| 661 | <td style="vertical-align: top; width: 40px;">3.6<br> |
|---|
| 662 | </td> |
|---|
| 663 | <td style="vertical-align: top; width: 27px;">N<br> |
|---|
| 664 | </td> |
|---|
| 665 | <td style="vertical-align: top; width: 497px;">This |
|---|
| 666 | release contains further additions, changes, and bugfixes for the |
|---|
| 667 | SGI-ICE system, as well as important changes like the splitting of the |
|---|
| 668 | user-interface into single files, a complete revision of the |
|---|
| 669 | pre-compiling mechanism allowing to use different make-depositories |
|---|
| 670 | simultaneously, and a revision of the reading-mechanism from restart |
|---|
| 671 | files, which has no more restrictions in case that the previous and |
|---|
| 672 | current run have different domain/subdomain sizes. A new handling of |
|---|
| 673 | all kind of output messages is introduced with this release and will be |
|---|
| 674 | applied to all existing messages within the next releases.<br> |
|---|
| 675 | |
|---|
| [392] | 676 | <br> |
|---|
| [481] | 677 | |
|---|
| 678 | <span style="font-weight: bold;">New features are:</span><br style="font-weight: bold;"> |
|---|
| 679 | |
|---|
| 680 | <br> |
|---|
| 681 | |
|---|
| 682 | Restart runs on SGI-ICE are working. (mrun)<br> |
|---|
| 683 | |
|---|
| 684 | 2d-decomposition is default on SGI-ICE systems. (init_pegrid)<br> |
|---|
| 685 | |
|---|
| 686 | <br> |
|---|
| 687 | |
|---|
| 688 | Ocean-atmosphere coupling realized with MPI-1. mrun adjusted for this |
|---|
| 689 | case (-Y option). Check that PALM is called with <span style="font-family: Courier New,Courier,monospace;">mrun</span>-option |
|---|
| 690 | "<span style="font-style: italic;">-K parallel</span>" for coupling. |
|---|
| 691 | Adjustments in <span style="font-family: Courier New,Courier,monospace;">mrun</span>, <span style="font-family: Courier New,Courier,monospace;">mbuild</span>, and |
|---|
| 692 | <span style="font-family: Courier New,Courier,monospace;">subjob</span> |
|---|
| 693 | for lcxt4.<br> |
|---|
| 694 | |
|---|
| 695 | <br> |
|---|
| 696 | |
|---|
| 697 | DVRP arguments changed to single precision, mode pathlines added.<br> |
|---|
| 698 | |
|---|
| 699 | <br> |
|---|
| 700 | |
|---|
| 701 | User can add additional routines in files <span style="font-family: Courier New,Courier,monospace;">user_additional_routines</span>.<br> |
|---|
| 702 | |
|---|
| 703 | <br> |
|---|
| 704 | |
|---|
| 705 | User can check user parameters and deduce further quantities in <span style="font-family: Courier New,Courier,monospace;">user_check_parameters</span></td> |
|---|
| 706 | <td style="vertical-align: top; width: 189px;">check_for_restart, |
|---|
| 707 | check_parameters, data_output_dvrp, init_dvrp, init_pegrid, local_stop, |
|---|
| 708 | modules, package_parin, palm, surface_coupler, timestep<br> |
|---|
| 709 | |
|---|
| 710 | <br> |
|---|
| 711 | |
|---|
| 712 | Makefile, mbuild, mrun, subjob<br> |
|---|
| 713 | |
|---|
| 714 | <br> |
|---|
| 715 | |
|---|
| 716 | <span style="font-weight: bold;">new:</span><br> |
|---|
| 717 | |
|---|
| 718 | init_coupling, user_additional_routines, user_check_parameters</td> |
|---|
| 719 | </tr> |
|---|
| 720 | <tr> |
|---|
| 721 | <td style="vertical-align: top; width: 27px;"><br> |
|---|
| [392] | 722 | </td> |
|---|
| [481] | 723 | <td style="vertical-align: top; width: 38px;"><br> |
|---|
| 724 | </td> |
|---|
| 725 | <td style="vertical-align: top; width: 42px;"><br> |
|---|
| 726 | </td> |
|---|
| 727 | <td style="vertical-align: top; width: 40px;"><br> |
|---|
| 728 | </td> |
|---|
| 729 | <td style="vertical-align: top; width: 27px;">C<br> |
|---|
| 730 | </td> |
|---|
| 731 | <td style="vertical-align: top; width: 497px;">User interface has been |
|---|
| 732 | split into single files.<br> |
|---|
| 733 | |
|---|
| 734 | <br> |
|---|
| 735 | |
|---|
| 736 | Type of variables <span style="font-family: Courier New,Courier,monospace;">count</span> and <span style="font-family: Courier New,Courier,monospace;">count_rate</span> |
|---|
| 737 | changed to INTEGER(8) in order to avoid "out of range" problems, which |
|---|
| 738 | result in measured negative time intervals. (cpu_log, local_tremain, |
|---|
| 739 | local_tremain_ini)<br> |
|---|
| 740 | |
|---|
| 741 | <br> |
|---|
| 742 | |
|---|
| 743 | Reading mechanism from restart files |
|---|
| 744 | completely revised. The subdomain/total domain size is now allowed to |
|---|
| 745 | vary arbitrarily between the current and previous run. (read_var_list, |
|---|
| 746 | read_3d_binary, user_read_restart_data)<br> |
|---|
| 747 | |
|---|
| 748 | <br> |
|---|
| 749 | |
|---|
| 750 | Precompilation mechanism (<span style="font-family: Courier New,Courier,monospace;">mbuild</span>) |
|---|
| 751 | completely revised: the source-tarfile remains in the source directory (<span style="font-family: Courier New,Courier,monospace;">trunk/SOURCE</span>). |
|---|
| 752 | One depository per block (given in the config-file) is created. Always |
|---|
| 753 | all files from the source directory are copied to the respective |
|---|
| 754 | depository. No additional file checks are done any more (mrun version |
|---|
| 755 | 2.1). The depository name used contains the conditions given by <span style="font-family: Courier New,Courier,monospace;">mrun</span>-option |
|---|
| 756 | "<span style="font-family: Courier New,Courier,monospace;">-K</span>". |
|---|
| 757 | (mbuild, mrun)<br> |
|---|
| 758 | |
|---|
| 759 | <br> |
|---|
| 760 | |
|---|
| 761 | Output messages (including required stop of execution) can now be |
|---|
| 762 | handled with the new subroutine <span style="font-family: Courier New,Courier,monospace;">message</span>. |
|---|
| 763 | All output messages will be replaced by this routine step by step |
|---|
| 764 | within the next revisions.<br> |
|---|
| 765 | |
|---|
| 766 | <br> |
|---|
| 767 | |
|---|
| 768 | Neumann boundary condition at <span style="font-family: Courier New,Courier,monospace;">k=nzb</span> is |
|---|
| 769 | explicitly set for better reading, although this has been already done |
|---|
| 770 | in boundary_conds. (advec_s_bc)<br> |
|---|
| 771 | |
|---|
| 772 | <br> |
|---|
| 773 | |
|---|
| 774 | Origin of the xy-coordinate system shifted from the center of the first |
|---|
| 775 | grid cell (indices <span style="font-family: Courier New,Courier,monospace;">i=0</span>, <span style="font-family: Courier New,Courier,monospace;">j=0</span>) to the |
|---|
| 776 | south-left corner of this cell. (netcdf) Topography definition |
|---|
| 777 | according to new user parameter <span style="font-family: Courier New,Courier,monospace;">topography_grid_convention</span>. |
|---|
| 778 | (init_grid, modules, user_header, user_init_grid, user_parin)</td> |
|---|
| 779 | <td style="vertical-align: top; width: 189px;">advec_s_bc, |
|---|
| 780 | check_parameters, cpu_log, init_grid, local_stop, local_tremain, |
|---|
| 781 | local_tremain_ini, modules, netcdf, read_3d_binary, read_var_list, |
|---|
| 782 | user_read_restart_data<br> |
|---|
| 783 | |
|---|
| 784 | <br> |
|---|
| 785 | |
|---|
| 786 | Makefile, mbuild, mrun<br> |
|---|
| 787 | |
|---|
| 788 | <br> |
|---|
| 789 | |
|---|
| 790 | <span style="font-weight: bold;">removed:</span><br> |
|---|
| 791 | |
|---|
| 792 | user_interface<br> |
|---|
| 793 | |
|---|
| 794 | <br> |
|---|
| 795 | |
|---|
| 796 | <span style="font-weight: bold;">new:</span><br> |
|---|
| 797 | |
|---|
| 798 | message, |
|---|
| 799 | user_3d_data_averaging, user_actions, user_advec_particles, |
|---|
| 800 | user_check_data_output, user_check_data_output_pr, user_data_output_2d, |
|---|
| 801 | user_data_output_3d, user_data_output_dvrp, user_define_netcdf_grid, |
|---|
| 802 | user_dvrp_coltab, user_header, user_init_3d_model, user_init, |
|---|
| 803 | user_init_grid, user_init_particles, user_init_plant_canopy, |
|---|
| 804 | user_last_actions, user_parin, user_parin, user_particle_attributes, |
|---|
| 805 | user_read_restart_data, user_spectra, user_statistics</td> |
|---|
| [392] | 806 | </tr> |
|---|
| [481] | 807 | <tr> |
|---|
| 808 | <td style="vertical-align: top; width: 27px;"><br> |
|---|
| 809 | </td> |
|---|
| 810 | <td style="vertical-align: top; width: 38px;"><br> |
|---|
| 811 | </td> |
|---|
| 812 | <td style="vertical-align: top; width: 42px;"><br> |
|---|
| 813 | </td> |
|---|
| 814 | <td style="vertical-align: top; width: 40px;"><br> |
|---|
| 815 | </td> |
|---|
| [482] | 816 | <td style="vertical-align: top; width: 27px;">B<br> |
|---|
| [481] | 817 | </td> |
|---|
| 818 | <td style="vertical-align: top; width: 497px;">Bugfix: array <span style="font-family: Courier New,Courier,monospace;">d</span> is |
|---|
| 819 | reallocated in case that multigrid is used. (calc_spectra)<br> |
|---|
| [392] | 820 | |
|---|
| [481] | 821 | <br> |
|---|
| 822 | |
|---|
| 823 | Bugfixes for nonparallel execution. (check_for_restart, cpu_statistics, |
|---|
| 824 | inflow_turbulence, timestep)<br> |
|---|
| 825 | |
|---|
| 826 | <br> |
|---|
| 827 | |
|---|
| 828 | Size of <span style="font-family: Courier New,Courier,monospace;">pf3d |
|---|
| 829 | </span>adjusted to the required output size (1 gridpoint less, |
|---|
| 830 | along all three dimensions), because output of a subset of the data (<span style="font-family: Courier New,Courier,monospace;">pf3d(nxa:nxe...)</span>) |
|---|
| 831 | in the NF90_PUT_VAR statement caused segmentation fault with the INTEL |
|---|
| 832 | compiler. (combine_plot_fields)<br> |
|---|
| 833 | |
|---|
| 834 | <br> |
|---|
| 835 | |
|---|
| 836 | Bugfix: error in <span style="font-family: Courier New,Courier,monospace;">zu</span> index in |
|---|
| 837 | case of <span style="font-family: Courier New,Courier,monospace;">section_xy=-1</span>. |
|---|
| 838 | (header)<br> |
|---|
| 839 | |
|---|
| 840 | <br> |
|---|
| 841 | |
|---|
| 842 | Inconsistency |
|---|
| 843 | removed: as the thermal stratification is not taken into account for |
|---|
| 844 | the evaluation of the wall fluxes at vertical walls, the eddy viscosity |
|---|
| 845 | <span style="font-family: Courier New,Courier,monospace;">km</span> |
|---|
| 846 | must not be used for the evaluation of the velocity gradients <span style="font-family: Courier New,Courier,monospace;">dudy</span>, <span style="font-family: Courier New,Courier,monospace;">dwdy</span>, <span style="font-family: Courier New,Courier,monospace;">dvdx</span> and <span style="font-family: Courier New,Courier,monospace;">dwdx</span>. |
|---|
| 847 | (production_e)<br> |
|---|
| 848 | |
|---|
| 849 | <br> |
|---|
| 850 | |
|---|
| 851 | Bugfix in calculating <span style="font-family: Courier New,Courier,monospace;">k</span> index in |
|---|
| 852 | case of oceans runs. (sort_particles)<br> |
|---|
| 853 | |
|---|
| 854 | <br> |
|---|
| 855 | |
|---|
| 856 | Bugfix: no output of particle concentration and radius unless particles |
|---|
| 857 | have been started. (data_output_2d)<br> |
|---|
| 858 | |
|---|
| 859 | <br> |
|---|
| 860 | |
|---|
| 861 | Bugfix: reading of <span style="font-family: Courier New,Courier,monospace;">spectrum_x|y</span> |
|---|
| 862 | from restart files ignored if total numbers of grid points do not |
|---|
| 863 | match. (read_3d_binary)<br> |
|---|
| 864 | |
|---|
| 865 | <br> |
|---|
| 866 | |
|---|
| 867 | Bugfix: abort in case that absolute temperature is below zero. |
|---|
| 868 | (init_cloud_physics)</td> |
|---|
| 869 | <td style="vertical-align: top; width: 189px;">advec_particles, |
|---|
| 870 | calc_spectra, check_for_restart, cpu_statistics, data_output_2d, |
|---|
| 871 | header, inflow_turbulence, init_cloud_physics, production_e, |
|---|
| 872 | read_3d_binary, timestep<br> |
|---|
| 873 | |
|---|
| 874 | <br> |
|---|
| 875 | |
|---|
| 876 | combine_plot_fields</td> |
|---|
| 877 | </tr> |
|---|
| 878 | </tbody> |
|---|
| 879 | </table> |
|---|
| 880 | <b><blink></blink></b><br> |
|---|
| [62] | 881 | <br> |
|---|
| 882 | <h2><a name="Kapitel2.0"></a>2.0 How |
|---|
| 883 | to change the model source |
|---|
| [392] | 884 | code</h2>Currently |
|---|
| 885 | (Sep 23rd, 2009), only selected users are allowed to commit changes to |
|---|
| 886 | the repository. Some rules about how to check in modified code will be |
|---|
| 887 | added later.<br> |
|---|
| 888 | <h2> |
|---|
| [62] | 889 | <a name="Kapitel3.0"></a>3.0 Description |
|---|
| 890 | of selected parts of the |
|---|
| [5] | 891 | model source code</h2> |
|---|
| [392] | 892 | <b></b> |
|---|
| 893 | <p>This section gives links to existing descriptions about special |
|---|
| 894 | parts of the PALM code, e.g. special numerical methods that are used |
|---|
| 895 | (including mathematical background). Documents are stored in single |
|---|
| 896 | files with different formats. Some are in German (sorry for the |
|---|
| 897 | inconvenience). |
|---|
| [62] | 898 | <br> |
|---|
| 899 | </p><ul> <li> <a name="UPS"></a>upstream-spline |
|---|
| [481] | 900 | advection scheme ( <a href="methods/upstream_spline/ups.pdf">.pdf</a> |
|---|
| 901 | , only in German)</li> <li> <a name="Cloud_physics"></a>cloud |
|---|
| 902 | physics module ( <a href="methods/cloud_physics/cloud_physics.pdf">.pdf</a> |
|---|
| 903 | )</li> |
|---|
| 904 | <li> <a name="Bitkompression"></a>data compression with <a href="methods/bit_compression/bit_compression.html">bit shifting method</a> (only in German)<br> |
|---|
| 905 | </li> |
|---|
| 906 | <li>Runge-Kutta |
|---|
| 907 | time integration scheme ( <a href="methods/runge_kutta/runge_kutta.pdf">.pdf</a> |
|---|
| 908 | , only in German)</li> |
|---|
| 909 | |
|---|
| [62] | 910 | </ul><br> |
|---|
| [5] | 911 | |
|---|
| [62] | 912 | <br> |
|---|
| 913 | <br> |
|---|
| 914 | <br> |
|---|
| 915 | <br> |
|---|
| 916 | <br> |
|---|
| 917 | <br> |
|---|
| 918 | <br> |
|---|
| [139] | 919 | </body></html> |
|---|
