<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 fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Dear Sandeep,
<div><br>
</div>
<div>You need to declare ucbx(lx1,ly1,lz1,lelt), etc.</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF845899" 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> Monday, January 16, 2017 5:33 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Rank mismatch in array reference in qvol<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Hello All,</div>
<div class="gmail_default" style="font-size:small">I would like to solve <br>
</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">dA/dt = Laplacian A + U cross (curl A) </div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">'A' is a combination of three scalar arrays, T, PS1 and PS2.<br>
</div>
<div class="gmail_default" style="font-size:small">I modified the following to switch off the advection for 'A' in the .rea file<br>
</div>
<div class="gmail_default" style="font-size:small"> 2.0000000E+00 p23 NPSCAL<br>
</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default">
<div class="gmail_default"> T IFFLOW<br>
</div>
</div>
</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default"> T IFHEAT</div>
<div class="gmail_default"> T IFTRAN</div>
<div class="gmail_default"> T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)</div>
<div class="gmail_default"> F F T T T F F F F F F F IFTMSH (IF mesh for this field is T mesh)</div>
</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">Followed by fluid, T, PS1, and PS2 bcs.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">In userchk: </div>
<div class="gmail_default">
<div class="gmail_default" style="font-size:small"> parameter (lt=lx1*ly1*lz1*lelt)</div>
<div class="gmail_default" style="font-size:small"> common /myvort/ b1(lt,3),work1(lt),work2(lt)</div>
<div class="gmail_default" style="font-size:small"> common /vpot/ ucbx(lt),ucby(lt),ucbz(lt)</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default"> </div>
<div class="gmail_default"> m = nx1*ny1*nz1*nelv</div>
<div> call comp_vort3(b1,work1,work2,t(1,1,1,1,1), t(1,1,1,1,2),t(1,1,1,1,3))<br>
</div>
</div>
<div class="gmail_default">
<div class="gmail_default"> call vcross(ucbx,ucby,ucbz,vx,vy,vz,b1(1,1),b1(1,2),b1(1,3),m)<br>
</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default">In userq(ix,iy,iz,ieg)</div>
<div class="gmail_default">
<div class="gmail_default"> parameter (lt=lx1*ly1*lz1*lelt)</div>
<div class="gmail_default"> common /vpot/ ucbx(lt),ucby(lt),ucbz(lt)</div>
<div class="gmail_default"> integer e,f,eg<br>
</div>
<div class="gmail_default"> e = gllel(eg)</div>
<div class="gmail_default"> if (ifield.eq.2) then<br>
</div>
<div class="gmail_default"> qvol = ucbx(ix,iy,iz,ieg)</div>
<div class="gmail_default"> elseif (ifield.eq.3) then<br>
</div>
<div class="gmail_default"> qvol = ucby(ix,iy,iz,ieg)</div>
<div class="gmail_default"> elseif (ifield.eq.4) then<br>
</div>
<div class="gmail_default"> qvol = ucbz(ix,iy,iz,ieg)</div>
<div class="gmail_default"> endif<br>
</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default">When I compile i get and error<br>
</div>
<div class="gmail_default">
<div class="gmail_default"> qvol = ucbx(ix,iy,iz,ieg) </div>
<div class="gmail_default"> 1</div>
<div class="gmail_default">Error: Rank mismatch in array reference at (1) (4/1)</div>
<div class="gmail_default">/home/sandeep_reddy/work/Dynamo/GP/vectpot2/conj_ht.f:187.20:</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default"> qvol = ucby(ix,iy,iz,ieg) </div>
<div class="gmail_default"> 1</div>
<div class="gmail_default">Error: Rank mismatch in array reference at (1) (4/1)</div>
<div class="gmail_default">/home/sandeep_reddy/work/Dynamo/GP/vectpot2/conj_ht.f:190.20:</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default"> qvol = ucbz(ix,iy,iz,ieg) </div>
<div class="gmail_default"> 1</div>
<div class="gmail_default">Error: Rank mismatch in array reference at (1) (4/1)</div>
<div><br>
</div>
<div>Could someone kindly help me with this?</div>
<div>Thanking you.</div>
<div><br>
</div>
<div>Your sincerely,</div>
<div>Sandeep</div>
</div>
</div>
<div style="font-size:small"><br>
</div>
</div>
<div style="font-size:small"><br>
</div>
</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>