<div dir="ltr">Hello,<div><br></div><div>I turned on uservp through param(30) in rea file and found that it messes up my output from drag_calc(). Shown below are outputs with and w/o uservp [ Repo : 560 ] .</div><div>Is there any specific flag for calculating integral quantities with uservp ?</div>
<div><br></div><div>With Uservp :</div><div>-----------------</div><div><br></div><div><div> Calculating eddy visosity</div><div> ediff max :-> 1.41399999999999989E-003</div><div> 5 2.2500E-02 3.0084E+00 9.7050E-01 3.0084E+00 9.7050E-01 umx</div>
<div> 5 2.25000E-02 1.50000E+00 3.46825E-01 6.48743E+00 6.48743E+00 cdiv</div><div> 5 1377 0.000E+00 6.316E+00 0.000E+00 4.436E-01 6.487E+00 1.203E-01divmnmx</div><div> 5 2.2500000E-02 0.0000000E+00 0.0000000E+00 6.5247004E-17 drgxt*****</div>
<div> 5 2.2500000E-02 0.0000000E+00 0.0000000E+00 1.0521227E-20 drgyt*****</div><div> 5 2.2500000E-02 0.0000000E+00 0.0000000E+00 -1.0486441E-20 drgzt*****</div></div><div><br></div><div>Without Uservp :</div>
<div>---------------------</div><div><br></div><div><div><div> ediff max :-> 1.41399999999999989E-003</div><div> 5 2.2500E-02 3.0107E+00 9.7141E-01 3.0107E+00 9.7141E-01 umx</div><div> 5 2.25000E-02 1.50000E+00 3.46825E-01 6.48743E+00 6.48743E+00 cdiv</div>
<div> 5 1377 0.000E+00 6.316E+00 0.000E+00 4.436E-01 6.487E+00 1.203E-01divmnmx</div><div> 5 2.2500000E-02 1.3228757E+00 1.0841687E+00 2.3870700E-01 dragx 1</div><div> 5 2.2500000E-02 -1.2363700E-04 -1.2361744E-04 -1.9553962E-08 dragy 1</div>
<div> 5 2.2500000E-02 1.6762005E-15 1.5786015E-15 9.7598988E-17 dragz 1</div><div>check: 6.28249446E+00 3.43007686E-16 0.00000000E+00 5</div><div> 5 2.2500000E-02 1.3228757E+00 1.0841687E+00 2.3870700E-01 drgxt 1</div>
<div> 5 2.2500000E-02 -1.2363700E-04 -1.2361744E-04 -1.9553962E-08 drgyt 1</div><div> 5 2.2500000E-02 1.6762005E-15 1.5786015E-15 9.7598988E-17 drgzt 1</div></div></div><div><br></div><div>Please let me know if rea/usr file is needed. Thanks.</div>
<div><br></div><div>Regards</div><div>Shriram<br>
</div></div>