[Nek5000-users] History points
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri May 7 16:03:27 CDT 2010
Hello Stefan,
It seems that this (very useful!) tool has not made it into the
repository. I tested out the code snippet you provided and am passing
on my results. It seems perhaps that it does not work with multiple
processes.
Cheers,
Frank
call userchk
reading hpts.in
found 2 points
reading hpts.in
found 2 points
reading hpts.in
found 2 points
reading hpts.in
found 2 points
reading hpts.in
found 2 points
reading hpts.in
found 2 points
reading hpts.in
found 2 points
reading hpts.in
found 2 points
forrtl: File exists
forrtl: File exists
forrtl: File exists
forrtl: File exists
forrtl: File exists
forrtl: severe (10): cannot overwrite existing file, unit 50,
file /data/fmuldoo/engineering-marangoni-flows/sim2/sim1/sim1/sim7/AR-.66-RE-1800-PR-04.0-dt-5.0e-5-tim
e-.0005-grid-2-Intel/hpts.out
Image PC Routine Line
Source
libintlc.so.5 00002B60A8F4CB1D Unknown Unknown
Unknown
libintlc.so.5 00002B60A8F4B625 Unknown Unknown
Unknown
libifcore.so.5 00002B60A8808C65 Unknown Unknown
Unknown
libifcore.so.5 00002B60A8779DC5 Unknown Unknown
Unknown
libifcore.so.5 00002B60A877955F Unknown Unknown
Unknown
libifcore.so.5 00002B60A87956D7 Unknown Unknown
Unknown
nek5000 000000000040470B Unknown Unknown
Unknown
nek5000 0000000000406D75 Unknown Unknown
Unknown
nek5000 00000000004059B5 Unknown Unknown
Unknown
nek5000 00000000004036AC Unknown forrtl: severe
(10): cannot overwrite existing file, unit 50,
file /data/fmuldoo/engineering-marangoni-flows
/sim2/sim1/sim1/sim7/AR-.66-RE-1800-PR-04.0-dt-5.0e-5-time-.0005-grid-2-Intel/hpts.out
Image PC Routine Line
Source
libintlc.so.5 00002ABBBCCE9B1D Unknown Unknown
Unknown
libintlc.so.5 00002ABBBCCE8625 Unknown Unknown
Unknown
libifcore.so.5 00002ABBBC5A5C65 Unknown Unknown
Unknown
libifcore.so.5 00002ABBBC516DC5 Unknown Unknown
Unknown
libifcore.so.5 00002ABBBC51655F Unknown Unknown
Unknown
libifcore.so.5 00002ABBBC5326D7 Unknown Unknown
Unknown
nek5000 000000000040470B Unknown Unknown
Unknown
nek5000 0000000000406D75 Unknown Unknown
Unknown
nek5000 00000000004059B5 Unknown Unknown
Unknown
nek5000 00000000004036AC Unknown forrtl: severe
(10): cannot overwrite existing file, unit 50,
file /data/fmuldoo/engineering-marangoni-flows
/sim2/sim1/sim1/sim7/AR-.66-RE-1800-PR-04.0-dt-5.0e-5-time-.0005-grid-2-Intel/hpts.out
Image PC Routine Line
Source
libintlc.so.5 00002B33F93D8B1D Unknown Unknown
Unknown
libintlc.so.5 00002B33F93D7625 Unknown Unknown
Unknown
libifcore.so.5 00002B33F8C94C65 Unknown Unknown
Unknown
libifcore.so.5 00002B33F8C05DC5 Unknown Unknown
Unknown
libifcore.so.5 00002B33F8C0555F Unknown Unknown
Unknown
libifcore.so.5 00002B33F8C216D7 Unknown Unknown
Unknown
nek5000 000000000040470B Unknown Unknown
Unknown
nek5000 0000000000406D75 Unknown Unknown
Unknown
nek5000 00000000004059B5 Unknown Unknown
Unknown
nek5000 00000000004036AC Unknown Unknown
Unknown
libc.so.6 00000030B041D994 Unknown Unknown
Unknown
nek5000 00000000004035A9 Unknown Unknown
Unknown
On Thu, 2010-05-06 at 15:49 +0200, nek5000-users at lists.mcs.anl.gov
wrote:
> Hi Guys,
>
> yes the history point module needs a major overhaul. I don't like the way it works!
> In fact we could do a much better job using our interpolation tool.
>
> I just coded something up (without testing the code snipplet) very quicky but I think I should work out of the box.
>
> How does it work:
>
> - create a 'hpts.in' file
> line 1: number of points
> line 2: X,Y,Z coordinates of a point 1
> ....
> line n+1: X,Y,Z coordinates of point n
>
> - set lpart >= n in the SIZE file
>
> - call hpts() in userchk()
>
> The results will be written to 'hpts.out' .
>
> Enjoy
> -Stefan
>
>
> ----------------------------------------------------------------------------------------------------------------------
>
> subroutine hpts
> c
> c dump velocity, temperature and ps-scalars for a given list
> c of points (read from file hpts.in) into a file (hpts.out)
> c
> INCLUDE 'SIZE'
> INCLUDE 'TOTAL'
>
> parameter(nmax=lpart,nfldmax=ldim+ldimt)
> parameter(mi=4,mr=1+2*ldim+nfldmax)
> real rTL(mr,nmax)
> integer iTL(mi,nmax)
> common /itlcb/ iTL
> common /rtlcb/ rTL
>
> common /outtmp / wrk(lx1,ly1,lz1,lelt,nfldmax)
>
> character*80 fname
> integer icalld,npoints,icount
> save icalld,npoints,icount
> data icalld /0/
> data npoints /0/
> data icount /0/
>
> nxyz = nx1*ny1*nz1
> ntot = nxyz*nelt
>
> if(icalld.eq.0) then
> icalld = 1
> npoints = 0
>
> write(6,*) 'reading hpts.in'
> open(50,file='hpts.in',status='old')
> read(50,*) npoints
> write(6,*) 'found ', npoints, ' points'
> do i = 1,npoints
> read(50,*) (rTL(1+j,i),j=1,ndim)
> enddo
> close(50)
>
> bb_tol = 1e-8 ! bounding-box test tolerance
> call intpts_setup(bb_tol)
>
> open(50,file='hpts.out',status='new')
> write(50,'(A)') '# time vx vy [vz] T PS1 PS2 ...'
> endif
>
> nflds = nfield + ndim-1 ! number of fields you want to interpolate
>
> ! pack working array
> call copy(wrk(1,1,1,1,1),vx,ntot)
> call copy(wrk(1,1,1,1,2),vy,ntot)
> if(if3d) call copy(wrk(1,1,1,1,2),vz,ntot)
> do i = 1,nfield-1
> call copy(wrk(1,1,1,1,ndim+i),T(1,1,1,1,i),ntot)
> enddo
>
> ! interpolate
> call intpts(wrk,nflds,iTL,mi,rTL,mr,npoints,nmax)
>
> ! write interpolation results to file
> if(nid.eq.0) then
> do ip = 1,npoints
> write(50,'(1p20E15.7)') time,
> & (rTL(1+2*ndim+ifld,ip), ifld=1,nflds)
> enddo
> endif
>
> return
> end
>
>
>
> On May 6, 2010, at 3:24 PM, nek5000-users at lists.mcs.anl.gov wrote:
>
> > P.S.
> >
> > Yes, Frank, the exact coordinates of the history points are written only when the code rans.
> >
> > But if for history points you choose the collocation points on the boundary of a spectral element -- those that have I,J,K equal to either 1 or lx1 -- than you can find out their coordinates from .rea ASCII file after the line
> >
> > *** MESH DATA ***
> >
> > where coordinates of the element edges are written.
> >
> > Aleks
> >
> >
> >
> >
> >
> > ----- Original Message -----
> > From: nek5000-users at lists.mcs.anl.gov
> > To: nek5000-users at lists.mcs.anl.gov
> > Sent: Wednesday, May 5, 2010 4:45:11 PM GMT -06:00 US/Canada Central
> > Subject: Re: [Nek5000-users] History points
> >
> > Hi Aleks,
> >
> > My question is what is the logfile? Perhaps you mean the output that
> > NEK writes to the screen? If so, then one would know the location of
> > the point at which data is being dumped only after entering the i,j,k
> > and element # of the point in the *.rea file....
> >
> > Cheers,
> > Frank
> >
> > On Wed, 2010-05-05 at 23:37 +0200, nek5000-users at lists.mcs.anl.gov
> > wrote:
> >> Thanks Aleks.
> >> Concerning grep, it is a linux command that searches for the string
> >> History in the logfile.
> >> _______________________________________________
> >> Nek5000-users mailing list
> >> Nek5000-users at lists.mcs.anl.gov
> >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> > --
> > 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
> > _______________________________________________
> > Nek5000-users mailing list
> > Nek5000-users at lists.mcs.anl.gov
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
--
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
More information about the Nek5000-users
mailing list