<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
  </head>
  <body text="#000000" bgcolor="#FFFFFF">
    <p>Hi Neks,<br>
      <br>
      I might have fixed my proj_ortho error issue. At least it is
      running now without errors on a small test case on my local
      machine. However, I am very much unsure if what I was changing is
      correct.<br>
      Can you please confirm that my adjustments are reasonable?<br>
      <br>
      In conduct.f when hsolve is called, I have adjusted the argument
      list like this:</p>
    <p>-      napproxt(1) = laxtt  ! Fix this... pff 10/10/15<br>
      +      napproxt(1,ifield-1) = laxtt  ! Fix this... pff 10/10/15<br>
       <br>
                  call hsolve  (name4t,TA,TB,H1,H2 <br>
            $                   ,tmask(1,1,1,1,ifield-1)<br>
            $                   ,tmult(1,1,1,1,ifield-1)<br>
            $                   ,imesh,tolht(ifield),nmxh,1<br>
      -     $                   ,approxt,napproxt,bintm1)<br>
      +     $                  
      ,approxt(1,1,ifield-1),napproxt(1,ifield-1)<br>
      +     $                   ,bintm1)<br>
      <br>
      such that approxt and napproxt refers to each ifield, ie.
      temperature, passive scalar 1, passive scalar 2, ... similar to
      velx, vely, velz.</p>
    <p><br>
    </p>
    <p>Accordingly, I have adjusted ORTHOT:</p>
    <p>-      integer         napproxt(2)<br>
      +<br>
      +C Add another dimension to (n)approxt for each passive scalar and
      temp.<br>
      +      integer         napproxt(2,ldimt)<br>
             common /trthoi/ napproxt<br>
       <br>
      -      real            approxt(ktott,0:laxtt)<br>
      +      real            approxt(ktott,0:laxtt,ldimt)<br>
      <br>
    </p>
    <p>and induct.f</p>
    <p>-      napproxt(1) = laxtt<br>
      +      napproxt(1,ifield-1) = laxtt<br>
      <br>
    </p>
    <p><br>
    </p>
    <pre wrap="">
</pre>
    <br>
    <div class="moz-cite-prefix">On 09/14/2017 06:29 PM,
      <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a> wrote:<br>
    </div>
    <blockquote type="cite"
      cite="mid:mailman.780.1505406553.10190.nek5000-users@lists.mcs.anl.gov">
      <pre wrap="">Date: Thu, 14 Sep 2017 15:34:39 +0000
From: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">nek5000-users@lists.mcs.anl.gov</a>
To: <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">"nek5000-users@lists.mcs.anl.gov"</a>
        <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><nek5000-users@lists.mcs.anl.gov></a>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
Message-ID:
        <a class="moz-txt-link-rfc2396E" href="mailto:mailman.770.1505403282.10190.nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><mailman.770.1505403282.10190.nek5000-users@lists.mcs.anl.gov></a>
Content-Type: text/plain; charset="windows-1252"


Hi Steffen,


Thanks for bringing this to our attention.


Adam and I have been looking into this... (no doubt, it's something

stupid done by me awhile back).


Hope to have it resolved by next week.


Best,


Paul


________________________________
From: Nek5000-users <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" moz-do-not-send="true"><nek5000-users-bounces@lists.mcs.anl.gov></a> on behalf of <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">nek5000-users@lists.mcs.anl.gov</a> <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><nek5000-users@lists.mcs.anl.gov></a>
Sent: Thursday, September 14, 2017 8:18:27 AM
To: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">nek5000-users@lists.mcs.anl.gov</a>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow

Hi,

I have made some more tests to find the cause of this issue.
Attached is a laminar pipe flow simulation with temperature and one
passive scalar. The original setup is from Paul's "nusselt.tgz", where
he answered another question of mine.

I found that I can reproduce the same error I see in my turbulent
simulations when there is an additional passive scalar with a different
type of boundary condition, i.e. Neumann boundary condition at the pipe
wall for temperature and Dirichlet boundary condition for the passive
scalar.

Can someone please check my setup and tell me why this causes the
projection to fail?


I hope to save some computational resources by using the projection.
However, if it is not possible, I will keep param(94) = 0.

Kind Regards,
Steffen Straub</pre>
    </blockquote>
    <br>
  </body>
</html>