#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<malloc.h>
#include <string.h>

int main()
{
	int r,y,m,rv,yv,mv;
	double sv,fv;

	float dv,wv; 
	char x0_fn[80],x_fn[180],out_fn[80]; 
	FILE *x0_ptr,*x_ptr,*out_ptr;		

//sprintf(out_fn,"../flux_xd/ffuel_flux_xd60m50.txt");
sprintf(out_fn,"../flux_xd/fireflux_xd60m50.txt");
    	if ((out_ptr=fopen(out_fn, "wt")) == NULL)
    	{
      		printf("\n Unable to open outp file,  exitting ...\n");
	  	exit(0);
    	}

	sprintf(x0_fn,"../fluxdata/co2fire_flux60x50.txt");
//	sprintf(x0_fn,"../fluxdata/co2ff_flux60x50.txt");
        if ((x0_ptr=fopen(x0_fn, "rt")) == NULL)
        {
        printf("\n Unable to open file0,  exitting program ...\n\n");
        exit(0);
        }

	//fv=-29.4+8;   // -29.4 permil Andres1990:28.4,Tans1981:27.4from Ciasi1995

	fv=-26.5+8;   //  -26.5 permil for fire

        for (y=2000;y<=2004;y++)
        {
        sv=0;

        for (m=1;m<=12;m++)
        {
		for (r=1;r<50;r++)
		{
		fscanf(x0_ptr,"%f ",&dv);	
//printf("%f\t%f\n",dv,fv);
		wv=(float)(dv*fv);    
		fprintf(out_ptr,"%f\t",wv);
		sv=sv+wv;
		}	

		fscanf(x0_ptr,"%f\n",&dv);
		wv=(float)(dv*fv);
                fprintf(out_ptr,"%f\n",wv);
                sv=sv+wv; 
		
	 }  //end of m
	printf("%d\t%f\n",y,sv);
	}
	fprintf(out_ptr,"%f\n",fv);

	fclose(x0_ptr);
   	fclose(out_ptr);

return 1;
}