02 AUG 2014 19 00 00
 JD =   2456872.29 MJD =   56871.7909
 P-angle =   11.4595449 B0 =   5.89171409
  Distance =   1.01472944 Radius =   945.710235

Additional limb height at radio:
limb=1.+66.6/960.0*(freq^(-0.60))=1.0374 at 2.8 GHz

FWHM FOV = 45 arcmin/f_GHz
= 16.1 arcmin @ 2.8 GHz
so need 8 arcmin separation of fields
try 2-3-3-2

cx=8.0*cos(indgen(101)/100.*2.*!pi)
cy=8.0*sin(indgen(101)/100.*2.*!pi)
sx=945.7/60.*cos(indgen(101)/100.*2.*!pi)*1.0374
sy=945.7/60.*sin(indgen(101)/100.*2.*!pi)*1.0374

plot,xx,yy,xrange=[-20,20],yrange=[-20,20],/nodata
oplot,sx,sy

; 2-3-3-2
xx=[-4,4,-9,0,9,-9,0,9,-4,4]
yy=[11,11,5,4,5,-5,-4,-5,-11,-11]
for i=0,9 do oplot,xx[i]+cx,yy[i]+cy,linestyle=2

; no, elongate at the equator, truncate at poles
; xx=[11,11,5,4,5,-5,-4,-5,-11,-11]
; yy=[-4,4,-10,0,10,-10,0,10,-4,4]
xx=[10.0,10.0,4.5,4,4.5,-4.5,-4,-4.5,-10.0,-10.0]
yy=[-4,4,-10.5,0,10.5,-10.5,0,10.5,-4,4]
for i=0,9 do oplot,xx[i]+cx,yy[i]+cy,linestyle=2

; version with one at disk center instead of two straddling - no
xx=[11,11,5,0,5,-5,-5,-11,-11]
yy=[-4,4,-10,0,10,-10,10,-4,4]
for i=0,8 do oplot,xx[i]+cx,yy[i]+cy,linestyle=2

; 2-4-4-2
xx=[-4,4,-12,-4,4,12,-12,-4,4,12,-4,4]
yy=[12,12,4,4,4,4,-4,-4,-4,-4,-12,-12]
for i=0,11 do oplot,xx[i]+cx,yy[i]+cy,linestyle=2

;-------------------------------------------------------------------------

for i=0,9 do print,sqrt((xx-xx[i])^2+(yy-yy[i])^2),format='(10(F4.1,1x))'
for i=0,9 do begin & vals=sqrt((xx-xx[i])^2+(yy-yy[i])^2) & print,vals[sort(vals)],format='(10(F4.1,1x))'

;-------------------------------------------------------------------------

; first draft, see below
xx=[11,11,5,4,5,-5,-4,-5,-11,-11]
yy=[-4,4,-10,0,10,-10,0,10,-4,4]
tref='02 AUG 2014 20 00 00'
for i=0,9 do begin & nsew_to_lb,tref,xx[i]*60.,yy[i]*60.,b,l & if (b ge 0.) then LAS='N' else LAS='S' & if (l ge 0.) then LOS='W' else LOS='E' & print,LAS+string(abs(b),format='(F04.1)')+LOS+string(abs(l),format='(F04.1)') & end

S10.5W45.0
N18.7W47.3
S33.8W22.4
N05.7W14.7
N44.6W26.3
S33.8E22.4
N05.7E14.7
N44.6E26.3
S10.5E45.0
N18.7E47.3

bb=[10.5,18.7,-33.8,5.7,44.6,-33.8,5.7,44.6,-10.5,18.7]
ll=[45.0,47.3,22.4,14.7,26.3,-22.4,-14.7,-26.3,-45.0,-47.3]

reg=['S10.5W45.0','N18.7W47.3','S33.8W22.4','N05.7W14.7','N44.6W26.3','S33.8E22.4','N05.7E14.7','N44.6E26.3','S10.5E45.0','N18.7E47.3']

for i=0,9 do generate_fields,'horizons_140802.ephem',reg[i],'20:00:00',ar1,order=3

rfl=['N05E14_140802.radec','N44E26_140802.radec','S33E22_140802.radec','N05W14_140802.radec','N44W26_140802.radec','S33W22_140802.radec','N18E47_140802.radec','S10E45_140802.radec','N18W47_140802.radec','S10W45_140802.radec']

plot_fields,'horizons_140802.ephem',rfl[i],2.8,pang=11.48

.run ../plot_fields_mosaic
plot_fields,'horizons_140802.ephem','_140802.radec',2.8

;-------------------------------------------------------------------------

; final
xx=[10.,10.,5,4,5,-5,-4,-5,-10.,-10.]             
yy=[-4,4,-10.5,0,10.5,-10.5,0,10.5,-4,4]
tref='02 AUG 2014 20 00 00'
for i=0,9 do begin & nsew_to_lb,tref,xx[i]*60.,yy[i]*60.,b,l & if (b ge 0.) then LAS='N' else LAS='S' & if (l ge 0.) then LOS='W' else LOS='E' & print,LAS+string(abs(b),format='(F04.1)')+LOS+string(abs(l),format='(F04.1)') & end

reg=['S10.2W40.0','N19.1W42.0','S36.2W23.1','N05.7W14.7','N46.9W27.6','S36.2E23.1','N05.7E14.7','N46.9E27.6','S10.2E40.0','N19.1E42.0']

bb=[-10.2,19.1,-36.2,5.7,46.9,-36.2,5.7,46.9,-10.2,19.1]
ll=[40.0,42.0,23.1,14.7,27.6,-23.1,-14.7,-27.6,-40.0,-42.0]

for i=0,9 do generate_fields,'horizons_140802.ephem',reg[i],'20:00:00',ar1,order=3

rfl=['N05E14_140802.radec','','N44E26_140802.radec','','S33E22_140802.radec','N05W14_140802.radec','','N44W26_140802.radec','','S33W22_140802.radec','N18E47_140802.radec','','S10E45_140802.radec','N18W47_140802.radec','','S10W45_140802.radec']

plot_fields,'horizons_140802.ephem',rfl[i],2.8,pang=11.48

.run ../plot_fields_mosaic
plot_fields_mosaic,'horizons_140802.ephem','_140802.radec',2.8,sunrad=1.0374*945.7/60.,pang=11.48