pro plotflux_lc

;------------------------------------------------------------------------------
; This program plots CO2, H2O, Temperature and Momentum fluxes from the 
; .sflx data file produced by Dayflux.pro.
;
; Created by Dana Carrington and Bruce Cook - August, 1999.
; Multiple day option added December 1999 (BCook)
; Modified by W wang for Lost Creek, Sep 20 2000
;------------------------------------------------------------------------------

print, 'Enter the desired date for flux plots (yyddd):'
startdate = STRING(5)
read, startdate

print, 'Enter number of days to plot (do not overlap Jan 1):'
days=1.
read,days

;extract file name information
yy=strtrim(strmid(startdate,0,2),2)
startday=fix(strmid(startdate,2,3))

;start loop for multiple days
for z=1,days do begin

ddd=startday+(z-1)
if (ddd lt 10) then filename=yy+'00'+strtrim(string(ddd),2) else $
	if (ddd lt 100) then filename=yy+'0'+strtrim(string(ddd),2) else filename=yy+strtrim(string(ddd),2)

cmd = "load, '/davis/s1/wang/d2wwg/lcreek/output/flx/'+filename+'.sflx', flx"
r = execute(cmd)

set_plot, 'ps'
fname = '/davis/s1/wang/d2wwg/lcreek/output/plots/flux/'+filename+'flux.ps'
device, filename = fname, /landscape 

!p.multi=[0,2,2,0,0]

labels=['0','4','8','12','16','20','24']

Plot, flx(1,*), title='CO2 FLUX', ytitle='WC(umol/(s*m^2))', $
      xticks=6,xtickname=labels,xminor=4,xstyle=1,min_value=-998.,ystyle=1

xyouts,.23,1.02,'LOST CREEK FLUX PLOT FOR '+filename, /normal,$
     charsize = 2,  charthick=2

plot, flx(2,*), title='LATENT HEAT FLUX', ytitle='WQ(W/m^2)', $
      xticks=6,xtickname=labels,xminor=4,xstyle=1,min_value=-998.

plot, flx(3,*), title='SENSIBLE HEAT FLUX', ytitle='WT(W/m^2)', $
      xticks=6,xtickname=labels,xminor=4,xstyle=1,min_value=-998.

plot, flx(4,*), title='MOMENTUM FLUX', ytitle='SQRT(WU**2+WV**2) (Pa)', $
      xticks=6,xtickname=labels,xminor=4,xstyle=1,min_value=-998.

cplot
device, /close
!p.multi = 0

endfor; end multiple day loop

return
end