PRO hw_2
; Homework 2

; Create a contour plot that displays the highest average
; temperature in Florida. There will be 3 plots. One for
; January, one for July, and the amount that July is warmer 
; than January.

; Collect the station data from at least these cities:
; Tallahassee, Jacksonville, Orlando, Key West, St. Petersburg,
; Miami, and Pensacola. Other Cities can be included.

; This goes on, but I don't have it all.

; initialize the arrays for lon lat and all the data

lon=FINDGEN(16)*1.0
lat=lon
data_jan=lon
data_jul=lon
data_dif=lon

;Load in all of the city Data. I used 16 cities in alphabetical order 
; 
; 
count=0
     ; Apalachicola
lon[count]=-85-1./60
lat[count]=29+44/60.
data_jan[count]=62.4
data_jul[count]=89.8
data_dif[count]=27.4

count=count+1
     ; Crestview
lon[count]=-86-31./60
lat[count]=30+47/60.
data_jan[count]=61.9
data_jul[count]=92.2
data_dif[count]=30.3



count=count+1
     ; Daytona Beach
lon[count]=-81-3./60
lat[count]=29+11/60.
data_jan[count]=69.7
data_jul[count]=91.
data_dif[count]=21.3



count=count+1
     ; Fort Myers
lon[count]=-81-52/60.
lat[count]=26+35/60.
data_jan[count]=75.3
data_jul[count]=91.7
data_dif[count]=16.4



count=count+1
     ; Gainesville
lon[count]=-82-17./60
lat[count]=29+42/60.
data_jan[count]=66.2
data_jul[count]=90.9
data_dif[count]=24.7



count=count+1
     ; Jacksonville
lon[count]=-81-42./60
lat[count]=30+30/60.
data_jan[count]=64.2
data_jul[count]=90.8
data_dif[count]=26.6



count=count+1
     ; Key West
lon[count]=-81-45./60
lat[count]=24+33/60.
data_jan[count]=75.3
data_jul[count]=90.5
data_dif[count]=18.8



count=count+1
     ; Melbourne
lon[count]=-80-39./60
lat[count]=28+6/60.
data_jan[count]=71.7
data_jul[count]=90.5
data_dif[count]=18.8



count=count+1
     ; Miami
lon[count]=-80-18./60
lat[count]=25+49/60.
data_jan[count]=76.5
data_jul[count]=90.9
data_dif[count]=14.4



count=count+1
     ; Orlando/Sanford
lon[count]=-81-16./60
lat[count]=28+48/60.
data_jan[count]=70.4
data_jul[count]=91.9
data_dif[count]=21.5



count=count+1
     ; Pensacola
lon[count]=-87-11./60
lat[count]=30+29/60.
data_jan[count]=61.2
data_jul[count]=90.2
data_dif[count]=29.5



count=count+1
     ; Saint Petersburg
lon[count]=-82-38./60
lat[count]=27+46/60.
data_jan[count]=69.3
data_jul[count]=90.2
data_dif[count]=20.9



count=count+1
     ; Tallahassee
lon[count]=-84-21./60
lat[count]=30+24/60.
data_jan[count]=63.8
data_jul[count]=92.0
data_dif[count]=28.2



count=count+1
     ; Tampa
lon[count]=-82-32./60
lat[count]=27+58/60.
data_jan[count]=70.1
data_jul[count]=89.7
data_dif[count]=19.6



count=count+1
     ; Vero Beach
lon[count]=-80-25./60
lat[count]=27+39/60.
data_jan[count]=73.3
data_jul[count]=90.4
data_dif[count]=17.1



count=count+1
     ; West Palm Beach
lon[count]=-80-6./60
lat[count]=26+41/60.
data_jan[count]=75.1
data_jul[count]=90.1
data_dif[count]=15.0

print, data_jan, data_jul

