Abrir - RDU

Apéndices 301
found=0
do i=1,np1
if (lam13(i).gt.desde) then
found=1
exit
endif
enddo
if (found.eq.0) then
write(*,*) ”desder13 not found”
write(52,*) ”desder13 not found”
stop
endif
nd13=i-1
!busca el ultimo pto. q es posterior al hasta, ”hastar13”
found=0
do j=nd13,np1
if (lam13(j).gt.hasta) then
found=1
exit
endif
enddo
if (found.eq.0) then
write(*,*) ”hastar13 not found”
write(52,*) ”hastar13 not found”
stop
endif
nh13=j
!spec2 sera el resampleado del spec1, va desde 2500-10500
en lambda
!igual q el R13, p comparar
!ahora resamplea usando los lambda1,2 mas proximos al
punto
do j=nd13+1,nh13-1
x=lam13(j)
x1=0.
x2=0.
do i=1,2e5
if (spec1(i,1).gt.x) then
nx2=i
x2=spec1(i,1)
exit
endif
enddo
210 continue
if (spec1(i,1).lt.x) then
nx1=i
x1=spec1(i,1)
else
i=i-1
goto 210
endif
if (x1.lt.1.) then
write(*,*) ”problems searching x1 for lamda x=”,x
write(52,*) ”problems searching x1 for lamda x=”,x
stop
endif
if (x2.lt.1.) then
write(*,*) ”problems searching x2 for lamda x=”,x
write(52,*) ”problems searching x2 for lamda x=”,x
stop
endif
!write(*,*) x1,x,x2
!tenemos x1,x2, ahora calcula las alturas respectivas y1,y2
y1=spec1(nx1,2)
y2=spec1(nx2,2)
call interp1d(x2,x1,y2,y1,x,px)
spec4(j)=px
enddo
spec4(nd13)=1
spec4(nh13)=1!los extremos del intervalo
!spec4 con el mismo sampleo y mismo rango 2500-10500
!que el archivo R13
end
!********************************************************
subroutine velrad(desde,hasta,spec1,vrestim,spec5)
!esta sub intenta determinar la vrad de una estrella
!usando el minimo chicuad
!entrada=desde,hasta,spec1 (todos SIN corregir x vrad)
!salida=vrestim (spec5 es temporal, no es de salida)
common /comunes/
spec2,lam13,vrot,points,ndesder13,nhastar13,
np1,wspec,file
real*8 spec2(2e5),lam13(2e5),vrot,vrad
integer points
integer ndesder13,nhastar13,np1,wspec
character file*80
save /comunes/
real*8 spec1(2e5,2),spec5(2e5,2),spec6(2e5)
!1 es el observado, 5 es correg. x vrad,6 es resampleado
real*8 sp1(2e5),sp2(2e5),sp3(2e5)!3 espectros para 3 temps.
!diferentes: 5000, 7500, 10000
real*8 vradi,frad,desde,hasta
real*8 desde1,hasta1
integer nd,nh,nd1,nh1
real*8 chi1,chi2,chi3,minchi1,minchi2,minchi3
real*8 minvr1,minvr2,minvr3
real*8 vrestim
open(54,file=”sp1.dat”,status=. o ld”)!lee los 3 espectros
open(55,file=”sp2.dat”,status=. o ld”)
open(56,file=”sp3.dat”,status=. o ld”)
i=1
do
read(54,*,end=221) sp1(i)
read(55,*,end=221) sp2(i)
read(56,*,end=221) sp3(i)
i=i+1
enddo
221close(54)