Abrir - RDU

Apéndices 300
write(fs2211,992) path,dir2211,f2211
write(fs1121,992) path,dir1121,f1121
write(fs2121,992) path,dir2121,f2121
write(fs1221,992) path,dir1221,f1221
write(fs2221,992) path,dir2221,f2221
write(fs1112,992) path,dir1112,f1112
write(fs2112,992) path,dir2112,f2112
write(fs1212,992) path,dir1212,f1212
write(fs2212,992) path,dir2212,f2212
write(fs1122,992) path,dir1122,f1122
write(fs2122,992) path,dir2122,f2122
write(fs1222,992) path,dir1222,f1222
write(fs2222,992) path,dir2222,f2222
open(30,file=fs1111,status=. o ld”)
open(31,file=fs2111,status=. o ld”)
open(32,file=fs1211,status=. o ld”)
open(33,file=fs2211,status=. o ld”)
open(34,file=fs1121,status=. o ld”)
open(35,file=fs2121,status=. o ld”)
open(36,file=fs1221,status=. o ld”)
open(37,file=fs2221,status=. o ld”)
open(38,file=fs1112,status=. o ld”)
open(39,file=fs2112,status=. o ld”)
open(40,file=fs1212,status=. o ld”)
open(41,file=fs2212,status=. o ld”)
open(42,file=fs1122,status=. o ld”)
open(43,file=fs2122,status=. o ld”)
open(44,file=fs1222,status=. o ld”)
open(45,file=fs2222,status=. o ld”)
do !j=1,np1
do i=1,16
read(29+i,*,end=1010) flux(i)
enddo
!tenemos la funcion en cada uno de los 16 puntos
!ahora la interpolacion 4d
call interp4d(x,vert,flux,p) !p es el valor interpolado
!ahora deberiamos guardar el par lon,p en alguna variable...
speckur(pos)=p
pos=pos+1
enddo
1010 close(30)
close(31)
close(32)
close(33)
close(34)
close(35)
close(36)
close(37)
close(38)
close(39)
close(40)
close(41)
close(42)
close(43)
close(44)
close(45)
!tenemos en speckur(1:pos) el espectro interpolado
if (wspec.eq.1) then
write(*,*) ”Writting (minimmum-chicuad) fitted spectra
to disk”
write(52,*) ”Writting (minimmum-chicuad) fitted spectra
to disk”
fout1=”dataout/min.chicuad.//file
fout2=”dataout/resampled.//file
open(48,file=fout1)
open(51,file=fout2)
write(48,*) teff,logg,feh,xita,vrot
write(48,*)
do i=1,pos
write(48,*) lam13(i),speckur(i),speckur(i)-spec2(i) !lam, flux,
difference with observed
enddo
do i=ndesder13, nhastar13
write(51,*) lam13(i),spec2(i)!lam, flux
enddo
close(48)
close(51)
endif
!speckur(1:points) tiene el mismo formato lamda que specobs
!entonces podemos comparar ambos espectros
call comspec1(spec2,speckur,ndesder13,nhastar13,chicuad)
125 funk=chicuad
write(*,*) ”funk called with ”,teff,logg,feh,xita,-¿”,chicuad
write(52,*) ”funk called with ”,teff,logg,feh,xita,-¿”,chicuad
end
!******************************************************
subroutine resample(desde,hasta,spec1,spec4,nd13,nh13)
!lleva al formato de la grilla lambda de munari
!entrada=desde,hasta,spec1
!deben ser conocidas (en common)= lam13, np1
!salida=spec4,nd13,nh13
real found
real*8 desde,hasta,spec1(2e5,2),spec4(2e5)
real*8 x2,x1,y2,y1,x,px
integer nx1,nx2,nd13,nh13
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/
!busca el 1er pto en R13 q es anterior al desde: ”desder13”