c****************************************************************** SUBROUTINE corr_bol_flower(teff,dteff,bc,d_bc,bmv,d_bmv) c subroutine PUBLIC du module mod_photo c determination de Bc (corr. bol) et B-V en fonction de Teff c suivant la table 3 de Flower 1996, ApJ 469, 355 c sorties c bmv, dbmv: B-V, erreur c bc, d_bc : corr. bol, erreur c entrees c teff, dteff: Teff, erreur c Auteur: P. Morel, Laboratoire Lagrange, O.C.A., CESAM2k c----------------------------------------------------------------------- USE mod_kind USE mod_numerique, ONLY : bsp1dn IMPLICIT NONE c m: ordre des splines INTEGER, PARAMETER :: pt=216, m=4 REAL (kind=dp), intent(in):: teff, dteff REAL (kind=dp), intent(out):: bc, d_bc, bmv, d_bmv REAL (kind=dp), PARAMETER :: dmb_sol=-0.01d0 REAL (kind=dp), SAVE, DIMENSION(2,pt) :: b REAL (kind=dp), SAVE, DIMENSION(pt+m) :: tt REAL (kind=dp), SAVE, DIMENSION(pt) :: t REAL (kind=dp), DIMENSION(2) :: fx, dfxdx INTEGER, SAVE :: l, knot c integer :: i LOGICAL, SAVE :: init=.TRUE. c----------------------------------------------------------------------- IF(init)THEN init=.FALSE. c initialisation des corrections bolometriques Bc, on utilise t comme Vt t= (/ 1-5.565d0,-5.004d0,-4.544d0,-4.152d0,-3.805d0,-3.502d0,-3.239d0,-3.007d0,-2.803d0,-2.620d0, 2-2.460d0,-2.318d0,-2.191d0,-2.078d0,-1.978d0,-1.885d0,-1.802d0,-1.728d0,-1.658d0,-1.595d0, 3-1.539d0,-1.486d0,-1.497d0,-1.393d0,-1.350d0,-1.312d0,-1.276d0,-1.241d0,-1.209d0,-1.178d0, 4-1.150d0,-1.122d0,-1.096d0,-1.071d0,-1.047d0,-1.024d0,-1.002d0,-0.980d0,-0.960d0,-0.939d0, 5-0.920d0,-0.901d0,-0.883d0,-0.865d0,-0.848d0,-0.831d0,-0.814d0,-0.798d0,-0.782d0,-0.766d0, 6-0.752d0,-0.736d0,-0.722d0,-0.707d0,-0.694d0,-0.679d0,-0.666d0,-0.652d0,-0.640d0,-0.626d0, 7-0.614d0,-0.601d0,-0.588d0,-0.576d0,-0.564d0,-0.552d0,-0.540d0,-0.528d0,-0.517d0,-0.505d0, 8-0.494d0,-0.482d0,-0.471d0,-0.459d0,-0.448d0,-0.437d0,-0.426d0,-0.414d0,-0.403d0,-0.392d0, 9-0.382d0,-0.372d0,-0.361d0,-0.352d0,-0.342d0,-0.332d0,-0.322d0,-0.313d0,-0.304d0,-0.295d0, 1-0.285d0,-0.277d0,-0.268d0,-0.259d0,-0.250d0,-0.242d0,-0.233d0,-0.225d0,-0.216d0,-0.208d0, 2-0.200d0,-0.192d0,-0.184d0,-0.176d0,-0.169d0,-0.161d0,-0.153d0,-0.146d0,-0.139d0,-0.132d0, 3-0.124d0,-0.117d0,-0.111d0,-0.104d0,-0.098d0,-0.091d0,-0.085d0,-0.079d0,-0.073d0,-0.067d0, 4-0.061d0,-0.055d0,-0.050d0,-0.045d0,-0.039d0,-0.035d0,-0.030d0,-0.025d0,-0.021d0,-0.016d0, 5-0.012d0,-0.008d0,-0.004d0,-0.001d0, 0.003d0, 0.006d0, 0.009d0, 0.012d0, 0.015d0, 0.018d0, 6 0.020d0, 0.022d0, 0.025d0, 0.026d0, 0.028d0, 0.030d0, 0.031d0, 0.032d0, 0.033d0, 0.034d0, 7 0.034d0, 0.035d0, 0.035d0, 0.035d0, 0.035d0, 0.035d0, 0.034d0, 0.034d0, 0.033d0, 0.033d0, 8 0.032d0, 0.031d0, 0.029d0, 0.028d0, 0.026d0, 0.024d0, 0.022d0, 0.019d0, 0.015d0, 0.010d0, 9 0.004d0,-0.004d0,-0.013d0,-0.024d0,-0.036d0,-0.050d0,-0.067d0,-0.085d0,-0.106d0,-0.129d0, 1-0.155d0,-0.184d0,-0.216d0,-0.252d0,-0.290d0,-0.334d0,-0.381d0,-0.432d0,-0.488d0,-0.549d0, 2-0.614d0,-0.684d0,-0.759d0,-0.839d0,-0.923d0,-1.013d0,-1.107d0,-1.205d0,-1.307d0,-1.414d0, 3-1.525d0,-1.641d0,-1.762d0,-1.891d0,-2.028d0,-2.177d0,-2.341d0,-2.523d0,-2.730d0,-2.966d0, 4-3.234d0,-3.534d0,-3.861d0,-4.197d0,-4.506d0,-4.720d0 /) c dmb_sol : correction de façon à ce que la correction bolometrique du c soleil soit -0.07 ce qui lui donne une magnitude absolue de 4.75 c avec comme zéro des magnitude absolues L=3.055W (R. Cayrel UAI 98) b(1,:)=t-dmb_sol !Bc c initialisation des B-V b(2,:)=(/ (1.8d0-(l-1)*0.01d0,l=1,pt) /) c initialisation des Teff t= (/ 1 2936.d0, 3000.d0, 3061.d0, 3119.d0, 3175.d0, 3228.d0, 3278.d0, 3326.d0, 3372.d0, 3415.d0, 2 3457.d0, 3496.d0, 3533.d0, 3568.d0, 3601.d0, 3634.d0, 3664.d0, 3693.d0, 3721.d0, 3748.d0, 3 3773.d0, 3797.d0, 3821.d0, 3843.d0, 3865.d0, 3885.d0, 3906.d0, 3926.d0, 3945.d0, 3964.d0, 4 3982.d0, 4001.d0, 4018.d0, 4036.d0, 4053.d0, 4070.d0, 4086.d0, 4103.d0, 4120.d0, 4137.d0, 5 4153.d0, 4169.d0, 4186.d0, 4202.d0, 4217.d0, 4234.d0, 4251.d0, 4266.d0, 4283.d0, 4300.d0, 6 4316.d0, 4333.d0, 4349.d0, 4366.d0, 4382.d0, 4399.d0, 4415.d0, 4433.d0, 4449.d0, 4466.d0, 7 4483.d0, 4500.d0, 4518.d0, 4535.d0, 4553.d0, 4569.d0, 4587.d0, 4605.d0, 4622.d0, 4640.d0, 8 4658.d0, 4676.d0, 4694.d0, 4713.d0, 4731.d0, 4748.d0, 4767.d0, 4786.d0, 4806.d0, 4825.d0, 9 4843.d0, 4862.d0, 4883.d0, 4902.d0, 4922.d0, 4943.d0, 4963.d0, 4984.d0, 5004.d0, 5025.d0, 1 5047.d0, 5068.d0, 5090.d0, 5114.d0, 5136.d0, 5159.d0, 5183.d0, 5207.d0, 5231.d0, 5256.d0, 2 5282.d0, 5307.d0, 5333.d0, 5359.d0, 5386.d0, 5413.d0, 5442.d0, 5470.d0, 5499.d0, 5528.d0, 3 5559.d0, 5589.d0, 5620.d0, 5653.d0, 5684.d0, 5717.d0, 5751.d0, 5784.d0, 5819.d0, 5855.d0, 4 5891.d0, 5927.d0, 5964.d0, 6002.d0, 6040.d0, 6078.d0, 6117.d0, 6158.d0, 6198.d0, 6240.d0, 5 6282.d0, 6324.d0, 6366.d0, 6409.d0, 6453.d0, 6496.d0, 6541.d0, 6587.d0, 6632.d0, 6678.d0, 6 6725.d0, 6771.d0, 6820.d0, 6867.d0, 6916.d0, 6964.d0, 7014.d0, 7063.d0, 7113.d0, 7164.d0, 7 7216.d0, 7267.d0, 7319.d0, 7372.d0, 7426.d0, 7481.d0, 7535.d0, 7592.d0, 7648.d0, 7707.d0, 8 7766.d0, 7825.d0, 7886.d0, 7950.d0, 8016.d0, 8083.d0, 8152.d0, 8222.d0, 8296.d0, 8373.d0, 9 8454.d0, 8538.d0, 8625.d0, 8717.d0, 8814.d0, 8916.d0, 9026.d0, 9139.d0, 9262.d0, 9392.d0, 1 9530.d0, 9680.d0, 9840.d0,10011.d0,10195.d0,10396.d0,10612.d0,10846.d0,11099.d0,11376.d0, 211678.d0,12008.d0,12368.d0,12764.d0,13197.d0,13677.d0,14203.d0,14787.d0,15434.d0,16154.d0, 316958.d0,17852.d0,18858.d0,19989.d0,21261.d0,22703.d0,24339.d0,26199.d0,28333.d0,30789.d0, 433620.d0,36897.d0,40719.d0,45196.d0,50477.d0,56728.d0 /) CALL bsp1dn(2,b,t,tt,pt,m,knot,.FALSE.,t(1),l,fx,dfxdx) PRINT* ; WRITE(*,1) 1 FORMAT('Bc (corr. bol) et B-V en fonction de Teff',/, 1 'suivant la table 3 de Flower 1996, ApJ 469, 355') PRINT* ENDIF CALL bsp1dn(2,b,t,tt,pt,m,knot,.TRUE.,teff,l,fx,dfxdx) bc=fx(1) ; d_bc=dfxdx(1)*dteff ; bmv=fx(2) ; d_bmv=dfxdx(2)*dteff RETURN END SUBROUTINE corr_bol_flower