#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,alf2oc;

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

sprintf(out_fn,"../flux_xd/ocnc13_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/oceanco2_oa_flux50r60m_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/ocnc13_disequ_11r12m.txt");
	if ((x_ptr=fopen(x_fn, "rt")) == NULL)
	{
	printf("\n Unable to open file2,  exitting ...\n");
	exit(0);
	}

	Ra = 0.011147;
	alf2oc = 0.998;

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

        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)  wv=0;
		else 
		wv=(float)(dv*dis[m][r-39]);    
		fprintf(out_ptr,"%9.7f ",wv);
		sv=sv+wv;
		}	

		fscanf(x0_ptr,"%f\n",&dv);
		wv=(float)(dv*dis[m][r-39]);
                fprintf(out_ptr,"%9.7f\n",wv);
                sv=sv+wv; 
		
	}  //end of m
	printf("%d\t%f\n",y,sv);

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

return 1;
}