Hi Alessandro,

fixed this - seems to work from my end. Let me know if you have any further 
issues.

r5206

just one question - both "rho" and "density" appear in the file - what is the 
difference? (it may be important for calculating the rendering weights)

Daniel

> On 6 Dec 2014, at 12:27 am, aatrani <[email protected]> wrote:
> 
> Thank you Daniel,
> 
> Yes, every output of the SPH codes of amuse also contain the smoothing length 
> of each particle. It's labelled "h_smooth".
> 
> I have attached two sample hdf files generated by the SPH codes bundled into 
> amuse, Fi and Gadget2.
> 
> Kind regards,
> 
> Alessandro
> 
> On Thursday, 4 December 2014 00:04:57 UTC+1, [email protected] wrote:
> Hi Alessandro, 
> 
> Thanks - I fixed the seg fault: 
> r5202 
> 
> However, this occurred because the smoothing length has not actually been 
> read from the file. This means none of the rendering options in splash will 
> work. 
> 
> Do you know if the file contains the smoothing length? If so, please post a 
> sample file so I can make sure the amuse read detects the smoothing length 
> correctly. 
> 
> If not, I am thinking of writing a splash utility to compute SPH-style 
> smoothing lengths for general point-like data. Would this be useful? 
> 
> Daniel 
> 
> > On 4 Dec 2014, at 1:25 am, aatrani <[email protected]> wrote: 
> > 
> >  Hello everybody, 
> > 
> > I'm trying to use SPLASH for the first time. 
> > 
> > I want to read an AMUSE hdf5 particle set generated by GADGET2. 
> > 
> > However, when I try to read the particle set with amsplash-hdf5 I get this 
> > error: 
> > 
> > splash.defaults not found: using default settings 
> >  reading AMUSE HDF5 format 
> > ----------------------- gd_t00000.hdf ----------------------- 
> >  npart:      107660 ncolumns:          19 
> > > allocating memory: parts =     107660 steps =      1 cols =   19 
> >  converting GADGET smoothing length on gas particles to usual SPH 
> > definition (x 0.5) 
> > At line 232 of file ../src/read_data_amuse_hdf5.f90 
> > Fortran runtime error: Index '0' of dimension 2 of array 'dat' outside of 
> > expected range (1:19) 
> > 
> > I commented the following lines in read_data_amuse_hdf5.f90: 
> > 
> >   if (required(ih) .and. size(dat(1,:,:)).ge.ih .and. 
> > npartoftype(1,i).gt.0) then 
> >      print "(a)",' converting GADGET smoothing length on gas particles to 
> > usual SPH definition (x 0.5)' 
> >      dat(1:npartoftype(1,i),ih,i) = 0.5*dat(1:npartoftype(1,i),ih,i) 
> >   endif 
> > 
> > After recompiling, amsplash-hdf5 seems to read the hfd5 file correctly. 
> > 
> > Do you know what could be the issue? 
> > 
> > Kind regards, 
> > 
> > Alessandro 
> >   
> > 
> > -- 
> > You received this message because you are subscribed to the Google Groups 
> > "SPLASH users forum" group. 
> > To unsubscribe from this group and stop receiving emails from it, send an 
> > email to [email protected]. 
> > To post to this group, send email to [email protected]. 
> > Visit this group at http://groups.google.com/group/splash-users. 
> > For more options, visit https://groups.google.com/d/optout. 
> 
> 
> -- 
> You received this message because you are subscribed to the Google Groups 
> "SPLASH users forum" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to [email protected].
> To post to this group, send email to [email protected].
> Visit this group at http://groups.google.com/group/splash-users.
> For more options, visit https://groups.google.com/d/optout.
> <gadget2_out.hdf><fi_out.hdf>

-- 
You received this message because you are subscribed to the Google Groups 
"SPLASH users forum" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/splash-users.
For more options, visit https://groups.google.com/d/optout.

Reply via email to