Abrir - RDU

Apéndices 299
teff =zx(1)
logg =zx(2)
feh =zx(3)
xita =zx(4)
x(1)=teff
x(2)=logg
x(3)=feh
x(4)=xita
if ((teff.gt.40000.).or.(teff.lt.3500.)) then
write(*,*) ”Teff ”,teff,”exceeding grid limits”
write(52,*) ”Teff ”,teff,”exceeding grid limits”
chicuad=1e7
goto 125
endif
if ((logg.gt.5.).or.(logg.lt.0.)) then
write(*,*) ”logg ”,logg,”exceeding grid limits”
write(52,*) ”logg ”,logg,”exceeding grid limits”
chicuad=1e7
goto 125
endif
if ((feh.gt.0.5).or.(feh.lt.-2.5)) then
write(*,*) ”feh ”,feh,”exceeding grid limits”
write(52,*) ”feh ”,feh,”exceeding grid limits”
chicuad=1e7
goto 125
endif
if ((xita.gt.4.0).or.(xita.lt.1.0)) then
write(*,*) ”Xita ”,xita,”exceeding grid limits”
write(52,*) ”Xita ”,xita,”exceeding grid limits”
chicuad=1e7
goto 125
endif
!valores de los 4 parametros mas cercanos en la grilla
!i.e. los que usara para realizar la interpolacion
call teffx(teff,teff1,teff2)
do r=0.0,4.5,0.5 !para la gravedad
if ((logg.ge.r).and.(logg.le.r+0.5)) then
logg1=r
logg2=r+0.5
exit
endif
enddo
call loggx(logg,logg1,logg2,teff,exceed)
if (exceed.eq.1) then
write(*,*) ”Logg ”,logg,”exceeding grid limits”
write(52,*) ”Logg ”,logg,”exceeding grid limits”
chicuad=1e7
goto 125
endif
do r=-2.5,0.0,0.5 !para la metalicidad
if ((feh.ge.r).and.(feh.le.r+0.5)) then
feh1=r
feh2=r+0.5
exit
endif
enddo
if ((xita.ge.1.0).and.(xita.lt.2.0)) then!para la veloc.microturb.
xita1=1.0
xita2=2.0
endif
if ((xita.ge.2.0).and.(xita.le.4.0)) then
xita1=2.0
xita2=4.0
endif
xita1=2
xita2=2
vert(1)=teff2
vert(2)=teff1
vert(3)=logg2
vert(4)=logg1
vert(5)=feh2
vert(6)=feh1
vert(7)=xita2
vert(8)=xita1
call vrotx(teff1,vrot,vrot1)
call vrotx(teff2,vrot,vrot2)
call concat1(teff1,logg1,feh1,vrot1,xita1,alpha,f1111,dir1111)
!formato del tipo T09740L25Fp05X02
call concat1(teff2,logg1,feh1,vrot2,xita1,alpha,f2111,dir2111)
call concat1(teff1,logg2,feh1,vrot1,xita1,alpha,f1211,dir1211)
call concat1(teff2,logg2,feh1,vrot2,xita1,alpha,f2211,dir2211)
call concat1(teff1,logg1,feh2,vrot1,xita1,alpha,f1121,dir1121)
call concat1(teff2,logg1,feh2,vrot2,xita1,alpha,f2121,dir2121)
call concat1(teff1,logg2,feh2,vrot1,xita1,alpha,f1221,dir1221)
call concat1(teff2,logg2,feh2,vrot2,xita1,alpha,f2221,dir2221)
call concat1(teff1,logg1,feh1,vrot1,xita2,alpha,f1112,dir1112)
call concat1(teff2,logg1,feh1,vrot2,xita2,alpha,f2112,dir2112)
call concat1(teff1,logg2,feh1,vrot1,xita2,alpha,f1212,dir1212)
call concat1(teff2,logg2,feh1,vrot2,xita2,alpha,f2212,dir2212)
call concat1(teff1,logg1,feh2,vrot1,xita2,alpha,f1122,dir1122)
call concat1(teff2,logg1,feh2,vrot2,xita2,alpha,f2122,dir2122)
call concat1(teff1,logg2,feh2,vrot1,xita2,alpha,f1222,dir1222)
call concat1(teff2,logg2,feh2,vrot2,xita2,alpha,f2222,dir2222)
pos=1
!abre los 16 archivos para interpolar
992 format(a12,a8,a31)
write(fs1111,992) path,dir1111,f1111
write(fs2111,992) path,dir2111,f2111
write(fs1211,992) path,dir1211,f1211