;$Id: po32.pro,v 1.8 2008/12/18 09:50:10 brandenb Exp $
@$HOME/tex/hubbard/memory/idl/fit_respfunct
@$HOME/tex/hubbard/memory/idl/realspace_respfunct
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=23,yoffset=3
  !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
end
;
;  mv idl.ps ~/tex/hubbard/memory/fig/forced_po32.ps
;
file='o32.dat' & urms0=.12 & kf=3.
a=rtable(file,6,head=1)
lam=reform(a(0,*))
om=reform(a(1,*))
;
om_uk=om/(urms0*kf)
alp0=urms0/3.
etat0=urms0/(3.*kf)
;
alp_re=reform(a(2,*))/alp0
alp_im=reform(a(3,*))/alp0
eta_re=reform(a(4,*))/etat0
eta_im=reform(a(5,*))/etat0
;
!p.charsize=1.8
!p.multi=[0,1,2]
!p.multi=[0,1,4] ;XX
!x.range=[.03,70]
!y.margin=[3.6,.8]
tilde='!9!s!aA!n!r!6'
circ_sym,1.4,1
fo='(a,2x,3f6.2)'
col=0
;
;  alpha
;
!y.range=[-.1,.65]
!x.title='!6'
!y.title=tilde+'!7a!6(!7x!6) / !7a!6!d0!n'
plot_oi,om_uk,-alp_re,ps=8
oplot,om_uk,-alp_im,li=2,ps=4
oplot,!x.range,!x.range*0.,li=0
;
;  automatic fit (blue is the initial guess)
;
fitparam=[1.0,2.0,.6]
ooo=grange(.01,50.,30,/log)
funct,ooo,fitparam,fff,pder
;oplot,ooo,float(fff),col=55
;oplot,ooo,imaginary(fff),li=2,col=55
fit_respfunct,om_uk,complex(-alp_re,-alp_im),fitparam,/oplo,col=col
print,'fitparam=',fitparam,fo=fo
;
;  phase (to be commented out)
;
plot_oi,om_uk,arg(complex(-alp_re,-alp_im))/!pi,ps=-8,yr=[-.1,1.1]*.5
;
;  etat
;
!y.range=[-.1,.7]
!x.title='!7x!6 / !8u!6!drms!n!8k!6!df!n'
!y.title=tilde+'!7g!6!dt!n(!7x!6) / !7g!6!dt0!n'
plot_oi,om_uk,eta_re,ps=8
oplot,om_uk,eta_im,li=2,ps=4
oplot,!x.range,!x.range*0.,li=0
;
;  automatic fit (blue is the initial guess)
;
fitparam=[.48,1.40,.39]
ooo=grange(.01,50.,30,/log)
funct,ooo,fitparam,fff,pder
;oplot,ooo,float(fff),col=55
;oplot,ooo,imaginary(fff),li=2,col=55
fit_respfunct,om_uk,complex(eta_re,eta_im),fitparam,/oplo,col=col
print,'fitparam=',fitparam,fo=fo
;
;  phase (to be commented out)
;
plot_oi,om_uk,arg(complex(eta_re,eta_im))/!pi,ps=-8,yr=[-.1,1.1]*.5
;
;realspace_respfunct,om_uk,complex(eta_re,eta_im)
;
END