<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>