#include "global_var.h" 

void main()
{    
	 
	 filehandle8=fopen("resul.dat","w");
	 default_setting();
     user_could_reset_some_setting_here();
	 generate_grid();
     user_initial_grid_value();
	 for (iter=0;iter<last;iter++)  
	 {  
		     user_correct_density();
             iter1=0;
			 do{ 
				 user_bound_condition_setup();
				 setup2();
				 printf("Time step=%d, iter=%d, smax=%9.4f\n",iter,iter1,smax);
				 iter1++;
       		   } while (iter1<10);
   		     user_correct_density();
			 user_output();
             time=time+dt;
     
	 }
 
	fclose(filehandle8);  
}