<div dir="ltr"><div><div>Hi Illias,<br><br></div>Thanks for the answer and clarification.<br><br></div>Sincerely,<br>Tanmoy<br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>Best Regards,<br></div>Tanmoy<br></div></div></div>
<br><div class="gmail_quote">On Mon, Jan 26, 2015 at 7:28 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">Hi Tanmoy,
thanks for the answer</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<br>
</p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">int_tp
is working perfectly for me now, the only issue is that it works only
for ASCII input files.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">For
my 3D simulation each text file occupies 100mb and I need more than
1000 time slices for analysis,</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">so
it becomes a real issue. </font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<br>
</p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">1.
I have not understood the purpose of you script,</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">to
make file.list I just execute a command like</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">ls
-1 *0.f* > file.list</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">What
was annoying for me is the manual confirmations after each file in
list.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">I
got rid of that by moving the part of int_tp.f file:</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
write(6,*) 'Input File Format: 0 - ASCII,
1 - binary(fld),</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
$ 2 - binary(0.f0000?) ? '</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
read (5,*) ib</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
ifbyte=.false.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
iffbin=.false.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
if(ib.ne.0) ifbyte = .true.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">---
if(ib.eq.2) iffbin = .true.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">out
of the cycle.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">But
what I can't figure out is how to make it working properly with
binary input format</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"> </font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">2.
If I am correct ifreguo = .true. allows you write on a
uniform mesh inside each element</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">instead
of GLL points. But I have a nonuniform mesh in terms of elements --</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">the
density of elements is greater near the boundaries.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">So
uniform mesh inside elements will not give me global uniform mesh.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<br>
</p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">3.
For Fourier analysis I use either fftw (<a href="http://www.fftw.org" target="_blank">www.fftw.org</a>)</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">or
routines from matlab, which are based on fftw.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">Really
i need to filter spacial waves on different time frequency, and for
that I need both time and space analysis (Hilbert transform),</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">which
can be also done with direct and inverse fft.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">From
results of int_tp I can construct a vector or square matrix (in 2D
case) in matlab and apply fft to it.</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<br>
</p>
<p style="margin-bottom:0in;line-height:100%">
<font color="#222222"><font face="arial, sans-serif"><font style="font-size:12pt">Best
regards, Ilias</font></font></font></p>
<p style="margin-bottom:0in;line-height:100%">
<br>
</p>
<p style="margin-bottom:0in;line-height:100%"><br>
</p></div><div><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi Ilias to circumvent, writing too mainy files in the "file.list"</div><div><div class="h5"><br>*--------------------------------------------------------------------------------*<br> program lst1<br> implicit none <br><br> integer, parameter :: N = 140<br> character*80,dimension(N) :: filename<br> integer :: i,j,k<br><br> <br> open(unit = 10,file='file.list',action='write')<br> write(10,"(I3)") N <br> do i = 1,N<br> <br> if(i .lt. 100)then<br> write(filename(i),'(A,I5.5)') 'turbChannel0.f',i<br> else<br> write(filename(i),'(A,I5.5)') 'turbChannel0.f',i<br> endif<br> enddo<br><br> do i = 1,N<br> write(10,*) trim(filename(i))<br> enddo<br><br><br> close(10)<br> stop<br> end program lst1<br>*--------------------------------------------------------------------------------*<br></div></div></div><br></div><div><div class="h5">For various post-processing format I use this script to write multiple files in file-list.<br><br><br></div></div></div><div><div class="h5"><div> 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.<br><br></div><div>Sincerely,<br>Tanmoy<br></div></div></div></div><div><div class="h5"><div class="gmail_extra"><br clear="all"><div><div><div dir="ltr"><div>Best Regards,<br></div>Tanmoy<br></div></div></div>
<br><div class="gmail_quote">On Sat, Jan 24, 2015 at 11:08 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>For time-space Fourier analysis I need data on a uniform mesh.</div><div>So I have used int_tp.</div><div><br></div><div>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.</div><div><br></div><div>So I tried to make uniform interpolation to the most simple tutorial examples, namely Rayleigh-Bendard ray_nn case.</div><div><br></div><div>Again if I apply int_tp to binary *0.f* of *fld*, then the result is as if the some space patches where switched,</div><div>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.</div><div><br></div><div>Is it a common problem, or I miss something? (input files in text format work for me but occupie too much space.)</div><div><br></div><div>Also, for each file in file.list int_tp asks for input file format.</div><div>But to make time-space analysis I need more that 1000 files in file.list, and</div><div>it is tiresome to type an answer that number of times.</div><div>It could be better to choose input file format just once at the beginning?</div><div><br></div><div>Thank you in advance, Ilias</div><div><br></div><div><br></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></div></div></blockquote></div><br></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>