Abrir - RDU
Abrir - RDU
Abrir - RDU
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Apéndices 302<br />
close(55)<br />
close(56)<br />
write(*,*) ”Estimating radial velocity...”<br />
write(52,*) ”Estimating radial velocity...”<br />
do i=1,2e5<br />
spec5(i,2)=spec1(i,2)<br />
enddo<br />
minchi1=1e10<br />
minchi2=1e10<br />
minchi3=1e10<br />
do vradi=-100.,100.,2.<br />
frad=( (-vradi/299800.)+1)**-1<br />
do i=1,2e5<br />
spec5(i,1)=spec1(i,1)*frad<br />
enddo<br />
desde1=desde*frad<br />
hasta1=hasta*frad<br />
!lo resamplea para comparar con munari, spec6 es el resampleado<br />
call resample(desde1,hasta1,spec5,spec6,nd1,nh1)<br />
call comspec1(spec6,sp1,nd1,nh1,chi1)<br />
call comspec1(spec6,sp2,nd1,nh1,chi2)<br />
call comspec1(spec6,sp3,nd1,nh1,chi3)<br />
if (chi1.lt.minchi1) then<br />
minchi1=chi1<br />
minvr1=vradi<br />
endif<br />
if (chi2.lt.minchi2) then<br />
minchi2=chi2<br />
minvr2=vradi<br />
endif<br />
if (chi3.lt.minchi3) then<br />
minchi3=chi3<br />
minvr3=vradi<br />
endif<br />
if (vradi/10.eq.int(vradi/10)) then<br />
write(* ,*) vradi,chi1,chi2,chi3<br />
write(52,*) vradi,chi1,chi2,chi3<br />
endif<br />
enddo<br />
if ((chi1.le.chi2).and.(chi1.le.chi3)) then<br />
!elige la q tiene menor chicuad entre las 3<br />
vrestim=minvr1<br />
endif<br />
if ((chi2.le.chi1).and.(chi2.le.chi3)) then<br />
vrestim=minvr2<br />
endif<br />
if ((chi3.le.chi1).and.(chi3.le.chi2)) then<br />
vrestim=minvr3<br />
endif<br />
minchi2=1e10<br />
minchi3=1e10<br />
do vradi=vrestim-2.,vrestim+2.,0.1<br />
frad=( (-vradi/299800.)+1)**-1<br />
do i=1,2e5<br />
spec5(i,1)=spec1(i,1)*frad<br />
enddo<br />
desde1=desde*frad<br />
hasta1=hasta*frad<br />
call resample(desde1,hasta1,spec5,spec6,nd1,nh1)<br />
call comspec1(spec6,sp1,nd1,nh1,chi1)<br />
call comspec1(spec6,sp2,nd1,nh1,chi2)<br />
call comspec1(spec6,sp3,nd1,nh1,chi3)<br />
if (chi1.lt.minchi1) then<br />
minchi1=chi1<br />
minvr1=vradi<br />
endif<br />
if (chi2.lt.minchi2) then<br />
minchi2=chi2<br />
minvr2=vradi<br />
endif<br />
if (chi3.lt.minchi3) then<br />
minchi3=chi3<br />
minvr3=vradi<br />
endif<br />
write(* ,*) vradi,chi1,chi2,chi3<br />
write(52,*) vradi,chi1,chi2,chi3<br />
enddo<br />
if ((chi1.le.chi2).and.(chi1.le.chi3)) then<br />
vrestim=minvr1<br />
endif<br />
if ((chi2.le.chi1).and.(chi2.le.chi3)) then<br />
vrestim=minvr2<br />
endif<br />
if ((chi3.le.chi1).and.(chi3.le.chi2)) then<br />
vrestim=minvr3<br />
endif<br />
write(*,* ) ”3 vrad values: ”,minvr1,minvr2,minvr3<br />
write(52,*) ”3 vrad values: ”,minvr1,minvr2,minvr3<br />
write(*,* ) ”Final Vrad estimated =”,vrestim<br />
write(52,*) ”Final Vrad estimated =”,vrestim<br />
end<br />
write(* ,*) ”Increasing presicion in Vrad estimation...”<br />
write(52,*) ”Increasing presicion in Vrad estimation...”<br />
minchi1=1e10