[Nek5000-users] uniform interpolation

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

Hi Tanmoy, thanks for the answer


 int_tp is working perfectly for me now, the only issue is that it works
only for ASCII input files.

For my 3D simulation each text file occupies 100mb and I need more than
1000 time slices for analysis,

so it becomes a real issue.


 1.  I have not understood the purpose of you script,

to make file.list I just execute a command like

ls -1  *0.f* > file.list

What was annoying for me is the manual confirmations after each file in
list.

I got rid of that by moving the part of int_tp.f file:

---         write(6,*) 'Input File Format: 0 - ASCII, 1 - binary(fld),

---     $   2 - binary(0.f0000?) ? '

---         read (5,*) ib

---         ifbyte=.false.

---         iffbin=.false.

---         if(ib.ne.0) ifbyte = .true.

---         if(ib.eq.2) iffbin = .true.

out  of the cycle.

But what I can't figure out is how to make it working  properly with binary
input format



2. If I am correct ifreguo = .true.  allows you write on a uniform mesh
inside each element

instead of GLL points. But I have a nonuniform mesh in terms of elements --

the density of elements is greater near the boundaries.

So uniform mesh inside elements will not give me global uniform mesh.


 3. For Fourier analysis I use either fftw (www.fftw.org)

or routines from matlab, which are based on fftw.

Really i need to filter spacial waves on different time frequency, and for
that I need both time and space analysis (Hilbert transform),

which can be also done with direct and inverse fft.

>From results of int_tp I can construct a vector or square matrix (in 2D
case) in matlab and apply fft to it.


 Best regards, Ilias





Hi Ilias to circumvent, writing too mainy files in the "file.list"
>
>
> *--------------------------------------------------------------------------------*
>  program lst1
>  implicit none
>
>  integer, parameter :: N = 140
>  character*80,dimension(N) :: filename
>  integer :: i,j,k
>
>
>    open(unit = 10,file='file.list',action='write')
>    write(10,"(I3)") N
>    do i = 1,N
>
>    if(i .lt. 100)then
>    write(filename(i),'(A,I5.5)') 'turbChannel0.f',i
>    else
>    write(filename(i),'(A,I5.5)') 'turbChannel0.f',i
>    endif
>    enddo
>
>    do i = 1,N
>      write(10,*) trim(filename(i))
>    enddo
>
>
>    close(10)
>    stop
>    end program lst1
>
> *--------------------------------------------------------------------------------*
>
> For various post-processing format I use this script to write multiple
> files in file-list.
>
>
>  I was previously using ifreguo = .true. in post-processing format of nek
> for uniform interpolation. But, how are you performing FFT after
> interpolating to uniform mesh. Are you using a library, could you
> elucidate/ describe a little more on the process, on how you do that? I was
> trying to do that as well.
>
> Sincerely,
> Tanmoy
>
> Best Regards,
> Tanmoy
>
> On Sat, Jan 24, 2015 at 11:08 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>>
>> For time-space Fourier analysis I need data on a uniform mesh.
>> So I have used int_tp.
>>
>> When I applied int_tp to binary *0.f*  files (param(66) = 6) I have got a
>> strange pictures which reflect the solution, but as if some elements were
>> switched before interpolation.
>>
>> So I tried to make uniform interpolation to the most simple tutorial
>> examples, namely Rayleigh-Bendard ray_nn case.
>>
>> Again if I apply int_tp to binary *0.f* of *fld*, then the result is as
>> if the some space patches where switched,
>> but if int_tp is applied ot ASCII fld files  (param(66) = 6) then int_tp
>> works and saves the result correctly to both text and binary format.
>>
>> Is it a common problem, or I miss something? (input files in text format
>> work for me but occupie too much space.)
>>
>> Also, for each file in file.list int_tp asks for input file format.
>> But to make time-space analysis I need more that 1000 files in file.list,
>> and
>> it is tiresome to type an answer that number of times.
>> It could be better to choose input file format just once at the beginning?
>>
>> Thank you in advance, Ilias
>>
>>
>>
>> _______________________________________________
>> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20150126/aafa1790/attachment.html>