ods _all_ close;
  ods listing;
OPTIONS MPRINT SYMBOLGEN MLOGIC NOCENTER nodate nonumber ;
ods pdf file="K:\sci\sci chesapeake bay program projects\CURRENT\ANNUAL_REPORT\GRAPHICS\Report\try1.pdf"  pdftoc=1;

LIBNAME SAS7BDAT "K:\sci\sci chesapeake bay program projects\CBP_DATABASE\CBPPHYTO\";
DATA ONE; SET SAS7BDAT.vapi1985_2014final;
IF LAYER="AP";
PROC SORT; BY STATION DATE;



%MACRO STAT(STAT,DIR,OUTSTAT);

%MACRO PLOTVAR (VAR,LABEL1,LABEL2,TITLE,OUTVAR,i);

GOPTIONS DEVICE=GIF DISPLAY NOPROMPT 
ROTATE=portrait GSFMODE=REPLACE  GSFNAME=&OUTSTAT&OUTVAR
HSIZE=5.07 IN VSIZE = 3.75
IN HORIGIN = 1 IN VORIGIN = 1 IN
HPOS=150 VPOS=70;
FILENAME &OUTSTAT&OUTVAR "K:\sci\sci chesapeake bay program projects\CURRENT\ANNUAL_REPORT\GRAPHICS\Report\&DIR\&OUTSTAT&OUTVAR..GIF";

DATA TWO; SET ONE;
YEAR=YEAR(DATE);
MONTH=MONTH(DATE);
DAY=DAY(DATE);

IF MONTH=2 THEN DD=DAY/28;

IF MONTH IN (1,3,5,7,8,10,12) THEN DO;
DD=DAY/31;
END;

IF MONTH IN (4,6,9,11) THEN DO;
DD=DAY/30;
END;

MM=((MONTH-1)+DD)/12;
PDATE=YEAR+MM;
IF STATION="&STAT";
PICO_C=PICO*1.5;

data two; set two;
IF &VAR=. THEN DELETE;
MERGEME=1;


PROC MEANS DATA=TWO NOPRINT;
VAR &VAR;
OUTPUT OUT=CHK STD=ST&VAR;

DATA CHK; SET CHK;
ST&VAR=3*(ST&VAR);
MERGEME=1;

PROC SORT; BY MERGEME;

DATA TWO; MERGE TWO CHK; BY MERGEME;

ods proclabel= "Figure &i..    Plot of above pycnocline &LABEL2 against time at station &TITLE for the period of 1985 through 2014." ;
PROC GPLOT;
AXIS1 VALUE=(H=1.5 FONT='Calibri' COLOR=BLACK)
         LABEL=(A=90 H=2 FONT='Calibri' COLOR=BLACK "&LABEL1")
         MAJOR=(H=2 COLOR=BLACK)
         MINOR=(NUMBER=9 H=1);
AXIS2 VALUE=(H=1.5 FONT='Calibri' COLOR=BLACK)
         LABEL=(H=2 FONT='Calibri' COLOR=BLACK 'Date')
         MAJOR=(H=2 COLOR=BLACK)
         MINOR=(H=1 COLOR=BLACK NUMBER=1)
         ORDER=(1989 TO 2015 by 2);
SYMBOL1 COLOR=BLACK V=STAR;
TITLE1 H=4 F='Calibri' C=BLACK "AP - &TITLE";

footnote1 justify=left C=BLACK  height=11pt font='calibri' bold wrap  "Figure" &i ".    Plot of above pycnocline &LABEL2 against time at station &TITLE for the period of 1985 through 2014." ;

PLOT &VAR*PDATE=1 / VAXIS=AXIS1  HAXIS=AXIS2 NOLEGEND;
RUN;

%MEND;

%PLOTVAR(TOTABUND,Total Abundance (#/L),Total Abundance,&STAT,01TAB,1);
%PLOTVAR(RATIO,Ratio of Biomass to Abundance,Ratio of Biomass to Abundance,&STAT,03RAT,2);
%PLOTVAR(TOTBMASS,Total Biomass (pg C/L),Total Biomass,&STAT,02TBM,3);
%PLOTVAR(DINOFL_C,Dinoflagellate Biomass (pg C/L),Dinoflagellate Biomass ,&STAT,06DIN,4);
%PLOTVAR(TOTABUND,Total Abundance (ind/L),Total Abundance,&STAT,01TAB,5);
%PLOTVAR(DIATOM_C,Diatom Biomass (pg C/L),Diatom Biomass,&STAT,05DIA,6);
%PLOTVAR(CYANO_C,Cyanobacteria Biomass (pg C/L),Cyanobacteria Biomass ,&STAT,07CYB,7);
%PLOTVAR(CHLORO_C,Chlorophyte Biomass (pg C/L),Chlorophyte Biomass,&STAT,08CHC,8);
%PLOTVAR(PICO_C,Picoplankton Biomass (ind/L),Picoplankton Biomass,&STAT,09PCA,9);
%PLOTVAR(CRYPTO_C,Cryptophyte Biomass (pg C/L),Cryptophyte Biomass,&STAT,10CRC,10);
%PLOTVAR(CYANO,Cyanobacteria Abundance (ind/L),Cyanobacteria Abundance,&STAT,11CYA,11);
%PLOTVAR(MARGALEF,Margalef Diversity Index,Margalef Diversity Index,&STAT,04MAR,12);




%MEND;

%STAT(SBE5,SBE5,SE5);

run;
quit;
ods pdf close;