[Nek5000-users] History points

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Fri May 7 19:06:13 CDT 2010


Frank,
I have not submitted it yet. Yes your right, some changes are needed to run it in parallel.

Stefan


On May 7, 2010, at 11:03 PM, nek5000-users at lists.mcs.anl.gov wrote:

> 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
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users




More information about the Nek5000-users mailing list