<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p><br>
</p>
<p>Hi,</p>
<p><br>
</p>
<p>Have you tested your routine on plane Poiseiulle flow,</p>
<p>u=1-y^2 ?</p>
<p><br>
</p>
<p>Paul</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Friday, June 23, 2017 8:12:37 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Help: compute dissipation</font>
<div> </div>
</div>
<div>Hi NEK users,
<div class=""><br class="">
</div>
<div class="">I am trying to calculate dissipation for fully developed turbulent channel flow. I tried the code below, but result I got from this totally disagrees with Vreman’s 2014 paper. x-axis is streamwise, y-axis is wall-normal and z-axis is spanwise.
 Any errors in the code? Any thoughts?</div>
<div class=""><br class="">
</div>
<div class=""><!--[if gte mso 9]><xml>
 <o:OfficeDocumentSettings>
  <o:AllowPNG/>
 </o:OfficeDocumentSettings>
</xml><![endif]--><!--[if gte mso 9]><xml>
 <w:WordDocument>
  <w:View>Normal</w:View>
  <w:Zoom>0</w:Zoom>
  <w:TrackMoves/>
  <w:TrackFormatting/>
  <w:PunctuationKerning/>
  <w:ValidateAgainstSchemas/>
  <w:SaveIfXMLInvalid>false</w:SaveIfXMLInvalid>
  <w:IgnoreMixedContent>false</w:IgnoreMixedContent>
  <w:AlwaysShowPlaceholderText>false</w:AlwaysShowPlaceholderText>
  <w:DoNotPromoteQF/>
  <w:LidThemeOther>EN-GB</w:LidThemeOther>
  <w:LidThemeAsian>X-NONE</w:LidThemeAsian>
  <w:LidThemeComplexScript>X-NONE</w:LidThemeComplexScript>
  <w:Compatibility>
   <w:BreakWrappedTables/>
   <w:SnapToGridInCell/>
   <w:WrapTextWithPunct/>
   <w:UseAsianBreakRules/>
   <w:DontGrowAutofit/>
   <w:SplitPgBreakAndParaMark/>
   <w:EnableOpenTypeKerning/>
   <w:DontFlipMirrorIndents/>
   <w:OverrideTableStyleHps/>
  </w:Compatibility>
  <m:mathPr>
   <m:mathFont m:val="Cambria Math"/>
   <m:brkBin m:val="before"/>
   <m:brkBinSub m:val="--"/>
   <m:smallFrac m:val="off"/>
   <m:dispDef/>
   <m:lMargin m:val="0"/>
   <m:rMargin m:val="0"/>
   <m:defJc m:val="centerGroup"/>
   <m:wrapIndent m:val="1440"/>
   <m:intLim m:val="subSup"/>
   <m:naryLim m:val="undOvr"/>
  </m:mathPr></w:WordDocument>
