<div dir="ltr"><div dir="ltr">Hi Ryan,</div><div dir="ltr"><br></div><div>In the file MagneticResistivity_interface.F90, does updating MagneticResistivity_fullState(...) to also accept your xx, yy arguments have any effect?<br></div><div><br></div><div>Or, would it be appropriate to /not/ pass xx, yy to this particular call to MagneticResistivity(...) ? I haven't used the MagneticResistivity unit before, so don't know what the right solution is.</div><div dir="ltr"><div><br></div><div>call MagneticResistivity(xCenter(i),yCenter(j),U(:,i,j,k), magResist(i,j,k))</div><div><br></div><div>sends in U(:,i,j,k), an array of NUNK_VARS reals, but the interface code you copied seems to expect temp, dens, and xn (real, real, array of reals).</div><div>The interface also includes "MagneticResistivity_fullState", which does take solnVec(NUNK_VARS), so it looks like MagneticResistivity(...) is a polymorphic function.</div><div><br></div><div>Best,</div><div><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr">Aaron</div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 20, 2020 at 10:26 AM Ryan Farber <<a href="mailto:rjfarber@umich.edu" target="_blank">rjfarber@umich.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi FLASH-users,<div><br></div><div>I expect there is an embarrassingly simple solution to my problem but since I'm at my wit's end:</div><div><br></div><div>I successfully modified the following files to make the magnetic resistivity eta = eta(x) for a test problem. However, when I do eta = eta(x,y), compilation fails with the following error in hy_uhd_getFaceFlux.F90</div><div><br></div><div>Files modified:</div><div> MagneticResistivity_interface.F90</div><div> MagneticResistivity.F90</div><div> Diffuse_computeDt.F90</div><div> hy_uhd_unsplitUpdate.F90</div><div> hy_uhd_getFaceFlux.F90</div><div><br></div><div>Error message:</div><div>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">hy_uhd_getFaceFlux.F90(333): error #6285: There is no matching specific subroutine for this generic subroutine call. <span> </span>[MAGNETICRESISTIVITY]</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>call MagneticResistivity(xCenter(i),yCenter(j),U(:,i,j,k), &</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">----------------------^</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">compilation aborted for hy_uhd_getFaceFlux.F90 (code 1)</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">make: *** [hy_uhd_getFaceFlux.o] Error 1</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)">
</p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">make: *** Waiting for unfinished jobs....</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p></div><div>Excerpt from getFaceFlux:</div><div>
<p style="margin:0px;font:11px Menlo;color:rgb(99,226,255)"><span style="font-variant-ligatures:no-common-ligatures">#if defined(FLASH_USM_MHD)<span> </span></span></p><p style="margin:0px;font:11px Menlo"><span style="font-variant-ligatures:no-common-ligatures"><font color="#000000"> </font></span><span style="color:rgb(235,239,24);font-variant-ligatures:no-common-ligatures">if</span><span style="color:rgb(0,0,0);font-variant-ligatures:no-common-ligatures"> (hy_use</span><span style="color:rgb(0,0,0);font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">Mag</span><span style="color:rgb(0,0,0);font-variant-ligatures:no-common-ligatures">neticResistivity) </span><span style="color:rgb(235,239,24);font-variant-ligatures:no-common-ligatures">then</span></p><p style="margin:0px;font:11px Menlo;color:rgb(50,244,241)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)"><span> </span></span><span style="font-variant-ligatures:no-common-ligatures">!! Get magnetic viscosity</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span></span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(50,244,241)"><b>call</b></span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">Mag</span><span style="font-variant-ligatures:no-common-ligatures">neticResistivity(xCenter(i),yCenter(j),U(:,i,j,k), </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(254,215,216)">&</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>magResist(i,j,k))</span></p><p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span></span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(235,239,24)">endif</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(99,226,255)">
</p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(99,226,255)"><span style="font-variant-ligatures:no-common-ligatures">#endif</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(99,226,255)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p></div><div>Excerpt from MagneticResistivity_interface:</div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span></span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(99,226,255)">interface</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">Mag</span><span style="font-variant-ligatures:no-common-ligatures">neticResistivity</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);min-height:13px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span></span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(99,226,255)">subroutine</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">Mag</span><span style="font-variant-ligatures:no-common-ligatures">neticResistivity(xx,yy,temp,dens,xn,magResist)</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(128,255,167)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>real</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">, </span><span style="font-variant-ligatures:no-common-ligatures">intent</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">(</span><span style="font-variant-ligatures:no-common-ligatures">IN</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">) :: xx,yy</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(128,255,167)"><span> </span>real</span><span style="font-variant-ligatures:no-common-ligatures">, </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(128,255,167)">intent</span><span style="font-variant-ligatures:no-common-ligatures">(</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(128,255,167)">IN</span><span style="font-variant-ligatures:no-common-ligatures">)<span> </span>:: temp, dens</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(128,255,167)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>real</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">, </span><span style="font-variant-ligatures:no-common-ligatures">intent</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">(</span><span style="font-variant-ligatures:no-common-ligatures">IN</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">), </span><span style="font-variant-ligatures:no-common-ligatures">dimension</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">(NSPECIES)<span> </span>:: xn</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(128,255,167)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>real</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">, </span><span style="font-variant-ligatures:no-common-ligatures">intent</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">(</span><span style="font-variant-ligatures:no-common-ligatures">OUT</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)">) :: magResist</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span></span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(99,226,255)">end subroutine</span><span style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">Mag</span><span style="font-variant-ligatures:no-common-ligatures">neticResistivity</span></p></div><div><br></div><div>I'm flummoxed since all I did was add "yCenter(j)" to the call of MagneticResistivity (and accordingly modified its arguments in the subroutine and in the interface).</div><div><br></div><div>If you can point out what silly thing I'm missing, that would be great!</div><div><br></div><div>Best,<br clear="all"><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr">--------<div>Ryan Farber<br>University of Michigan</div><div>PhD Candidate</div></div></div></div></div></div></div></div></div>
</blockquote></div></div>