Hi Ali,
Seems like it has read nfiles=1 from the header which explains your behaviour.
So either this is set incorrectly or other information in the header is somehow
corrupted.
If you setenv SPLASH_DEBUG yes then you should get some debugging output
indicating this:
...
DEBUG: ntoti = xxx ntotall = xxx nfiles = xxx
...
this is from reading the first line in the file:
read(iunit,iostat=ierr) npartoftypei(1:6),massoftypei(1:6),timetemp,ztemp, &
iFlagSfr,iFlagFeedback,Nall(1:6),iFlagCool,nfiles
So I don't think it's related to the SFR not being dumped.
I would be happy to investigate further, if you can post the files somewhere
useful (e.g. dropbox, on web).
Daniel
On 20 May 2014, at 10:20 pm, Ali Snedden <[email protected]> wrote:
> Also, I have setenv GSPLASH_IGNORE_IFLAGCOOL 'YES' enabled. Thanks!
>
> ~ali
>
>
> On Tue, May 20, 2014 at 8:18 AM, Ali Snedden <[email protected]> wrote:
> Hello,
>
> Below is the Splash output. Yes I do get an error reading Column 11 (I think
> b/c I didn't write the SFR. Thank you for your time.
>
> Best Regards,
> ~ali
>
>
> [asnedden@crcfe01 100Mpc]$ gsplash lcdm/snapshot_000
> _ _
> (_) _ _ _ _ (_)_
> _ (_) ___ _ __ | | __ _ ___| |__ (_) _ (_)
> _ (_) _ / __| '_ \| |/ _` / __| '_ \ _ (_)
> (_) _ (_) \__ \ |_) | | (_| \__ \ | | | _ (_) _
> (_) _ |___/ .__/|_|\__,_|___/_| |_| (_) _ (_)
> (_) (_)|_| (_) (_) (_)(_) (_)(_) (_)(_)
>
> ( B | y ) ( D | a | n | i | e | l ) ( P | r | i | c | e )
>
> ( v2.2.2 [10th May 2013] Copyright (C) 2005-2013 )
>
> * SPLASH comes with ABSOLUTELY NO WARRANTY.
> This is free software; and you are welcome to redistribute it
> under certain conditions (see LICENSE file for details). *
>
> Comments, bugs, suggestions and queries to: [email protected]
> Check for updates at: http://users.monash.edu.au/~dprice/splash
> Please cite Price (2007), PASA, 24, 159-173 (arXiv:0709.0832) if you
> use SPLASH for scientific work and if you plot something beautiful,
> why not send me a copy for the gallery?
>
> splash.defaults: file not found: using program settings
> Using default M4 cubic spline (2h) kernel
>
> ----------------------- lcdm/snapshot_000.0 -----------------------
> >> reading default Gadget format <<
> star formation flag on: assuming star formation rate dumped
> time : 9.09090896686154948E-002
> z (redshift) : 10.00 (set GSPLASH_USE_Z=yes to use in legend)
> Npart (by type) : 16792216 16790451 0 0
> 389 0
> Mass (by type) : 0.0000000000000000 4.72310291121784340E-002
> 0.0000000000000000 0.0000000000000000 0.0000000000000000
> 0.0000000000000000
> N_gas : 16792216
> N_total : 33583056
> N data columns : 11
> > allocating memory: parts = 33583056 steps = 1 cols = 11
> positions 33583056
> velocities 33583056
> particle masses 16792605
> gas properties 16792216
> ERROR READING PARTICLE DATA from column 11
> converting GADGET smoothing length on gas particles to usual SPH definition
> (x 0.5)
>
> *** ERRORS detected during data read: data will be corrupted
> Please REPORT this and/or fix your file ***
>
> (set GSPLASH_IGNORE_ERRORS=yes to skip this message)
> > Press any key to bravely proceed anyway <
>
>
> Simulations employ 48.6 neighbours,
> corresponding to h = 1.13*(m/rho)^(1/3) in 3D
>
> >> last step ntot = 33583056
> time for data read = 21.13 s
>
> splash.units not found
>
> setting plot limits: steps 1-> 1 cols 1-> 11
> plot limits set
>
> splash.limits not found
>
> You may choose from a delectable sample of plots
> -------------------------------------------------------
> 1) x 7) particle mass
> 2) y 8) u
> 3) z 9) density
> 4) v\dx 10) h
> 5) v\dy 11) column 11
> 6) v\dz
> -------------------------------------------------------
> 12) multiplot [ 4 ] m) set multiplot
> -------------------------------------------------------
> d(ata) p(age) o(pts) l(imits) le(g)end h(elp)
> r(ender) v(ector) x(sec/rotate) s,S(ave) q(uit)
> -------------------------------------------------------
> Please enter your selection now (y axis or option):2
> (x axis) (default=1): 1
> (render) (0=none) ([0:11], default=0): 9
> (vector plot) (0=none, 4=v) ([0:4], default=0):
> Graphics device/type (? to see list, default /xw):
> using colour scheme red
> ----- t = 0.09, dump # 1 ------------------
> auto-selecting device resolution = 537 x 537
> projecting from particles to pixels...
> (use accelerated rendering option to make this faster)
> ( 10% - 1679222 particles done)
> ( 20% - 3358444 particles done)
> ( 30% - 5037665 particles done)
> ( 40% - 6716887 particles done)
> ( 50% - 8396108 particles done)
> ( 60% - 10075330 particles done)
> ( 70% - 11754552 particles done)
> ( 80% - 13433773 particles done)
> ( 90% - 15112995 particles done)
> (100% - 16792216 particles done)
> completed in 8.35 s
> adapting render limits
> x min, max = 3.05878669E-02 99999.906
> y min, max = 0.0000000 99999.992
> column density min, max = 0.0000000 4.04501846E-03
> rendering... 537 x 538 = 288906 pixels
>
> -------------- interactive mode commands -----------------------------
> SPACE BAR (or n) : skip to next timestep/file
> 0,1,2,3..9 and click : go forward/back n timesteps (back=r.click)
> left click (or A) : zoom/select
> right click (or X or b) : previous timestep
> shift+left click : IRREGULAR particle selection
> left click on colour bar : change rendering limits
> +/- : zoom IN/OUT (_ for out by 20%)
> a : (a)dapt plot limits (inside box, over axes or colour bar)
> l : (l)og / unlog axis (over x/y axis or colour bar)
> o/C : re-centre plot on (o)rigin/(C)ursor position
> backspace: delete annotation (over axes, legend, title or colour bar)
> r : (r)eplot current plot
> R : (R)eset/remove all range restrictions
> p/c : label closest (p)article/plot (c)ircle of interaction
> t : t)rack closest particle/turn tracking off (coord plots only)
> g : plot a line and find its g)radient
> ctrl-t : add text annotation at current position
> G/T/H : move le(G)end, (T)itle or (H) vector legend to current position
> m/M/i : change colour map (m=next,M=previous,i=invert) (rendered plots
> only)
> v/V/w : decrease/increase/adapt arrow size on vector plots (Z for x10)
> k/K : decrease/increase opacity on opacity-rendered plots (Z for x10)
> e/E : use current frame/settings as end point to animation sequence
> , . < > : rotate about z axis by +(-) 15,30 degrees (coord plots only)
> [ ] { } : rotate about x axis by +/- 15,30 degrees (coord plots only)
> / \ ? | : rotate about y axis by +/- 15,30 degrees (coord plots only)
> x : take cross section (coord plots only)
> u/U/d/D : move cross section/perspective pos. up/down (towards/away from
> observer)
> s : (s)ave current settings for all steps
> q/Q/esc : (q)uit plotting
> z/Z(oom) : timestepping, zoom and limits-changing multiplied by 10
> ----------------------------------------------------------------------
> log10(xmin = 0): min set to 9.99999996E-13
> ----- t = 0.09, dump # 1 ------------------
> auto-selecting device resolution = 537 x 537
> projecting from particles to pixels...
> (use accelerated rendering option to make this faster)
> ( 10% - 1679222 particles done)
> ( 20% - 3358444 particles done)
> ( 30% - 5037665 particles done)
> ( 40% - 6716887 particles done)
> ( 50% - 8396108 particles done)
> ( 60% - 10075330 particles done)
> ( 70% - 11754552 particles done)
> ( 80% - 13433773 particles done)
> ( 90% - 15112995 particles done)
> (100% - 16792216 particles done)
> completed in 8.35 s
> x min, max = 3.05878669E-02 99999.906
> y min, max = 0.0000000 99999.992
> log column density min, max = -12.000000 -2.3930795
> rendering... 537 x 538 = 288906 pixels
> entering interactive mode...press h in plot window for help
>
>
>
>
>
> On Mon, May 19, 2014 at 10:45 PM, Daniel Price <[email protected]>
> wrote:
> Hi Ali,
>
> Do you get any error messages? What is the log output?
>
> Daniel
>
> On 20 May 2014, at 10:10 am, Ali Snedden <[email protected]> wrote:
>
> > Hello,
> >
> > I am trying to use splash on several Gadget2 simulations in format "1".
> > Star formation is enabled, but the star formation rate is not saved to the
> > file (to save disk space).
> >
> > I can plot single files correctly, albeit with the errors complaining about
> > being unable to read the 11th column. However when I try to plot a single
> > snapshot composed of multiple files (i.e. snapshot_000.0, snapshot_000.1
> > ... snapshot_000.11), it only displays snapshot_000.0 and none of the
> > others. It appears that in the file header num_files is written correctly.
> >
> > Question:
> > Is there a way around this? Or will I need to write something in the Star
> > formation column?
> >
> > Attached is an example of what is happening. It is snapshot_000.0 of 12
> > snapshots. I tried getting it to work by entering "gsplash snapshot_000".
> > On files with star formation written, this successfully plots the all the
> > files in the snapshot.
> >
> > Thank you so much for your time.
> >
> > Best Regards,
> > Ali
> > <snapshot_000.0 of 12.png>
>
>
>
--
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.
