PRO noahplots, basename1, basename2, nsoil, title=title, name=name

  IF n_elements(nsoil) EQ 0 THEN nsoil=8

  IF n_elements(basename1) gt 100 then begin
     data1=basename1
  endif else $
    j=load_cols(basename1, data1)
  
  IF n_elements(basename2) gt 100 then begin
     data2=basename2
  endif else $
    j=load_cols(basename2, data2)

  

  xr=[223,240]
  yr=[285,325]

  day=lindgen(n_elements(data1[0,*]))/48.

  if n_elements(day) lt 10000 then day=day+220

  i=2
  plot, day, data1[i,*], xr=xr, /xs, yr=yr, /ys, $
    line=2, title=title+' Skin Temperature', $
    ytitle='Skin Temperature (K)', xtitle='Day'
  oplot, day, data2[i,*], color=1
  oplot, day, data1[i,*], color=3
  
  i=8
  plot, day, data1[i,*], xr=xr, yr=yr, /xs, /ys, $
    line=2, title=title+' Temperature', $
    ytitle='Temperature (K)', xtitle='Day'
  oplot, day, data2[i,*], color=1
  oplot, day, data1[i,*], color=3

  i=8+nsoil
  plot, day, data1[i,*], xr=xr, yr=[0,0.5], /xs, /ys, $
    line=2, title=title+' Soil Moisture and Precip', $
    ytitle='Moisture Content (cm^3/cm^3)', xtitle='Day'
  oplot, day, data2[i,*], color=1
  oplot, day, data1[i,*], color=3
  oplot, day, data2[i+1,*], color=1, line=2
  oplot, day, data1[i+1,*], color=3, line=2
  oplot, day, data2[i+2,*], color=1, line=1
  oplot, day, data1[i+2,*], color=3, line=1
  oplot, day, data2[i+3,*], color=1, line=1
  oplot, day, data1[i+3,*], color=3, line=1
  oplot, day, data2[i+4,*], color=1, line=1
  oplot, day, data1[i+4,*], color=3, line=1
  oplot, day, data2[i+5,*], color=1, line=1
  oplot, day, data1[i+5,*], color=3, line=1
  oplot, day, data2[i+6,*], color=1, line=1
  oplot, day, data1[i+6,*], color=3, line=1
;precip
  oplot, day, data1[3,*]*100

  i=6
  plot, day, data1[i,*], xr=xr, /xs, yr=[-100,500], /ys, $
    line=2, title=title+' Sensible Heat Flux', $
    ytitle='Sensible Heat (W/m^2)', xtitle='Day'
  oplot, day, data2[i,*], color=1
  oplot, day, data1[i,*], color=3
  
  i=7
  plot, day, data1[i,*], xr=xr, /xs, yr=[-100,500], /ys, $
    line=2, title=title+' Latent Heat Flux', $
    ytitle='Latent Heat (W/m^2)', xtitle='Day'
  oplot, day, data2[i,*], color=1
  oplot, day, data1[i,*], color=3

  i=11+2*nsoil
  plot, day, data1[i,*], xr=xr, /xs, yr=[-100,500], /ys, $
    line=2, title=title+' Ground Heat Flux', $
    ytitle='Ground Heat Flux (W/m^2)', xtitle='Day'
  oplot, day, data2[i,*], color=1
  oplot, day, data1[i,*], color=3
  
  
;  device, /close

;  !p.multi=oldp
;  set_plot, olddevice

end