; Plots image and overlays contours in a .ps file
pro posplot, imag, nl_xy, nl_info, coords, edpt, fn, nl_xy_2=nl_xy_2, nl_info_2=nl_info_2, coords_2=coords_2, nl_xy_3=nl_xy_3, nl_info_3=nl_info_3, coords_3=coords_3, no_lines=no_lines, t_obs=t_obs

if (n_elements(no_lines) eq 0) then no_lines=0
sz=size(imag)

if n_elements(fn) eq 0 then begin
print, "filename:"
fn='' & read, fn
endif

set_plot,'ps'
device,filename=(fn+'.ps'),xsize=16,ysize=16*sz(2)/sz(1),xoffset=2,yoffset=5, bits=8,/color 
r=indgen(256) & g=r & b=g & g(253:255)=[0,255,0] & b(253:255)=[255,0,0] & r(253:255)=[0,0,255]
tvlct,r,g,b

tv,bytscl(imag,min=-500,max=500)<252,0,0,xsize=1,ysize=1,/norm
	; set ploting

contour,imag,pos=[0,0,1,1],xstyle=5,ystyle=5,/norm,/noerase,/nodata

	if (no_lines eq 0) then begin

loadct,8

;for i=0, n_elements(nl_info)-1 do begin
;if nl_info(i).typ lt 2 then if nl_info(i).typ eq 1 then s=[indgen(nl_info(i).n),0] else s=[indgen(nl_info(i).n)]
;if nl_info(i).typ ne 2 then BEGIN
;  oplot, nl_xy(0,nl_info(i).offset+s),nl_xy(1,nl_info(i).offset+s), color=0, thick=23, linestyle=0
;  oplot, nl_xy(0,nl_info(i).offset+s),nl_xy(1,nl_info(i).offset+s), color=150, thick=14,linestyle=2 
;ENDIF ELSE BEGIN
;  oplot, [nl_xy(0,nl_info(i).offset), nl_xy(0,nl_info(i).n)],[nl_xy(1,nl_info(i).offset), nl_xy(1,nl_info(i).n)], color=0, thick=23, linestyle=0
;  oplot, [nl_xy(0,nl_info(i).offset), nl_xy(0,nl_info(i).n)],[nl_xy(1,nl_info(i).offset), nl_xy(1,nl_info(i).n)], color=150, thick=14, linestyle=2
;ENDELSE
;endfor

loadct, 3
                                                                                
;if n_elements(nl_xy_2) ne 0 then begin
;for i=0, n_elements(nl_info_2)-1 do begin if nl_info_2(i).typ lt 2 then if nl_info_2(i).typ eq 1 then s=[indgen(nl_info_2(i).n),0] else s=[indgen(nl_info_2(i).n)]
;if nl_info_2(i).typ ne 2 then BEGIN
; oplot, nl_xy_2(0,nl_info_2(i).offset+s),nl_xy_2(1,nl_info_2(i).offset+s),color=0,thick=23,linestyle=0
; oplot,nl_xy_2(0,nl_info_2(i).offset+s),nl_xy_2(1,nl_info_2(i).offset+s),color=130,thick=2.5,psym=6
;ENDIF else BEGIN 
;  oplot, [nl_xy_2(0,nl_info_2(i).offset), nl_xy_2(0,nl_info_2(i).n)],[nl_xy_2(1,nl_info_2(i).offset), nl_xy_2(1,nl_info_2(i).n)], color=0, thick=23,linestyle=0
;  oplot, [nl_xy_2(0,nl_info_2(i).offset), nl_xy_2(0,nl_info_2(i).n)],[nl_xy_2(1,nl_info_2(i).offset), nl_xy_2(1,nl_info_2(i).n)], color=130, thick=2.5,psym=6
;ENDELSE
;endfor
;endif

loadct,1

if n_elements(nl_xy_3) ne 0 then begin
for i=0, n_elements(nl_info_3)-1 do begin if nl_info_3(i).typ lt 2 then if nl_info_3(i).typ eq 1 then s=[indgen(nl_info_3(i).n),0] else s=[indgen(nl_info_3(i).n)]
if nl_info_3(i).typ ne 2 then BEGIN
  oplot,nl_xy_3(0,nl_info_3(i).offset+s),nl_xy_3(1,nl_info_3(i).offset+s),color=0,thick=23,linestyle=0
  oplot,nl_xy_3(0,nl_info_3(i).offset+s),nl_xy_3(1,nl_info_3(i).offset+s),color=230,thick=16,linestyle=2
ENDIF else BEGIN 
  oplot, [nl_xy_3(0,nl_info_3(i).offset), nl_xy_3(0,nl_info_3(i).n)],[nl_xy_3(1,nl_info_3(i).offset), nl_xy_3(1,nl_info_3(i).n)], color=0, thick=23, linestyle=0
  oplot, [nl_xy_3(0,nl_info_3(i).offset), nl_xy_3(0,nl_info_3(i).n)],[nl_xy_3(1,nl_info_3(i).offset), nl_xy_3(1,nl_info_3(i).n)], color=230, thick=16, linestyle=2

ENDELSE
endfor
endif

;for i=0, n_elements(edpt)-1 do begin
;for j=0, n_elements(where(edpt(i).pt))-1 do begin
;m=where(edpt(i).pt)
;oplot, [coords(0,edpt(i).pt(m(j)))],[coords(1,edpt(i).pt(m(j)))],symsize=.5,psym=2, thick=2, color=255
;endfor
;endfor

	endif

xyouts,sz(1)-100, 30, t_obs, color=0, charsize=3, charthick=4.5, /device

device,/close

print, format='("neutral line output to ",A)',(fn+'.ps')

end