Changeset 90


Ignore:
Timestamp:
May 30, 2007 9:18:47 AM (14 years ago)
Author:
raasch
Message:

New:
---
Calculation and output of user-defined profiles. New &userpar parameters data_output_pr_user and max_pr_user.

check_parameters, flow_statistics, modules, parin, read_var_list, user_interface, write_var_list

Changed:


Division through dt_3d replaced by multiplication of the inverse. For performance optimisation, this is done in the loop calculating the divergence instead of using a seperate loop. (pres.f90) var_hom and var_sum renamed pr_palm.

data_output_profiles, flow_statistics, init_3d_model, modules, parin, pres, read_var_list, run_control, time_integration

Errors:


Bugfix: work_fft*_vec removed from some PRIVATE-declarations (poisfft).

Bugfix: field_chr renamed field_char (user_interface).

Bugfix: output of use_upstream_for_tke (header).

header, poisfft, user_interface

Location:
palm/trunk
Files:
18 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_3.5.1.html

    r62 r90  
    11<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
    22<html><head>
    3 <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM
    4 chapter 3.5.1</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;11451548"> <meta name="CHANGED" content="20041112;14401801"> <meta name="KEYWORDS" content="parallel LES model"> <style>
     3<meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM chapter 3.5.1</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;11451548"> <meta name="CHANGED" content="20041112;14401801"> <meta name="KEYWORDS" content="parallel LES model"> <style>
    54<!--
    65@page { size: 21cm 29.7cm }
     
    204203quantities can be computed here. Also, additional time series
    205204quantities can be calculated. They have to be defined before in routine
    206 <a href="#user_init"><span style="font-family: Courier New,Courier,monospace;">user_init</span></a>.
    207 The routine is called once for each defined statistic region (see <a href="chapter_4.1.html#statistic_region">statistic_region</a>).</p><p></p><p>The
     205<a style="font-family: monospace;" href="#user_init">user_init</a>.
     206The routine is called once for each defined statistic region (see <a href="chapter_4.1.html#statistic_region">statistic_region</a>).</p><p>The
    208207routine contains some simple examples (as comment lines) in order to
    209 demonstrate how to use it.</p><p><br>Profiles have to
    210 be written to
    211 the array <tt><font style="font-size: 10pt;" size="2">sums</font></tt>.
    212 Please use profile numbers starting from 70 upward (thus e.g. <tt><font style="font-size: 10pt;" size="2">sums (&hellip;,
    213 60) =</font></tt><font style="font-size: 10pt;" size="2">&hellip;)</font>. These additional
     208demonstrate how to use it.<br><br>Profile sums have to
     209be stored in
     210the array <tt><font style="font-size: 10pt;" size="2">sums_l</font></tt>.&nbsp;The profile identification number (second index of array <span style="font-family: monospace;">sums_l</span>) must be within the range [<span style="font-family: monospace;"> pr_palm+1 , pr_palm+max_pr_user</span> ], where <span style="font-family: monospace;">pr_palm</span> is an internal parameter and&nbsp; <span style="font-family: Courier New,Courier,monospace;">max_pr_user</span> is the number of&nbsp;user-defined profiles as given by parameter <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.3.html#data_output_pr_user">data_output_pr_user</a> in the respective PALM run.</p><p> These additional
    214211profiles
    215 are temporally averaged in <tt><font style="font-size: 10pt;" size="2">flow_statistics</font></tt>.</p> </td>
     212are further processed (e.g. temporally averaged), like the standard profiles, in <tt><font style="font-size: 10pt;" size="2">flow_statistics</font></tt>.</p> </td>
    216213</tr> <tr><td style="vertical-align: top;">End
    217214of <span style="font-family: monospace;">check_parameters.f90</span>
    218 where output quantities are checked</td><td style="vertical-align: top;"><span style="font-family: monospace;"><a name="user_check_data_output"></a>user_check_data_<br>output(
    219 variable, unit )</span></td><td style="vertical-align: top;">Here the <span style="font-family: monospace;">unit</span>(s) of
    220 user-defined output quantities (see <a href="chapter_4.3.html#data_output_user">data_output_user</a>)
     215where 2d/3d output quantities are checked</td><td style="vertical-align: top;"><span style="font-family: monospace;"><a name="user_check_data_output"></a>user_check_data_<br>output(
     216variable, unit )</span></td><td style="vertical-align: top;">Here the physical <span style="font-family: monospace;">unit</span>(s) of
     217user-defined 2d/3d output quantities (see <a href="chapter_4.3.html#data_output_user">data_output_user</a>)
    221218have to be set. Input parameter <span style="font-family: monospace;">variable</span>
    222219contains the string identifier of the respective variable. For those
     220variables not recognized by the user, the parameter <span style="font-family: monospace;">unit</span> is set to
     221"illegal", which tells the calling routine that the output variable is
     222not defined and leads to a program abort.<br><br>See <a href="chapter_3.5.4.html">chapter 3.5.4</a> about
     223creating user-defined output quantities.</td></tr><tr><td style="vertical-align: top;">End
     224of <span style="font-family: monospace;">check_parameters.f90</span>
     225where vertical profile output quantities are checked</td><td style="vertical-align: top;"><a name="user_check_data_output_pr"></a><br><span style="font-family: monospace;">user_check_data_</span><br style="font-family: monospace;"><span style="font-family: monospace;">output_pr(</span><br style="font-family: monospace;"><span style="font-family: monospace;">variable, var_count, unit )</span></td><td style="vertical-align: top;">Here the physical <span style="font-family: monospace;">unit</span>(s) and vertical grid (u- or w-grid) of
     226user-defined vertical profile output quantities (see <a href="chapter_4.3.html#data_output_pr_user">data_output_pr_user</a>)
     227have to be set. Input parameter <span style="font-family: monospace;">variable</span>
     228contains the string identifier of the respective variable. Parameter<span style="font-family: monospace;"> var_count </span>contains the internal profile number and must not be changed by the user. For those
    223229variables not recognized by the user, the parameter <span style="font-family: monospace;">unit</span> is set to
    224230"illegal", which tells the calling routine that the output variable is
  • palm/trunk/DOC/tec/technical_documentation.html

    r83 r90  
    1212<br><table nosave="" cellpadding="0" cellspacing="0"> <caption>&nbsp; <br> </caption><tbody>
    1313</tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Current
    14 model version:</b></td> <td><span style="font-weight: bold;">3.2a</span></td> </tr>
     14model version:</b></td> <td><span style="font-weight: bold;">3.2b</span></td> </tr>
    1515<tr nosave=""> <td nosave=""><b>Last
    1616change of this document</b>:&nbsp;</td> <td nosave=""><b>$Id$</b></td> </tr>
     
    23082308header, init_dvrp, init_particles, init_1d_model, init_dvrp,
    23092309init_pegrid, local_getenv, local_system, local_tremain,
    2310 local_tremain_ini, modules, palm, parin, run_control<br><br><span style="font-weight: bold;">new:</span><br>local_flush<br><br><span style="font-weight: bold;">changed scripts:</span><br>mbuild, mrun, .mrun.config.default<br></td></tr>
     2310local_tremain_ini, modules, palm, parin, run_control<br><br><span style="font-weight: bold;">new:</span><br>local_flush<br><br><span style="font-weight: bold;">changed scripts:</span><br>mbuild, mrun, .mrun.config.default<br></td></tr><tr><td style="vertical-align: top;">30/05/07</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.2b</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">Calculation and output of user-defined profiles: New&nbsp;userpar-parameter <span style="font-family: Courier New,Courier,monospace;">data_output_pr_user</span>. In the user-interface, routine <span style="font-family: Courier New,Courier,monospace;">user_check_data_output_pr</span> has been added and <span style="font-family: Courier New,Courier,monospace;">user_statistics</span> has been extended.<br></td><td style="vertical-align: top;">check_parameters, flow_statistics, modules, parin, read_var_list, user_interface, write_var_list</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;">C</td><td style="vertical-align: top;">In <span style="font-family: Courier New,Courier,monospace;">pres.f90</span>, division through <span style="font-family: Courier New,Courier,monospace;">dt_3d</span>
     2311replaced by multiplication of the inverse. For performance
     2312optimisation, this is done in the loop calculating the divergence
     2313instead of using a seperate loop.<br><br>Variables <span style="font-family: Courier New,Courier,monospace;">var_hom</span> and <span style="font-family: Courier New,Courier,monospace;">var_sum</span> are both renamed <span style="font-family: Courier New,Courier,monospace;">pr_palm</span>.</td><td style="vertical-align: top;">data_output_profiles, flow_statistics, init_3d_model, modules, parin, pres, read_var_list, run_control, time_integration</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: <span style="font-family: Courier New,Courier,monospace;">work_fft*_vec</span> removed from some PRIVATE-declarations (<span style="font-family: Courier New,Courier,monospace;">poisfft</span>).<br><br>Bugfix: <span style="font-family: Courier New,Courier,monospace;">field_chr</span> renamed <span style="font-family: Courier New,Courier,monospace;">field_char</span> (<span style="font-family: Courier New,Courier,monospace;">user_interface</span>).<br><br>Bugfix: output of <span style="font-family: Courier New,Courier,monospace;">use_upstream_for_tke</span> (<span style="font-family: Courier New,Courier,monospace;">header</span>).</td><td style="vertical-align: top;">header, poisfft, user_interface</td></tr>
    23112314</tbody>
    23122315</table>&nbsp;<b><blink>Attention:</blink></b>
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r87 r90  
    11New:
    22---
    3 Calculation and output of user-defined profiles. New &userpar parameters data_output_pr_user and max_pr_user.
    4 
    5 check_parameters, flow_statistics, modules, parin, read_var_list, user_interface, write_var_list
    63
    74
    85Changed:
    96-------
    10 Division through dt_3d replaced by multiplication of the inverse. For performance optimisation, this is done in the loop calculating the divergence instead of using a seperate loop. (pres.f90) var_hom and var_sum renamed pr_palm.
    11 
    12 data_output_profiles, flow_statistics, init_3d_model, modules, parin, pres, read_var_list, run_control, time_integration
    137
    148
    159Errors:
    1610------
    17 Bugfix: work_fft*_vec removed from some PRIVATE-declarations (poisfft).
    18 
    19 Bugfix: field_chr renamed field_char (user_interface).
    20 
    21 Bugfix: output of use_upstream_for_tke (header).
    22 
    23 header, poisfft, user_interface
  • palm/trunk/SOURCE/check_parameters.f90

    r89 r90  
    44! Actual revisions:
    55! -----------------
    6 ! Check for user-defined profiles
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 89 2007-05-25 12:08:31Z raasch
     13! Check for user-defined profiles
    1114!
    1215! 75 2007-03-22 09:54:05Z raasch
  • palm/trunk/SOURCE/data_output_profiles.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! var_hom renamed pr_palm
     6!
    77!
    88! Former revisions:
     
    1010! $Id$
    1111! RCS Log replace by Id keyword, revision history cleaned up
     12!
     13! 87 2007-05-22 15:46:47Z raasch
     14! var_hom renamed pr_palm
    1215!
    1316! Revision 1.18  2006/08/16 14:27:04  raasch
  • palm/trunk/SOURCE/flow_statistics.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
    613! Two more arguments added to user_statistics, which is now also called for
    714! user-defined profiles,
    815! var_hom and var_sum renamed pr_palm
    9 !
    10 ! Former revisions:
    11 ! -----------------
    12 ! $Id$
    1316!
    1417! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/header.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: output of use_upstream_for_tke
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
     13! Bugfix: output of use_upstream_for_tke
    1114!
    1215! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/init_3d_model.f90

    r87 r90  
    77! Actual revisions:
    88! -----------------
    9 ! var_hom and var_sum renamed pr_palm
     9!
    1010!
    1111! Former revisions:
    1212! -----------------
    1313! $Id$
     14!
     15! 87 2007-05-22 15:46:47Z raasch
     16! var_hom and var_sum renamed pr_palm
    1417!
    1518! 75 2007-03-22 09:54:05Z raasch
  • palm/trunk/SOURCE/modules.f90

    r89 r90  
    55! Actual revisions:
    66! -----------------
    7 ! +data_output_pr_user, max_pr_user, size of data_output_pr, dopr_index,
    8 ! dopr_initial_index and dopr_unit enlarged,
    9 ! var_hom and var_sum renamed pr_palm
     7!
    108!
    119! Former revisions:
    1210! -----------------
    1311! $Id$
     12!
     13! 89 2007-05-25 12:08:31Z raasch
     14! +data_output_pr_user, max_pr_user, size of data_output_pr, dopr_index,
     15! dopr_initial_index and dopr_unit enlarged,
     16! var_hom and var_sum renamed pr_palm
    1417!
    1518! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/palm.f90

    r83 r90  
    6363    INTEGER           ::  i, run_description_header_i(80)
    6464
    65     version = 'PALM 3.2a'
     65    version = 'PALM 3.2b'
    6666
    6767#if defined( __parallel )
  • palm/trunk/SOURCE/parin.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! Size of hom increased by the maximum number of user-defined profiles,
    7 ! var_hom renamed pr_palm
     6!
    87!
    98! Former revisions:
    109! -----------------
    1110! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
     13! Size of hom increased by the maximum number of user-defined profiles,
     14! var_hom renamed pr_palm
    1215!
    1316! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/poisfft.f90

    r85 r90  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: work_fft*_vec removed from some PRIVATE-declarations
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 85 2007-05-11 09:35:14Z raasch
     13! Bugfix: work_fft*_vec removed from some PRIVATE-declarations
    1114!
    1215! 76 2007-03-29 00:58:32Z raasch
  • palm/trunk/SOURCE/pres.f90

    r85 r90  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 85 2007-05-11 09:35:14Z raasch
    613! Division through dt_3d replaced by multiplication of the inverse.
    714! For performance optimisation, this is done in the loop calculating the
    815! divergence instead of using a seperate loop.
    9 !
    10 ! Former revisions:
    11 ! -----------------
    12 ! $Id$
    1316!
    1417! 75 2007-03-22 09:54:05Z raasch
  • palm/trunk/SOURCE/read_var_list.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! +max_pr_user (version 3.1), var_hom renamed pr_palm
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
     13! +max_pr_user (version 3.1), var_hom renamed pr_palm
    1114!
    1215! 75 2007-03-22 09:54:05Z raasch
  • palm/trunk/SOURCE/run_control.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! var_hom renamed pr_palm
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
     13! var_hom renamed pr_palm
    1114!
    1215! 82 2007-04-16 15:40:52Z raasch
  • palm/trunk/SOURCE/time_integration.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! var_hom renamed pr_palm
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
     13! var_hom renamed pr_palm
    1114!
    1215! 75 2007-03-22 09:54:05Z raasch
  • palm/trunk/SOURCE/user_interface.f90

    r89 r90  
    44! Actual revisions:
    55! -----------------
     6!
     7!
     8! Former revisions:
     9! -----------------
     10! $Id$
     11!
     12! 89 2007-05-25 12:08:31Z raasch
    613! Calculation and output of user-defined profiles: new routine
    714! user_check_data_output_pr, +data_output_pr_user, max_pr_user in userpar,
    815! routine user_statistics has got two more arguments
    916! Bugfix: field_chr renamed field_char
    10 !
    11 ! Former revisions:
    12 ! -----------------
    13 ! $Id$
    1417!
    1518! 60 2007-03-11 11:50:04Z raasch
  • palm/trunk/SOURCE/write_var_list.f90

    r87 r90  
    44! Actual revisions:
    55! -----------------
    6 ! +max_pr_user (version 3.1)
     6!
    77!
    88! Former revisions:
    99! -----------------
    1010! $Id$
     11!
     12! 87 2007-05-22 15:46:47Z raasch
     13! +max_pr_user (version 3.1)
    1114!
    1215! 75 2007-03-22 09:54:05Z raasch
Note: See TracChangeset for help on using the changeset viewer.