#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,alpha_ph;
	double ra = 0.011147;
	float dv,wv,bigdlt; 
	char x0_fn[80],x_fn[180],out_fn[80]; 
	FILE *x0_ptr,*x_ptr,*out_ptr;		

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

        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(x_ptr,"%d  %d  %d  %f\n",&yv,&mv,&rv,&bigdlt);
		    
		  wv=(float)(dv*bigdlt);    
		 }
		 else  wv=0;
 
		  fprintf(out_ptr,"%9.6f ",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,"-18.0\n");

	fclose(x0_ptr);
   	fclose(out_ptr);

return 1;
}