[Nek5000-users] Uservp and drag_calc routine
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Sat Aug 14 18:49:10 CDT 2010
Hello,
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 ] .
Is there any specific flag for calculating integral quantities with uservp ?
With Uservp :
-----------------
Calculating eddy visosity
ediff max :-> 1.41399999999999989E-003
5 2.2500E-02 3.0084E+00 9.7050E-01 3.0084E+00 9.7050E-01 umx
5 2.25000E-02 1.50000E+00 3.46825E-01 6.48743E+00 6.48743E+00
cdiv
5 1377 0.000E+00 6.316E+00 0.000E+00 4.436E-01 6.487E+00
1.203E-01divmnmx
5 2.2500000E-02 0.0000000E+00 0.0000000E+00 6.5247004E-17
drgxt*****
5 2.2500000E-02 0.0000000E+00 0.0000000E+00 1.0521227E-20
drgyt*****
5 2.2500000E-02 0.0000000E+00 0.0000000E+00 -1.0486441E-20
drgzt*****
Without Uservp :
---------------------
ediff max :-> 1.41399999999999989E-003
5 2.2500E-02 3.0107E+00 9.7141E-01 3.0107E+00 9.7141E-01 umx
5 2.25000E-02 1.50000E+00 3.46825E-01 6.48743E+00 6.48743E+00
cdiv
5 1377 0.000E+00 6.316E+00 0.000E+00 4.436E-01 6.487E+00
1.203E-01divmnmx
5 2.2500000E-02 1.3228757E+00 1.0841687E+00 2.3870700E-01 dragx
1
5 2.2500000E-02 -1.2363700E-04 -1.2361744E-04 -1.9553962E-08 dragy
1
5 2.2500000E-02 1.6762005E-15 1.5786015E-15 9.7598988E-17 dragz
1
check: 6.28249446E+00 3.43007686E-16 0.00000000E+00 5
5 2.2500000E-02 1.3228757E+00 1.0841687E+00 2.3870700E-01 drgxt
1
5 2.2500000E-02 -1.2363700E-04 -1.2361744E-04 -1.9553962E-08 drgyt
1
5 2.2500000E-02 1.6762005E-15 1.5786015E-15 9.7598988E-17 drgzt
1
Please let me know if rea/usr file is needed. Thanks.
Regards
Shriram
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100814/e1ae5652/attachment.html>
More information about the Nek5000-users
mailing list