Hi Diana,
I would try stepping through the file read in Gadget2 in a debugger (for
convenience: don't run parallel) (or spamming print statements), and
checking whether everything (especially the header) is read as it
should. Changing the line (read_ic.c)
printf("too many particles\n");
by
printf("too many particles Numpart, n_for_this_task, All.MaxPart: %d %d
%d\n", Numpart, n_for_this_task, All.MaxPart);
might already give you a better clue as to what is going wrong.
You could also try whether reading the file works using other programs
(e.g. splash, hyplot).
Cheers,
Sander
Diana Juncher wrote:
> Hi all,
>
> Up until now I have always used Absoft's fortran-compiler to create my
> IC-files for Gadget2 and everything has worked just fine.
> Unfortunately, I don't have access to that compiler anymore, so I am
> instead giving gfortran a go. This has resulted in some complications,
> though. I use the exact same source code as I always have, but now I
> get the following error:
>
> ============================================================
> This is Gadget, version `2.0'.
>
> Running on 4 processors.
>
> Allocated 30 MByte communication buffer per processor.
>
> Communication buffer has room for 413910 particles in gravity computation
> Communication buffer has room for 131072 particles in density computation
> Communication buffer has room for 103477 particles in hydro computation
> Communication buffer has room for 106274 particles in domain decomposition
>
> Hubble (internal units) = 3.24078e-18
> G (internal units) = 1
> UnitMass_in_g = 1
> UnitTime_in_s = 1
> UnitVelocity_in_cm_per_s = 1
> UnitDensity_in_cgs = 1
> UnitEnergy_in_cgs = 1
>
> too many particles
> task 1: endrun called with an error level of 1313
>
>
> too many particles
> task 2: endrun called with an error level of 1313
>
> too many particles
> task 3: endrun called with an error level of 1313
>
>
>
> reading file `om1.dat' on task=0 (contains 19000 particles.)
> distributing this file to tasks 0-3
> --------------------------------------------------------------------------
> MPI_ABORT was invoked on rank 3 in communicator MPI_COMM_WORLD
> with errorcode 1313.
>
> NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
> You may or may not see output from other processes, depending on
> exactly when Open MPI kills them.
> --------------------------------------------------------------------------
> Type 0 (gas): 10000 (tot= 0000010000) masstab=1e-05
> Type 1 (halo): 9000 (tot= 0000009000) masstab=0.0001
> Type 2 (disk): 0 (tot= 0000000000) masstab=0
> Type 3 (bulge): 0 (tot= 0000000000) masstab=0
> Type 4 (stars): 0 (tot= 0000000000) masstab=0
> Type 5 (bndry): 0 (tot= 0000000000) masstab=0
>
> too many particles
> task 0: endrun called with an error level of 1313
>
> --------------------------------------------------------------------------
> mpirun has exited due to process rank 0 with PID 1835 on
> node node486.steno.dcsc.ku.dk <http://node486.steno.dcsc.ku.dk>
> exiting without calling "finalize". This may
> have caused other processes in the application to be
> terminated by signals sent by mpirun (as reported here).
> --------------------------------------------------------------------------
> ============================================================
>
> I am guessing it is a compiler problem, since I haven't changed the
> source code nor re-installed Gadget-2.
>
> The gfortran configuration is:
> >gfortran -v
> Configured with: ../configure --prefix=/usr --mandir=/usr/share/man
> --infodir=/usr/share/info --enable-shared --enable-threads=posix
> --enable-checking=release --with-system-zlib --enable-__cxa_atexit
> --disable-libunwind-exceptions --enable-libgcj-multifile
> --enable-languages=c,c++,objc,obj-c++,java,fortran,ada
> --enable-java-awt=gtk --disable-dssi --enable-plugin
> --with-java-home=/usr/lib/jvm/java-1.4.2-gcj-1.4.2.0/jre
> --with-cpu=generic --host=x86_64-redhat-linux
> Thread model: posix
> gcc version 4.1.2 20080704 (Red Hat 4.1.2-48)
>
> When I compile I use the command:
> >gfortran -frecord-marker=4 program.f90
>
> Without the frecord-marker I get an "incorrect header format" error.
>
> Thank you for your time,
> Diana
>
>
>
>
--
--------------------------------------------------
| Sander Valcke |
| Sterrenkundig Observatorium, Universiteit Gent |
| Krijgslaan 281 S9, B-9000 Gent, Belgium |
| phone +32-9-264-4796 | fax +32-9-264-4989 |
| http://sourceforge.net/apps/wordpress/hyplot/ |
--------------------------------------------------
Lawful Dungeon Master -- and they're MY laws!
--------------------------------------------------
Disclaimer:
By sending an email to ANY of my addresses you are agreeing that:
1. I am by definition, "the intended recipient"
2. All information in the email is mine to do with as I see fit and
make such financial profit, political mileage, or good joke as
it lends itself to.
3. I may take the contents as representing the views of your company.
4. This overrides any disclaimer or statement of confidentiality
that may be included in your message.
Received on 2010-12-08 16:23:27