[Nek5000-users] CFL

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

Hi Frank,

Presently, the max cfl is written -- the /10 is because
of the p-format in the following format stmt, e.g.,

    1pe12.4

is like e-format,  0.1234e07, save that it shifts the
significant digits to the left to give 1.2345e06, which
provides one more digit in the same space.  (A trick I
learned from the linpack/lapack source.)

Oddly, the "1p" has the effect of shifting all reals that
appear to the right of the first instance of the p-format
to the left, even if they are printed in std. "f-format".
The "/10" unwinds this shift.

If you want the max in each direction you could make a
copy of the routine that computes it, rename it,
instrument it, add it to your .usr file and call it
from userchk --- this way you'll continue to be repo-
compatible but have full access to the information you
desire.   I usually call it blahx if the original was blah.

Paul



On Thu, 15 Apr 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Hello all,

Is it possible to have the maximum CFL # (even better the maximum in
each direction x,y,z) written out to the screen?  I see the line in
"drive2.f":
WRITE (6,100) ISTEP,TIME,DT,COURNO/10,TTIME,TTIME_STP

It seems this is the maximum CFL #, but I don't understand the "/10".

Cheers,
Frank

-- 
Frank Herbert Muldoon, Ph.D. Mechanical Engineering
Technische Universität Wien (Technical University of Vienna)
Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
Mechanics and Heat Transfer)
Resselgasse 3
1040 Wien
Tel: +4315880132232
Fax: +4315880132299 
Cell:+436765203470
fmuldoo (skype)
http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users