stn=['AAF', 'CEW', 'DAB', 'FMY', 'GNV', 'JAX', 'EYW', 'MLB', 'MIA', 'ORL', 'PNS',$
	 'SPG', 'TLH', 'TPA', 'VRB', 'PBI']
 ;    aaf=[29, 44, 85, 01, 62.4, 89.8, 27.4]	; Apalachicola
 ;    cew=[30, 47, 86, 31, 61.9, 92.2, 30.3]	; Crestview
 ;    dab=[29, 11, 81, 03, 69.7, 91.0, 21.3]	; Daytona Beach
 ;    fmy=[26, 35, 81, 52, 75.3, 91.7, 16.4]	; Fort Myers
 ;    gnv=[29, 42, 82, 17, 66.2, 90.9, 24.7]	; Gainesville
 ;    jax=[30, 30, 81, 42, 64.2, 90.8, 26.6]	; Jacksonville
 ;    eyw=[24, 33, 81, 45, 75.3, 89.4, 14.1]	; Key West
 ;    mlb=[28, 06, 80, 39, 71.7, 90.5, 18.8]	; Melbourne
 ;    mia=[25, 49, 80, 18, 76.5, 90.9, 14.4]	; Miami
 ;    orl=[28, 48, 81, 16, 70.4, 91.9, 21.5]	; Orlando
 ;    pns=[30, 29, 87, 11, 61.2, 90.7, 29.5]	; Pensacola
 ;    spg=[27, 46, 82, 38, 69.3, 90.2, 20.9]	; St. Petersburg
 ;    tlh=[30, 24, 84, 21, 63.8, 92.0, 28.2]		; Tallahassee
 ;    tpa=[27, 58, 82, 32, 70.1, 89.7, 19.6]	; Tampa
 ;    vrb=[27, 39, 80, 25, 73.3, 90.4, 17.1]	; Vero Beach
;     pbi=[26, 41, 80, 06, 75.1, 90.1, 15.0]	; West Palm Beach

     DEVICE, DECOMPOSED=0 ; I want to see something useful
     ; It is a numerical color arrangement for the color table assigned next.
     LOADCT, 39 ; Loads the Rainbow color table.

     Window, 1, Title='MEAN JANUARY HIGH TEMPERATURE FOR 16 FLORIDA CITIES(1971-2000)', $
	XSize=450, YSize=400, XPOS=0, YPOS=900

     !P.BACKGROUND=255 ; Forces the background to be white.
     !P.COLOR=0 ; Forces default draw color to black.

TRIANGULATE, lon, lat, TRIANGLES, SPHERE=s, FVALUE=data_jan, /DEGREES

gdata=TRIGRID(data_jan, SPHERE=s, [0.5,0.5], [MIN(lon), MIN(lat), $
	MAX(lon), MAX(lat)], /DEGREES, XGRID=glon, YGRID=glat)

MAP_SET, /ISOTROPIC, /MERCATOR, LIMIT=[MIN(glat), MIN(glon), $
	MAX(glat), MAX(glon)]

cint = (MAX(gdata) - MIN(gdata))/10.
levels = MIN(gdata) + cint*FINDGEN(10)
labels = 1 + INTARR(10)
CONTOUR, gdata, glon, glat, LEVELS = levels, C_LABELS = labels,/FILL,/OVERPLOT
CONTOUR, gdata, glon, glat, LEVELS = levels, C_LABELS = labels,/OVERPLOT

MAP_CONTINENTS, COLOR= 123
MAP_GRID, LATDEL=2, LONDEL=2, LATLAB=MIN(glon)+0.5, LONLAB=MIN(glat) $
	+0.5, COLOR= 123, GLINESTYLE=2
;
;
;
; This will allow me to plot the station ID's and a symbol for each
;
;
FOR i= 0, count DO BEGIN
     PLOTS, lon[i], lat[i], PSYM=2, SYMSIZE=1.5, COLOR=0
     XYOUTS, lon[i], lat[i], stn[i], COLOR = 0, CHARSIZE=1.5, $
	 CHARTHICK=1.5
