<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body class="" style="word-wrap:break-word" fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Dear JC,
<div><br>
</div>
<div>Yes --- just pull a copy of the common block out of torq_calc in navier5.f:</div>
<div><br>
</div>
<div><br>
</div>
<div>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
      common /ctorq/ dragx(0:maxobj),dragpx(0:maxobj),dragvx(0:maxobj)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , dragy(0:maxobj),dragpy(0:maxobj),dragvy(0:maxobj)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , dragz(0:maxobj),dragpz(0:maxobj),dragvz(0:maxobj)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
c</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , torqx(0:maxobj),torqpx(0:maxobj),torqvx(0:maxobj)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , torqy(0:maxobj),torqpy(0:maxobj),torqvy(0:maxobj)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , torqz(0:maxobj),torqpz(0:maxobj),torqvz(0:maxobj)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
c</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , dpdx_mean,dpdy_mean,dpdz_mean</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
     $             , dgtq(3,4)</p>
<p style="margin-right: 0px; margin-left: 0px; line-height: normal; font-family: Courier;">
c</p>
<div><br>
</div>
<div><br>
</div>
<div>There you have access to all components, including breakdown into pressure and viscous</div>
<div>contributions.  </div>
<div><br>
</div>
<div>If you have only one object (or surface),</div>
<div><br>
</div>
<div>      dragx(0) = force in x direction  (typ., = drag)</div>
<div>      dragy(0) = force in y direction (typ., = lift)</div>
<div><br>
</div>
<div>If you have multiple objects/surfaces, then dragx(k)=force in x direction on object k, k=1,...</div>
<div>and dragx(0) = sum_k=1,.. dragx(k)</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF303273" style="direction: ltr;"><font face="Tahoma" size="2" 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> Thursday, March 30, 2017 3:33 AM<br>
<b>To:</b> Nek 5000<br>
<b>Subject:</b> [Nek5000-users] Accessing the value of drag and lift coefficients in usrchk<br>
</font><br>
</div>
<div></div>
<div>Hi Nek,
<div class=""><br class="">
</div>
<div class="">I was wondering whether it was possible to access the values of the drag and lift coefficient directly into usrchk?</div>
<div class="">I know how to compute them but they only get printed in the logfile. I would need them for a control problem.</div>
<div class=""><br class="">
</div>
<div class="">Thanks a lot</div>
<div class=""><br class="">
<div class="">
<div class="" style="word-wrap:break-word">
<div style="color:rgb(0,0,0); font-family:Helvetica; font-size:12px; font-style:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">
Jean-Christophe B. Loiseau</div>
<div style="color:rgb(0,0,0); font-family:Helvetica; font-size:12px; font-style:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">
Maître de conférence</div>
<div style="color:rgb(0,0,0); font-family:Helvetica; font-size:12px; font-style:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">
DynFluid, Arts et Métiers ParisTech</div>
<div style="color:rgb(0,0,0); font-family:Helvetica; font-size:12px; font-style:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">
<br class="">
</div>
<div style="color:rgb(0,0,0); font-family:Helvetica; font-size:12px; font-style:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">
<a href="https://sites.google.com/site/loiseaujc/" class="" target="_blank">Homepage</a></div>
</div>
</div>
<br class="">
</div>
</div>
</div>
</div>
</div>
</body>
</html>