#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,Ra,alf2ph;

	float dv,wv,fv,dis[40],rat; 
	char x0_fn[80],x_fn[180],x1_fn[180],out_fn[80]; 
	FILE *x0_ptr,*x_ptr,*x1_ptr,*out_ptr;		

sprintf(out_fn,"../flux_xd/bioc13_flux_xd60m50r_pt2.txt");
    	if ((out_ptr=fopen(out_fn, "wt")) == NULL)
    	{
      		printf("\n Unable to open outp file,  exitting ...\n");
	  	exit(0);
    	}
	
sprintf(x0_fn,"../fluxdata/beps_HR_50r60m_gt.txt");
        if ((x0_ptr=fopen(x0_fn, "rt")) == NULL)
        {
        printf("\n Unable to open file0,  exitting program ...\n\n");
        exit(0);
        }

	sprintf(x_fn,"../fluxdata/disequ_land39rgn_permil042.txt");
	if ((x_ptr=fopen(x_fn, "rt")) == NULL)
	{
	printf("\n Unable to open file2,  exitting ...\n");
	exit(0);
	}

sprintf(x1_fn,"../fluxdata/bigDlt_land_39reg60m.txt");
        if ((x1_ptr=fopen(x1_fn, "rt")) == NULL)
        {
        printf("\n Unable to open file3,  exitting ...\n");
        exit(0);
        }

//	for (m=1;m<=12;m++)
          for (r=1;r<40;r++)
          {
	   fscanf(x_ptr,"%d %f %f\n",&rv,&fv,&dv);
	   dis[r]=fv;   // using permil
	  }
	fclose(x_ptr);

	alf2ph = 0.981;
	Ra = 0.011147;

        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);	
		if (r<40) fscanf(x1_ptr,"%d %d %d %f\n",&yv,&mv,&rv,&rat);

		if (r>=40)  wv=0;
		else 
		{
		
		wv = dv*dis[r];    
		}
		fprintf(out_ptr,"%9.7f ",wv);
		sv=sv+wv;
		}	

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

	}
	fprintf(out_ptr,"0.42\n");
	fclose(x0_ptr);
   	fclose(out_ptr);

return 1;
}