<- [[:programming]]

====== base flow and fluctuations? ======

//Emmanuel Germaine, 2009-02-10 07:47 //

I was working with the channel flow code V2.1.2. Last week, I have installed the latest version available on the wiki. 1.3.3.
It seems that this version does not generate the same results than the previous that I used.
If  I save the velocity field using :

  FlowField u(Nx,Ny,Nz,3,Lx,Lz,a,b);
  ...
   for(int nz=0;nz<u.Nz();++nz){
     for(int ny=0;ny<u.Ny();++ny){
       for(int nx=0;nx<u.Nx();++nx){
         fileName <<  x(nx) << ' ' << y(ny) << ' ' << z(nz) <<' ' << ' ';
         fileName <<  u(nx,ny,nz,0) << ' ' << u(nx,ny,nz,1) << ' '<< u(nx,ny,nz,2) <<  endl ;
       }
     }
   }

The version v1.2.2 saved an instantaneous field. The new version saves a velocity field that looks centered on zero (there is no mean). I would like to work with the instantaneous velocity field instead of its fluctuating part. Do you know if there is a way to do this ?

thanks!

~~DISCUSSION~~