Home


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

Changeset 197 for palm/trunk


Ignore:
Timestamp:
Sep 16, 2008 3:29:03 PM (16 years ago)
Author:
raasch
Message:

further adjustments for SGI and other small changes

Location:
palm/trunk
Files:
1 added
21 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/DOC/app/chapter_3.4.html

    r134 r197  
    88 
    99  <title>PALM chapter 3.4</title></head>
    10 
    1110<body>
    1211
    … …  
    22122211accuracy has to be
    22132212indicated
    2214 with the help of <a href="chapter_4.2.html#do3d_precision">do3d_precision</a>.
     2213with the help of <a href="chapter_4.2.html#do3d_comp_prec">do3d_comp_prec</a>.
    22152214In favorable cases the file size can be reduced by
    22162215a factor of up to 5. For the visualization of compressed data a special

  • palm/trunk/DOC/app/chapter_3.8.html

    r115 r197  
    11<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
    2 <html>
    3 <head>
     2<html><head>
    43
    54  <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252">
    6   <title>PALM chapter 3.8</title>
     5 
    76
    87
    9   <meta name="GENERATOR" content="StarOffice 7 (Win32)">
     8  <title>PALM chapter 3.8</title><meta name="GENERATOR" content="StarOffice 7 (Win32)">
    109
    1110  <meta name="AUTHOR" content="Siegfried Raasch">
    … …  
    2120@page { size: 21cm 29.7cm }
    2221-->
    23   </style>
    24 </head>
    25 
     22  </style></head>
    2623
    2724<body style="direction: ltr;" lang="en-US">
    … …  
    126123<p style="line-height: 100%;"><br>
    127124
    128 <font color="#000080"><font color="#000080"><a href="chapter_3.7.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" src="left.gif" name="Grafik1"></font></a><a href="index.html"><font color="#000080"><img src="up.gif" name="Grafik2" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_4.0.html"><font color="#000080"><img src="right.gif" name="Grafik3" align="bottom" border="2" height="32" width="32"></font></a></font></font></p>
     125<font color="#000080"><font color="#000080"><a href="chapter_3.7.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" src="left.gif" name="Grafik1"></font></a><a href="index.html"><font color="#000080"><img src="up.gif" name="Grafik2" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_3.9.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" src="right.gif" name="Grafik3"></font></a></font></font></p>
    129126
    130127<p style="line-height: 100%;"><i>Last
    131128change:&nbsp;</i> $Id$</p>
    132129
    133 </body>
    134 </html>
     130</body></html>

  • palm/trunk/DOC/app/chapter_4.1.html

    r166 r197  
    55105510
    55115511
    5512     <tr>
     5512    <tr><td style="vertical-align: top;"><p><a name="humidity"></a><b>humidity</b></p></td><td style="vertical-align: top;">L</td><td style="vertical-align: top;"><i>.F.</i></td><td style="vertical-align: top;"><p>Parameter to
     5513switch on the prognostic equation for specific
     5514humidity q.<br>
     5515
     5516
     5517
     5518
     5519
     5520
     5521 </p>
     5522
     5523
     5524
     5525
     5526
     5527
     5528 
     5529     
     5530     
     5531     
     5532     
     5533     
     5534     
     5535      <p>The initial vertical
     5536profile of q can be set via parameters <a href="chapter_4.1.html#q_surface">q_surface</a>, <a href="chapter_4.1.html#q_vertical_gradient">q_vertical_gradient</a>
     5537and <a href="chapter_4.1.html#q_vertical_gradient_level">q_vertical_gradient_level</a>.&nbsp;
     5538Boundary conditions can be set via <a href="chapter_4.1.html#q_surface_initial_change">q_surface_initial_change</a>
     5539and <a href="chapter_4.1.html#surface_waterflux">surface_waterflux</a>.<br>
     5540
     5541
     5542
     5543
     5544
     5545
     5546
     5547      </p>
     5548
     5549
     5550
     5551
     5552
     5553
     5554
     5555If the condensation scheme is switched on (<a href="chapter_4.1.html#cloud_physics">cloud_physics</a>
     5556= .TRUE.), q becomes the total liquid water content (sum of specific
     5557humidity and liquid water content).</td></tr><tr><td style="vertical-align: top;"><span style="font-weight: bold;"><a name="inflow_damping_height"></a>inflow_damping_height</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">from precursor run</span></td><td style="vertical-align: top;">Height below which the turbulence signal is used for turbulence recycling (in m).<br><br>In case of a turbulent inflow (see <a href="chapter_4.1.html#turbulent_inflow">turbulent_inflow</a>),
     5558this parameter defines the vertical thickness of the turbulent layer up
     5559to which the turbulence extracted at the recycling plane (see <a href="chapter_4.1.html#recycling_width">recycling_width</a>)
     5560shall be imposed to the inflow. Above this level the turbulence signal
     5561is linearly damped to zero. The transition range within which the
     5562signal falls to zero is given by the parameter <a href="chapter_4.1.html#inflow_damping_width">inflow_damping_width</a>.<br><br>By default, this height is set as the height of the convective boundary layer as calculated from a precursor run. See <a href="chapter_3.9.html">chapter 3.9</a> about proper settings for getting this CBL height from a precursor run. </td></tr><tr><td style="vertical-align: top;"><span style="font-weight: bold;"><a name="inflow_damping_width"></a>inflow_damping_width</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">0.1 * <a href="chapter_4.1.html#inflow_damping_height">inflow_damping</a></span><a href="chapter_4.1.html#inflow_damping_height"><br style="font-style: italic;"><span style="font-style: italic;">_height</span></a></td><td style="vertical-align: top;">Transition range within which the turbulance signal is damped to zero (in m).<br><br>See <a href="chapter_4.1.html#inflow_damping_height">inflow_damping_height</a> for explanation.</td></tr><tr>
    55135563
    55145564
    … …  
    60226072must be switched of in the equation of motion&nbsp; for w (this
    60236073requires the user to comment out the call of <span style="font-family: monospace;">buoyancy</span> in the
    6024 source code of <span style="font-family: monospace;">prognostic_equations.f90</span>).</p>
    6025 
    6026 
    6027 
    6028 
    6029 
    6030 
    6031 
    6032      
    6033      
    6034      
    6035      
    6036      
    6037      
    6038       </ul>
    6039 
    6040 
    6041 
    6042 
    6043 
    6044 
    6045  
    6046      
    6047      
    6048      
    6049      
    6050      
    6051      
    6052       <p style="font-style: normal;">Values may be
     6074source code of <span style="font-family: monospace;">prognostic_equations.f90</span>).</p></ul>
     6075
     6076
     6077
     6078
     6079
     6080
     6081 
     6082     
     6083     
     6084     
     6085     
     6086     
     6087     
     6088      <p style="font-style: italic;">'read_data_for_recycling'</p><p style="font-style: normal; margin-left: 40px;">Here,
     60893d-data from a precursor run are read by the initial (main) run. The
     6090precursor run is allowed to have a smaller domain along x and y
     6091compared with the main run. Also, different numbers of processors can
     6092be used for these two runs. Limitations are that the precursor run must
     6093use cyclic horizontal boundary conditions and that the subdomains of
     6094the main run must not be larger than the subdomains of the precursor
     6095run. If the total domain of the main run is larger than that of the precursor
     6096run, the domain is filled by cyclic repetition&nbsp;of the (cyclic)
     6097precursor data. This initialization method is recommended if a
     6098turbulent inflow is used (see <a href="chapter_4.1.html#turbulent_inflow">turbulent_inflow</a>). 3d-data must be made available to the run by activating an appropriate file connection statement for local file BININ. See <a href="chapter_3.9.html">chapter 3.9</a> for more details, where usage of a turbulent inflow is explained. </p><p style="font-style: normal;">Values may be
    60536099combined, e.g. <b>initializing_actions</b> = <span style="font-style: italic;">'set_constant_profiles
    60546100initialize_vortex'</span>, but the values of <span style="font-style: italic;">'set_constant_profiles'</span>,
    … …  
    60696115     
    60706116     
    6071       <p style="font-style: italic;"> </p>
     6117     
    60726118
    60736119
    … …  
    67496795
    67506796
    6751  <tr>
    6752 
    6753 
    6754 
    6755 
    6756 
    6757 
    6758  <td style="vertical-align: top;">
    6759      
    6760      
    6761      
    6762      
    6763      
    6764      
    6765       <p><a name="humidity"></a><b>humidity</b></p>
    6766 
    6767 
    6768 
    6769 
    6770 
    6771 
    6772 
    6773       </td>
    6774 
    6775 
    6776 
    6777 
    6778 
    6779 
    6780  <td style="vertical-align: top;">L</td>
    6781 
    6782 
    6783 
    6784 
    6785 
    6786 
    6787 
    6788       <td style="vertical-align: top;"><i>.F.</i></td>
    6789 
    6790 
    6791 
    6792 
    6793 
    6794 
    6795 
    6796       <td style="vertical-align: top;">
    6797      
    6798      
    6799      
    6800      
    6801      
    6802      
    6803       <p>Parameter to
    6804 switch on the prognostic equation for specific
    6805 humidity q.<br>
    6806 
    6807 
    6808 
    6809 
    6810 
    6811 
    6812  </p>
    6813 
    6814 
    6815 
    6816 
    6817 
    6818 
    6819  
    6820      
    6821      
    6822      
    6823      
    6824      
    6825      
    6826       <p>The initial vertical
    6827 profile of q can be set via parameters <a href="chapter_4.1.html#q_surface">q_surface</a>, <a href="chapter_4.1.html#q_vertical_gradient">q_vertical_gradient</a>
    6828 and <a href="chapter_4.1.html#q_vertical_gradient_level">q_vertical_gradient_level</a>.&nbsp;
    6829 Boundary conditions can be set via <a href="chapter_4.1.html#q_surface_initial_change">q_surface_initial_change</a>
    6830 and <a href="chapter_4.1.html#surface_waterflux">surface_waterflux</a>.<br>
    6831 
    6832 
    6833 
    6834 
    6835 
    6836 
    6837 
    6838       </p>
    6839 
    6840 
    6841 
    6842 
    6843 
    6844 
    6845 
    6846 If the condensation scheme is switched on (<a href="chapter_4.1.html#cloud_physics">cloud_physics</a>
    6847 = .TRUE.), q becomes the total liquid water content (sum of specific
    6848 humidity and liquid water content).</td>
    6849 
    6850 
    6851 
    6852 
    6853 
    6854 
    6855  </tr>
     6797 
    68566798
    68576799
    … …  
    78097751
    78107752
    7811     <tr>
    7812 
    7813 
    7814 
    7815 
    7816 
    7817 
    7818  <td style="vertical-align: top;">
    7819      
    7820      
    7821      
    7822      
    7823      
    7824      
    7825       <p><a name="npex"></a><b>npex</b></p>
    7826 
    7827 
    7828 
    7829 
    7830 
    7831 
    7832  </td>
    7833 
    7834 
    7835 
    7836 
    7837 
    7838 
    7839 
    7840       <td style="vertical-align: top;">I</td>
    7841 
    7842 
    7843 
    7844 
    7845 
    7846 
    7847  <td style="vertical-align: top;"><br>
    7848 
    7849 
    7850 
    7851 
    7852 
    7853 
    7854  </td>
    7855 
    7856 
    7857 
    7858 
    7859 
    7860 
    7861  <td style="vertical-align: top;">
    7862      
    7863      
    7864      
    7865      
    7866      
    7867      
    7868       <p>Number of processors
    7869 along x-direction of the virtual
    7870 processor
    7871 net.&nbsp; </p>
    7872 
    7873 
    7874 
    7875 
    7876 
    7877 
    7878  
    7879      
    7880      
    7881      
    7882      
    7883      
    7884      
    7885       <p>For parallel runs, the total
    7886 number of processors to be used
    7887 is given by
    7888 the <span style="font-weight: bold;">mrun</span>
    7889 option <a href="http://www.muk.uni-hannover.de/software/mrun_beschreibung.html#Opt-X">-X</a>.
    7890 By default, depending on the type of the parallel computer, PALM
    7891 generates a 1d processor
    7892 net (domain decomposition along x, <span style="font-weight: bold;">npey</span>
    7893 = <span style="font-style: italic;">1</span>) or a
    7894 2d-net (this is
    7895 favored on machines with fast communication network). In case of a
    7896 2d-net, it is tried to make it more or less square-shaped. If, for
    7897 example, 16 processors are assigned (-X 16), a 4 * 4 processor net is
    7898 generated (<span style="font-weight: bold;">npex</span>
    7899 = <span style="font-style: italic;">4</span>, <span style="font-weight: bold;">npey</span>
    7900 = <span style="font-style: italic;">4</span>).
    7901 This choice is optimal for square total domains (<a href="#nx">nx</a>
    7902 = <a href="#ny">ny</a>),
    7903 since then the number of ghost points at the lateral boundarys of
    7904 the subdomains is minimal. If <span style="font-weight: bold;">nx</span>
    7905 nd <span style="font-weight: bold;">ny</span>
    7906 differ extremely, the
    7907 processor net should be manually adjusted using adequate values for <span style="font-weight: bold;">npex</span> and <span style="font-weight: bold;">npey</span>.&nbsp; </p>
    7908 
    7909 
    7910 
    7911 
    7912 
    7913 
    7914 
    7915      
    7916      
    7917      
    7918      
    7919      
    7920      
    7921       <p><b>Important:</b> The value of <span style="font-weight: bold;">npex</span> * <span style="font-weight: bold;">npey</span> must exactly
    7922 correspond to the
    7923 value assigned by the <span style="font-weight: bold;">mrun</span>-option
    7924       <tt>-X</tt>.
    7925 Otherwise the model run will abort with a corresponding error
    7926 message.&nbsp; <br>
    7927 
    7928 
    7929 
    7930 
    7931 
    7932 
    7933 
    7934 Additionally, the specification of <span style="font-weight: bold;">npex</span>
    7935 and <span style="font-weight: bold;">npey</span>
    7936 may of course
    7937 override the default setting for the domain decomposition (1d or 2d)
    7938 which may have a significant (negative) effect on the code performance.
    7939       </p>
    7940 
    7941 
    7942 
    7943 
    7944 
    7945 
    7946  </td>
    7947 
    7948 
    7949 
    7950 
    7951 
    7952 
    7953  </tr>
    7954 
    7955 
    7956 
    7957 
    7958 
    7959 
    7960  <tr>
    7961 
    7962 
    7963 
    7964 
    7965 
    7966 
    7967  <td style="vertical-align: top;">
    7968      
    7969      
    7970      
    7971      
    7972      
    7973      
    7974       <p><a name="npey"></a><b>npey</b></p>
    7975 
    7976 
    7977 
    7978 
    7979 
    7980 
    7981 
    7982       </td>
    7983 
    7984 
    7985 
    7986 
    7987 
    7988 
    7989  <td style="vertical-align: top;">I</td>
    7990 
    7991 
    7992 
    7993 
    7994 
    7995 
    7996 
    7997       <td style="vertical-align: top;"><br>
    7998 
    7999 
    8000 
    8001 
    8002 
    8003 
    8004  </td>
    8005 
    8006 
    8007 
    8008 
    8009 
    8010 
    8011  <td style="vertical-align: top;">
    8012      
    8013      
    8014      
    8015      
    8016      
    8017      
    8018       <p>Number of processors
    8019 along y-direction of the virtual
    8020 processor
    8021 net.&nbsp; </p>
    8022 
    8023 
    8024 
    8025 
    8026 
    8027 
    8028  
    8029      
    8030      
    8031      
    8032      
    8033      
    8034      
    8035       <p>For further information see <a href="#npex">npex</a>.</p>
    8036 
    8037 
    8038 
    8039 
    8040 
    8041 
    8042  </td>
    8043 
    8044 
    8045 
    8046 
    8047 
    8048 
    8049  </tr>
     7753   
     7754
     7755
     7756
     7757
     7758
     7759
     7760 
    80507761
    80517762
    … …  
    1153011241
    1153111242
    11532  <tr>
     11243 <tr><td style="vertical-align: top;"><span style="font-weight: bold;"><a name="recycling_width"></a>recycling_width</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">0.1 * <a href="chapter_4.1.html#nx">nx</a> * <a href="chapter_4.1.html#dx">dx</a></span></td><td style="vertical-align: top;">Distance of the recycling plane from the inflow boundary (in m).<br><br>This
     11244parameter sets the horizontal extension (along the direction of the
     11245main flow) of the so-called recycling domain which is used to generate
     11246a turbulent inflow (see <a href="chapter_4.1.html#turbulent_inflow">turbulent_inflow</a>). <span style="font-weight: bold;">recycling_width</span> must be larger than the grid spacing (dx) and smaller than the length of the total domain (nx * dx).</td></tr><tr>
    1153311247
    1153411248
    … …  
    1437614090= <a href="#scalar_advec">scalar_advec</a>
    1437714091= '<i>pw-scheme'</i>, <a href="chapter_4.2.html#psolver">psolver</a>
    14378 = <i>'poisfft'</i> or '<i>poisfft_hybrid'</i>,
     14092/= <i>'sor</i><i>'</i>,
    1437914093      <i>&nbsp;</i><a href="#alpha_surface">alpha_surface</a>
    14380 = 0.0, <a href="#bc_lr">bc_lr</a> = <a href="#bc_ns">bc_ns</a> = <span style="font-style: italic;">'cyclic'</span>,&nbsp;<a style="" href="#galilei_transformation">galilei_transformation</a>
     14094= 0.0,<span style="font-style: italic;"></span>&nbsp;<a style="" href="#galilei_transformation">galilei_transformation</a>
    1438114095= <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>,&nbsp;&nbsp;<a href="#humidity">humidity</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
    1438214096
    … …  
    1495214666
    1495314667
    14954     <tr>
     14668    <tr><td style="vertical-align: top;"><a name="turbulent_inflow"></a><span style="font-weight: bold;">turbulent_inflow</span></td><td style="vertical-align: top;">L</td><td style="vertical-align: top;"><span style="font-style: italic;">.F.</span></td><td style="vertical-align: top;">Generates a turbulent inflow at side boundaries using a turbulence recycling method.<br><br>Turbulent inflow is realized using the turbulence recycling method from Lund et al. (1998, J. Comp. Phys., <span style="font-weight: bold;">140</span>, 233-258) modified by Kataoka and Mizuno (2002, Wind and Structures, <span style="font-weight: bold;">5</span>, 379-392).<br><br>A turbulent inflow requires Dirichlet conditions at the respective inflow boundary. <span style="font-weight: bold;">So far, a turbulent inflow is realized from the left (west) side only, i.e. </span><a style="font-weight: bold;" href="chapter_4.1.html#bc_lr">bc_lr</a><span style="font-weight: bold;">&nbsp;=</span><span style="font-style: italic; font-weight: bold;"> 'dirichlet/radiation'</span><span style="font-weight: bold;"> is required!</span><br><br>The initial (quasi-stationary) turbulence field should be generated by a precursor run and used by setting <a href="chapter_4.1.html#initializing_actions">initializing_actions</a> =<span style="font-style: italic;"> 'read_data_for_recycling'</span>.<br><br>The distance of the recycling plane from the inflow boundary can be set with parameter <a href="chapter_4.1.html#recycling_width">recycling_width</a>.
     14669The heigth above ground above which the turbulence signal is not used
     14670for recycling and the width of the layer within&nbsp;the magnitude of
     14671the turbulence signal is damped from 100% to 0% can be set with
     14672parameters <a href="chapter_4.1.html#inflow_damping_height">inflow_damping_height</a> and <a href="chapter_4.1.html#inflow_damping_width">inflow_damping_width</a>.<br><br>The detailed setup for a turbulent inflow is described in <a href="chapter_3.9.html">chapter 3.9</a>.</td></tr><tr>
    1495514673
    1495614674

  • palm/trunk/DOC/app/chapter_4.2.html

    r189 r197  
    25752575
    25762576
    2577  <tr>
     2577 <tr><td align="undefined" valign="undefined"><span style="font-style: italic; color: rgb(255, 0, 0);">p</span></td><td align="undefined" valign="undefined">Perturbation pressure (in Pa)</td></tr><tr>
    25782578
    25792579
    … …  
    37173717assigned. This
    37183718yields a loss of accuracy, but the file size is clearly reduced. The
    3719 parameter <a href="chapter_4.2.html#do3d_precision">do3d_precision</a>
     3719parameter <a href="chapter_4.2.html#do3d_comp_prec">do3d_comp_prec</a>
    37203720can be used to separately define the number of significant digits for
    37213721each quantity.<br>
    … …  
    37473747     
    37483748     
    3749       <p><a name="do3d_precision"></a><b>do3d_precision</b></p>
     3749      <p><a name="do3d_comp_prec"></a><b>do3d_comp_prec</b></p>
    37503750
    37513751
    … …  
    38173817     
    38183818     
    3819       <p>The default assignment is <span style="font-weight: bold;">do3d_precision</span>
     3819      <p>The default assignment is <span style="font-weight: bold;">do3d_comp_prec</span>
    38203820= <span style="font-style: italic;">'u2'</span>, <span style="font-style: italic;">'v2'</span>, <span style="font-style: italic;">'w2'</span>, <span style="font-style: italic;">'p5'</span>, <span style="font-style: italic;">'pt2'</span>.</p>
    38213821
    … …  
    48174817
    48184818            <td style="vertical-align: top;">Subgrid-scale sensible
    4819 heat flux near the surface (horizontal
    4820 average)
    4821 between z = 0 and z = z<sub>p</sub> = zu(1) (there it
    4822 corresponds to
    4823 the total heat flux) (in K m/s).</td>
     4819heat flux at k=0 (horizontal
     4820average), constant within Prandtl-layer (in K m/s).</td>
    48244821
    48254822
    … …  
    48564853
    48574854
    4858  <tr>
     4855 <tr><td align="undefined" valign="undefined"><span style="font-style: italic;">w"u"0</span></td><td align="undefined" valign="undefined">Subgrid-scale momentum flux (u-component) at k=0 (horizontal average), constant within Prandtl-layer (in m<sup>2</sup>/s<sup>2</sup>).</td></tr><tr><td align="undefined" valign="undefined"><span style="font-style: italic;">w"v"0</span></td><td align="undefined" valign="undefined">Subgrid-scale momentum flux (v-component) at k=0 (horizontal average), constant within Prandtl-layer (in m<sup>2</sup>/s<sup>2</sup>).</td></tr><tr><td align="undefined" valign="undefined"><span style="font-style: italic;">w"q"0</span></td><td align="undefined" valign="undefined">Subgrid-scale
     4856humidity flux at k=0 (horizontal average), constant within
     4857Prandtl-layer (in kg/kg m/s). Zero values are output if humidity is not
     4858used.</td></tr><tr>
    48594859
    48604860
    … …  
    48864886
    48874887
    4888  <tr>
    4889 
    4890 
    4891 
    4892             <td style="vertical-align: top; font-style: italic;">splptx</td>
    4893 
    4894 
    4895 
    4896             <td style="vertical-align: top;">Percentage of grid points
    4897 using upstream scheme along x with
    4898 upstream-spline advection switched on.</td>
    4899 
    4900 
    4901  </tr>
    4902 
    4903 
    4904  <tr>
    4905 
    4906 
    4907 
    4908             <td style="vertical-align: top; font-style: italic;">splpty</td>
    4909 
    4910 
    4911 
    4912             <td style="vertical-align: top;">Percentage of grid points
    4913 using upstream scheme along y with
    4914 upstream-spline
    4915 advection switched on.</td>
    4916 
    4917 
    4918  </tr>
    4919 
    4920 
    4921  <tr>
    4922 
    4923 
    4924  <td style="vertical-align: top; font-style: italic;">splptz</td>
    4925 
    4926 
    4927 
    4928             <td style="vertical-align: top;">Percentage of grid points
    4929 using upstream scheme along z with
    4930 upstream-spline
    4931 advection switched on.<br>
    4932 
    4933 
    4934  </td>
    4935 
    4936 
    4937  </tr>
     4888 
     4889
     4890
     4891 
     4892
     4893
     4894 
    49384895
    49394896
    … …  
    56405597It is only possible to gather data from a level larger than the one
    56415598determined automatically. A test run may be neccessary to determine
    5642 this level.</td>
     5599this level.<br><br>Setting of <span style="font-weight: bold;">mg_switch_to_pe0_level</span> = <span style="font-style: italic;">-1</span> prevents that data are collected on PE0 at all, i.e. coarsening of grids is limited by the subdomains.</td>
    56435600
    56445601
    … …  
    58485805
    58495806
    5850  <tr>
     5807 <tr><td style="vertical-align: top;"><a href="npex"><span style="font-weight: bold;"></span></a><span style="font-weight: bold;"><a name="npex"></a>npex</span></td><td style="vertical-align: top;">I</td><td style="vertical-align: top;"></td><td style="vertical-align: top;"><p>Number of processors
     5808along x-direction of the virtual
     5809processor
     5810net.&nbsp; </p>
     5811
     5812
     5813
     5814
     5815
     5816
     5817 
     5818     
     5819     
     5820     
     5821     
     5822     
     5823     
     5824      <p>For parallel runs, the total
     5825number of processors to be used
     5826is given by
     5827the <span style="font-weight: bold;">mrun</span>
     5828option <a href="http://www.muk.uni-hannover.de/software/mrun_beschreibung.html#Opt-X">-X</a>.
     5829By default, depending on the type of the parallel computer, PALM
     5830generates a 1d processor
     5831net (domain decomposition along x, <span style="font-weight: bold;">npey</span>
     5832= <span style="font-style: italic;">1</span>) or a
     58332d-net (this is
     5834favored on machines with fast communication network and/or large number of processors (&gt;256)). In case of a
     58352d-net, it is tried to make it more or less square-shaped. If, for
     5836example, 16 processors are assigned (-X 16), a 4 * 4 processor net is
     5837generated (<span style="font-weight: bold;">npex</span>
     5838= <span style="font-style: italic;">4</span>, <span style="font-weight: bold;">npey</span>
     5839= <span style="font-style: italic;">4</span>).
     5840This choice is optimal for square total domains (<a href="chapter_4.1.html#nx">nx</a>
     5841= <a href="chapter_4.1.html#ny">ny</a>),
     5842since then the number of ghost points at the lateral boundarys of
     5843the subdomains reaches a minimum. If <span style="font-weight: bold;">nx</span>
     5844and <span style="font-weight: bold;">ny</span>
     5845differ extremely, the
     5846processor net should be manually adjusted using adequate values for <span style="font-weight: bold;">npex</span> and <span style="font-weight: bold;">npey</span>.&nbsp; </p>
     5847
     5848
     5849
     5850
     5851
     5852
     5853
     5854     
     5855     
     5856     
     5857     
     5858     
     5859     
     5860      <p><b>Important:</b> The value of <span style="font-weight: bold;">npex</span> * <span style="font-weight: bold;">npey</span> must exactly match the
     5861value assigned by the <span style="font-weight: bold;">mrun</span>-option
     5862      <tt>-X</tt>.
     5863Otherwise the model run will abort with a corresponding error
     5864message.&nbsp; <br>
     5865
     5866
     5867
     5868
     5869
     5870
     5871
     5872Additionally, the specification of <span style="font-weight: bold;">npex</span>
     5873and <span style="font-weight: bold;">npey</span>
     5874may of course
     5875override the default setting for the domain decomposition (1d or 2d)
     5876which may have a significant (negative) effect on the code performance.
     5877      </p></td></tr><tr><td style="vertical-align: top;"><span style="font-weight: bold;"><a name="npey"></a>npey</span></td><td style="vertical-align: top;">I</td><td style="vertical-align: top;"></td><td style="vertical-align: top;"><p>Number of processors
     5878along y-direction of the virtual
     5879processor
     5880net.&nbsp; </p>
     5881
     5882
     5883
     5884
     5885
     5886
     5887 
     5888     
     5889     
     5890     
     5891     
     5892     
     5893     
     5894      <p>For further information see <a href="chapter_4.2.html#npex">npex</a>.</p></td></tr><tr>
    58515895
    58525896

  • palm/trunk/DOC/app/chapter_4.4.1.html

    r97 r197  
    151151virtual processor topology (grid) of 8*1 (x*y) is used. (<span style="font-weight: bold;">Note:</span> the user may
    152152adjust this
    153 default domain decomposition with the help of the parameters <a href="chapter_4.1.html#npex">npex</a>
    154 and <a href="chapter_4.1.html#npey">npey</a>).
     153default domain decomposition with the help of the parameters <a href="chapter_4.2.html#npex">npex</a>
     154and <a href="chapter_4.2.html#npey">npey</a>).
    155155</p><p style="line-height: 100%;">Provided that the
    156156parameters

  • palm/trunk/DOC/app/chapter_4.4.2.html

    r97 r197  
    155155virtual processor topology (grid) of 8*1 (x*y) is used. (<span style="font-weight: bold;">Note:</span> the user may
    156156adjust this
    157 default domain decomposition with the help of the parameters <a href="chapter_4.1.html#npex">npex</a>
    158 and <a href="chapter_4.1.html#npey">npey</a>).
     157default domain decomposition with the help of the parameters <a href="chapter_4.2.html#npex">npex</a>
     158and <a href="chapter_4.2.html#npey">npey</a>).
    159159</p><p style="line-height: 100%;">Provided that the
    160160parameters

  • palm/trunk/DOC/app/chapter_4.5.6.html

    r62 r197  
    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 4.5.6</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20041029;13171001"> <meta name="CHANGED" content="20041117;12013945"> <meta name="KEYWORDS" content="parallel LES model"> <style>
     3<meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM chapter 4.5.6</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20041029;13171001"> <meta name="CHANGED" content="20041117;12013945"> <meta name="KEYWORDS" content="parallel LES model"> <style>
    54<!--
    65@page { size: 21cm 29.7cm }
    76-->
    87</style></head>
    9 
    108<body style="direction: ltr;" lang="en-US"><h4 style="line-height: 100%;"><font size="4">4.5.6
    119Plots of compressed
    … …  
    2220= <i>.TRUE.</i> and by assigning the desired accuracy
    2321(number of
    24 significant digits) via <a href="chapter_4.2.html#do3d_precision">do3d_precision</a>.
     22significant digits) via <a href="chapter_4.2.html#do3d_comp_prec">do3d_comp_prec</a>.
    2523A separate accuracy must be assigned for each individual quantity.
    2624Using uncompressed data, each grid point value is

  • palm/trunk/DOC/app/chapter_4.6.html

    r189 r197  
    34443444
    34453445
    3446  <td><a href="chapter_4.2.html#do3d_precision"><span style="font-weight: bold;">do3d_precision</span></a></td>
     3446 <td><a href="chapter_4.2.html#do3d_comp_prec"><span style="font-weight: bold;">do3d_comp_prec</span></a></td>
    34473447
    34483448
    … …  
    64086408
    64096409
    6410     <tr>
     6410    <tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#inflow_damping_height"><span style="font-weight: bold;">inflow_damping_height</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">from precursor run</span></td><td align="undefined" valign="undefined">Height below which the turbulence signal is used for turbulence recycling (in m).</td></tr><tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#inflow_damping_width"><span style="font-weight: bold;">inflow_damping_width</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">0.1 * <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#inflow_damping_height">inflow_damping</a></span><a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#inflow_damping_height"><span style="font-style: italic;">_height</span></a></td><td align="undefined" valign="undefined">Transition range within which the turbulance signal is damped to zero (in m).</td></tr><tr>
    64116411
    64126412
    … …  
    77047704     
    77057705     
    7706       <p><a href="chapter_4.1.html#npex"><b>npex</b></a></p>
     7706      <p><a href="chapter_4.2.html#npex"><b>npex</b></a></p>
    77077707
    77087708
    … …  
    77177717     
    77187718     
    7719       <p>I</p>
     7719      <p>R</p>
    77207720
    77217721
    … …  
    77817781     
    77827782     
    7783       <p><a href="chapter_4.1.html#npey"><b>npey</b></a></p>
     7783      <p><a href="chapter_4.2.html#npey"><b>npey</b></a></p>
    77847784
    77857785
    … …  
    77947794     
    77957795     
    7796       <p>I</p>
     7796      <p>R</p>
    77977797
    77987798
    … …  
    1165011650
    1165111651
    11652  <tr>
     11652 <tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#recycling_width"><span style="font-weight: bold;">recycling_width</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">0.1 * <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#nx">nx</a> * <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.1.html#dx">dx</a></span></td><td align="undefined" valign="undefined">Distance of the recycling plane from the inflow boundary (in m).<br></td></tr><tr>
    1165311653
    1165411654
    … …  
    1390513905
    1390613906
    13907     <tr>
     13907    <tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#turbulent_inflow"><span style="font-weight: bold;">turbulent_inflow</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">L</td><td align="undefined" valign="undefined"><span style="font-style: italic;">.F.</span></td><td align="undefined" valign="undefined">Generates a turbulent inflow at side boundaries using a turbulence recycling method.</td></tr><tr>
    1390813908
    1390913909

  • palm/trunk/DOC/app/index.html

    r140 r197  
    160160 <u3:LatentStyles DefLockedState="false" LatentStyleCount="156">  </u3:LatentStyles>
    161161</xml><![endif]--></head>
    162 
    163162<body style="" lang="DE" link="blue" vlink="blue">
    164163
    … …  
    177176<p><span style="font-family: Thorndale;" lang="EN-GB">Current
    178177model
    179 version:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;3.4a<br>
     178version:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;3.5<br>
    180179For date of last change see bottom line of each page. <o:p></o:p></span></p>
    181180

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

    r139 r197  
    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.4a</span></td> </tr>
     14model version:</b></td> <td><span style="font-weight: bold;">3.5</span></td> </tr>
    1515<tr nosave=""> <td nosave=""><b>Last
    1616change of this document</b>:&nbsp;</td> <td nosave=""><b>$Id$</b></td> </tr>
    … …  
    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>
     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;"></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>
    23542354</tbody>
    23552355</table>&nbsp;<b><blink>Attention:</blink></b>

  • palm/trunk/SCRIPTS/mrun

    r196 r197  
    22492249          [[ $? != 0 ]]  &&  compile_error=true
    22502250          rm -rf  $TEMPDIR_COMPILE
     2251       elif [[ $localhost = lcsgib ]]
     2252       then
     2253          printf "      compiler is called via ssh on \"bicegate\" \n"
     2254          ssh  130.73.232.102  -l $usern  ". /usr/share/modules/init/bash; module load mpt; cd $TEMPDIR; make -f Makefile PROG=a.out  F90=$compiler_name  COPT=\"$cpp_options\"  F90FLAGS=\"$fopts\"  LDFLAGS=\"$lopts\" "
     2255          [[ ! -f a.out ]]  &&  compile_error=true
     2256          continue   # ANDERENFALLS IST STATUS=1, FALLS A.OUT VORHANDEN
     2257       elif [[ $localhost = lcsgih ]]
     2258       then
     2259          printf "      compiler is called via ssh on \"hicegate\" \n"
     2260          ssh  130.75.4.102  -l $usern  ". /usr/share/modules/init/bash; module load mpt; cd $TEMPDIR; make -f Makefile PROG=a.out  F90=$compiler_name  COPT=\"$cpp_options\"  F90FLAGS=\"$fopts\"  LDFLAGS=\"$lopts\" 2>&1 "
     2261#          ssh  130.75.4.102  -l $usern  ". /usr/share/modules/init/bash; module load mvapich2; cd $TEMPDIR; make -f Makefile PROG=a.out  F90=$compiler_name  COPT=\"$cpp_options\"  F90FLAGS=\"$fopts\"  LDFLAGS=\"$lopts\" 2>&1 "
     2262#          print  ". /usr/share/modules/init/bash; module load mpt; cd $TEMPDIR; make -f Makefile PROG=a.out  F90=$compiler_name  COPT=\"$cpp_options\"  F90FLAGS=\"$fopts\"  LDFLAGS=\"$lopts\" " | ssh  ${usern}@130.75.4.102  2>&1
     2263          [[ ! -f a.out ]]  &&  compile_error=true
     2264          continue   # ANDERENFALLS IST STATUS=1, FALLS A.OUT VORHANDEN
    22512265       else
    22522266          [[ $localhost = lctit ]]  &&  export LM_LICENSE_FILE=27050@tggls
    … …  
    22602274
    22612275  # WORKAROUND: REMOVE IF CONSTRUCT LATER, BUT KEEP THE EXIT!
    2262           if [[ $localhost != lcsgib  &&  $localhost != lcsgih ]]
    2263           then
     2276#          if [[ $localhost != lcsgib  &&  $localhost != lcsgih ]]
     2277#          then
    22642278             exit
    2265           else
    2266              locat=normal
    2267           fi
     2279#          else
     2280#             locat=normal
     2281#          fi
    22682282       else
    22692283          printf "$striche\n  *** compilation finished \n"
    … …  
    28392853                if [[ $host = lcsgih  ||  $host = lcsgib ]]
    28402854                then
    2841                    ulimit -s 30000   # A too large stack size causes problems
     2855#                   ulimit -s 300000   # A too large stack size causes problems
    28422856#                   export MV2_NUM_PORTS=2
    28432857#                   export MV2_CPU_MAPPING=0:1:2:3
    2844 #                   mpiexec    ./palm  $ROPTS
    2845                    mpiexec_mpt -np $ii   ./palm  $ROPTS
     2858#                   mpiexec    ./a.out  $ROPTS
     2859                   mpiexec_mpt -np $ii   ./a.out  $ROPTS
    28462860                elif [[ $host = lcxt4 ]]
    28472861                then
    … …  
    28592873                printf "\n\n"
    28602874
    2861                 if [[ $host = lcxt4 ]]
     2875                if [[ $host = lcsgih  ||  $host = lcsgib ]]
     2876                then
     2877
     2878                   mpiexec  -n $iii  a.out  $ROPTS  <  runfile_atmos &
     2879                   mpiexec  -n $iii  a.out  $ROPTS  <  runfile_ocean &
     2880#                   head  -n $iii  $PBS_NODEFILE  >  nodefile_atmos
     2881#                   echo "--- nodefile_atmos:"
     2882#                   cat nodefile_atmos
     2883#                   tail  -n $iii  $PBS_NODEFILE  >  nodefile_ocean
     2884#                   echo "--- nodefile_ocean:"
     2885#                   cat nodefile_ocean
     2886#                   export  PBS_NODEFILE=${PWD}/nodefile_atmos
     2887#                   mpiexec_mpt -np $iii   ./a.out  $ROPTS  <  runfile_atmos &
     2888#                   export  PBS_NODEFILE=${PWD}/nodefile_ocean
     2889#                   mpiexec_mpt -np $iii   ./a.out  $ROPTS  <  runfile_ocean &
     2890
     2891
     2892                elif [[ $host = lcxt4 ]]
    28622893                then
    28632894                   aprun  -n $iii  -N $tasks_per_node  a.out < runfile_atmos  $ROPTS  &

  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r192 r197  
    2828Particle advection adjusted for ocean runs.
    2929
     30Multigrid levels are limited by subdomains if mg_switch_to_pe0_level = -1.
     31
    3032User-defined spectra.
     33
     34Output of q*2 profile added.
    3135
    3236Argument -c introduced to most of the subjob calls, which allows the user to
    … …  
    9195(combine_plot_fields)
    9296
    93 advec_s_ups, advec_u_ups, advec_v_ups, advec_w_ups, calc_spectra, check_open, check_parameters, cpu_statistics, data_output_spectra, 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
     97Time coordinate t=0 stored on netcdf-file only if an output is required for
     98this time for at least one of the profiles
     99
     100Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0) added
     101(flow_statistics)
     102
     103advec_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
    94104
    95105
    … …  
    148158Bugfix: misplaced #endif directives (combine_plot_fields)
    149159
    150 calc_spectra, check_parameters, 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
     160data are collected from PE0 in an ordered sequence which seems to avoid hanging of processes on SGI-ICE (cpu_statistics)
    151161
     162calc_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/check_parameters.f90

    r181 r197  
    1010! Leaf area density (LAD) explicitly set to its surface value at k=0
    1111! Case of reading data for recycling included in initializing_actions,
    12 ! check of turbulent_inflow and calculation of recycling_plane
     12! check of turbulent_inflow and calculation of recycling_plane.
     13! q*2 profile added
    1314!
    1415! Former revisions:
    … …  
    22562257             hom(:,2,69,:) = SPREAD( zu, 2, statistic_regions+1 )
    22572258
     2259          CASE ( 'q*2' )
     2260             IF ( .NOT. humidity )  THEN
     2261                IF ( myid == 0 )  THEN
     2262                   PRINT*, '+++ check_parameters:  data_output_pr = ', &
     2263                           data_output_pr(i),                          &
     2264                           '    is not implemented for humidity = FALSE'
     2265                ENDIF
     2266                CALL local_stop
     2267             ELSE
     2268                dopr_index(i) = 70
     2269                dopr_unit(i)  = 'kg2/kg2'
     2270                hom(:,2,70,:) = SPREAD( zu, 2, statistic_regions+1 )
     2271             ENDIF
    22582272
    22592273          CASE DEFAULT

  • palm/trunk/SOURCE/cpu_statistics.f90

    r181 r197  
    44! Actual revisions:
    55! -----------------
    6 ! Format adjustments in order to allow CPU# > 999
     6! Format adjustments in order to allow CPU# > 999,
     7! data are collected from PE0 in an ordered sequence which seems to avoid
     8! hanging of processes on SGI-ICE
    79!
    810! Former revisions:
    … …  
    6971       DO  i = 1, numprocs-1
    7072          CALL MPI_RECV( pe_max(1), SIZE( log_point ), MPI_REAL, &
    71                          MPI_ANY_SOURCE, MPI_ANY_TAG, comm2d, status, ierr )
     73                         i, i, comm2d, status, ierr )
    7274          sender = status(MPI_SOURCE)
    7375          pe_log_points(:,sender) = pe_max
    … …  
    8789!--       Calculate rms
    8890          DO  i = 0, numprocs-1
    89 !             IF ( log_point(iii)%place == 'run_control' )  THEN
    90 !                PRINT*, 'pe_rms=',pe_rms(iii),' plp=',pe_log_points(iii,i), &
    91 !                        ' lps=',log_point(iii)%sum
    92 !             ENDIF
    9391             pe_rms(iii) = pe_rms(iii) + ( &
    9492                                 pe_log_points(iii,i) - log_point(iii)%sum &
    … …  
    103101       ALLOCATE( pe_max( SIZE( log_point ) ) )
    104102       pe_max = log_point%sum
    105        CALL MPI_SEND( pe_max(1), SIZE( log_point ), MPI_REAL, 0, 0, comm2d, &
     103       CALL MPI_SEND( pe_max(1), SIZE( log_point ), MPI_REAL, 0, myid, comm2d, &
    106104                      ierr )
    107105#endif

  • palm/trunk/SOURCE/data_output_profiles.f90

    r90 r197  
    44! Actual revisions:
    55! -----------------
    6 !
     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
    78!
    89! Former revisions:
    … …  
    3940
    4041    INTEGER ::  i, id, ilc, ils, j, k, sr
     42    LOGICAL ::  output_for_t0
    4143    REAL    ::  uxma, uxmi
    4244
    … …  
    101103#if defined( __netcdf )
    102104!
    103 !--          Store initial time (t=0) to time axis         
    104              nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, (/ 0.0 /), &
    105                                      start = (/ 1 /), count = (/ 1 /) )
    106              IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 329 )
     105!--          Store initial time (t=0) to time axis, but only if an output
     106!--          is required for at least one of the profiles
     107             output_for_t0 = .FALSE.
     108             DO  i = 1, dopr_n
     109                IF ( dopr_initial_index(i) /= 0 )  THEN
     110                   nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr,  &
     111                                           (/ 0.0 /), start = (/ 1 /), &
     112                                           count = (/ 1 /) )
     113                   IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 329 )
     114                   output_for_t0 = .TRUE.
     115                   EXIT
     116                ENDIF
     117             ENDDO
    107118
    108119!
    … …  
    331342          ENDDO   ! Loop over dopr_n for initial profiles
    332343
    333           IF ( netcdf_output )  dopr_time_count = dopr_time_count + 1
     344          IF ( netcdf_output  .AND.  output_for_t0 )  THEN
     345             dopr_time_count = dopr_time_count + 1
     346          ENDIF
    334347
    335348       ENDIF   ! Initial profiles

  • palm/trunk/SOURCE/flow_statistics.f90

    r142 r197  
    44! Actual revisions:
    55! -----------------
     6! Spline timeseries splptx etc. removed, timeseries w'u', w'v', w'q' (k=0)
     7! added
    68! Bugfix: divide sums(k,8) (e) and sums(k,34) (e*) by ngp_2dh_s_inner(k,sr)
    79! (like other scalars)
    … …  
    109111       sums_l(nzb+9,pr_palm,0)  = sums_divold_l(sr)  ! old divergence from pres
    110112       sums_l(nzb+10,pr_palm,0) = sums_divnew_l(sr)  ! new divergence from pres
    111 !--    WARNING: next four lines still may have to be adjusted for OpenMP
    112        sums_l(nzb:nzb+2,pr_palm-1,0)    = sums_up_fraction_l(1,1:3,sr)! upstream
    113        sums_l(nzb+3:nzb+5,pr_palm-1,0)  = sums_up_fraction_l(2,1:3,sr)! parts
    114        sums_l(nzb+6:nzb+8,pr_palm-1,0)  = sums_up_fraction_l(3,1:3,sr)! from
    115        sums_l(nzb+9:nzb+11,pr_palm-1,0) = sums_up_fraction_l(4,1:3,sr)! spline
    116113
    117114!
    … …  
    380377             sums_l(nzb+3,pr_palm,tn) = sums_l(nzb+3,pr_palm,tn) + &
    381378                                        ts(j,i)   * rmask(j,i,sr)
     379             IF ( humidity )  THEN
     380                sums_l(nzb+12,pr_palm,tn) = sums_l(nzb+12,pr_palm,tn) + &
     381                                            qs(j,i)   * rmask(j,i,sr)
     382             ENDIF
    382383          ENDDO
    383384       ENDDO
    … …  
    807808!--    above the topography, they are being divided by ngp_2dh(sr)
    808809       sums(nzb:nzb+3,pr_palm)    = sums(nzb:nzb+3,pr_palm)    / &
     810                                    ngp_2dh(sr)
     811       sums(nzb+12,pr_palm)       = sums(nzb+12,pr_palm)       / &    ! qs
    809812                                    ngp_2dh(sr)
    810813!--    eges, e*
    … …  
    882885       hom(:,1,68,sr) = sums(:,68)     ! w*p*
    883886       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
     887       hom(:,1,70,sr) = sums(:,70)     ! q*2
    884888
    885889       hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)
    … …  
    9991003       ts_value(17,sr) = hom(nzb,1,4,sr)            ! pt(0)
    10001004       ts_value(18,sr) = hom(nzb+1,1,4,sr)          ! pt(zp)
    1001        ts_value(19,sr) = hom(nzb+9,1,pr_palm-1,sr)  ! splptx
    1002        ts_value(20,sr) = hom(nzb+10,1,pr_palm-1,sr) ! splpty
    1003        ts_value(21,sr) = hom(nzb+11,1,pr_palm-1,sr) ! splptz
     1005       ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)    ! u'w'    at k=0
     1006       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
     1007       ts_value(21,sr) = hom(nzb+12,1,pr_palm,sr)   ! w'q'    at k=0
     1008
    10041009       IF ( ts_value(5,sr) /= 0.0 )  THEN
    10051010          ts_value(22,sr) = ts_value(4,sr)**2 / &

  • palm/trunk/SOURCE/header.f90

    r189 r197  
    44! Actual revisions:
    55! -----------------
     6! TEST: mg_switch_to_pe0_level = -1!!!!!!!!!
    67! allow 100 spectra levels instead of 10 for consistency with
    78! define_netcdf_header
    … …  
    117118    ELSEIF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 )  THEN
    118119       run_classification = '3D - run without 1D - prerun'
    119     ELSEIF ( INDEX(initializing_actions, 'set_1d-model_profiles') /= 0 ) THEN
     120    ELSEIF ( INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN
    120121       run_classification = '3D - run with 1D - prerun'
     122    ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 )  THEN
     123       run_classification = '3D - run initialized by user'
    121124    ELSE
    122125       PRINT*,'+++ header:  unknown action(s): ',initializing_actions
    … …  
    184187          WRITE ( io, 136 )  nxr_mg(1)-nxl_mg(1)+1, nyn_mg(1)-nys_mg(1)+1, &
    185188                             nzt_mg(1)
    186        ELSE
     189       ELSEIF (  mg_switch_to_pe0_level /= -1 )  THEN
    187190          WRITE ( io, 137 )  mg_switch_to_pe0_level,            &
    188191                             mg_loc_ind(2,0)-mg_loc_ind(1,0)+1, &

  • palm/trunk/SOURCE/init_pegrid.f90

    r181 r197  
    55! -----------------
    66! ATTENTION: nnz_x undefined problem still has to be solved!!!!!!!!
     7! multigrid levels are limited by subdomains if mg_switch_to_pe0_level = -1,
    78! nz is used instead nnz for calculating mg-levels
    89! Collect on PE0 horizontal index bounds from all other PEs,
    … …  
    719720!--    Find out, if the total domain allows more levels. These additional
    720721!--    levels are processed on PE0 only.
    721        IF ( numprocs > 1 )  THEN
     722       IF ( numprocs > 1  .AND.  mg_switch_to_pe0_level /= -1 )  THEN
    722723          IF ( mg_levels_z > MIN( mg_levels_x, mg_levels_y ) )  THEN
    723724             mg_switch_to_pe0_level_l = maximum_grid_level

  • palm/trunk/SOURCE/modules.f90

    r189 r197  
    788788             'vmax   ', 'wmax   ', 'div_new', 'div_old', 'z_i_wpt', 'z_i_pt ', &
    789789             'w*     ', 'w"pt"0 ', 'w"pt"  ', 'wpt    ', 'pt(0)  ', 'pt(zp) ', &
    790              'splptx ', 'splpty ', 'splptz ', 'mo_L   ',                       &
     790             'w"u"0  ', 'w"v"0  ', 'w"q"0 ', 'mo_L   ',                       &
    791791             ( 'unknown', i9 = 1, 78) /)
    792792
    … …  
    795795             'm/s    ', 'm/s    ', 's-1    ', 's-1    ', 'm      ', 'm      ', &
    796796             'm/s    ', 'K m/s  ', 'K m/s  ', 'k m/s  ', 'K      ', 'K      ', &
    797              '%      ', '%      ', '%      ', 'm      ',                       &
     797             'm2/s2  ', 'm2/s2  ', 'kg m/s ', 'm      ',                       &
    798798             ( 'unknown', i9 = 1, 78 ) /)
    799799

  • palm/trunk/SOURCE/palm.f90

    r164 r197  
    6666    INTEGER           ::  i, run_description_header_i(80)
    6767
    68     version = 'PALM 3.4a'
     68    version = 'PALM 3.5'
    6969
    7070#if defined( __parallel )

  • palm/trunk/SOURCE/parin.f90

    r153 r197  
    88! +inflow_damping_height, inflow_damping_width, recycling_width,
    99! turbulent_inflow in inipar
     10! -skip_time_dosp in d3par,
    1011! Allocation of hom_sum moved from init_3d_model to here,
    1112! npex, npey moved from inipar to d3par, setting of myid_char_14 removed,
    … …  
    109110             surface_heatflux, surface_pressure, surface_scalarflux, &
    110111             surface_waterflux, s_surface, s_surface_initial_change, &
    111              s_vertical_gradient, s_vertical_gradient_level, top_heatflux, &
    112              top_momentumflux_u, top_momentumflux_v, top_salinityflux, &
    113              timestep_scheme, topography, turbulent_inflow, ug_surface, &
     112             s_vertical_gradient, s_vertical_gradient_level, timestep_scheme,
     113             topography, top_heatflux, top_momentumflux_u, top_momentumflux_v, &
     114             top_salinityflux, turbulent_inflow, ug_surface, &
    114115             ug_vertical_gradient, ug_vertical_gradient_level, ups_limit_e, &
    115116             ups_limit_pt, ups_limit_u, ups_limit_v, ups_limit_w, &
    … …  
    143144                       section_yz, skip_time_data_output, &
    144145                       skip_time_data_output_av, skip_time_dopr, &
    145                        skip_time_dosp, skip_time_do2d_xy, skip_time_do2d_xz, &
     146                       skip_time_do2d_xy, skip_time_do2d_xz, &
    146147                       skip_time_do2d_yz, skip_time_do3d, &
    147148                       termination_time_needed, use_prior_plot1d_parameters, &
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |