LOCAL void print_ses_tri_solution_data(
FILE *file,
Front *fr,
Wave *wv,
const char *tri_header,
const char *table_name,
size_t num_plots,
float (**plot_fns)(float*,Front*,POINTER,COMPONENT,Locstate),
SESAME_EOS *seos)
{
RECT_GRID *rgr;
float area;
static const char *FORMAT =
"\n stop_time = %-10g stop_step = %-10d\n";
if (file != stdout)
{
record_print_version(file);
print_title_for_sesame(file,seos);
}
rgr = fr->rect_grid;
(void) foutput(file);
(void) fprintf(file,"\t\t\tINITIAL DATA:\n\n\n");
fprint_rectangular_grid(file,rgr);
(void) fprintf(file,FORMAT,0.0,0);
area = (rgr->U[0] - rgr->L[0])*(rgr->U[1] - rgr->L[1]);
(void) fprintf(file,"\n\t\tComputational Area = %g\n",area);
(void) fprintf(file,"\n\t\tRemap Geometry: IDENTITY_REMAP\n");
(void) fprintf(file,"\n\t\tPrinting Interval: 1 mesh units\n\n\n\n\n");
(void) fprintf(file,"\n\n\n\n");
(void) foutput(file);
(void) fprintf(file,"%s INTERFACE\n",table_name);
(void) fprintf(file,"\n\n\n\n");
(void) foutput(file);
(void) fprintf(file,"\t\t\tFRONT DATA:\n");
(void) fprintf(file," \n\t\t\tFront Rectangular Grid:\n\n");
fprint_rectangular_grid(file,rgr);
(void) fprintf(file,"\n\t\t\tInterface Topological Grid:\n\n");
fprint_rectangular_grid(file,&topological_grid(fr->interf));
fprint_interface(file,fr->interf);
(void) fprintf(file,"\n\n\n\n");
(void) foutput(file);
(void) fprintf(file,"\t\t\tEND OF FRONT DATA:\n");
(void) foutput(file);
(void) fprintf(file,"\t\t\tSTATE DATA:\n");
(void) foutput(file);
(void) fprintf(file,"%s",tri_header);
(void) fprintf(file,"\n");
(void) fprintf(file,"#Point Source Data\n#0\n");
print_tri_soln(file,fr,wv,wave_tri_soln(wv),num_plots,plot_fns);
(void) fprintf(file,"END OF TRI_SOLN\n");
(void) fprintf(file,"\n\n\n\n");
(void) foutput(file);
(void) fprintf(file,"\t\t\tEND OF STATE DATA\n");
} /*end print_ses_tri_solution_data*/
EXPORT void print_rectangular_grid(
const RECT_GRID *grid)
{
fprint_rectangular_grid(stdout,grid);
} /*end print_rectantular_grid*/
