;program to make plots of radiation data
;****must read all data first with read_radiation.pro****


cd,'/psd3data/mshupe/mosaic/dataeval'

mons=['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec']
;Get the data subset (must build to be same length)
for m=10,12 do begin
  wht=where(mont eq m)
  wh3=where(mon3 eq m)
  wh4=where(mon4 eq m)
  wh5=where(mon5 eq m)
  crossplot,lwdt[wht],jtit[wht],lwd3[wh3],jti3[wh3],lwd4[wh4],jti4[wh4],lwd5[wh5],jti5[wh5],title='LWD, '+mons[m-1],unit='[W/m2]',titl1='TOW',titl2='A30',titl3='A40',titl4='A50',$
    /nomissing,outfile='LWD.'+mons[m-1]+'.compplot.png'
  crossplot,lwut[wht],jtit[wht],lwu3[wh3],jti3[wh3],lwu4[wh4],jti4[wh4],lwu5[wh5],jti5[wh5],title='LWU, '+mons[m-1],unit='[W/m2]',titl1='TOW',titl2='A30',titl3='A40',titl4='A50',$
    /nomissing,outfile='LWU.'+mons[m-1]+'.compplot.png'
  crossplot,swdt[wht],jtit[wht],swd3[wh3],jti3[wh3],swd4[wh4],jti4[wh4],swd5[wh5],jti5[wh5],title='SWD, '+mons[m-1],unit='[W/m2]',titl1='TOW',titl2='A30',titl3='A40',titl4='A50',$
    /nomissing,outfile='SWD.'+mons[m-1]+'.compplot.png'
  crossplot,swut[wht],jtit[wht],swu3[wh3],jti3[wh3],swu4[wh4],jti4[wh4],swu5[wh5],jti5[wh5],title='SWU, '+mons[m-1],unit='[W/m2]',titl1='TOW',titl2='A30',titl3='A40',titl4='A50',$
    /nomissing,outfile='SWU.'+mons[m-1]+'.compplot.png'
endfor

end