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

int main()
{
	int r,m,yr,r1,yr1,m1;
	float Rc1,Rc2;
	double alpha_ph,sv;
	double ra = 0.011147;
	float bigdlt;	
	char x_fn[80],in_fn[80]; 
	FILE *x_ptr,*in_ptr;		

sprintf(x_fn,"../factor4resp50reg_Tans93_b8o16.txt");
    if ((x_ptr=fopen(x_fn, "wt")) == NULL)
    {
      printf("\n Unable to open outp file,  exitting program ...\n\n");
	  exit(0);
    }

	/** open land  file **/
sprintf(in_fn,"../../flux2conc/fluxdata/bigDlt_land_39reg60m.txt");
	if ((in_ptr=fopen(in_fn, "r")) == NULL)
	{
	  printf("\n Unable to open inp file,  exitting program ...\n\n");
	  exit(0);
	}

	printf("processing... \n");
	
	 for (yr=2000;yr<=2004;yr++)
        {
 	  for (m=1;m<=12;m++)
	  {
		sv=0;
		for (r=1;r<=39;r++)
		{
		 fscanf(in_ptr,"%d\t%d\t%d\t%f\n",&yr1,&m1,&r1,&bigdlt);
		Rc1=bigdlt-8;
		fprintf(x_ptr,"%d\t%d\t%d\t%8.6f\n",yr,m,r,Rc1);
		sv=sv+Rc1;
		}
		
	//	for (r=1;r<=39;r++) 
	//	 fprintf(x_ptr,"%d\t%d\t%d\t%8.6f\n",yr,m,r,sv/39);

		Rc2 = -2-16;
		for (r=40;r<=50;r++)
		fprintf(x_ptr,"%d\t%d\t%d\t%8.6f\n",yr,m,r,Rc2);
		
	  }   // end of m

	}  // end of yr

	fclose(in_ptr);
    fclose(x_ptr);
return 1;

}