Hi Nek's,<br><br>I have some troubles with the load_fld subroutine. It seems to load the files I want, however it sets everything (vx, vy, vz) to zero...<br>Here is what I have in my usrchk routine:<br><br> n = nx1*ny1*nz1*nelt<br>
ifto = .true.<br><br> if(nid.EQ.0) then<br> open(unit=199,file='file.list',form='formatted',status='old')<br> read(199,'(A80)') (filename(i),i=1,ncli+1)<br> close(199)<br>
endif<br><br> call bcast(filename,(ncli+1)*80)<br><br> do i = 1,ncli+1<br> call load_fld(filename(i))<br> if(nid.EQ.0) write(*,*) 'File', filename(i), 'loaded.'<br> call lambda2(t(1,1,1,1,1))<br>
call outpost(vx,vy,vz,pr,t,'PRT')<br> enddo<br><br>And here is what I get in my logfile for the first file for instance:<br><br> nsteps=0 -> skip time loop<br> running solver in post processing mode<br>
<br> call userchk<br> set initial conditions<br> nekuic (1) for ifld 1<br> call nekuic for vel<br> xyz min 0.0000 0.0000 0.0000<br> uvwpt min 0.99960E-20 0.99960E-20 0.99960E-20 0.0000 0.99960E-20<br>
PS min 0.0000 0.0000 0.99000E+22<br> xyz max 1.0000 1.0000 1.0000<br> uvwpt max 0.80000E-19 0.80000E-19 0.80000E-19 0.0000 0.80000E-19<br> PS max 0.0000 0.0000 -0.99000E+22<br>
done :: set initial conditions<br> <br> File prtcav0.f00001 loaded.<br><br>Have I been doing something wrong using the load_fld routine?<br><br>Regards,<br><br><br><div class="gmail_quote">On 12 March 2012 13:00, <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"><a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> a écrit :<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi everyone, i am trying to post process some data but it seems load_fld doesn't load the velocity field correctly.<br>
My binary files are default .f00* files (param(66)=param(67)=0)<br>
I want to recalculate the vortex size (aa in my code) in a 2D case.<br>
So here is what i want my userchk() to do<br>
Open one by one .f00* files,<br>
Calculate the vorticity field from velocity field<br>
Restrain the domain to avoid Boundary conditions<br>
Calculate the circulation<br>
Calculate The vorticity barycenter coordinates<br>
Calculate The vorticity Size (aa)<br>
Output to "Rayon_Circulation_domaine_<u></u>entier" file the quantity sqrt(aa)<br>
Does load_fld() work only with .fld files ?? Because with this i got only negative values for the vorticity size and the circulation value is not good<br>
B.regards<br>
Can<br>
<br>
<br>
<br>
My post-precessing part in my userchk() is like<br>
<br>
! read file-list<br>
if (nid.eq.0) then<br>
open(unit=199,file='file.list'<u></u>,form='formatted',status='old'<u></u>)<br>
read(199,*) nfiles<br>
read(199,'(A80)') (filename(j),j=1,nfiles)<br>
close(199)<br>
endif<br>
call bcast(nfiles,isize)<br>
call bcast(filename,nfiles*80) do j = 1,nfiles<br>
call load_fld(filename(j))<br>
call comp_vort3(vort,work1,work2,<u></u>vx,vy,vz)<br>
do i=1, ntot<br>
if((abs(xm1(i,1,1,1))<14.9) .and. (abs(ym1(i,1,1,1))<14.9))<br>
& then x_new(i)=xm1(i,1,1,1)<br>
y_new(i)=ym1(i,1,1,1)<br>
rId(i)=1.0<br>
<br>
else<br>
<br>
x_new(i)=0.0<br>
y_new(i)=0.0<br>
rId(i)=0.0<br>
<br>
endif<br>
<br>
enddo<br>
circ=glsc3(rId,bm1,vort(1,1),<u></u>ntot)<br>
x_c=glsc3(x_new,bm1,vort(1,1),<u></u>ntot)/circ ! Calculation of vorticity barrycenter<br>
y_c=glsc3(y_new,bm1,vort(1,1),<u></u>ntot)/circ<br>
<br>
do i=1, ntot<br>
if((abs(xm1(i,1,1,1))<14.9) .and. (abs(ym1(i,1,1,1))<14.9))<br>
& then rr(i)=(xm1(i,1,1,1)-x_c(1))**2 + (ym1(i,1,1,1)<br>
& -y_c(1))**2<br>
else<br>
rr(i)=0.0<br>
endif<br>
enddo<br>
<br>
aa=glsc3(rr,bm1,vort(1,1),<u></u>ntot)/circ ! Calculation of vorticity size<br>
print*,'aa=',aa<br>
if (nid.eq.0) then<br>
open(UNIT=1,FILE='Rayon_<u></u>Circulation_domaine_entier',<u></u>STATUS='unknown') !Output to 'Rayon_Circulation_domaine_<u></u>entier' file<br>
write(1,'(1p20E15.7)') time,sqrt(aa),circ<br>
endif<br>
enddo<br>
!we are done<br>
call exitt<br>
c ==============================<u></u>==============================<u></u>==============================<u></u>======================= <br>
<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote></div></div>
Well sorry for this, load_fld() works great with .f00* files<br>
Le problem was related to my grid in the domain which was not tight enough for x>10<br>
Thx anyways<div class="HOEnZb"><div class="h5"><br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Jean-Christophe<br>