To obtain a beam pattern IDL> beam=norh_beam(index) CR This is a pixels array. To display this pattern (contour of half power beam width [HPBW]) IDL> ocontour,beam,levels=[0.5] CR