</xml><![endif]--><!--[if gte mso 9]><xml>
 <w:LatentStyles DefLockedState="false" DefUnhideWhenUsed="false"
  DefSemiHidden="false" DefQFormat="false" DefPriority="99"
  LatentStyleCount="382">
  <w:LsdException Locked="false" Priority="0" QFormat="true" Name="Normal"/>
  <w:LsdException Locked="false" Priority="9" QFormat="true" Name="heading 1"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 2"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 3"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 4"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 5"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 6"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 7"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 8"/>
  <w:LsdException Locked="false" Priority="9" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="heading 9"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 6"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 7"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 8"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index 9"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 1"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 2"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 3"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 4"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 5"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 6"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 7"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 8"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" Name="toc 9"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Normal Indent"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="footnote text"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="annotation text"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="header"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="footer"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="index heading"/>
  <w:LsdException Locked="false" Priority="35" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="caption"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="table of figures"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="envelope address"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="envelope return"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="footnote reference"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="annotation reference"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="line number"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="page number"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="endnote reference"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="endnote text"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="table of authorities"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="macro"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="toa heading"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Bullet"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Number"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Bullet 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Bullet 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Bullet 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Bullet 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Number 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Number 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Number 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Number 5"/>
  <w:LsdException Locked="false" Priority="10" QFormat="true" Name="Title"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Closing"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Signature"/>
  <w:LsdException Locked="false" Priority="1" SemiHidden="true"
   UnhideWhenUsed="true" Name="Default Paragraph Font"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text Indent"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Continue"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Continue 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Continue 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Continue 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="List Continue 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Message Header"/>
  <w:LsdException Locked="false" Priority="11" QFormat="true" Name="Subtitle"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Salutation"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Date"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text First Indent"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text First Indent 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Heading"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text Indent 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Body Text Indent 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Block Text"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Hyperlink"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="FollowedHyperlink"/>
  <w:LsdException Locked="false" Priority="22" QFormat="true" Name="Strong"/>
  <w:LsdException Locked="false" Priority="20" QFormat="true" Name="Emphasis"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Document Map"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Plain Text"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="E-mail Signature"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Top of Form"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Bottom of Form"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Normal (Web)"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Acronym"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Address"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Cite"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Code"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Definition"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Keyboard"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Preformatted"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Sample"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Typewriter"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="HTML Variable"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Normal Table"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="annotation subject"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="No List"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Outline List 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Outline List 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Outline List 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Simple 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Simple 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Simple 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Classic 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Classic 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Classic 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Classic 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Colorful 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Colorful 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Colorful 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Columns 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Columns 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Columns 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Columns 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Columns 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 6"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 7"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Grid 8"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 6"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 7"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table List 8"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table 3D effects 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table 3D effects 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table 3D effects 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Contemporary"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Elegant"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Professional"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Subtle 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Subtle 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Web 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Web 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Web 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Balloon Text"/>
  <w:LsdException Locked="false" Priority="39" Name="Table Grid"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Table Theme"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 1"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 2"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 3"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 4"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 5"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 6"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 7"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 8"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Note Level 9"/>
  <w:LsdException Locked="false" SemiHidden="true" Name="Placeholder Text"/>
  <w:LsdException Locked="false" Priority="1" QFormat="true" Name="No Spacing"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 1"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List Accent 1"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 1"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 1"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 1"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 1"/>
  <w:LsdException Locked="false" SemiHidden="true" Name="Revision"/>
  <w:LsdException Locked="false" Priority="34" QFormat="true"
   Name="List Paragraph"/>
  <w:LsdException Locked="false" Priority="29" QFormat="true" Name="Quote"/>
  <w:LsdException Locked="false" Priority="30" QFormat="true"
   Name="Intense Quote"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 1"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 1"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 1"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 1"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List Accent 1"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 1"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 1"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 1"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 2"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List Accent 2"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 2"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 2"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 2"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 2"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 2"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 2"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 2"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 2"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List Accent 2"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 2"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 2"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 2"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 3"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List Accent 3"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 3"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 3"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 3"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 3"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 3"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 3"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 3"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 3"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List Accent 3"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 3"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 3"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 3"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 4"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List Accent 4"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 4"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 4"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 4"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 4"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 4"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 4"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 4"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 4"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List Accent 4"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 4"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 4"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 4"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 5"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List Accent 5"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 5"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 5"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 5"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 5"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 5"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 5"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 5"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 5"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List Accent 5"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 5"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 5"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 5"/>
  <w:LsdException Locked="false" Priority="60" Name="Light Shading Accent 6"/>
  <w:LsdException Locked="false" Priority="61" Name="Light List Accent 6"/>
  <w:LsdException Locked="false" Priority="62" Name="Light Grid Accent 6"/>
  <w:LsdException Locked="false" Priority="63" Name="Medium Shading 1 Accent 6"/>
  <w:LsdException Locked="false" Priority="64" Name="Medium Shading 2 Accent 6"/>
  <w:LsdException Locked="false" Priority="65" Name="Medium List 1 Accent 6"/>
  <w:LsdException Locked="false" Priority="66" Name="Medium List 2 Accent 6"/>
  <w:LsdException Locked="false" Priority="67" Name="Medium Grid 1 Accent 6"/>
  <w:LsdException Locked="false" Priority="68" Name="Medium Grid 2 Accent 6"/>
  <w:LsdException Locked="false" Priority="69" Name="Medium Grid 3 Accent 6"/>
  <w:LsdException Locked="false" Priority="70" Name="Dark List Accent 6"/>
  <w:LsdException Locked="false" Priority="71" Name="Colorful Shading Accent 6"/>
  <w:LsdException Locked="false" Priority="72" Name="Colorful List Accent 6"/>
  <w:LsdException Locked="false" Priority="73" Name="Colorful Grid Accent 6"/>
  <w:LsdException Locked="false" Priority="19" QFormat="true"
   Name="Subtle Emphasis"/>
  <w:LsdException Locked="false" Priority="21" QFormat="true"
   Name="Intense Emphasis"/>
  <w:LsdException Locked="false" Priority="31" QFormat="true"
   Name="Subtle Reference"/>
  <w:LsdException Locked="false" Priority="32" QFormat="true"
   Name="Intense Reference"/>
  <w:LsdException Locked="false" Priority="33" QFormat="true" Name="Book Title"/>
  <w:LsdException Locked="false" Priority="37" SemiHidden="true"
   UnhideWhenUsed="true" Name="Bibliography"/>
  <w:LsdException Locked="false" Priority="39" SemiHidden="true"
   UnhideWhenUsed="true" QFormat="true" Name="TOC Heading"/>
  <w:LsdException Locked="false" Priority="41" Name="Plain Table 1"/>
  <w:LsdException Locked="false" Priority="42" Name="Plain Table 2"/>
  <w:LsdException Locked="false" Priority="43" Name="Plain Table 3"/>
  <w:LsdException Locked="false" Priority="44" Name="Plain Table 4"/>
  <w:LsdException Locked="false" Priority="45" Name="Plain Table 5"/>
  <w:LsdException Locked="false" Priority="40" Name="Grid Table Light"/>
  <w:LsdException Locked="false" Priority="46" Name="Grid Table 1 Light"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark"/>
  <w:LsdException Locked="false" Priority="51" Name="Grid Table 6 Colorful"/>
  <w:LsdException Locked="false" Priority="52" Name="Grid Table 7 Colorful"/>
  <w:LsdException Locked="false" Priority="46"
   Name="Grid Table 1 Light Accent 1"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 1"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 1"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 1"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 1"/>
  <w:LsdException Locked="false" Priority="51"
   Name="Grid Table 6 Colorful Accent 1"/>
  <w:LsdException Locked="false" Priority="52"
   Name="Grid Table 7 Colorful Accent 1"/>
  <w:LsdException Locked="false" Priority="46"
   Name="Grid Table 1 Light Accent 2"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 2"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 2"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 2"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 2"/>
  <w:LsdException Locked="false" Priority="51"
   Name="Grid Table 6 Colorful Accent 2"/>
  <w:LsdException Locked="false" Priority="52"
   Name="Grid Table 7 Colorful Accent 2"/>
  <w:LsdException Locked="false" Priority="46"
   Name="Grid Table 1 Light Accent 3"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 3"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 3"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 3"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 3"/>
  <w:LsdException Locked="false" Priority="51"
   Name="Grid Table 6 Colorful Accent 3"/>
  <w:LsdException Locked="false" Priority="52"
   Name="Grid Table 7 Colorful Accent 3"/>
  <w:LsdException Locked="false" Priority="46"
   Name="Grid Table 1 Light Accent 4"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 4"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 4"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 4"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 4"/>
  <w:LsdException Locked="false" Priority="51"
   Name="Grid Table 6 Colorful Accent 4"/>
  <w:LsdException Locked="false" Priority="52"
   Name="Grid Table 7 Colorful Accent 4"/>
  <w:LsdException Locked="false" Priority="46"
   Name="Grid Table 1 Light Accent 5"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 5"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 5"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 5"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 5"/>
  <w:LsdException Locked="false" Priority="51"
   Name="Grid Table 6 Colorful Accent 5"/>
  <w:LsdException Locked="false" Priority="52"
   Name="Grid Table 7 Colorful Accent 5"/>
  <w:LsdException Locked="false" Priority="46"
   Name="Grid Table 1 Light Accent 6"/>
  <w:LsdException Locked="false" Priority="47" Name="Grid Table 2 Accent 6"/>
  <w:LsdException Locked="false" Priority="48" Name="Grid Table 3 Accent 6"/>
  <w:LsdException Locked="false" Priority="49" Name="Grid Table 4 Accent 6"/>
  <w:LsdException Locked="false" Priority="50" Name="Grid Table 5 Dark Accent 6"/>
  <w:LsdException Locked="false" Priority="51"
   Name="Grid Table 6 Colorful Accent 6"/>
  <w:LsdException Locked="false" Priority="52"
   Name="Grid Table 7 Colorful Accent 6"/>
  <w:LsdException Locked="false" Priority="46" Name="List Table 1 Light"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark"/>
  <w:LsdException Locked="false" Priority="51" Name="List Table 6 Colorful"/>
  <w:LsdException Locked="false" Priority="52" Name="List Table 7 Colorful"/>
  <w:LsdException Locked="false" Priority="46"
   Name="List Table 1 Light Accent 1"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 1"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 1"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 1"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 1"/>
  <w:LsdException Locked="false" Priority="51"
   Name="List Table 6 Colorful Accent 1"/>
  <w:LsdException Locked="false" Priority="52"
   Name="List Table 7 Colorful Accent 1"/>
  <w:LsdException Locked="false" Priority="46"
   Name="List Table 1 Light Accent 2"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 2"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 2"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 2"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 2"/>
  <w:LsdException Locked="false" Priority="51"
   Name="List Table 6 Colorful Accent 2"/>
  <w:LsdException Locked="false" Priority="52"
   Name="List Table 7 Colorful Accent 2"/>
  <w:LsdException Locked="false" Priority="46"
   Name="List Table 1 Light Accent 3"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 3"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 3"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 3"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 3"/>
  <w:LsdException Locked="false" Priority="51"
   Name="List Table 6 Colorful Accent 3"/>
  <w:LsdException Locked="false" Priority="52"
   Name="List Table 7 Colorful Accent 3"/>
  <w:LsdException Locked="false" Priority="46"
   Name="List Table 1 Light Accent 4"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 4"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 4"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 4"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 4"/>
  <w:LsdException Locked="false" Priority="51"
   Name="List Table 6 Colorful Accent 4"/>
  <w:LsdException Locked="false" Priority="52"
   Name="List Table 7 Colorful Accent 4"/>
  <w:LsdException Locked="false" Priority="46"
   Name="List Table 1 Light Accent 5"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 5"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 5"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 5"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 5"/>
  <w:LsdException Locked="false" Priority="51"
   Name="List Table 6 Colorful Accent 5"/>
  <w:LsdException Locked="false" Priority="52"
   Name="List Table 7 Colorful Accent 5"/>
  <w:LsdException Locked="false" Priority="46"
   Name="List Table 1 Light Accent 6"/>
  <w:LsdException Locked="false" Priority="47" Name="List Table 2 Accent 6"/>
  <w:LsdException Locked="false" Priority="48" Name="List Table 3 Accent 6"/>
  <w:LsdException Locked="false" Priority="49" Name="List Table 4 Accent 6"/>
  <w:LsdException Locked="false" Priority="50" Name="List Table 5 Dark Accent 6"/>
  <w:LsdException Locked="false" Priority="51"
   Name="List Table 6 Colorful Accent 6"/>
  <w:LsdException Locked="false" Priority="52"
   Name="List Table 7 Colorful Accent 6"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Mention"/>
  <w:LsdException Locked="false" SemiHidden="true" UnhideWhenUsed="true"
   Name="Smart Hyperlink"/>
 </w:LatentStyles>
