FUNCTION fit_simul_dainty_wind_function,p

common func_simul_dainty_wind_xy, x, y1, y2, y3, y4

;The function we are trying to fit is
;y=((d+(k0/kl)^(1/(b-4))*D)/po))^(b-2)
;
;p(0)=ko
;p(1)=wind
;p(2)=kl
;p(3)=kt
;p(4)=b
;p(5)=po

kl=(p(4)-3)*p(3)
a=1./(p(4)-4)
b=p(4)-2
koxl=float(p(0)/kl)
koyl=float(p(0)/kl))
koxt=float(p(0)/p(3))
koyt=float(p(0)/p(3))


;xtrans
;For this use
;kt to fit xt
;k0 to fit xo
xtrans_tot=total((y1-p(3)*(((x+(koxt^a)*.17)/p(5))^b*$
	          (5e-7/(x+(koxt^a)*.17))^2))^2)
;xlong
;For this use
;w*kl to fit xl
;w*k0 to fit x0
xlong_tot=total((y2-p(1)*kl*(((x+(koxl^a)*.17)/p(5))^b*$
                  (5e-7/(x+(koxl^a)*.17))^2))^2)
;ytrans
;For this use
;w*kt to fit yt
;w*k0 to fit y0
ytrans_tot=total((y3-p(1)*p(3)*(((x+(koyt^a)*.17)/p(5))^b*$
                  (5e-7/(x+(koyt^a)*.17))^2))^2)

;ylong
;For this use
;kl to fit yl
;k0 to fit y0
ylong_tot=total((y4-kl*(((x+(koyl^a)*.17)/p(5))^b*$
                  (5e-7/(x+(koyl^a)*.17))^2))^2)

print,((ylong_tot+ytrans_tot+xlong_tot+xtrans_tot))
return,((ylong_tot+ytrans_tot+xlong_tot+xtrans_tot))

end