Abrir - RDU

Apéndices 297
enddo
y(2)=funk(zx)
do i=1,4
zx(i)=p2(i)
enddo
y(3)=funk(zx)
do i=1,4
zx(i)=p3(i)
enddo
y(4)=funk(zx)
do i=1,4
zx(i)=p4(i)
enddo
y(5)=funk(zx)
write(*,*) ”Starting chicuad values =”
write(*,*) y(1),y(2),y(3),y(4),y(5)
write(52,*) ”Starting chicuad values =”
write(52,*) y(1),y(2),y(3),y(4),y(5)
!llamada a la subrutina downhill para obtener el minimo
de chicuad...
write(*,*)
write(*,*) ”Entering downhill iteration...”
write(52,*)
write(52,*) ”Entering downhill iteration...”
ftol=0.0005
call amoeba(p,y,5,4,4,ftol,funk,iter)
write(*,*)
write(*,*) ”Downhill convergence before ”,iter,”steps.”
write(52,*) ”Final solution: Teff logg [Fe/H] xita (vrot
vrad)”
write(52,*) p(1,1),p(1,2),p(1,3),p(1,4),”(”,vrot,vrad,”)”
write(52,*)
write(52,*) ”Downhill convergence before ”,iter,”steps.”
write(52,*) ”Final solution: Teff logg [Fe/H] xita (vrot
vrad)”
write(52,*) p(1,1),p(1,2),p(1,3),p(1,4),”(”,vrot,vrad,”)”
write(50,123) file,p(1,1),p(1,2),p(1,3),p(1,4),
+vrot,vrad,”(”,teff0,logg0,feh0,xita0,”)”
123 format(a30,6f10.2,2x,a1,4f10.2,a1)
do i=1,4 !vuelve a llamar funk p escrbir el arch.salida
zx(i)=p(1,i)
enddo
wspec=1
chicuad=funk(zx) !...aqui es dde lo escribe
end
!***********************************************
subroutine force1(mint,minl,minf,minx)
real t1,l1,f1,x1
real ex,oldt1
real min,mint,minl,minf,minx
real*8 spec3(2e5),chicuad,vrot1
character fname1*31,fname2*51,path*12,dir*8
character*31 munari(68200),a20(2e5)*20
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/
!lee nombre archivos de la grilla
open(3,file=”munari.filenames.txt”,status=. o ld”)
i=1
do
read(3,’(a31)’,end=240) munari(i)
i=i+1
enddo
240close(3)
!el bucle principal
write(*,*) ”Writting some examples of grid chicuad-values...”
write(*,*) ”teff logg feh xita vsini –¿chicuad”
write(52,*) ”Writting some examples of grid chicuad-values...”
write(52,*) ”teff logg feh xita vsini –¿chicuad”
min=1e7
oldt1=0
do t1=3500,40000,250
do l1=0,5,0.5
do f1=-2.5,0.5, 0.5
x1=2.
call vrotx(t1,vrot,vrot1)
call concat1(t1,l1,f1,vrot1,x1,”S”,fname1,dir)
path=/dos/munari/”
fname2=path//dir//fname1
!se fija si existe el nombre del archivo
ex=0
do i=1,68200
if (fname1.eq.munari(i)) then
ex=1
exit
endif
enddo
if (ex.eq.0) then
cycle
endif
open(46,file=fname2,status=. o ld”)
do i=1,np1-1
read(46,’(a20)’) a20(i)