PRO jjj
read_voy,"$vgrDATA/v2/ha/v2_jwb_",90,95
n=n_elements(tt)
yr=dd(n-1).yr
tt = (dd.yr ) + (dd.day + dd.hr/24.0)/366.5.
pr= dd.df * ( dd.vr * dd.vr * dd.rad * dd.rad ) * 1.67e-6
plot,tt,pr,psym=3,xstyle=1,ystyle=1
n=n-1
day=fltarr(n)
pra=fltarr(n)
pr=0
j=0
nn=0
for i=1 , n-1 do begin
	if ( dd(i).day ne dd(i-1).day ) then begin
		pra = 1.67e-6 * pr / nn
		if ( dd(i-1).yr mod 4 ) eq 0 then yl = 366 else yl = 365
		day(j) = dd(i-1).yr + dd(i-1).day / yl
		j = j+1
		pr = 0
		nn = 0
	endif
	pr = pr + dd(i).df * dd(i).vr * dd(i).vr * dd(i).rad * dd(i).rad
	nn = nn + dd(i).num
endfor
ii=where( pra ne 0.)
day=day(ii)
pra = pra(ii)
plot,day,pra,psym=3,xstyle=1,ystyle=1
end ; jjj