bin2simple.c

/* $Header$
 * $Log$
 * Revision 1.2  2003/06/24 23:40:41  trq
 * Added ascii2std program.
 *
 * Revision 1.1  2001/05/17 19:33:34  trq
 * Print out tipsy files in a "simple" format.
 *
 * Revision 1.1.1.1  1995/01/20  20:02:45  trq
 * Initial revision
 *
 */
/*
 */

#include <stdio.h>
#include <malloc.h>
#include <assert.h>
#include "tipsydefs.h"

int
main()
{
    struct gas_particle *gas_particles, *gp, *lastgp;
    struct dark_particle *dark_particles, *dp, *lastdp;
    struct star_particle *star_particles, *sp, *lastsp;
    struct dump header;
    
    while(1) 
      {
	if(fread((char *)&header,sizeof(header),1,stdin) == 0)
	  break;
	if(header.nsph != 0) {
	    gas_particles = (struct gas_particle *)
				malloc(header.nsph*sizeof(*gas_particles));
	    if(gas_particles == NULL) {
		printf("<sorry, no memory for gas particles, master>\n") ;
		return -1;
	    }
	}
	if(header.ndark != 0) {
	    dark_particles = (struct dark_particle *)
				malloc(header.ndark*sizeof(*dark_particles));
	    if(dark_particles == NULL) {
		printf("<sorry, no memory for dark particles, master>\n") ;
		return -1;
	    }
	}
	if(header.nstar != 0) {
	    star_particles = (struct star_particle *)
				malloc(header.nstar*sizeof(*star_particles));
	    if(star_particles == NULL) {
		printf("<sorry, no memory for star particles, master>\n") ;
		return -1;
	    }
	}
    
	printf("%d %d %d %d\n" ,header.nbodies, header.nsph,
		header.ndark, header.nstar) ;

	fread((char *)gas_particles,sizeof(struct gas_particle),
			 header.nsph,stdin) ;

	lastgp = gas_particles + header.nsph ;
	for(gp=gas_particles; gp < lastgp;  gp++){
	    fprintf(stdout,"%g %g %g %g %g %g %g %g %g\n",gp->mass,
		    gp->pos[0], gp->pos[1], gp->pos[2],
		    gp->vel[0], gp->vel[1], gp->vel[2], gp->rho, gp->temp);
	}
	fread((char *)dark_particles,sizeof(struct dark_particle),
			 header.ndark,stdin) ;

	lastdp = dark_particles + header.ndark ;
	for(dp=dark_particles; dp < lastdp;  dp++){
	    fprintf(stdout,"%g %g %g %g %g %g %g\n",dp->mass,
		    dp->pos[0], dp->pos[1], dp->pos[2],
		    dp->vel[0], dp->vel[1], dp->vel[2]);
	}
	fread((char *)star_particles,sizeof(struct star_particle),
			 header.nstar,stdin) ;

	lastsp = star_particles + header.nstar ;
	for(sp=star_particles; sp < lastsp;  sp++){
	    fprintf(stdout,"%g %g %g %g %g %g %g\n",sp->mass,
		    sp->pos[0], sp->pos[1], sp->pos[2],
		    sp->vel[0], sp->vel[1], sp->vel[2]);
	}
	fprintf(stderr,"read time %f\n",header.time) ;

	if(header.nsph != 0) {
	  free(gas_particles);
	}
	if(header.ndark != 0) {
	  free(dark_particles);
	}
	if(header.nstar != 0) {
	  free(star_particles);
	}
    }
    return(0);
}