<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 saw that there was a new commit on Github "import Paul's
      projection fix"  on October, 4th so I tried the newest Nek5000
      version with my setup. <br>
      Projection seems to be running fine if only temperature is solved
      additionally to pressure and velocity. However, when I want to
      solve for an additional passive scalar the projection fails (for
      both unperturbed initial conditions as well as chaotic fields).<br>
      <br>
      I tried to restrict the projection only for velocity and
      temperature in hsolve subroutine by adding the following code
      snippet:<br>
    </p>
    <p>...<br>
      if (ifstdh) then<br>
          call hmholtz(name,u,r,h1,h2,vmk,vml,imsh,tol,maxit,isd)<br>
      ! My restriction for using projection only on ifield 1 and 2 but
      not on 3 (first passive scalar)<br>
      elseif (ifield.gt.2) then<br>
          call hmholtz(name,u,r,h1,h2,vmk,vml,imsh,tol,maxit,isd)<br>
      else<br>
      ...<br>
      <br>
      <br>
      If I restrict the projection in such a way, I do not get errors
      for the projection of the temperature field but on the other hand
      no projection for passive scalar(s) is done either. <br>
      I have tested this in a simple test case which is uploaded here:<br>
    </p>
    <p><a class="moz-txt-link-freetext" href="https://bwsyncandshare.kit.edu/dl/fiVn2UB9QqWZN7awPsJHJSVF">https://bwsyncandshare.kit.edu/dl/fiVn2UB9QqWZN7awPsJHJSVF</a><br>
      <br>
      <br>
      I believe there is a problem when more than one temperature field
      is present and one activates projection (p94 != 0). Unfortunately,
      my knowledge of the code is not sufficient enough to really dig
      deeper and solve this problem.<br>
      The adjustments from my previous message (to extend approxt and
      napproxt) only postponed the projection errors to later timesteps.<br>
    </p>
    <p>Maybe you can point me in the right direction?<br>
      <br>
      Best Regards,<br>
      Steffen Straub<br>
      <br>
      <br>
      <br>
      <br>
    </p>
    <br>
    <div class="moz-cite-prefix">On 09/28/2017 07:00 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.23.1506618014.75224.nek5000-users@lists.mcs.anl.gov">
      <pre wrap="">Message: 1
Date: Thu, 28 Sep 2017 14:40:27 +0200
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>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
Message-ID:
        <a class="moz-txt-link-rfc2396E" href="mailto:mailman.558.1506602432.6921.nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><mailman.558.1506602432.6921.nek5000-users@lists.mcs.anl.gov></a>
Content-Type: text/plain; charset="utf-8"; Format="flowed"

Hi Neks,

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.
Can you please confirm that my adjustments are reasonable?

In conduct.f when hsolve is called, I have adjusted the argument list 
like this:

-????? napproxt(1) = laxtt? ! Fix this... pff 10/10/15
+????? napproxt(1,ifield-1) = laxtt? ! Fix this... pff 10/10/15

 ??????????? call hsolve? (name4t,TA,TB,H1,H2
 ????? $?????????????????? ,tmask(1,1,1,1,ifield-1)
 ????? $?????????????????? ,tmult(1,1,1,1,ifield-1)
 ????? $?????????????????? ,imesh,tolht(ifield),nmxh,1
-???? $?????????????????? ,approxt,napproxt,bintm1)
+???? $ ,approxt(1,1,ifield-1),napproxt(1,ifield-1)
+???? $?????????????????? ,bintm1)

such that approxt and napproxt refers to each ifield, ie. temperature, 
passive scalar 1, passive scalar 2, ... similar to velx, vely, velz.


Accordingly, I have adjusted ORTHOT:

-????? integer???????? napproxt(2)
+
+C Add another dimension to (n)approxt for each passive scalar and temp.
+????? integer???????? napproxt(2,ldimt)
 ?????? common <i class="moz-txt-slash"><span class="moz-txt-tag">/</span>trthoi<span class="moz-txt-tag">/</span></i> napproxt

-????? real??????????? approxt(ktott,0:laxtt)
+????? real??????????? approxt(ktott,0:laxtt,ldimt)

and induct.f

-????? napproxt(1) = laxtt
+????? napproxt(1,ifield-1) = laxtt





Message: 2
Date: Thu, 28 Sep 2017 14:04:58 +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.565.1506607685.6921.nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><mailman.565.1506607685.6921.nek5000-users@lists.mcs.anl.gov></a>
Content-Type: text/plain; charset="us-ascii"


Dear Steffen,


I believe this issue may have been resolved in an update to the git repo

last week.


Best,


Paul</pre>
    </blockquote>
    <br>
  </body>
</html>