*==================================================================
OPTIONS LS=132 PS=32400;
*==================================================================
* SAS CODE TO PROCESS GRIDDED CO2 FOSSIL-FUEL EMISSION DATA,
*  GRID CELL LOCATIONS, COUNTRY NAMES AND COUNTRY IDENTIFICATION 
*  CODES  FOR THE YEARS 1950, 1960, 1970, 1980 AND 1990
* WHERE:
* GEIAID = 1000*JGRID+IGRID
* LAT(ITUDE) = (JGRID-91)+0.5 (DEGREES)
* LONG(ITUDE) = (IGRID-181)+0.5 (DEGREES)
* FF = CO2 EMISSION FROM FOSSIL FUELS 
*      (1000 METRIC TONS C / GRID CELL / YEAR)
* UNID = UNITED NATIONS COUNTRY CODE
* GNAME =  UNITED NATIONS COUNTRY NAME
* IGIS = NASA-GISS COUNTRY NAME
* IJCODE = NASA-GISS COUNTRY/PROVINCE-REGION CODE
* CIJNAME = NASA-GISS COUNTRY/PROVINCE-REGION NAME
*
*==================================================================
 * 1950;
 DATA NEW5;
 INFILE 'GRIDALL.50';
 INPUT @1 GEIAID $CHAR6. @9 LAT 6.1 @17 LONG 6.1 @25 FF5 12.6 
       @39 UNID 3. @44 GNAME $CHAR42. 
       @88 IGIS 6. @95 IJCODE 6. @102 CIJNAME $CHAR14.;
 RUN;
 ;
 PROC SORT DATA=NEW5;
 BY UNID;
 PROC MEANS NOPRINT;
 BY UNID;
 ID GNAME IGIS;
 VAR FF5;
 OUTPUT OUT=FINAL SUM=SFF;
 PROC PRINT;
 VAR SFF GNAME IGIS UNID;
 RUN;
 ;
 * 1960;
 DATA NEW6;
 INFILE 'GRIDALL.60';
 INPUT @1 GEIAID $CHAR6. @9 LAT 6.1 @17 LONG 6.1 @25 FF6 12.6
       @39 UNID 3. @44 GNAME $CHAR42. 
       @88 IGIS 6. @95 IJCODE 6. @102 CIJNAME $CHAR14.;
 RUN;
 PROC SORT DATA=NEW6;
 BY UNID;
 PROC MEANS NOPRINT;
 BY UNID;
 ID GNAME IGIS;
 VAR FF6;
 OUTPUT OUT=FINAL SUM=SFF;
 PROC PRINT;
 VAR SFF GNAME IGIS UNID;
 RUN;
 ;
 * 1970;
 DATA NEW7;
 INFILE 'GRIDALL.70';
 INPUT @1 GEIAID $CHAR6. @9 LAT 6.1 @17 LONG 6.1 @25 FF7 12.6
       @39 UNID 3. @44 GNAME $CHAR42. 
       @88 IGIS 6. @95 IJCODE 6. @102 CIJNAME $CHAR14.;
 RUN;
 PROC SORT DATA=NEW7;
 BY UNID;
 PROC MEANS NOPRINT;
 BY UNID;
 ID GNAME IGIS;
 VAR FF7;
 OUTPUT OUT=FINAL SUM=SFF;
 PROC PRINT;
 VAR SFF GNAME IGIS UNID;
 RUN;
 ;
 * 1980;
 DATA NEW8;
 INFILE 'GRIDALL.80';
 INPUT @1 GEIAID $CHAR6. @9 LAT 6.1 @17 LONG 6.1 @25 FF8 12.6
       @39 UNID 3. @44 GNAME $CHAR42. 
       @88 IGIS 6. @95 IJCODE 6. @102 CIJNAME $CHAR14.;
 RUN;
 PROC SORT DATA=NEW8;
 BY UNID;
 PROC MEANS NOPRINT;
 BY UNID;
 ID GNAME IGIS;
 VAR FF8;
 OUTPUT OUT=FINAL SUM=SFF;
 PROC PRINT;
 VAR SFF GNAME IGIS UNID;
 RUN;
 ;
 * 1990;
 DATA NEW9;
 INFILE 'GRIDALL.90';
 INPUT @1 GEIAID $CHAR6. @9 LAT 6.1 @17 LONG 6.1 @25 FF9 12.6
       @39 UNID 3. @44 GNAME $CHAR42. 
       @88 IGIS 6. @95 IJCODE 6. @102 CIJNAME $CHAR14.;
 RUN;
 PROC SORT DATA=NEW9;
 BY UNID;
 PROC MEANS NOPRINT;
 BY UNID;
 ID GNAME IGIS;
 VAR FF9;
 OUTPUT OUT=FINAL SUM=SFF;
 PROC PRINT;
 VAR SFF GNAME IGIS UNID;
 RUN;
 ;
*==================================================================