[Nek5000-users] uniform interpolation

Mon Jan 26 10:48:32 CST 2015

Hi Illias,

Thanks for the answer and clarification.

Sincerely,
Tanmoy

Best Regards,
Tanmoy

On Mon, Jan 26, 2015 at 7:28 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

>  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
>>
>>
>
> _______________________________________________
> 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/6406482e/attachment-0001.html>