[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