</xml><![endif]--><!--[if gte mso 10]>
<style>
 /* Style Definitions */
table.MsoNormalTable
        {mso-style-name:"Table Normal";
        mso-tstyle-rowband-size:0;
        mso-tstyle-colband-size:0;
        mso-style-noshow:yes;
        mso-style-priority:99;
        mso-style-parent:"";
        mso-padding-alt:0in 5.4pt 0in 5.4pt;
        mso-para-margin-top:0in;
        mso-para-margin-right:0in;
        mso-para-margin-bottom:8.0pt;
        mso-para-margin-left:0in;
        line-height:107%;
        mso-pagination:widow-orphan;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;
        mso-ascii-font-family:Calibri;
        mso-ascii-theme-font:minor-latin;
        mso-hansi-font-family:Calibri;
        mso-hansi-theme-font:minor-latin;
        mso-fareast-language:EN-US;}
</style>
<![endif]--><!--StartFragment-->
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     real diss_sqr(lx1,ly1,lz1, lelt)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      real s11<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      real s12<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      real s13<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      real s22<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      real s23<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      real s33<o:p class=""></o:p></span></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">  common /plane/  uavg_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  vavg_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  wavg_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  urms_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  vrms_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  wrms_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  uvms_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  diss_sqr_pl(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  yy(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  w1(ly1*lely),w2(ly1*lely)<o:p class=""></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Consolas",sans-serif;mso-bidi-font-family:
Consolas" class="">     $             ,  ffx_avg, dragx_avg<o:p class=""></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Consolas",sans-serif;mso-bidi-font-family:
Consolas" class=""> </span></p>
<p class="MsoNormal"><span style="font-family:"Consolas",sans-serif;mso-bidi-font-family:
Consolas" class=""> </span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">parameter (lr=lx1*ly1*lz1)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      common /scruz/  ur(lr),us(lr),ut(lr)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  vr(lr),vs(lr),vt(lr)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $             ,  wr(lr),ws(lr),wt(lr)<o:p class=""></o:p></span></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">common /avg/    uavg(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  vavg(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  wavg(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  urms(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  vrms(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  wrms(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  uvms(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     &             ,  dissavg(lx1,ly1,lz1,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">common /scrns/  sij (lx1*ly1*lz1,6,lelv)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class=""> </span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">if(icalld.eq.0) then<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(uavg,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(urms,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(vrms,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(wrms,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(uvms,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(dissavg,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(wavg,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(vavg,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">        call rzero(diss_sqr,n)<o:p class=""></o:p></span></p>
<p class="MsoNormal"><o:p class=""> </o:p> </p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">call comp_sij(sij,6,vx,vy,vz,ur,us,ut,vr,vs,vt,wr,ws,wt)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class=""> </span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family: Consolas, sans-serif;" class="">       </span><span style="font-family: Consolas, sans-serif;" class="">do e=1,nelt</span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">          do i = 1,nx1*ny1*nz1<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">            s11=sij(i,1,e)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">            s12=sij(i,4,e)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">            s22=sij(i,2,e)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">            s13=sij(i,6,e)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">            s23=sij(i,5,e)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">            s33=sij(i,3,e)<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      
<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">          diss_sqr(i,1,1,e) = (s11*s11) +<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $                       (s22*s22) +
<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $                       (s33*s33) +<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $                       (s12*s12)*2 + (s23*s23)*2 +<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">     $                       (s13*s13)*2<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class=""> </span><span style="font-family: Consolas, sans-serif;" class="">             </span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">         enddo<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">         enddo<o:p class=""></o:p></span></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<p class="MsoNormal"><span style="font-family: Consolas, sans-serif;" class="">call avg1(dissavg,diss_sqr,alpha,beta,n,'diss',ifverbose)</span></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<div class=""><span style="font-family: Consolas, sans-serif;" class="">call planar_average_s(diss_sqr_pl,dissavg,w1,w2)</span></div>
<p class="MsoNormal"><o:p class=""> </o:p> </p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">do i=1,ny1*nely/2<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           uavg_pl(i) = 0.5 * (uavg_pl(i) + uavg_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           vavg_pl(i) = 0.5 * (vavg_pl(i) + vavg_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           wavg_pl(i) = 0.5 * (wavg_pl(i) + wavg_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           diss_sqr_pl(i) = 0.5 * (diss_sqr_pl(i) + diss_sqr_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           urms_pl(i) = 0.5 * (urms_pl(i) + urms_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           vrms_pl(i) = 0.5 * (vrms_pl(i) + vrms_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">           wrms_pl(i) = 0.5 * (wrms_pl(i) + wrms_pl(m-i+1))<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class=""> </span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">       enddo<o:p class=""></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:
normal;tab-stops:48.0pt 96.0pt 2.0in 192.0pt 240.0pt 4.0in 336.0pt 384.0pt 6.0in 480.0pt 528.0pt 8.0in 624.0pt 672.0pt 10.0in 768.0pt 816.0pt 12.0in 912.0pt 960.0pt 14.0in 1056.0pt 1104.0pt 16.0in 1200.0pt 1248.0pt 18.0in 1344.0pt 1392.0pt 20.0in 1488.0pt 1536.0pt;
mso-layout-grid-align:none;text-autospace:none">
<span style="font-family:"Consolas",sans-serif;
mso-bidi-font-family:Consolas" class="">      <o:p class=""></o:p></span></p>
<p class="MsoNormal"><o:p class=""> </o:p></p>
<!--EndFragment--></div>
<div class=""><br class="">
</div>
<div class="">Many thanks.</div>
<div class="">
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
</div>
</div>
</body>
</html>