Abrir - RDU

Apéndices 298
enddo
msj=0
do i=1,np1-1
if (index(a20(i),”nan”).gt.0) then
spec3(i)=1
if (msj.ne.1) then
msj=1
write(*,*) ”munari problems at lambda ”,lam13(i)
write(52,*) ”munari problems at lambda ”,lam13(i)
endif
cycle
endif
read(a20(i),*) spec3(i)
enddo
close(46)
call comspec1(spec2,spec3,ndesder13,nhastar13,chicuad)
if (t1.ne.oldt1) then
write(*,270) t1,l1,f1,x1,vrot1,-¿”,chicuad
oldt1=t1
endif
270format(f7.0,3f5.1,f8.1,a5,f10.1)
if (chicuad.lt.min) then
min=chicuad
mint=t1
minl=l1
minf=f1
minx=x1
endif
220enddo
enddo
enddo
!enddo
write(*,*) ”minimo chicuad ”,mint,minl,minf,minx,-¿”,min
write(52,*) ”minimo chicuad ”,mint,minl,minf,minx,-¿”,min character *51 fs1121,fs2121,fs1221,fs2221
end
character *51 fs1112,fs2112,fs1212,fs2212
character *51 fs1122,fs2122,fs1222,fs2222
!***********************************************
subroutine comspec1(cspec1,cspec2,ndesde,nhasta,chicuad)
!compara 2 espectros y calcula chicuad entre ellos
real*8 cspec1(2e5),cspec2(2e5)
real*8 chicuad,int2,int3
integer ndesde,nhasta
chicuad=0
do i=ndesde,nhasta
int2=cspec1(i)
int3=cspec2(i)
if ((int2.eq.9999).or.(int3.eq.9999).or.(int2.gt.2).or.
(int3.gt.2))then
int2=0
int3=0
endif
chicuad=chicuad+(int3-int2)**2
enddo
end
!*********************************************
!funcion que calcula chicuad para un espectro con cierta
Teff,logg,Feh,xita
!tiene que hacer dos cosas: 1ero, generar un espectro virtual
para los
!valores dados de los 4 parametros, y segundo comparar
con el esp.observado
!para obtener chicuad
real*8 function funk(zx)
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 teff,logg,feh,xita
real teff1,logg1,feh1,xita1
real teff2,logg2,feh2,xita2
character *31 f1111,f2111,f1211,f2211 !nombre archiv.espect.
a interpolar
character *31 f1121,f2121,f1221,f2221
character *31 f1112,f2112,f1212,f2212
character *31 f1122,f2122,f1222,f2222
character *51 fs1111,fs2111,fs1211,fs2211 !nombre archivos
path incluido
character *8 dir1111,dir2111,dir1211,dir2211 !nombre de
los directorios
character *8 dir1121,dir2121,dir1221,dir2221
character *8 dir1112,dir2112,dir1212,dir2212
character *8 dir1122,dir2122,dir1222,dir2222
real*8 zx(4),p
real*8 lam,lon,flux(16),x(4),vert(8),speckur(2e5)
integer ilam,pos,found,start,exceed
character a200*200,path*12,fout1*80,fout2*80
real*8 chicuad,desde,hasta,delta
real*8 vrot1,vrot2
!real vrot
character alpha*1,dir*8
!vrot=0.
alpha=”S”
path=/dos/munari/”