Module datagraph_XXX ! USE MTD_STRUCT ! Implicit None ! Character(len = 2), parameter :: Symb = 'XX' Integer, parameter :: Z_IONI = 0 Character(len = 4) :: cpage = 'DA4L' Integer, parameter :: lenterm = 1 Character(len = 32) :: continuum = '0.00' ! eV Real, dimension(lenterm) :: SLP = (/100/) Real, dimension(lenterm) :: iterm = (/ 1/) character(len = 7), dimension(lenterm) :: cterm = (/"[1$P[o$"/) ! Position et taille du graphique Integer, parameter :: posx = 250, posy = 1980 Integer, parameter :: lenx = 2600, leny = 1800 ! Abscisses du continuum Real, dimension(1), parameter :: abscontinuum = (/ 1 /) ! Noms des niveaux Integer, parameter :: levelsize = 1 character(len = 5),dimension(levelsize) :: levelname = (/'4s'/) ! Position des niveaux Real, dimension(levelsize, 2) :: levelpos = reshape(source=(/ 0.4, 0.05 /), shape=(/levelsize, 2/), order=(/2, 1/)) ! ! SUPPLEMENTARY LEVELS ! Integer, parameter :: Klevels = 1 Type(LEVELS), dimension(Klevels), save :: S_level_sup=(/LEVELS('spdfghiklm', -1, 1.797D+308, 1.797D+308, 1.797D+308, 9., -1)/) ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! Integer, parameter :: App1_N = 1 Real, dimension(App1_N), parameter :: App1_lambda = (/ 4226.73 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ 150 /) ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia Integer, dimension(App1_N) :: GNi = (/ 1 /) Integer, dimension(App1_N) :: GNj = (/ 3 /) ! End Module datagraph_XXX ! !*********************************************************************** Module datagraph_BoI ! USE MTD_STRUCT ! Implicit None ! Character(len = 2), parameter :: Symb = 'Bo' Integer, parameter :: Z_IONI = 1 Character(len = 4) :: cpage = 'DA4L' Integer, parameter :: lenterm = 4 Character(len = 32) :: continuum = '5.14' ! eV Real, dimension(lenterm) :: SLP = (/200,211,220,231/) Real, dimension(lenterm) :: iterm = (/1,2,3,4/) character(len = 7), dimension(lenterm) :: cterm = (/"[2$S ","[2$P[o$","[2$D ","[2$F[o$"/) ! Position et taille du graphique Integer, parameter :: posx = 250, posy = 1980 Integer, parameter :: lenx = 2600, leny = 1800 ! Abscisses du continuum Real, dimension(1), parameter :: abscontinuum = (/ 1 /) ! Noms des niveaux Integer, parameter :: levelsize = 8 character(len = 3),dimension(levelsize) :: levelname = (/'3s ','11s','3p ','10p','3d ','10d','4f ','10f'/) ! Position des niveaux Real, dimension(levelsize, 2) :: levelpos = reshape(source=(/ 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, & 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.05,0.4, 0.05/), & shape=(/levelsize, 2/), & order=(/2, 1/) ) ! ! SUPPLEMENTARY LEVELS ! Integer, parameter :: Klevels = 1 Type(LEVELS), dimension(Klevels), save :: S_level_sup=(/LEVELS('spdfghiklm', -1, 1.797D+308, 1.797D+308, 1.797D+308, 9., -1)/) ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! Integer, parameter :: App1_N = 1 Real, dimension(App1_N), parameter :: App1_lambda = (/ 0.0 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ 0.0 /) ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia Integer, dimension(App1_N) :: GNi = (/ 1 /) Integer, dimension(App1_N) :: GNj = (/ 3 /) End Module datagraph_BoI ! !*********************************************************************** Module datagraph_AlI ! USE MTD_STRUCT ! Implicit None ! Character(len = 2), parameter :: Symb = 'Al' Integer, parameter :: Z_IONI = 1 Character(len = 4) :: cpage = 'DA4L' Integer, parameter :: lenterm = 4 Character(len = 32) :: continuum = '5.14' ! eV Real, dimension(lenterm) :: SLP = (/200,211,220,231/) Real, dimension(lenterm) :: iterm = (/1,2,3,4/) character(len = 7), dimension(lenterm) :: cterm = (/"[2$S ","[2$P[o$","[2$D ","[2$F[o$"/) ! Position et taille du graphique Integer, parameter :: posx = 250, posy = 1980 Integer, parameter :: lenx = 2600, leny = 1800 ! Abscisses du continuum Real, dimension(1), parameter :: abscontinuum = (/ 1 /) ! Noms des niveaux Integer, parameter :: levelsize = 8 character(len = 3),dimension(levelsize) :: levelname = (/'3s ','11s','3p ','10p','3d ','10d','4f ','10f'/) ! Position des niveaux Real, dimension(levelsize, 2) :: levelpos = reshape(source=(/ 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.4, & 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.05,0.4, 0.05/), & shape=(/levelsize, 2/), & order=(/2, 1/) ) ! ! SUPPLEMENTARY LEVELS ! Integer, parameter :: Klevels = 1 Type(LEVELS), dimension(Klevels), save :: S_level_sup=(/LEVELS('spdfghiklm', -1, 1.797D+308, 1.797D+308, 1.797D+308, 9., -1)/) ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! Integer, parameter :: App1_N = 1 Real, dimension(App1_N), parameter :: App1_lambda = (/ 0.0 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ 0.0 /) ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia Integer, dimension(App1_N) :: GNi = (/ 1 /) Integer, dimension(App1_N) :: GNj = (/ 3 /) End Module datagraph_AlI ! !*********************************************************************** ! Module datagraph_NaI ! USE MTD_STRUCT ! Implicit None ! Character(len = 2), parameter :: Symb = 'Na' Integer, parameter :: Z_IONI = 1 Character(len = 4) :: cpage = 'DA4L' Integer, parameter :: lenterm = 6 Character(len = 32) :: continuum = '5.14' ! eV Real, dimension(lenterm) :: SLP = (/200,211,220,231,240,251/) Real, dimension(lenterm) :: iterm = (/1,2,3,4,5,6/) character(len = 7), dimension(lenterm) :: cterm = (/"[2$S ","[2$P[o$","[2$D ","[2$F[o$", "[2$G ", "[2$H[o$"/) ! Position et taille du graphique Integer, parameter :: posx = 250, posy = 1980 Integer, parameter :: lenx = 2600, leny = 1800 ! Abscisses du continuum Real, dimension(1), parameter :: abscontinuum = (/ 1 /) ! Noms des niveaux Integer, parameter :: levelsize = 8 character(len = 3),dimension(levelsize) :: levelname = (/'3s ','11s','3p ','10p','3d ','10d','4f ','10f'/) ! Position des niveaux Real, dimension(levelsize, 2) :: levelpos = reshape(source=(/ 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.4, & 0.05, 0.4, 0.05, 0.4, 0.05, 0.4, 0.05,0.4, 0.05/), & shape=(/levelsize, 2/), & order=(/2, 1/) ) ! ! SUPPLEMENTARY LEVELS ! Integer, parameter :: Klevels = 1 Type(LEVELS), dimension(Klevels), save :: S_level_sup=(/LEVELS('spdfghiklm', -1, 1.797D+308, 1.797D+308, 1.797D+308, 9., -1)/) ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! Integer, parameter :: App1_N = 1 Real, dimension(App1_N), parameter :: App1_lambda = (/ 0.0 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ 0.0 /) ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia Integer, dimension(App1_N) :: GNi = (/ 1 /) Integer, dimension(App1_N) :: GNj = (/ 3 /) End Module datagraph_NaI ! !*********************************************************************** ! Module datagraph_MgI ! USE MTD_STRUCT Implicit None ! Character(len = 2), parameter :: Symb = 'Mg' Integer, parameter :: Z_IONI = 1 Character(len = 4) :: cpage = 'DA4L' Integer,parameter :: lenterm = 16 Character(len = 32) :: continuum = '7.65' Real, dimension(lenterm) :: SLP = (/100,111,120,131,140,300,310,311,320,331,340,351,360,371,380,391/) Real, dimension(lenterm) :: iterm = (/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16/) Character(len = 10), dimension(lenterm) :: cterm = & (/"[1$S ","[1$P[o$","[1$D ","[1$F[o$","[1$G ",& ! Système de multiplicité 1 "[3$S ","[3$P ","[3$P[o$","[3$D ","[3$F[o$","[3$G ",& ! Système de multiplicité 3 "H[o$ ","I ","K[o$ ","L ","M[o$ "/) ! Système de multiplicité dégénéré ! Position et taille du graphique Integer, parameter :: posx = 250, posy = 1980 Integer, parameter :: lenx = 2600, leny = 1800 ! Abscisses du continuum Real, dimension(18), parameter :: abscontinuum = (/0.3,1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,16.5/) ! Noms des niveaux Integer,parameter :: levelsize = 32 Character(len= 5),dimension(levelsize):: levelname = & (/'3s ','4s ','10s ','3p ','10p ','3d ','10d ','4f ','10f ','5g ','10g ',& '4s ','10s ','3p[2$','3p ','4p ','10p ','3d ','10d ','4f ','10f ','5g ','10g ',& '6h ','10h ','7i ','10i ','8k ','10k ','9l ','10l ','10m '/) ! Position des niveaux Real,dimension(levelsize,2) :: levelpos = reshape(source=& (/ 0.5,0.1, 0.5,5.5, 0.42,7.6, 1.5,4.5, 1.4,7.6, 2.4,5.9, 2.4,7.6, 3.5,6.95, 3.4,7.6, 4.5,7.18, 4.4,7.6,& 5.5,5.2, 5.4,7.6, 6.45,7.3, 7.5,2.85, 7.5,6.05, 7.4,7.6, 8.5,6.05, 8.4,7.6, 9.5,6.9, 9.4,7.6, 10.5,7.2, 10.4,7.6,& 11.5,7.3, 11.4,7.6, 12.5,7.4, 12.4,7.6, 13.5,7.45, 13.4,7.6, 14.58,7.47, 14.48,7.63, 15.3,7.63/),& shape=(/levelsize, 2/),& order=(/2, 1/)) ! ! SUPPLEMENTARY LEVELS ! Integer, parameter :: Klevels = 31 Type(LEVELS), dimension(Klevels), parameter :: S_level_sup = (/ & ! ! Ajout de niveaux supplémentaires pour être consistent ! LEVELS('3s.10p',311,60303.620_8,5.4952702E-01,7.4766955E+00,1.0,63,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3s.10p',311,60303.620_8,5.4952702E-01,7.4766955E+00,3.0,63,"STOPB"),& ! avec levée de dégénérescence LEVELS('3s.10p',311,60303.620_8,5.4952702E-01,7.4766955E+00,5.0,63,"STOPB"),& ! avec levée de dégénérescence LEVELS('3s.9g',140,60312.317_8,5.4960627E-01,7.4777738E+00,9.0,67,"SNIST"),& ! En(NIST) 3s.9g dégénéré LEVELS('3s.9g',340,60312.317_8,5.4960627E-01,7.4777738E+00,7.0,68,"SNIST"),& ! En(NIST) 3s.9g dégénéré LEVELS('3s.9g',340,60312.317_8,5.4960627E-01,7.4777738E+00,9.0,68,"SNIST"),& ! En(NIST) 3s.9g dégénéré LEVELS('3s.9g',340,60312.317_8,5.4960627E-01,7.4777738E+00,11.0,68,"SNIST"),& ! En(NIST) 3s.9g dégénéré LEVELS('3s.10g',140,60572.020_8,5.5197286E-01,7.5099729E+00,9.0,77,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3s.10g',340,60572.020_8,5.5197286E-01,7.5099729E+00,7.0,78,"STOPB"),& ! avec levée de dégénérescence LEVELS('3s.10g',340,60572.020_8,5.5197286E-01,7.5099729E+00,9.0,78,"STOPB"),& ! avec levée de dégénérescence LEVELS('3s.10g',340,60572.020_8,5.5197286E-01,7.5099729E+00,11.0,78,"STOPB"),& ! avec levée de dégénérescence LEVELS('3s.5g',140,57262.760_8,5.2181666E-01,7.0998173E+00,9.0,27,"NIST"),& ! LEVELS('3s.5g',340,57262.760_8,5.2181666E-01,7.0998173E+00,7.0,26,"NIST"),& ! NIST but put here for compatibility LEVELS('3s.5g',340,57262.760_8,5.2181666E-01,7.0998173E+00,9.0,26,"NIST"),& ! LEVELS('3s.5g',340,57262.760_8,5.2181666E-01,7.0998173E+00,11.0,26,"SNIST"),& ! Suppl perso avec En(NIST) LEVELS('3s.8g',340,59950.860_8,5.4631243E-01,7.4331069E+00,7.0,57,"SNIST"),& ! Suppl perso avec En(NIST) ! ! Ajout de niveaux supplémentaires présents dans NIST mais pas dans la TOPBASE ! LEVELS('3s.6h',351,58618.942_8,0.d0,7.2678221D+00,44.0,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3s.7h',351,59428.853_8,0.d0,7.3682382D+00,44.0,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3s.8h',351,59954.523_8,0.d0,7.4334130D+00,44.0,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3s.9h',351,60314.902_8,0.d0,7.4780943D+00,44.0,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3s.7i',360,59430.517_8,0.d0,7.3684445D+00,52.0,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3s.9i',360,60315.722_8,0.d0,7.4781960D+00,52.0,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3s.9k',371,60316.041_8,0.d0,7.4782355D+00,60.0,0,"SNIST"),& ! NIST sans TOPBASE ! ! Ajout de niveaux supplémentaires calculé par la théorie de la polarisation de Chang et Noyes (1983) ! LEVELS('3s.10h',351,60572.63_8,0.d0,7.5100485D+00,44.0,0,"CN83"),& LEVELS('3s.8i',360,59955.66_8,0.d0,7.4336780D+00,52.0,0,"CN83"),& LEVELS('3s.10i',360,60573.25_8,0.d0,7.5101253D+00,52.0,0,"CN83"),& LEVELS('3s.8k',371,59956.08_8,0.d0,7.4336060D+00,60.0,0,"CN83"),& LEVELS('3s.10k',371,60573.48_8,0.d0,7.5101539D+00,60.0,0,"CN83"), & LEVELS('3s.9l',380,60316.15_8,0.d0,7.4782490D+00,68.0,0,"CN83"),& LEVELS('3s.10l',380,60573.57_8,0.d0,7.5101650D+00,68.0,0,"CN83"),& LEVELS('3s.10m',391,60573.62_8,0.d0,7.5101712D+00,76.0,0,"CN83") /) ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 139 Type(LINES), dimension(Klines), parameter :: S_line_sup = (/ & LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.8g', 340, 59950.871_8, 4.0, 5274.1504_8, 1896.0400, -0.697, 4.140E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.8g', 340, 59950.871_8, 5.0, 5274.1221_8, 1896.0500, -0.582, 4.420E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.8g', 340, 59950.871_8, 4.0, 5274.0947_8, 1896.0599, -1.873, 2.760E+04, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.7g', 340, 59423.535_8, 3.0, 4746.9854_8, 2106.6001, -3.358, 9.410E+02, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.7g', 340, 59423.535_8, 3.0, 4746.9854_8, 2106.6001, -1.541, 6.180E+04, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.7g', 340, 59423.535_8, 4.0, 4746.9854_8, 2106.6001, -1.541, 4.810E+04, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.7g', 340, 59423.535_8, 4.0, 4746.9854_8, 2106.6001, -0.366, 7.190E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.652_8, 2.0, '3s.7g', 340, 59423.535_8, 3.0, 4746.8818_8, 2106.6460, -0.483, 7.050E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.7g', 340, 59423.535_8, 5.0, 4746.7832_8, 2106.6899, -0.251, 7.670E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.6g', 340, 58610.797_8, 3.0, 3934.2195_8, 2541.8000, -1.069, 1.260E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.6g', 340, 58610.797_8, 4.0, 3934.2195_8, 2541.8000, 0.105, 1.460E+06, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.652_8, 2.0, '3s.6g', 340, 58610.797_8, 3.0, 3934.1406_8, 2541.8511, -0.011, 1.440E+06, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.6g', 340, 58610.797_8, 3.0, 3934.0649_8, 2541.8999, -2.887, 1.910E+03, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.6g', 340, 58610.797_8, 4.0, 3934.0649_8, 2541.8999, -1.069, 9.780E+04, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.6g', 340, 58610.797_8, 5.0, 3934.0403_8, 2541.9160, 0.221, 1.560E+06, 0.00), & !+ NIST LINES('3s.5f', 131, 57204.164_8, 3.0, '3s.8g', 340, 59950.871_8, 4.0, 2746.7246_8, 3640.7000, -0.200, 3.530E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.266_8, 3.0, '3s.8g', 340, 59950.871_8, 4.0, 2746.6489_8, 3640.8000, -0.229, 3.300E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.8g', 340, 59950.871_8, 5.0, 2746.5737_8, 3640.8999, -0.113, 3.520E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.8g', 340, 59950.871_8, 4.0, 2746.5737_8, 3640.8999, -1.404, 2.200E+04, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.8f', 131, 59935.371_8, 3.0, 2672.6533_8, 3741.6001, -3.391, 2.760E+02, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.8f', 131, 59935.371_8, 3.0, 2672.6533_8, 3741.6001, -2.190, 4.400E+03, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.8f', 331, 59935.371_8, 4.0, 2672.6533_8, 3741.6001, -5.209, 3.270E+00, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.8f', 331, 59935.371_8, 4.0, 2672.6533_8, 3741.6001, -3.392, 2.150E+02, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.8f', 331, 59935.371_8, 2.0, 2672.6533_8, 3741.6001, -2.334, 4.420E+03, 0.00), & !+ NIST LINES('3s.5g', 140, 57262.762_8, 4.0, '3s.8f', 131, 59935.371_8, 3.0, 2672.6533_8, 3741.6001, -2.216, 4.130E+03, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.5g', 340, 57262.762_8, 4.0, 2586.1177_8, 3866.8000, 0.921, 4.130E+06, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.652_8, 2.0, '3s.5g', 340, 57262.762_8, 3.0, 2586.1057_8, 3866.8179, 0.803, 4.050E+06, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.699_8, 3.0, '3s.5g', 340, 57262.762_8, 3.0, 2586.0588_8, 3866.8879, -0.254, 3.550E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.5g', 340, 57262.762_8, 4.0, 2586.0508_8, 3866.8999, -0.254, 2.760E+05, 0.00), & !+ NIST LINES('3s.4f', 331, 54676.754_8, 4.0, '3s.5g', 340, 57262.762_8, 3.0, 2586.0046_8, 3866.9692, -2.071, 5.410E+03, 0.00), & !+ NIST LINES('3s.5f', 131, 57204.164_8, 3.0, '3s.7g', 340, 59423.535_8, 4.0, 2219.4109_8, 4505.7002, 0.209, 5.900E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.227_8, 2.0, '3s.7g', 340, 59423.535_8, 3.0, 2219.3125_8, 4505.8999, 0.063, 5.420E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.266_8, 3.0, '3s.7g', 340, 59423.535_8, 3.0, 2219.3125_8, 4505.8999, -0.994, 4.760E+04, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.266_8, 3.0, '3s.7g', 340, 59423.535_8, 4.0, 2219.3125_8, 4505.8999, 0.182, 5.540E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.7g', 340, 59423.535_8, 3.0, 2219.2632_8, 4506.0000, -2.813, 7.220E+02, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.7g', 340, 59423.535_8, 4.0, 2219.2632_8, 4506.0000, -0.993, 3.710E+04, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.7g', 340, 59423.535_8, 5.0, 2219.2632_8, 4506.0000, 0.297, 5.910E+05, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.7f', 331, 59400.762_8, 4.0, 2138.0313_8, 4677.2002, -2.925, 4.020E+02, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.7f', 131, 59400.762_8, 3.0, 2138.0313_8, 4677.2002, -2.927, 5.150E+02, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.7f', 331, 59400.762_8, 4.0, 2138.0313_8, 4677.2002, -4.743, 6.120E+00, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.7f', 131, 59400.762_8, 3.0, 2138.0313_8, 4677.2002, -1.724, 8.230E+03, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.7f', 331, 59400.762_8, 2.0, 2138.0313_8, 4677.2002, -1.869, 8.240E+03, 0.00), & !+ NIST LINES('3s.5g', 140, 57262.762_8, 4.0, '3s.7f', 131, 59400.762_8, 3.0, 2138.0313_8, 4677.2002, -1.751, 7.720E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.10f', 131, 60562.637_8, 3.0, 1951.8671_8, 5123.2998, -2.093, 2.930E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.10f', 331, 60562.637_8, 4.0, 1951.8671_8, 5123.2998, -5.114, 2.170E+00, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.10f', 331, 60562.637_8, 2.0, 1951.8671_8, 5123.2998, -2.238, 2.940E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.10f', 331, 60562.637_8, 4.0, 1951.8671_8, 5123.2998, -3.296, 1.430E+02, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.10f', 131, 60562.637_8, 3.0, 1951.8671_8, 5123.2998, -3.296, 1.840E+02, 0.00), & !+ NIST LINES('3s.6g', 140, 58610.797_8, 4.0, '3s.10f', 131, 60562.637_8, 3.0, 1951.8671_8, 5123.2998, -2.120, 2.750E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 5.0, '3s.10f', 331, 60562.637_8, 4.0, 1951.8671_8, 5123.2998, -2.005, 2.790E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.9f', 131, 60301.281_8, 3.0, 1690.5027_8, 5915.3999, -1.795, 4.360E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.9f', 331, 60301.281_8, 2.0, 1690.5027_8, 5915.3999, -1.912, 4.670E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.9f', 331, 60301.281_8, 4.0, 1690.5027_8, 5915.3999, -2.970, 2.270E+02, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.9f', 131, 60301.281_8, 3.0, 1690.5027_8, 5915.3999, -2.970, 2.920E+02, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.9f', 331, 60301.281_8, 4.0, 1690.5027_8, 5915.3999, -4.788, 3.450E+00, 0.00), & !+ NIST LINES('3s.6g', 140, 58610.797_8, 4.0, '3s.9f', 131, 60301.281_8, 3.0, 1690.5027_8, 5915.3999, -1.767, 4.660E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 5.0, '3s.9f', 331, 60301.281_8, 4.0, 1690.5027_8, 5915.3999, -1.680, 4.430E+03, 0.00), & !+ NIST LINES('3s.5f', 131, 57204.164_8, 3.0, '3s.6g', 340, 58610.797_8, 4.0, 1406.6478_8, 7109.1001, 0.866, 1.080E+06, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.227_8, 2.0, '3s.6g', 340, 58610.797_8, 3.0, 1406.5687_8, 7109.5000, 0.724, 9.990E+05, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.266_8, 3.0, '3s.6g', 340, 58610.797_8, 3.0, 1406.5291_8, 7109.7002, -0.335, 8.710E+04, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.266_8, 3.0, '3s.6g', 340, 58610.797_8, 4.0, 1406.5291_8, 7109.7002, 0.840, 1.010E+06, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.6g', 340, 58610.797_8, 3.0, 1406.5094_8, 7109.7998, -2.153, 1.330E+03, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.6g', 340, 58610.797_8, 4.0, 1406.5094_8, 7109.7998, -0.336, 6.770E+04, 0.00), & !+ NIST LINES('3s.5f', 331, 57204.305_8, 4.0, '3s.6g', 340, 58610.797_8, 5.0, 1406.5094_8, 7109.7998, 0.954, 1.080E+06, 0.00), & !+ NIST LINES('3s.6f', 131, 58575.477_8, 3.0, '3s.8g', 340, 59950.871_8, 4.0, 1375.4022_8, 7270.6001, 0.218, 2.320E+05, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 4.0, '3s.8g', 340, 59950.871_8, 4.0, 1375.3456_8, 7270.8999, -0.985, 1.450E+04, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 3.0, '3s.8g', 340, 59950.871_8, 4.0, 1375.3456_8, 7270.8999, 0.191, 2.180E+05, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 4.0, '3s.8g', 340, 59950.871_8, 5.0, 1375.3456_8, 7270.8999, 0.306, 2.320E+05, 0.00), & !+ NIST LINES('3s.6g', 140, 58610.797_8, 4.0, '3s.8f', 131, 59935.371_8, 3.0, 1324.5911_8, 7549.5000, -1.307, 8.250E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.8f', 331, 59935.371_8, 4.0, 1324.5911_8, 7549.5000, -4.326, 6.130E+00, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.8f', 131, 59935.371_8, 3.0, 1324.5911_8, 7549.5000, -2.508, 5.180E+02, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.8f', 331, 59935.371_8, 4.0, 1324.5911_8, 7549.5000, -2.509, 4.030E+02, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.8f', 331, 59935.371_8, 2.0, 1324.5911_8, 7549.5000, -1.451, 8.290E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 5.0, '3s.8f', 331, 59935.371_8, 4.0, 1324.5911_8, 7549.5000, -1.217, 7.880E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.8f', 131, 59935.371_8, 3.0, 1324.5911_8, 7549.5000, -1.333, 7.760E+03, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.6f', 331, 58575.527_8, 4.0, 1312.7838_8, 7617.3999, -3.945, 1.450E+01, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.6f', 331, 58575.527_8, 3.0, 1312.7838_8, 7617.3999, -2.126, 1.230E+03, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.6f', 331, 58575.527_8, 4.0, 1312.7838_8, 7617.3999, -2.127, 9.540E+02, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 3.0, '3s.6f', 331, 58575.527_8, 2.0, 1312.7838_8, 7617.3999, -1.069, 1.960E+04, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.6f', 331, 58575.527_8, 3.0, 1312.7838_8, 7617.3999, -0.952, 1.830E+04, 0.00), & !+ NIST LINES('3s.5g', 340, 57262.762_8, 4.0, '3s.6f', 131, 58575.477_8, 3.0, 1312.7322_8, 7617.7002, -0.925, 1.950E+04, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 5.0, '3s.10f', 331, 60562.637_8, 4.0, 1139.1079_8, 8778.7998, -1.405, 3.790E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 4.0, '3s.10f', 131, 60562.637_8, 3.0, 1139.1079_8, 8778.7998, -1.521, 3.730E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.10f', 331, 60562.637_8, 2.0, 1139.1079_8, 8778.7998, -1.638, 3.990E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 4.0, '3s.10f', 331, 60562.637_8, 4.0, 1139.1079_8, 8778.7998, -2.696, 1.940E+02, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.10f', 131, 60562.637_8, 3.0, 1139.1079_8, 8778.7998, -2.696, 2.490E+02, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.10f', 331, 60562.637_8, 4.0, 1139.1079_8, 8778.7998, -4.512, 2.960E+00, 0.00), & !+ NIST LINES('3s.7g', 140, 59423.535_8, 4.0, '3s.10f', 131, 60562.637_8, 3.0, 1139.1079_8, 8778.7998, -1.493, 3.970E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 5.0, '3s.9f', 331, 60301.281_8, 4.0, 877.7473_8, 11392.7998, -0.954, 6.340E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 4.0, '3s.9f', 131, 60301.281_8, 3.0, 877.7473_8, 11392.7998, -1.071, 6.240E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.9f', 331, 60301.281_8, 2.0, 877.7473_8, 11392.7998, -1.188, 6.660E+03, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 4.0, '3s.9f', 331, 60301.281_8, 4.0, 877.7473_8, 11392.7998, -2.246, 3.240E+02, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.9f', 131, 60301.281_8, 3.0, 877.7473_8, 11392.7998, -2.246, 4.160E+02, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.9f', 331, 60301.281_8, 4.0, 877.7473_8, 11392.7998, -4.065, 4.920E+00, 0.00), & !+ NIST LINES('3s.7g', 140, 59423.535_8, 4.0, '3s.9f', 131, 60301.281_8, 3.0, 877.7473_8, 11392.7998, -1.041, 6.670E+03, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 4.0, '3s.7g', 340, 59423.535_8, 5.0, 848.0102_8, 11792.3105, 0.910, 3.540E+05, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 4.0, '3s.7g', 340, 59423.535_8, 4.0, 848.0102_8, 11792.3105, -0.381, 2.220E+04, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 4.0, '3s.7g', 340, 59423.535_8, 3.0, 848.0102_8, 11792.3105, -2.199, 4.340E+02, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 3.0, '3s.7g', 340, 59423.535_8, 4.0, 848.0102_8, 11792.3105, 0.794, 3.320E+05, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 3.0, '3s.7g', 340, 59423.535_8, 3.0, 848.0102_8, 11792.3105, -0.381, 2.850E+04, 0.00), & !+ NIST LINES('3s.6f', 331, 58575.527_8, 2.0, '3s.7g', 340, 59423.535_8, 3.0, 848.0102_8, 11792.3105, 0.677, 3.250E+05, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.7f', 131, 59400.762_8, 3.0, 789.9705_8, 12658.7002, -0.529, 1.760E+04, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.7f', 331, 59400.762_8, 4.0, 789.9705_8, 12658.7002, -3.550, 1.300E+01, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.7f', 131, 59400.762_8, 3.0, 789.9705_8, 12658.7002, -1.732, 1.100E+03, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 4.0, '3s.7f', 331, 59400.762_8, 4.0, 789.9705_8, 12658.7002, -1.732, 8.570E+02, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 3.0, '3s.7f', 331, 59400.762_8, 2.0, 789.9705_8, 12658.7002, -0.675, 1.760E+04, 0.00), & !+ NIST LINES('3s.6g', 340, 58610.797_8, 5.0, '3s.7f', 331, 59400.762_8, 4.0, 789.9705_8, 12658.7002, -0.441, 1.670E+04, 0.00), & !+ NIST LINES('3s.6g', 140, 58610.797_8, 4.0, '3s.7f', 131, 59400.762_8, 3.0, 789.9705_8, 12658.7002, -0.557, 1.650E+04, 0.00), & !+ NIST LINES('3s.8g', 140, 59950.871_8, 4.0, '3s.10f', 131, 60562.637_8, 3.0, 611.7704_8, 16346.0000, -0.855, 4.970E+03, 0.00), & !+ NIST LINES('3s.7f', 331, 59400.762_8, 4.0, '3s.8g', 340, 59950.871_8, 5.0, 550.1095_8, 18178.1992, 0.886, 1.410E+05, 0.00), & !+ NIST LINES('3s.7f', 131, 59400.762_8, 3.0, '3s.8g', 340, 59950.871_8, 4.0, 550.1095_8, 18178.1992, 0.770, 1.320E+05, 0.00), & !+ NIST LINES('3s.7f', 331, 59400.762_8, 4.0, '3s.8g', 340, 59950.871_8, 4.0, 550.1095_8, 18178.1992, -0.405, 8.820E+03, 0.00), & !+ NIST LINES('3s.7f', 131, 59400.762_8, 3.0, '3s.8g', 140, 59950.871_8, 4.0, 550.1095_8, 18178.1992, 0.800, 1.410E+05, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 5.0, '3s.8f', 331, 59935.371_8, 4.0, 511.8336_8, 19537.5996, -0.197, 1.230E+04, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 4.0, '3s.8f', 131, 59935.371_8, 3.0, 511.8336_8, 19537.5996, -0.312, 1.220E+04, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.8f', 331, 59935.371_8, 2.0, 511.8336_8, 19537.5996, -0.430, 1.300E+04, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 4.0, '3s.8f', 331, 59935.371_8, 4.0, 511.8336_8, 19537.5996, -1.487, 6.330E+02, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.8f', 131, 59935.371_8, 3.0, 511.8336_8, 19537.5996, -1.487, 8.130E+02, 0.00), & !+ NIST LINES('3s.7g', 340, 59423.535_8, 3.0, '3s.8f', 331, 59935.371_8, 4.0, 511.8336_8, 19537.5996, -3.305, 9.620E+00, 0.00), & !+ NIST LINES('3s.7g', 140, 59423.535_8, 4.0, '3s.8f', 131, 59935.371_8, 3.0, 511.8336_8, 19537.5996, -0.286, 1.290E+04, 0.00), & !+ NIST LINES('3s.8g', 340, 59950.871_8, 5.0, '3s.9f', 331, 60301.281_8, 4.0, 350.4100_8, 28538.0000, -0.025, 8.590E+03, 0.00), & !+ NIST LINES('3s.8g', 340, 59950.871_8, 4.0, '3s.9f', 131, 60301.281_8, 3.0, 350.4100_8, 28538.0000, -0.141, 8.460E+03, 0.00), & !+ NIST LINES('3s.8g', 340, 59950.871_8, 4.0, '3s.9f', 331, 60301.281_8, 4.0, 350.4100_8, 28538.0000, -1.316, 4.400E+02, 0.00), & !+ NIST + 123 transitions LINES('3s.6h', 351, 58618.942_8, 21.5, '3s.7i', 360, 59430.517_8, 25.5, 811.5750, 12321.7000, 1.950, 7.530E+05, 0.00), & !+ Green LINES('3s.7h', 351, 59428.853_8, 21.5, '3s.9i', 360, 60315.722_8, 25.5, 886.8690, 11275.6000, 1.120, 1.330E+05, 0.00), & !+ Green LINES('3s.7i', 360, 59430.517_8, 25.5, '3s.9k', 370, 60316.041_8, 29.5, 885.524, 11292.7000, 1.080, 1.050E+05, 0.00), & !+ Green LINES('3s.6g', 140, 58610.825_8, 4.0, '3s.7h', 351, 59428.853_8, 21.5, 818.028, 12224.5000, 1.128, 1.362E+05, 0.00), & !+ Green LINES('3s.6g', 340, 58610.805_8, 3.0, '3s.7h', 351, 59428.853_8, 21.5, 818.048, 12224.2000, 1.019, 1.060E+05, 0.00), & !+ Green LINES('3s.6g', 340, 58610.795_8, 4.0, '3s.7h', 351, 59428.853_8, 21.5, 818.058, 12224.1000, 1.128, 1.362E+05, 0.00), & !+ Green LINES('3s.6g', 340, 58610.815_8, 5.0, '3s.7h', 351, 59428.853_8, 21.5, 818.038, 12224.4000, 1.215, 1.664E+05, 0.00), & !+ Green LINES('3s.7i', 360, 59430.517_8, 25.5, '3s.8k', 371, 59956.080_8, 29.5, 525.563, 19027.2146, 2.086, 3.741E+05, 0.00), & !++ Green LINES('3s.8k', 371, 59956.080_8, 29.5, '3s.9l', 380, 60316.150_8, 33.5, 360.068, 27772.3776, 2.205, 2.040E+05, 0.00), & !+ Green LINES('3s.9l', 380, 60316.150_8, 33.5, '3s.10m', 391, 60573.620_8, 37.5, 257.470, 38839.4764, 2.310, 1.189E+05, 0.00), & !+ Green LINES('3s.5g', 140, 57262.760_8, 4.0, '3s.6h', 351, 58618.942_8, 21.5, 1356.182, 7373.6416, 1.183, 4.253E+05, 0.00), & !+ Green LINES('3s.5g', 340, 57262.760_8, 3.0, '3s.6h', 351, 58618.942_8, 21.5, 1356.182, 7373.6416, 1.074, 3.308E+05, 0.00), & !+ Green LINES('3s.5g', 340, 57262.760_8, 4.0, '3s.6h', 351, 58618.942_8, 21.5, 1356.182, 7373.6416, 1.183, 4.253E+05, 0.00), & !+ Green LINES('3s.5g', 340, 57262.760_8, 5.0, '3s.6h', 351, 58618.942_8, 21.5, 1356.182, 7373.6416, 1.271, 5.199E+05, 0.00), & !+ Green LINES('3s.3d', 320, 47957.027_8, 2.0, '3s.7f', 331, 59400.763_8, 3.0, 11443.736, 873.6006, -0.860, 1.720E+06, 0.00), & !+ KALD LINES('3s.3d', 320, 47957.045_8, 3.0, '3s.7f', 131, 59400.763_8, 3.0, 11443.718, 873.6019, -2.170, 8.440E+04, 0.00) /) !+ KALD ! ! Force de collisions effectives avec H de Barklem (communication privée 2011) ! 21 transitions bb 7 transitions bf integer, parameter :: BarkN = 55 integer, dimension(BarkN, 8), parameter :: BarkConf = reshape(source = (/ & 3,0,100,1,3,1,311,1,& 3,0,100,1,3,1,311,3,& 3,0,100,1,3,1,311,5,& 3,0,100,1,3,1,111,3,& 3,0,100,1,4,0,300,3,& 3,0,100,1,4,0,100,1,& 3,0,100,1,3,2,120,5,& 3,0,100,1,4,1,311,1,& 3,0,100,1,4,1,311,3,& 3,0,100,1,4,1,311,5,& 3,1,311,1,3,1,311,3,& 3,1,311,1,3,1,311,5,& 3,1,311,1,3,1,111,3,& 3,1,311,1,4,0,300,3,& 3,1,311,1,4,0,100,1,& 3,1,311,1,3,2,120,5,& 3,1,311,1,4,1,311,1,& 3,1,311,1,4,1,311,3,& 3,1,311,1,4,1,311,5,& 3,1,311,3,3,1,311,5,& 3,1,311,3,3,1,111,3,& 3,1,311,3,4,0,300,3,& 3,1,311,3,4,0,100,1,& 3,1,311,3,3,2,120,5,& 3,1,311,3,4,1,311,1,& 3,1,311,3,4,1,311,3,& 3,1,311,3,4,1,311,5,& 3,1,311,5,3,1,111,3,& 3,1,311,5,4,0,300,3,& 3,1,311,5,4,0,100,1,& 3,1,311,5,3,2,120,5,& 3,1,311,5,4,1,311,1,& 3,1,311,5,4,1,311,3,& 3,1,311,5,4,1,311,5,& 3,1,111,3,4,0,300,3,& 3,1,111,3,4,0,100,1,& 3,1,111,3,3,2,120,5,& 3,1,111,3,4,1,311,1,& 3,1,111,3,4,1,311,3,& 3,1,111,3,4,1,311,5,& 4,0,300,3,4,0,100,1,& 4,0,300,3,3,2,120,5,& 4,0,300,3,4,1,311,1,& 4,0,300,3,4,1,311,3,& 4,0,300,3,4,1,311,5,& 4,0,100,1,3,2,120,5,& 4,0,100,1,4,1,311,1,& 4,0,100,1,4,1,311,3,& 4,0,100,1,4,1,311,5,& 3,2,120,5,4,1,311,1,& 3,2,120,5,4,1,311,3,& 3,2,120,5,4,1,311,5,& 4,1,311,1,4,1,311,3,& 4,1,311,1,4,1,311,5,& 4,1,311,3,4,1,311,5 /), & shape = (/ BarkN, 8/), order = (/ 2, 1/)) Real, dimension(BarkN, 5), parameter :: Bark_Upsilon = reshape(source = (/ & 3.28E-07,5.93E-07,1.53E-06,3.68E-06,7.95E-06,& 9.83E-07,1.78E-06,4.59E-06,1.10E-05,2.38E-05,& 1.64E-06,2.96E-06,7.65E-06,1.84E-05,3.97E-05,& 3.66E-06,6.30E-06,9.89E-06,1.53E-05,2.33E-05,& 3.23E-05,4.22E-05,4.96E-05,6.22E-05,7.74E-05,& 2.78E-05,4.17E-05,4.58E-05,5.46E-05,7.15E-05,& 4.00E-05,3.28E-05,3.50E-05,4.29E-05,5.24E-05,& 3.54E-06,4.26E-06,4.87E-06,5.50E-06,6.26E-06,& 1.06E-05,1.28E-05,1.46E-05,1.65E-05,1.88E-05,& 1.77E-05,2.13E-05,2.43E-05,2.75E-05,3.13E-05,& 0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,& 0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,& 6.62E-03,9.47E-03,1.12E-02,1.26E-02,1.41E-02,& 2.39E-02,2.88E-02,2.79E-02,2.79E-02,2.84E-02,& 1.25E-02,1.78E-02,1.69E-02,1.61E-02,1.62E-02,& 1.41E-02,1.06E-02,9.21E-03,9.02E-03,9.08E-03,& 9.28E-04,1.09E-03,1.06E-03,1.00E-03,9.38E-04,& 2.78E-03,3.26E-03,3.17E-03,3.00E-03,2.81E-03,& 4.64E-03,5.44E-03,5.29E-03,5.00E-03,4.69E-03,& 0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,& 1.99E-02,2.84E-02,3.36E-02,3.79E-02,4.22E-02,& 7.16E-02,8.64E-02,8.38E-02,8.38E-02,8.53E-02,& 3.74E-02,5.35E-02,5.06E-02,4.82E-02,4.86E-02,& 4.22E-02,3.19E-02,2.76E-02,2.71E-02,2.72E-02,& 2.78E-03,3.26E-03,3.17E-03,3.00E-03,2.81E-03,& 8.35E-03,9.79E-03,9.52E-03,9.00E-03,8.44E-03,& 1.39E-02,1.63E-02,1.59E-02,1.50E-02,1.41E-02,& 3.31E-02,4.74E-02,5.60E-02,6.32E-02,7.04E-02,& 1.19E-01,1.44E-01,1.40E-01,1.40E-01,1.42E-01,& 6.23E-02,8.92E-02,8.43E-02,8.03E-02,8.09E-02,& 7.03E-02,5.31E-02,4.61E-02,4.51E-02,4.54E-02,& 4.64E-03,5.44E-03,5.29E-03,5.00E-03,4.69E-03,& 1.39E-02,1.63E-02,1.59E-02,1.50E-02,1.41E-02,& 2.32E-02,2.72E-02,2.64E-02,2.50E-02,2.34E-02,& 1.95E-01,3.30E-01,4.26E-01,5.34E-01,6.69E-01,& 8.40E-02,1.34E-01,1.49E-01,1.61E-01,1.77E-01,& 2.55E-02,3.80E-02,4.90E-02,6.15E-02,7.73E-02,& 2.68E-03,5.04E-03,6.12E-03,6.85E-03,7.68E-03,& 8.03E-03,1.51E-02,1.84E-02,2.05E-02,2.30E-02,& 1.34E-02,2.52E-02,3.06E-02,3.42E-02,3.84E-02,& 1.20E-01,2.24E-01,3.00E-01,3.65E-01,4.34E-01,& 2.88E-02,4.56E-02,6.72E-02,1.02E-01,1.44E-01,& 1.02E-02,1.50E-02,1.98E-02,2.42E-02,2.85E-02,& 3.07E-02,4.50E-02,5.94E-02,7.26E-02,8.54E-02,& 5.12E-02,7.51E-02,9.91E-02,1.21E-01,1.42E-01,& 3.05E-01,9.42E-01,1.16E+00,1.26E+00,1.33E+00,& 2.68E-03,1.27E-02,2.32E-02,3.23E-02,4.05E-02,& 8.03E-03,3.80E-02,6.95E-02,9.69E-02,1.22E-01,& 1.34E-02,6.33E-02,1.16E-01,1.62E-01,2.03E-01,& 1.67E-02,7.06E-02,1.35E-01,1.92E-01,2.43E-01,& 5.00E-02,2.12E-01,4.06E-01,5.75E-01,7.30E-01,& 8.33E-02,3.53E-01,6.77E-01,9.58E-01,1.22E+00,& 0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,& 0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,& 0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00 /), & shape = (/BarkN, 5/), order = (/2, 1/)) ! integer, parameter :: BarkNi = 11 integer, dimension(BarkNi, 4), parameter :: BarkiConf = reshape(source = (/ & 3,0,100,1,& 3,1,311,1,& 3,1,311,3,& 3,1,311,5,& 3,1,111,3,& 4,0,300,3,& 4,0,100,1,& 3,2,120,5,& 4,1,311,1,& 4,1,311,3,& 4,1,311,5 /),& shape = (/ BarkNi, 4/), order = (/ 2, 1/)) Real, dimension(BarkNi, 5), parameter :: Barki_Upsilon = reshape(source = (/ & 2.72E-05,4.73E-05,7.60E-05,1.08E-04,1.41E-04,& 1.91E-03,4.57E-03,8.45E-03,1.14E-02,1.38E-02,& 5.74E-03,1.37E-02,2.53E-02,3.42E-02,4.13E-02,& 9.56E-03,2.28E-02,4.22E-02,5.71E-02,6.88E-02,& 1.68E-01,1.75E-01,2.46E-01,3.42E-01,4.35E-01,& 1.03E+00,1.19E+00,1.41E+00,1.68E+00,1.99E+00,& 1.37E+01,1.69E+01,1.92E+01,2.10E+01,2.27E+01,& 5.59E+00,6.68E+00,7.38E+00,8.00E+00,8.57E+00,& 1.70E-01,2.10E-01,2.33E-01,2.52E-01,2.70E-01,& 5.11E-01,6.29E-01,6.98E-01,7.55E-01,8.10E-01,& 8.51E-01,1.05E+00,1.16E+00,1.26E+00,1.35E+00/),& shape = (/BarkNi, 5/), order = (/2, 1/)) ! ! Force de collisions effectives avec les électrons de Zatsarinny et al. (2009) ! integer, parameter :: ZatsN = 14 integer, dimension(ZatsN, 8), parameter :: ZatsConf = reshape(source = (/ & 3,0,100,1,3,1,311,1,& 3,0,100,1,3,1,311,3,& 3,0,100,1,3,1,311,5,& 3,0,100,1,3,1,111,3,& 3,0,100,1,4,0,300,3,& 3,0,100,1,4,0,100,1,& 3,0,100,1,3,2,120,5,& 3,0,100,1,4,1,311,1,& 3,0,100,1,4,1,311,3,& 3,0,100,1,4,1,311,5,& 3,0,100,1,3,2,320,3,& 3,0,100,1,3,2,320,5,& 3,0,100,1,3,2,320,7,& 3,0,100,1,4,1,111,3/),& shape = (/ ZatsN, 8/), order = (/ 2, 1/)) Real, dimension(ZatsN, 8), parameter :: Zats_Upsilon = reshape(source = (/ & 2.39E-02,6.99E-02,1.67E-01,2.61E-01,3.50E-01,4.34E-01,6.22E-01,7.84E-01,& 7.17E-02,2.10E-01,5.02E-01,7.83E-01,1.05E+00,1.30E+00,1.86E+00,2.35E+00,& 1.19E-01,3.50E-01,8.36E-01,1.31E+00,1.75E+00,2.17E+00,3.11E+00,3.92E+00,& 4.06E-02,1.03E-01,3.03E-01,5.97E-01,9.80E-01,1.45E+00,3.03E+00,5.16E+00,& 1.24E-02,4.27E-02,1.18E-01,1.92E-01,2.59E-01,3.21E-01,4.52E-01,5.62E-01,& 4.00E-03,1.43E-02,4.00E-02,6.83E-02,9.83E-02,1.30E-01,2.20E-01,3.29E-01,& 2.01E-02,4.97E-02,1.15E-01,1.94E-01,2.91E-01,4.06E-01,7.73E-01,1.24E+00,& 1.74E-03,4.26E-03,9.29E-03,1.40E-02,1.82E-02,2.22E-02,3.15E-02,4.03E-02,& 5.21E-03,1.28E-02,2.79E-02,4.19E-02,5.47E-02,6.67E-02,9.44E-02,1.21E-01,& 8.68E-03,2.13E-02,4.65E-02,6.98E-02,9.12E-02,1.11E-01,1.57E-01,2.01E-01,& 1.85E-03,5.53E-03,1.47E-02,2.51E-02,3.63E-02,4.81E-02,8.06E-02,1.17E-01,& 3.09E-03,9.22E-03,2.46E-02,4.19E-02,6.05E-02,8.02E-02,1.34E-01,1.95E-01,& 4.32E-03,1.29E-02,3.44E-02,5.87E-02,8.47E-02,1.12E-01,1.88E-01,2.73E-01,& 8.86E-03,2.34E-02,5.74E-02,9.64E-02,1.42E-01,1.97E-01,3.86E-01,6.59E-01/),& shape = (/ZatsN, 8/), order = (/2, 1/)) Integer, parameter :: App1_N = 19 Real, dimension(App1_N), parameter :: App1_lambda = (/ & 3829.36, & 4167.27, & 4571.10, & 4702.99, & 4730.03, & 5167.32, & 5172.68, & 5183.60, & 5528.41, & 5711.08, & 7657.60, & 8471.27, & 8715.26, & 8736.02, & 8806.75, & 8923.56, & 8997.15, & 10312.52, & 11828.18 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ & 100, & 100, & 35, & 100, & 50, & 100, & 100, & 100, & 100, & 50, & 35, & 35, & 35, & 35, & 100, & 50, & 35, & 35, & 150 /) ! vecteur d'entiers représentant les niveaux bas/hauts des transitions indexées pour Gaia Integer, dimension(App1_N) :: GNi = (/ 1, 2, 3, 4, 5, 5, 5, 5, 5, 5, 5, 7, 15, 6, 11, 9, 9, 12, 12/) Integer, dimension(App1_N) :: GNj = (/ 3, 6, 6, 6, 7, 8, 17, 18, 31, 32, 71, 29, 71, 24, 70, 73, 90, 79, 82/) ! End Module datagraph_MgI ! !*********************************************************************** ! module datagraph_MgII use MTD_STRUCT implicit none Character(len = 2), parameter :: Symb = 'Mg' Integer, parameter :: Z_IONI = 2 Character(len = 4) :: cpage = 'DA4L' ! DA4P integer,parameter :: lenterm = 6 character(len=32) :: continuum = '15.04' real,dimension(lenterm) :: SLP = (/200,211,220,231,240,251/) real,dimension(lenterm) :: iterm = (/1,2,3,4,5,6/) character(len=7),dimension(lenterm) :: cterm = (/"[2$S ","[2$P[o$","[2$D ","[2$F[o$","[2$G ","[2$H[o$"/) ! Position et taille du graphique integer, parameter :: posx = 250, posy = 1980 integer, parameter :: lenx = 2600, leny = 1800 !integer, parameter :: posx = 200, posy = 2850 !integer, parameter :: lenx = 1800, leny = 2600 ! Abscisses du continuum real, dimension(8), parameter :: abscontinuum = (/0.,1.,2.,3.,4.,5.,6.,7./) ! Noms des niveaux integer,parameter :: levelsize = 13 character(len= 5),dimension(levelsize):: levelname = & (/'3s ','4s ','10s','3p ','10p','3d ','10d','4f ','10f','5g ','10g','6h ','10h'/) ! Position des niveaux real,dimension(levelsize,2) :: levelpos = reshape(source=& (/ 0.5,0.2, 0.5,8.8, 0.5,14.5, 1.5,4.6, 1.5,14.6, 2.5,9.0, 2.5,14.6, 3.5,11.8, 3.5,14.6, & 4.5,13., 4.5,14.7, 5.5,13.6, 5.5,14.7/),& shape=(/levelsize,2/),& order=(/2,1/)) integer,parameter :: Klevels = 24 type(LEVELS), dimension(Klevels),save :: S_level_sup = (/ & ! Ajout de niveaux supplémentaires pour être consistent LEVELS('10p',211,116206.92_8,1.0589554D+00,1.4407821D+01,2.,0), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('10p',211,116206.92_8,1.0589554D+00,1.4407821D+01,4.,0),& ! avec levée de dégénérescence ! Ajout de niveaux supplémentaires présents dans NIST mais pas dans la TOPBASE LEVELS('2p6.5g',240,103705.660_8,9.4503551E-01,1.2858116E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.5g',240,103705.660_8,9.4503551E-01,1.2858116E+01,8.0,0),& ! NIST sans TOPBASE LEVELS('2p6.6g',240,109072.050_8,9.9393765E-01,1.3523477E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.6g',240,109072.050_8,9.9393765E-01,1.3523477E+01,8.0,0),& ! NIST sans TOPBASE LEVELS('2p6.7g',240,112307.790_8,1.0234239E+00,1.3924665E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.7g',240,112307.790_8,1.0234239E+00,1.3924665E+01,8.0,0),& ! NIST sans TOPBASE LEVELS('2p6.8g',240,114407.880_8,1.0425613E+00,1.4185048E+01,8.0,0),& ! NIST sans TOPBASE LEVELS('2p6.8g',240,114407.880_8,1.0425613E+00,1.4185048E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.9g',240,115847.670_8,1.0556816E+00,1.4363563E+01,8.0,0),& ! NIST sans TOPBASE LEVELS('2p6.9g',240,115847.670_8,1.0556816E+00,1.4363563E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.10g',240,116877.540_8,1.0650665E+00,1.4491253E+01,8.0,0),& ! NIST sans TOPBASE LEVELS('2p6.10g',240,116877.540_8,1.0650665E+00,1.4491253E+01,10.0,0),&! NIST sans TOPBASE LEVELS('2p6.6h',251, 109073.970_8,9.9395515E-01,1.3523715E+01,12.0,0),& ! NIST sans TOPBASE LEVELS('2p6.6h',251, 109073.970_8,9.9395515E-01,1.3523715E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.7h',251, 112309.060_8,1.0234355E+00,1.3924823E+01,12.0,0),& ! NIST sans TOPBASE LEVELS('2p6.7h',251, 112309.060_8,1.0234355E+00,1.3924823E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.8h',251, 114408.740_8,1.0425692E+00,1.4185155E+01,12.0,0),& ! NIST sans TOPBASE LEVELS('2p6.8h',251, 114408.740_8,1.0425692E+00,1.4185155E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.9h',251, 115848.280_8,1.0556872E+00,1.4363639E+01,12.0,0),& ! NIST sans TOPBASE LEVELS('2p6.9h',251, 115848.280_8,1.0556872E+00,1.4363639E+01,10.0,0),& ! NIST sans TOPBASE LEVELS('2p6.10h',251, 116878.040_8,1.0650711E+00,1.4491315E+01,12.0,0),& ! NIST sans TOPBASE LEVELS('2p6.10h',251, 116878.040_8,1.0650711E+00,1.4491315E+01,10.0,0) /) ! NIST sans TOPBASE ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! ! Force de collisions effectives de Sigut & Pradhan(J Phys B Sigut & Pradhan 1995) ! 36 transitions (17 permises, 19 interdites) integer, parameter :: SiguN = 36 integer, dimension(SiguN, 4), parameter :: SiguConf = reshape( & source = (/ 3, 0, 3, 1, & ! Transitions permises 3, 0, 4, 1, & 3, 0, 5, 1, & 3, 1, 4, 0, & 3, 1, 3, 2, & 3, 1, 5, 0, & 3, 1, 4, 2, & 4, 0, 4, 1, & 4, 0, 5, 1, & 3, 2, 4, 1, & 3, 2, 4, 3, & 3, 2, 5, 1, & 4, 1, 5, 0, & 4, 1, 4, 2, & 5, 0, 5, 1, & 4, 2, 4, 3, & 4, 2, 5, 1, & 3, 0, 4, 0, & ! Transitions interdites 3, 0, 3, 2, & 3, 0, 5, 0, & 3, 0, 4, 2, & 3, 0, 4, 3, & 3, 1, 4, 1, & 3, 1, 4, 3, & 3, 1, 5, 1, & 4, 0, 3, 2, & 4, 0, 5, 0, & 4, 0, 4, 2, & 4, 0, 4, 3, & 3, 2, 5, 0, & 3, 2, 4, 2, & 4, 1, 4, 3, & 4, 1, 5, 1, & 5, 0, 4, 2, & 5, 0, 4, 3, & 4, 3, 5, 1/), & shape = (/ SiguN, 4 /), & order = (/ 2, 1 /) ) Real, dimension(SiguN, 6), parameter :: Sigu_Upsilon = reshape( & source = (/ 15.1, 15.5, 15.9, 16.9, 18.6, 20.2, & ! Transitions permises .597, .705, .763, .790, .751, .715, & .311, .330, .340, .345, .346, .345, & 8.78, 7.18, 6.48, 5.82, 5.50, 5.62, & 30.8, 32.2, 33.1, 34.9, 38.1, 41.5, & 1.95, 1.62, 1.53, 1.45, 1.40, 1.39, & 6.00, 5.74, 5.62, 5.41, 5.20, 5.12, & 32.5, 38.4, 43.6, 54.6, 74.8, 93.4, & 2.49, 2.61, 2.70, 2.78, 2.81, 2.81, & 61.3, 68.1, 73.1, 82.5, 97.0, 10.9, & 64.7, 70.8, 76.2, 87.2, 104., 118., & 12.3, 12.4, 12.3, 11.9, 11.5, 11.4, & 14.8, 13.7, 14.0, 15.9, 21.4, 27.1, & 52.0, 53.4, 57.9, 71.9, 101., 128., & 53.1, 74.9, 96.7, 147., 235., 312., & 1720., 1840., 1980., 2280., 2720., 3010., & 129., 173., 214., 299., 417., 501., & 1.92, 1.60, 1.45, 1.30, 1.17, 1.11, & ! Transitions interdites 2.70, 2.82, 2.89, 2.99, 3.07, 3.10, & .463, .430, .425, .419, .413, .409, & .719, .733, .733, .718, .698, .687, & .590, .613, .618, .619, .614, .610, & 6.49, 7.39, 7.72, 7.75, 7.31, 6.99, & 5.18, 5.50, 5.67, 5.90, 6.06, 6.13, & 2.67, 2.76, 2.78, 2.75, 2.71, 2.68, & 16.3, 15.0, 15.1, 14.7, 13.5, 12.6, & 5.80, 5.31, 5.24, 5.16, 5.05, 4.99, & 9.92, 9.50, 9.54, 9.59, 9.58, 9.54, & 9.36, 9.51, 9.61, 9.83, 10.1, 10.2, & 10.4, 8.34, 7.70, 7.01, 6.48, 6.26, & 51.2, 45.2, 42.7, 39.8, 37.7, 36.9, & 36.7, 36.2, 36.3, 36.7, 37.1, 37.3, & 21.3, 22.6, 23.1, 23.4, 23.3, 23.2, & 32.4, 34.2, 33.8, 32.7, 31.8, 31.4, & 22.3, 23.1, 22.0, 19.9, 18.0, 17.2, & 83.0, 82.4, 80.6, 77.2, 74.4, 73.2/),& shape = (/ SiguN, 6/), & order = (/ 2, 1/) ) ! Integer, parameter :: App1_N = 1 Real, dimension(App1_N), parameter :: App1_lambda = (/ 0.0 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ 0.0 /) ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia Integer, dimension(App1_N) :: GNi = (/ 1 /) Integer, dimension(App1_N) :: GNj = (/ 3 /) ! end module datagraph_MgII ! !*********************************************************************** ! Module datagraph_CaI use MTD_STRUCT implicit none Character(len = 2), parameter :: Symb = 'Ca' integer, parameter :: Z_IONI = 1 character(len=4) :: cpage = 'DA4L' integer,parameter :: lenterm = 20 character(len=32) :: continuum = '6.11' !eV real,dimension(lenterm) :: SLP = (/100,111,120,121,131,140,151,& 300,310,311,320,321,330,331,340,351,& 360,371,380,391/) real,dimension(lenterm) :: iterm = (/1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20/) character(len=7),dimension(lenterm) :: cterm = (/ & "[1$S ","[1$P[o$","[1$D ","[1$D[o$","[1$F[o$","[1$G ","[1$H[o$",& ! Système de multiplicité 1 "[3$S ","[3$P ","[3$P[o$","[3$D ","[3$D[o$","[3$F ","[3$F[o$","[3$G ",& ! Système de multiplicité 3 "[3$H[o$","I ","K[o$ ","L ","M[o$ "/) ! Système de multiplicité dégénéré ! Position et taille du graphique integer, parameter :: posx = 250, posy = 1980 integer, parameter :: lenx = 2600, leny = 1800 !1800 ! Abscisses du continuum real, dimension(16), parameter :: abscontinuum = (/0.3,1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,20.5/) ! Noms des niveaux integer,parameter :: levelsize = 44 character(len= 5),dimension(levelsize):: levelname = & (/'4s ','4p[2$','10s ','4p ','10p ','3d ','4p[2$','3d5s ','10d ','3d4p ','3d4p ','4f ','10f ','5g ','3d[2$',& '10g ','6h ','10h ',& '5s ','10s ','4p[2$','3d[2$','4p ','3d4p ','10p ','3d ','3d5s ','10d ','3d4p ',& '3d[2$','3d4p ','4f ','10f ','5g ','10g ','6h ','10h ','7i ','10i ','8k ','10k ','9l ','10l ','10m '/) ! Position des niveaux real,dimension(levelsize,2) :: levelpos = reshape(source=& (/ 0.4,0.05, 0.4,5.3, 0.35,6., 1.4,2.95, 1.3,6., 2.4,2.8, 2.4,5.2, 3.3,5.95, 2.3,6.07, & 3.1,4.55, 4.1,5.15, 4.4,5.35, 4.4,6.07, 5.4,5.66, 5.32,5.95, 5.35,6.15, 6.4,5.8, 6.35,6.07,& 7.4,4.01, 7.4,6.0, 8.4,4.87, 8.4,6.08, 9.4,1.98, 9.1,4.97, 9.35,6.04, & 11.3,2.6, 11.3,5.95, 10.35,6.08, 11.1,4.83, 12.4,5.49, 14.3,4.54, 13.4,5.33, 13.4,6.07,& 14.4,5.65, 14.35,6.07, & 15.4,5.80, 15.35,6.07, 16.4,5.85, 16.35,6.07, 17.4,5.95, 17.35,6.07, 18.4,5.95, 18.35,6.07, 19.35,6.10 /),& shape=(/levelsize,2/),& order=(/2,1/)) integer,parameter :: Klevels = 48 !27 type(LEVELS), dimension(Klevels),save :: S_level_sup = (/ & ! Ajout de niveaux supplémentaires pour être consistent LEVELS('3p6.4s.6g',140,46246.090_8,4.2142538D-01,5.7337840D+00,9.,43,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.7g',140,47057.004_8,4.2881497D-01,5.8343246D+00,9.,56,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.8g',140,47583.458_8,4.3361237D-01,5.8995965D+00,9.,69,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.9g',140,47946.502_8,4.3692067D-01,5.9446082D+00,9.,84,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.9p',311,47084.000_8,4.2906098D-01,5.8376716D+00,1.,61,"SNIST"),& ! A la main à partir de NIST LEVELS('3p6.4s.10p',311,47603.750_8,4.3379729D-01,5.9021124D+00,1.,74,"SNIST"),& ! A la main à partir de NIST LEVELS('3p6.4s.5g',340,44874.860_8,4.0892981D-01,5.5638829D+00,7,33,"SNIST"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.5g',340,44874.860_8,4.0892981D-01,5.5638829D+00,9,33,"SNIST"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.6g',340,46246.221_8,4.2142657D-01,5.7338002D+00,7.,44,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.6g',340,46246.221_8,4.2142657D-01,5.7338002D+00,9.,44,"STOPB"),& ! + levée de dégénérescence LEVELS('3p6.4s.6g',340,46246.221_8,4.2142657D-01,5.7338002D+00,11.,44,"STOPB"),& ! + levée de dégénérescence LEVELS('3p6.4s.7g',340,47057.257_8,4.2881728D-01,5.8343559D+00,7.,57,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.7g',340,47057.257_8,4.2881728D-01,5.8343559D+00,9.,57,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.7g',340,47057.257_8,4.2881728D-01,5.8343559D+00,11.,57,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.8g',340,47584.314_8,4.3362017D-01,5.8997027D+00,7.,70,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.8g',340,47584.314_8,4.3362017D-01,5.8997027D+00,9.,70,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.8g',340,47584.314_8,4.3362017D-01,5.8997027D+00,11.,70,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.9g',340,47945.833_8,4.3691458D-01,5.9445253D+00,7.,83,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.9g',340,47945.833_8,4.3691458D-01,5.9445253D+00,9.,83,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.9g',340,47945.833_8,4.3691458D-01,5.9445253D+00,11.,83,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.10g',340,48204.462_8,4.3927138D-01,5.9765912D+00,7.,95,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.10g',340,48204.462_8,4.3927138D-01,5.9765912D+00,9.,95,"STOPB"),& ! NIST + levée de dégénérescence LEVELS('3p6.4s.10g',340,48204.462_8,4.3927138D-01,5.9765912D+00,11.,95,"STOPB"), &! NIST + levée de dégénérescence ! ! Ajout de niveaux supplémentaires présents dans NIST mais pas dans la TOPBASE ! LEVELS('3p6.3d.5s',120,47449.083_8,4.3238786D-01,5.8830522D+00,5.,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3p6.3d.5s',320,47456.452_8,4.3245501D-01,5.8839658D+00,3.,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3p6.3d.5s',320,47466.014_8,4.3254215D-01,5.8851514D+00,5.,0,"SNIST"),& ! NIST sans TOPBASE LEVELS('3p6.3d.5s',320,47475.915_8,4.3263237D-01,5.8863790D+00,7.,0,"SNIST"),& ! NIST sans TOPBASE ! ! AJOUT DU 31/05/2011 ! LEVELS('3p6.3d2',140,47759.433_8,4.3521598D-01,5.9214145D+00,9.,76),& ! [Eion(NIST) - En(TOPBASE)]cm-1 ! LEVELS('3p6.4s.6h',151,46257.415_8,4.2152858D-01,5.7351880D+00,11.,45,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.6h',351,46257.415_8,4.2152858D-01,5.7351880D+00,33.,46,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.7h',151,47066.058_8,4.2889748D-01,5.8354470D+00,33.,58,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.7h',351,47066.058_8,4.2889748D-01,5.8354470D+00,11.,59,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.8h',151,47590.965_8,4.3368078D-01,5.9005271D+00,33.,72,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.8h',351,47590.965_8,4.3368078D-01,5.9005271D+00,11.,73,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.9h',151,47950.870_8,4.3696048D-01,5.9451497D+00,11.,85,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.9h',351,47950.870_8,4.3696048D-01,5.9451497D+00,33.,86,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.10h',151,48208.325_8,4.3930658D-01,5.9770700D+00,11.,97,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.4s.10h',351,48208.325_8,4.3930658D-01,5.9770700D+00,33.,98,"STOPB"),& ! [Eion(NIST) - En(TOPBASE)]cm-1 ! ! Formule de polarisation AJOUT DU 31/05/2011 ! LEVELS('3p6.4s.7i',360, 47066.420_8,4.2890077D-01,5.8354918D+00,52.,0,"CN83"),& LEVELS('3p6.4s.8i',360, 47591.310_8,4.3368393D-01,5.9005699D+00,52.,0,"CN83"),& LEVELS('3p6.4s.9i',360, 47951.174_8,4.3696325D-01,5.9451874D+00,52.,0,"CN83"),& LEVELS('3p6.4s.10i',360,48208.582_8,4.3930893D-01,5.9771019D+00,52.,0,"CN83"),& LEVELS('3p6.4s.8k',371, 47591.320_8,4.3368401D-01,5.9005711D+00,60.,0,"CN83"),& LEVELS('3p6.4s.9k',371, 47951.181_8,4.3696331D-01,5.9451882D+00,60.,0,"CN83"),& LEVELS('3p6.4s.10k',371,48208.587_8,4.3930897D-01,5.9771025D+00,60.,0,"CN83"),& LEVELS('3p6.4s.9l',380, 47951.184_8,4.3696334D-01,5.9451886D+00,68.,0,"CN83"),& LEVELS('3p6.4s.10l',380,48208.590_8,4.3930900D-01,5.9771029D+00,68.,0,"CN83"),& LEVELS('3p6.4s.10m',391,48208.591_8,4.3930901D-01,5.9771030D+00,76.,0,"CN83") /) ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! ! Force de collisions effectives de Samson & Berrington (ADNDT Samson & Berrington 2001) ! 10 transitions integer, parameter :: SamsN = 10 integer, dimension(SamsN, 6), parameter :: SamsConf = reshape(source = (/ 4, 0, 100, 4, 1, 311, & 4, 0, 100, 3, 2, 320, & 4, 0, 100, 3, 1, 120, & 4, 0, 100, 4, 1, 111, & 4, 0, 100, 5, 0, 300, & 4, 0, 100, 5, 0, 100, & 4, 0, 100, 4, 1, 331, & 4, 0, 100, 4, 1, 121, & 4, 0, 100, 5, 1, 311, & 4, 0, 100, 5, 1, 111/), & shape = (/ SamsN, 6 /), & order = (/ 2, 1 /) ) Real, dimension(SamsN, 9), parameter :: Sams_Upsilon = reshape( & source = (/ 3.201, 4.183, 4.440, 4.437, 4.340, 4.211, 4.075, 3.943, 3.701, & 1.472, 1.723, 1.926, 2.118, 2.288, 2.433, 2.552, 2.649, 2.789, & .2312, .2972, .3471, .3982, .4536, .5124, .5730, .6342, .7553, & .1879, .3393, .5268, .7315, .9455, 1.165, 1.388, 1.613, 2.068, & .2065, .2867, .3051, .3048, .2990, .2917, .2844, .2776, .2661, & .4552, .5369, .5394, .5269, .5129, .5005, .4903, .4820, .4699, & .0747, .1745, .2761, .3622, .4300, .4819, .5211, .5506, .5889, & .0410, .0729, .0896, .0977, .1024, .1059, .1090, .1122, .1190, & .3644, .3888, .3979, .4002, .3996, .3974, .3945, .3913, .3846, & .1850, .1939, .1952, .1977, .2028, .2104, .2198, .2304, .2533 /), & shape = (/SamsN, 9/), & order = (/2, 1/)) ! Integer, parameter :: App1_N = 18 Real, dimension(App1_N), parameter :: App1_lambda = (/ & 4226.73, & 4425.44, & 4578.55, & 5512.98, & !5588.75, & 7 26 !5590.11, & 5 25 5867.56, & 6102.72, & 6122.22, & 6161.30, & 6162.17, & 6166.44, & 6169.04, & 6169.56, & !6439.07, & 7 15 6455.60, & !6471.66, & 7 13 !6493.78, & 5 12 !6499.65, & 6 12 6572.78, & !7148.15, & 8 14 7326.15, & 8525.72, & 8583.32, & 8633.95 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ & 150, & 30, & 20, & 15, & !15, & !15, & 10, & 25, & 40, & 10, & 40, & 15, & 10, & 15, & !15, & 10, & !15, & !20, & !20, & 20, & !25, & 15, & 10, & 10, & 10 /) Integer, dimension(20) :: GNi = & ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia (/ 1, 1, 2, 2, 3, 4, 5, 5, 6, 6, 6, 7, 9, 9, 9, 12, 13, 15, 14, 14 /) Integer, dimension(20) :: GNj = & (/ 3, 9, 10, 21, 10, 10, 16, 40, 14, 17, 18, 18, 20, 35, 38, 110, 111, 112, 109, 110/) ! End Module datagraph_CaI ! !*********************************************************************** ! Module datagraph_CaII use MTD_STRUCT implicit none Character(len = 2), parameter :: Symb = 'Ca' integer, parameter :: Z_IONI = 2 character(len=4) :: cpage = 'DA4L' !'DA4P' integer,parameter :: lenterm = 6 character(len=32) :: continuum = '11.87' real,dimension(lenterm) :: SLP = (/200,211,220,231,240,251/) real,dimension(lenterm) :: iterm = (/1,2,3,4,5,6/) character(len=7),dimension(lenterm) :: cterm = (/"[2$S ","[2$P[o$","[2$D ","[2$F[o$","[2$G ","[2$H[o$"/) ! Position et taille du graphique integer, parameter :: posx = 250, posy = 1980 integer, parameter :: lenx = 2600, leny = 1800 !integer, parameter :: posx = 200, posy = 2850 !integer, parameter :: lenx = 1800, leny = 2820 ! Abscisses du continuum real, dimension(8), parameter :: abscontinuum = (/0.3,1.,2.,3.,4.,5.,6.,6.5/) ! Noms des niveaux integer,parameter :: levelsize = 12 character(len= 5),dimension(levelsize):: levelname = & (/'4s ','10s','4p ','10p','3d ','10d','4f ','10f','5g ','10g','6h ','10h' /) ! Position des niveaux real,dimension(levelsize,2) :: levelpos = reshape(source=& (/ 1.2,0.1, 1.2,11.18, 2.2,3.24, 2.2,11.3, 3.2,1.79, 3.2,11.40,& 4.2,8.53, 4.2,11.48, 5.2,9.79, 5.2,11.48, 6.2,10.45, 6.2,11.48/),& shape=(/levelsize,2/),& order=(/2,1/)) integer,parameter :: Klevels = 16 type(LEVELS), dimension(Klevels),save :: S_level_sup = (/ & ! Ajout de niveaux supplémentaires pour être consistent LEVELS('3p6.7p',211,81657.10_8,7.4411424D-01,1.0124189D+01,2.,15), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.7p',211,81657.10_8,7.4411424D-01,1.0124189D+01,4.,15),& ! avec levée de dégénérescence LEVELS('3p6.8p',211,85625.94_8,7.8028094D-01,1.0616262D+01,2.,21), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.8p',211,85625.94_8,7.8028094D-01,1.0616262D+01,4.,21),& ! avec levée de dégénérescence LEVELS('3p6.9p',211,88124.53_8,8.0304984D-01,1.0926049D+01,2.,27), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.9p',211,88124.53_8,8.0304984D-01,1.0926049D+01,4.,27),& ! avec levée de dégénérescence LEVELS('3p6.10p',211,89799.71_8,8.1831541D-01,1.1133744D+01,2.,33), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.10p',211,89799.71_8,8.1831541D-01,1.1133744D+01,4.,33),& ! avec levée de dégénérescence LEVELS('3p6.10g',240,91362.27_8,8.3255424D-01,1.1327477D+01,8.,41), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.10g',240,91362.27_8,8.3255424D-01,1.1327477D+01,10.,41),& ! avec levée de dégénérescence LEVELS('3p6.6h',251,83558.85_8,7.6144424D-01,1.0359976D+01,10.,18), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.6h',251,83558.85_8,7.6144424D-01,1.0359976D+01,12.,18),& ! avec levée de dégénérescence LEVELS('3p6.7h',251,86793.72_8,7.9092254D-01,1.0761049D+01,10.,24), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.7h',251,86793.72_8,7.9092254D-01,1.0761049D+01,12.,24), & ! avec levée de dégénérescence LEVELS('3p6.9h',251,90332.75_8,8.2317254D-01,1.1199832D+01,10.,36), & ! [Eion(NIST) - En(TOPBASE)]cm-1 LEVELS('3p6.9h',251,90332.75_8,8.2317254D-01,1.1199832D+01,10.,36) /) ! avec levée de dégénérescence ! Force de collisions de Burgess et Tully (A&A Burgess, Chidichimo, Tully 1995) ! 21 transitions permises (10 permises 11 interdites) integer, parameter :: BurgN = 21 integer, dimension(BurgN, 4), parameter :: BurgConf = reshape(source = (/ 4, 0, 4, 1, & 4, 0, 5, 1, & 3, 2, 4, 1, & 3, 2, 5, 1, & 3, 2, 4, 3, & 4, 1, 5, 0, & 4, 1, 4, 2, & 5, 0, 5, 1, & 4, 2, 5, 1, & 4, 2, 4, 3, & 4, 0, 3, 2, & 4, 0, 5, 0, & 4, 0, 4, 2, & 4, 0, 4, 3, & 3, 2, 5, 0, & 3, 2, 4, 2, & 4, 1, 5, 1, & 4, 1, 4, 3, & 5, 0, 4, 2, & 5, 0, 4, 3, & 5, 1, 4, 3 /), & shape = (/ BurgN, 4 /), & order = (/ 2, 1 /) ) real, dimension(BurgN), parameter :: BurgC = (/ 2.9, 2.5, 5.6, 2.8, 2.0, 2.6, 3.3, 5.6, 7.0, 5.5, & 1.4, 1.5, 0.5, 2.0, 4.4, 2.3, 0.9, 6.5, 3.8, 4.5, 13./) real, dimension(BurgN), parameter :: BurgP1 = (/ 21.96, 1.247, 43.42, 3.137, 3.208, 7.667, 30.19, 53.82, 154.7, 178.9, & 2.560, 2.601, 4.398, 1.944, 0.9722, 30.18, 8.586, 11.53, 29.66, 12.80, 32.88/) real, dimension(BurgN), parameter :: BurgP2 = (/ 29.53, 1.604, 33.18, 2.375, 4.757, 5.636, 32.30, 85.39, 173.8, 264.8, & 2.954, 1.778, 4.335, 2.108, 0.7578, 25.24, 10.34, 19.37, 32.41, 15.41, 69.61/) real, dimension(BurgN), parameter :: BurgP3 = (/ 32.35, 1.225, 29.32, 1.516, 5.973, 7.140, 45.41, 119.5, 202.3, 330.2, & 3.762, 1.981, 4.153, 1.935, 0.6744, 22.07, 12.40, 21.89, 37.08, 16.04, 90.66/) real, dimension(BurgN), parameter :: BurgP4 = (/ 37.43, 0.5296, 27.86, 0.6059, 7.892, 11.49, 64.65, 152.6, 226.4, 403.4, & 4.339, 2.261, 4.551, 1.754, 0.5929, 19.80, 15.62, 24.57, 41.75, 16.02, 111.7/) real, dimension(BurgN), parameter :: BurgP5 = (/ 37.56, 6.531E-4, 24.18, 7.291E-2,12.41, 17.08, 79.09, 157.9,211.2, 424.1, & 5.023, 3.382, 6.024, 1.198, 0.4160, 16.97, 20.72, 22.18, 47.00, 14.88, 124.4/) ! Fraction des forces de collisions effectives pour les composantes d'un multiplet 26(+ 6 interdites) * K_special, gi, gj Integer, dimension(32, 3), parameter :: BurgFracIndex = RESHAPE(source = (/ 1, 2, 2, & 1, 2, 4, & 2, 2, 2, & 2, 2, 4, & 3, 4, 2, & 3, 4, 4, & 3, 6, 4, & 3, 6, 2, & ! transition interdite 4, 4, 2, & 4, 4, 4, & 4, 6, 4, & 4, 6, 2, & ! transition interdite 5, 4, 6, & 5, 6, 6, & 5, 6, 8, & 5, 4, 8, & ! transition interdite 6, 2, 2, & 6, 4, 2, & 7, 2, 4, & 7, 4, 4, & 7, 4, 6, & 7, 2, 6, & ! transition interdite 8, 2, 2, & 8, 2, 4, & 9, 4, 2, & 9, 4, 4, & 9, 6, 4, & 9, 6, 2, & ! transition interdite 10, 4, 6, & 10, 6, 6, & 10, 6, 8, & 10, 4, 8 /), & ! transition interdite SHAPE = (/ 32, 3 /), & ORDER = (/ 2, 1 /) ) ! Fractions basés sur la somme des gf pour chaque multiplet (source des gf : VALD) + ! même fraction pour la transition interdite issue du même niveau bas (que la transition permise) quand il y a lieu Real, dimension(32), parameter :: BurgFrac = (/30.8, 69.2, & ! 100% de Upsilon(K_special = 1) 2.2, 97.8, & ! 100% de Upsilon(K_special = 2) 41.3, 8.3, 50.4, 50.4, & ! 100% de Upsilon(K_special = 3) + 1 interdite 39.8, 8.7, 51.5, 51.5, & ! 100% de Upsilon(K_special = 4) + 1 interdite 50.0, 2.4, 47.6, 50.0, & ! 100% de Upsilon(K_special = 5) + 1 interdite 49.7, 50.3, & ! 100% de Upsilon(K_special = 6) 49.8, 5.0, 45.1, 49.8, & ! 100% de Upsilon(K_special = 7) + 1 interdite 33.2, 66.8, & ! 100% de Upsilon(K_special = 8) 39.7, 8.6, 51.7, 51.7, & ! 100% de Upsilon(K_special = 9) + 1 interdite 45.5, 2.2, 52.3, 45.5 /) ! 100% de Upsilon(K_special = 10) + 1 interdite ! ! SUPPLEMENTARY LINES ! Integer, parameter :: Klines = 1 Type(LINES), dimension(Klines), parameter :: S_line_sup = & (/ LINES('spdfghiklm', -1, 1.797D+308, 3.402E+38, & 'spdfghiklm', -1, 1.797D+308, 3.402E+38, & 1.797D+308, 3.402E+38, 3.402E+38, 3.402E+38, 3.402E+38 ) /) ! Transitions concernées par Application 1 Integer, parameter :: App1_N = 9 Real, dimension(App1_N), parameter :: App1_lambda = (/ & 3933.66, & 3968.47, & 8248.79, & 8498.02, & 8542.09, & 8662.14, & 8912.07, & 9854.76, & 11949.75 /) Real, dimension(App1_N), parameter :: App1_QMAX = (/ & 700, & 700, & 15, & 400, & 400, & 400, & 25, & 15, & 25 /) Integer, dimension(9) :: GNi = & ! vecteur d'entiers représentant les niveaux bas/hauts des transistions indexées pour Gaia (/ 1, 1, 10, 2, 3, 2, 8, 10, 6 /) Integer, dimension(9) :: GNj = & (/ 5, 4, 15, 5, 5, 4, 12, 13, 10 /) End Module datagraph_CaII ! !*********************************************************************** ! Module DATAGRAPH ! USE DATAGRAPH_NaI ! End Module DATAGRAPH