PRO rcp37,imin,imax         

; pspc,17   1-10
; on 16 bit monitors first execute: device,decomposed=0

col_ind9

!p.multi=0

rcp_name = [ $
'22583.rcp','22842.rcp','23042.rcp','23307.rcp','24066.rcp', $
'24961.rcp','25390.rcp','26084.rcp','26655.rcp','44823.rcp', $
'46039.rcp','46086.rcp' ]

init = 0
ist =7
sy = findgen(37) * (!PI*2/36.)

g1 = fltarr(37)
g2 = g1 & gav = g1
az = g1 & azav = az
el = g1 & elav = el
null = el
ich = strarr(37)
col = [2,3,4,5,6,7,8,2,3,4,5,6,7]

for i=imin,imax do begin

readcol,rcp_name(i),f='i,f,f,f,f,f', $
      ich,g1,g2,az,el,dum,ski=0,num=37

az(17) = 36.
el(17) = -42.

gav = gav + g1
azav = azav + az
elav = elav + el
print,(total(g1)+0.)/36.

    usersym,5.3*cos(sy),5.3*sin(sy)

    if init eq 0 then plot,az,el, title='MAMBO-37 rcp_2002may',$
      xrange=[-80,80],yrange=[-80,80],psym=8,/nodata 
        if init eq 0 then  xyouts,az+8,el,string(ich,format='(i2)')
    init = 1

    for j=0,36 do begin
        usersym,3*g1(j)*cos(sy),3*g2(j)*sin(sy)
        oplot,[az(j)],[el(j)],psym=8,color=col(i-imin)
    endfor

endfor

nn = imax - imin + 1
gav = gav / nn
gav = gav / ( total(gav)/36. )
elav = elav / nn
azav = azav / nn

    for j=0,36 do begin
        usersym,3*gav(j)*cos(sy),3*gav(j)*sin(sy)
        oplot,[azav(j)],[elav(j)],psym=8,color=0
        print,j+1,gav(j),gav(j),azav(j),elav(j),null(j), $
                format='(i4,2f10.3,3f10.3)'
    endfor


END

;==========================================================

PRO rcp117,imin,imax         ; pspc,17   1-10

col_ind9
!p.multi=0

rcp_name = [ $
 '14576.rcp','14687.rcp','15182.rcp','17532.rcp','17859.rcp', $
 '17965.rcp','20569.rcp','34104.rcp','36390.rcp','36526.rcp', $
 '44651.rcp']

rcp_1 = [ '14687.rcp','15182.rcp','17532.rcp','17965.rcp' ]
rcp_2 = [ '34104.rcp','36390.rcp','36526.rcp','44651.rcp']
rcp_3 = [ rcp_1, rcp_2 ]
rcp_4 = [ '120ave.rcp','14576.rcp','17859.rcp','20569.rcp']

init = 0
ist =7
sy = findgen(37) * (!PI*2/36.)

g1 = fltarr(117)
g2 = g1 & gav = g1
az = g1 & azav = az
el = g1 & elav = el
null = el
cnt = intarr(117)*0
ich = strarr(117)
;      1 2 3 4 5 6 7 8 9 0 1
col = [0,2,3,4,5,6,7,8,2,3,4,5]


for i=imin,imax do begin

readcol,rcp_name(i),f='i,f,f,f,f,f', $
      ich,g1,g2,az,el,dum,ski=0,num=117

az = az - az(52)
el = el - el(52)

good = where(g1 ge 0.1 and g1 le 2., igood)
bad  = where(g1 le 0.1 or  g1 gt 2., ibad)
print,rcp_name(i),igood,ibad,bad,format='(a10,100i4)'
g_good = g1(good)
g1(good) = g1(good) / mean( g_good )
if ibad ge 1 then g1(bad)  = 0.
if ibad ge 1 then az(bad)  = 100.
if ibad ge 1 then el(bad)  = -120.


;az(17) = 36. & el(17) = -42.

gav(good)  = gav(good) + g1(good)
azav(good) = azav(good) + az(good)
elav(good) = elav(good) + el(good)
cnt(good)  = cnt(good) + 1

;print,(total(g1)+0.)/36.

    if init eq 0 then plot,az,el, title='MAMBO-117 rcp_2001/2',$
      xrange=[-130,130],yrange=[-130,130],psym=3,/nodata 
        if init eq 0 then  xyouts,az+5,el,string(ich,format='(i3)')
    init = 1

    for j=0,116 do begin
        usersym,2*g1(j)*cos(sy),2*g1(j)*sin(sy)
        oplot,[az(j)],[el(j)],psym=8,color=col(i-imin)
    endfor

endfor

cnt(where(cnt eq 0)) = 1
nn = imax - imin + 1
gav = gav / cnt
;gav = gav / ( total(gav)/36. )
elav = elav / cnt
azav = azav / cnt

good = where(gav gt 0.1, igood)
bad  = where(gav le 0.1 or gav gt 2., ibad)

print,'mean=', mean( gav(good) )
gav = gav / mean( gav(good) )
gav(bad) = 1.0

    for j=0,116 do begin
        usersym,2*gav(j)*cos(sy),2*gav(j)*sin(sy)
        oplot,[azav(j)],[elav(j)],psym=8 ;,color=0
        print,j+1,gav(j),gav(j),azav(j),elav(j),null(j), $
                format='(i4,2f10.3,3f10.3)'
    endfor


END