I read below and have a question about <span style="font-family: monospace;"></span>macroscopic averages. If my system have a surface between two materials which have different lattice constant. How to set the awin?How to macroscopic it twice? <br><br>Thank you!<br><pre>1<br>filename<br>1.0<br>3000<br>3<br>2.59<br><br>The explanation is in the header of the file PP/average.f90:<br><br>  !      nfile        ! the number of 3D-FFT files<br>  ! for each file:<br>  !      filename     ! the name of the 3D-FFT file<br>  !      weight       ! the weight of the quantity in this file<br>  !      .<br>  !      .<br>  ! end<br>  !      npt          ! the number of points of the thick mesh<br>  !      idir         ! 1,2 or 3. It is the fixed index which defines<br>  !                   ! the planes of the planar average<br>  !      awin         ! the size of the window for macroscopic averages.<br>  !<br><br>With the example above you'll get a 3-column output file for the <br>quantity to average<br>V(u):<br>- column 1: coordinate u, defined over a mesh of npt=3000 points along <br>the direction<br>idir=3<br>- column 2: averaged quantity V<br>- column 3: macroscopically averaged quantity V', over a windows with size<br>awin=2.50 au<br><br>If you want to average a weighted sum of different quantities V1,...,VN<br>defined over the same FFT mesh, just set nfile to N and specify N file names<br>with the relative weight.<br><br>Regards,<br>Giovanni<br><br></pre><br>