Abrir - RDU

Apéndices 296
path=/dos/munari/”
fname=path//’LAMBDA R13.DAT’
open(1,file=fname,status=. o ld”)
np1=1
do
read(1,*,end=200) lam13(np1)
np1=np1+1
enddo
200close(1)
call resample(desde,hasta,spec1,spec2,ndesder13,nhastar13)
if (vrad.eq.9999.) then
call velrad(desde,hasta,spec1,vrestim,spec5)
vrad=vrestim
frad=( (-vrad/299800.)+1)**-1
desde=desde*frad
hasta=hasta*frad
do i=1,2e5
spec1(i,1)=spec1(i,1)*frad
enddo
call resample(desde,hasta,spec1,spec2,ndesder13,nhastar13)
endif
!4to vertice
p3(1)=p0(1)
p3(2)=p0(2)
p3(3)=p0(3)+0.2
p3(4)=p0(4)
!5to vertice
p4(1)=p0(1)
p4(2)=p0(2)
p4(3)=p0(3)
p4(4)=p0(4)!+1.0
!escribimos ahora la matriz p usando los 5 vertices
p(1,1)=p0(1)
p(1,2)=p0(2)
p(1,3)=p0(3)
p(1,4)=p0(4)
p(2,1)=p1(1)
p(2,2)=p1(2)
p(2,3)=p1(3)
p(2,4)=p1(4)
!call force1(teff0,logg0,feh0,xita0)
if (logg0.ge.5.0) then
logg0=4.5
endif
write(*,*)
write(*,*) ”Initial guess of 4d points: Teff logg [Fe/H] Xita
”
write(52,*)
write(52,*) ”Initial guess of 4d points: Teff logg [Fe/H]
Xita ”
+,teff0,logg0,feh0,xita0
!ajusta el valor inicial de los datos downhill
!comenzamos por los 5 vertices del simplex inicial
!1er vertice
p0(1)=teff0
p0(2)=logg0
p0(3)=feh0
p0(4)=xita0
!2do vertice
p1(1)=p0(1)+100.
p1(2)=p0(2)
p1(3)=p0(3)
p1(4)=p0(4)
!3er vertice
p2(1)=p0(1)
p2(2)=p0(2)+0.2
p2(3)=p0(3)
p2(4)=p0(4)
p(3,1)=p2(1)
p(3,2)=p2(2)
p(3,3)=p2(3)
p(3,4)=p2(4)
p(4,1)=p3(1)
p(4,2)=p3(2)
p(4,3)=p3(3)
p(4,4)=p3(4)
p(5,1)=p4(1)
p(5,2)=p4(2)
p(5,3)=p4(3)
p(5,4)=p4(4)
write(*,*) ”Initial matrix p”
write(52,*) ”Initial matrix p”
!obtener el valor de la funcion, valuada en los 5 vertices
del simplex
!write(*,*)
write(*,*) . O btaining chicuad values for 5 Simplex vertices...”
write(52,*) . O btaining chicuad values for 5 Simplex vertices...”
do i=1,4
zx(i)=p0(i)
enddo
y(1)=funk(zx)
do i=1,4
zx(i)=p1(i)