PRO jjj
@~vgr/idl/common.types
read_voy,"$vgrDATA/v2/ha/v2_jwb_",77,95
n=n_elements(tt)
yr=dd(n-1).yr
tt = dd.yr - 1900 + (dd.day + dd.hr/24.0)/366.
pr= dd.df * ( dd.vr * dd.vr * dd.rad * dd.rad ) * 1.67e-6
wset  ,0
plot,tt,pr,psym=3,xstyle=1,ystyle=1
help,tt,pr,dd
n=n-1
day=fltarr(n)
pra=fltarr(n)
pr=0
j=0
nn=0
i = 0l
idayo = dd(i).yr*366 + dd(i).day
for i=1l , n do begin
	if  ( i eq n ) or ( dd(i).yr*366 + dd(i).day ge idayo + 25 ) then begin
		pra(j) = 1.67e-6 * pr / nn
		if ( dd(i-1).yr mod 4 ) eq 0 then yl = 366.0 else yl = 365.0
		day(j) = dd(i-1).yr + dd(i-1).day / yl
		idayo = dd(i).yr*366 + dd(i).day
		j = j+1
		pr = 0
		nn = 0
	endif
	pr = pr + dd(i).num * 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)
day = day - 1900
pra = pra(ii)
help,day,pra,n
wset  ,1
plot,day,pra,psym=3,xstyle=1,ystyle=1
wset,2
plot,day,pra,psym=0,xstyle=1,ystyle=1
print,day(0),day(n_elements(day)-1)
print,pra(0),pra(n_elements(pra)-1)
end ; jjj