#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; 
	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_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/ocean_netflux50r60m_jon2gt.txt");
        if ((x0_ptr=fopen(x0_fn, "rt")) == NULL)
        {
        printf("\n Unable to open file0,  exitting program ...\n\n");
        exit(0);
        }

	alf2oc = -2;
	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);	
//		fscanf(x_ptr,"%d  %d  %d  %f\n",&yv,&mv,&rv,&fv);
//printf("%f\t%f\n",dv,fv);
		wv=(float)(dv*alf2oc);    
		fprintf(out_ptr,"%9.6f ",wv);
		sv=sv+wv;
		}	

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

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

return 1;
}