ENDFOR
;
;
; I will now generate Jul highs
;
;
;

     Window, 2, Title='MEAN JULY HIGH TEMPERATURE FOR 16 FLORIDA CITIES(1971-2000)', $
	XSize=450, YSize=400, XPOS=500, YPOS=900

TRIANGULATE, lon, lat, TRIANGLES, SPHERE=s, FVALUE=data_jul, /DEGREES

gdata=TRIGRID(data_jul, SPHERE=s, [0.5,0.5], [MIN(lon), MIN(lat), $
	MAX(lon), MAX(lat)], /DEGREES, XGRID=glon, YGRID=glat)

MAP_SET, /ISOTROPIC, /MERCATOR, LIMIT=[MIN(glat), MIN(glon), $
	MAX(glat), MAX(glon)]

cint = (MAX(gdata) - MIN(gdata))/10.
levels = MIN(gdata) + cint*FINDGEN(10)
labels = 1 + INTARR(10)
CONTOUR, gdata, glon, glat, LEVELS = levels, C_LABELS = labels,/FILL,/OVERPLOT
CONTOUR, gdata, glon, glat, LEVELS = levels, C_LABELS = labels,/OVERPLOT

MAP_CONTINENTS, COLOR= 123
MAP_GRID, LATDEL=2, LONDEL=2, LATLAB=MIN(glon)+0.5, LONLAB=MIN(glat) $
	+0.5, COLOR= 123, GLINESTYLE=2
;
;
;
; This will allow me to plot the station ID's and a symbol for each
;
;
FOR i= 0, count DO BEGIN
     PLOTS, lon[i], lat[i], PSYM=2, SYMSIZE=1.5, COLOR=0
     XYOUTS, lon[i], lat[i], stn[i], COLOR = 0, CHARSIZE=1.5, $
	 CHARTHICK=1.5
ENDFOR
;
;
; I will now generate the differences
;
;
;
     Window, 3, Title='DIFFERENCES IN MEAN HIGH TEMPERATURE FOR 16 FLORIDA CITIES(1971-2000)', $
	XSize=450, YSize=400, XPOS=0, YPOS=000

     !P.BACKGROUND=255 ; Forces the background to be white.
     !P.COLOR=0 ; Forces default draw color to black.

TRIANGULATE, lon, lat, TRIANGLES, SPHERE=s, FVALUE=data_dif, /DEGREES

gdata=TRIGRID(data_dif, SPHERE=s, [0.5,0.5], [MIN(lon), MIN(lat), $
	MAX(lon), MAX(lat)], /DEGREES, XGRID=glon, YGRID=glat)

MAP_SET, /ISOTROPIC, /MERCATOR, LIMIT=[MIN(glat), MIN(glon), $
	MAX(glat), MAX(glon)]

cint = (MAX(gdata) - MIN(gdata))/10.
levels = MIN(gdata) + cint*FINDGEN(10)
labels = 1 + INTARR(10)
CONTOUR, gdata, glon, glat, LEVELS = levels, C_LABELS = labels,/FILL,/OVERPLOT
CONTOUR, gdata, glon, glat, LEVELS = levels, C_LABELS = labels,/OVERPLOT

MAP_CONTINENTS, COLOR= 123
MAP_GRID, LATDEL=2, LONDEL=2, LATLAB=MIN(glon)+0.5, LONLAB=MIN(glat) $
	+0.5, COLOR= 123, GLINESTYLE=2
;
;
;
; This will allow me to plot the station ID's and a symbol for each
;
;
FOR i= 0, count DO BEGIN
     PLOTS, lon[i], lat[i], PSYM=2, SYMSIZE=1.5, COLOR=0
     XYOUTS, lon[i], lat[i], stn[i], COLOR = 0, CHARSIZE=1.5, $
	 CHARTHICK=1.5
ENDFOR






; 
; 
; 
; 
; 
END