[Nek5000-users] Uservp and drag_calc routine

[nek5000-users at lists.mcs.anl.gov]

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>