Hi,
The box size is *not* set in the file make_glass.pro; the value given there for "BoxSize" is irrelevant and can be anything non-zero (it will be rescaled automatically by N-GenIC).
The box size is set in the parameterfile given to N-GenIC, in the system of units specified there.
Volker
On Aug 3, 2012, at 3:26 PM, Agarwal, Shankar wrote:
>
> Hi,
>
> While generating ICs with N-GenIC, the boxsize needs to be set in Mpc/h or Kpc/h ?
>
> I reproduce below, the make_glass.pro file.
>
> Thanks
> Shankar
>
> ---------------------------------------
> fout = "dummy_glass.dat"
> N = 16L
> Ntot = N * N * N
>
> npart=lonarr(6)
> massarr=dblarr(6)
> time=0.0D
> redshift=0.0D
> flag_sfr=0L
> flag_feedback=0L
> npartall=lonarr(6)
> flag_cooling= 0L
> num_files= 1L
> BoxSize = 0.0D
>
> bytesleft=120
> la=intarr(bytesleft/2)
>
> BoxSize = 100.0D
> npart(1) = Ntot
> npartall(1) = Ntot
>
> pos= fltarr(3, Ntot)
>
> for i=0L, N-1 do begin
> for j=0L, N-1 do begin
> for k=0L, N-1 do begin
> pos(0, (i*N+j)*N+k) = (i+0.0)/N * BoxSize
> pos(1, (i*N+j)*N+k) = (j+0.0)/N * BoxSize
> pos(2, (i*N+j)*N+k) = (k+0.0)/N * BoxSize
> endfor
> endfor
> endfor
>
> openw,1,fout,/f77_unformatted
> writeu,1, npart,massarr,time,redshift,flag_sfr,flag_feedback,npartall,flag_cooling,num_files,BoxSize,la
> writeu,1, pos
> close,1
>
> end
>
>
>
>
> -----------------------------------------------------------
>
> If you wish to unsubscribe from this mailing, send mail to
> minimalist_at_MPA-Garching.MPG.de with a subject of: unsubscribe gadget-list
> A web-archive of this mailing list is available here:
> http://www.mpa-garching.mpg.de/gadget/gadget-list
Received on 2012-08-03 15:47:05