;   14-Mar-2014  DG
;      Renamed by appending ovsa_ to avoid potential conflict with the copies of the original
;      routines located in the ssw\hessi tree
pro ovsa_plot_vis2vis_simpl, vis_obs, vis_mdl

;big_tit ='I M A G R / V i s 2 v i s - s i m p l e'
big_tit ='Vis_obs (psym10) vs. Vis_model (line)'

sze=size(vis_obs)
npts=sze[4]
yobs=reform(vis_obs,npts)
ymdl=reform(vis_mdl,npts)
	xtit='u,v points'
	ytit='Visibility [s.f.u]'
	tit=['V_Real [s.f.u]','V_Imaginary [s.f.u]']
;------------------------window size/pos stuff
 xwt 	=(800.-400.)*float(npts-500)/500.
 xwt    =(xwt>400.)<790. 
 xw	= xwt-75
 ywt    =360 & yw  =120
 xpos	=[55,  55]+8
 ypos	=[200, 45] 
 tpos	=[xpos[0],ypos[0]+yw+20]

;------------------------ plot

window,3,xs=xwt, ys=ywt
device,decompo=0
!p.multi=[0,1,2]
!p.charsize=1.2
!p.font=-1

;-----upper panel = Vis_real

yy = float(yobs)
zz = float(ymdl)
xp0=xpos[0] & xp1=xp0+xw & yp0=ypos[0] & yp1=yp0+yw
 
  plot,yy,psym=10, pos=[xp0,yp0,xp1,yp1],/dev, $
	ytit=tit[0], xst=1
  oplot,zz,thick=2

;-----lower panel = Vis_Imaginary

yy = imaginary(yobs)
zz = imaginary(ymdl)
xp0=xpos[1] & xp1=xp0+xw & yp0=ypos[1] & yp1=yp0+yw
 
  plot,yy,psym=10, pos=[xp0,yp0,xp1,yp1],/dev, $
	xtit=xtit,ytit=tit[1],xst=1
  oplot,zz,thick=2

xyouts,tpos[0],tpos[1],big_tit,charsiz=2,charthic=2,/device

return
end