c**********************************************************************
c*
c* Einstein coefficients
c*
c* fgg mar 2010 (adapted from MALV get_a subroutine)
c*
c**********************************************************************
c                                                
! isotope  626		for this isotope the values are from mlp    
c
      real*8 a1_103, a1_102, a1_101, a1_023, a1_022, a1_021
c
  parameter (a1_103       = 1260.d0)
  parameter (a1_102       = 820.d0)
  parameter (a1_101       = 410.d0)
  parameter (a1_023       = 1220.d0)
  parameter (a1_022       = 814.d0)
  parameter (a1_021       = 406.d0)
c
c
      real*8 	a1_003, a1_002, a1_301, a1_221, a1_141, a1_061
c
  parameter (a1_003       = 1200.d0)
  parameter (a1_002       = 780.d0)
  parameter (a1_301       = 400.d0)
  parameter (a1_221       = 400.d0)
  parameter (a1_141       = 400.d0)
  parameter (a1_061       = 400.d0)
c
c
      real*8 a1_201, a1_121, a1_041, a1_220, a1_130, a1_050
c
  parameter (a1_201       = 424.d0)
  parameter (a1_121       = 418.d0)
  parameter (a1_041       = 418.d0)
  parameter (a1_220       = 1.43d0)
  parameter (a1_130       = 2.3d0)
  parameter (a1_050       = 7.15d0)
c
c
      real*8 a1_120, a1_55, a1_44, a1_300, a1_a43, a1_210, a1_200
c
  parameter (a1_120       = 2.08d0)
  parameter (a1_55        = (5.67d0+a1_120)/2.d0)
  parameter (a1_44        = 3.68d0)
  parameter (a1_300       = 4.3d0)
  parameter (a1_a43       = 2.44d0)
  parameter (a1_210       = a1_a43+a1_44)
  parameter (a1_200       = 4.53d0)
c
c
      real*8 a1_140, a1_060, a1_0510, a1_040,  a1_0420, a1_1220
c
  parameter (a1_140       = 2.27d0)
  parameter (a1_060       = 2.32d0)
  parameter (a1_0510      = a1_0510_0420+a1_0510_0400)
  parameter (a1_040       = 2.35d0)
  parameter (a1_0420      = 3.3d0)
  parameter (a1_1220      = 4.71d0)
c
c
      real*8 a1_1111_1110,a1_1111_0110, a1_0311_0310,a1_0311_0110
c
  parameter (a1_1111_1110 = 367.0d0)
  parameter (a1_1111_0110 = 17.0d0)
  parameter (a1_0311_0310 = 370.8d0)
  parameter (a1_0311_0110 = 9.65d0)
c
c
      real*8 a1_0111_010, a1_0201_0200, a1_0201_000, a1_0221_0220
c
  parameter (a1_0111_010  = 388.4d0)
  parameter (a1_0201_000  = 10.39d0)
  parameter (a1_0201_0200 = 383.756d0)
  parameter (a1_0221_0220 = 386.56d0)
c
c
      real*8 a1_1001_1000, a1_1001_000, a1_001_000,a1_010_000
c
  parameter (a1_1001_1000 = 381.222d0)
  parameter (a1_1001_000  = 16.68d0)
  parameter (a1_001_000   = 389.34d0)
  parameter (a1_010_000   = 1.3546d0)
c
c
      real*8 a1_0421_1220, a1_0421_0420, a1_0421_000
c
  parameter (a1_0421_000  = 0.66d0)
  parameter (a1_0421_0420 = 360.5d0)
  parameter (a1_0421_1220 = 0.129d0)
c
c
      real*8 a1_1221_0111, a1_1221_1220, a1_1221_0420
c
  parameter (a1_1221_0420 = 0.01227d0)
  parameter (a1_1221_1220 = 355.2d0)
  parameter (a1_1221_0111 = 0.0219d0)
c
c
      real*8 a1_0510_000, a1_0510_0420, a1_0510_0400
c
  parameter (a1_0510_000  = 4.d-4)
  parameter (a1_0510_0420 = 1.086d0)
  parameter (a1_0510_0400 = 2.406d0)
c
c
      real*8 a1_0510_1000, a1_0510_0220, a1_0510_0200
  parameter (a1_0510_1000 = 3.44d-5)
  parameter (a1_0510_0220 = 6.339d-4)
  parameter (a1_0510_0200 = 1.5d-3)
c
c
      real*8 a1_1310_000 , a1_1310_0420 , a1_1310_0400
c
  parameter (a1_1310_000  = 4.45d-4)
  parameter (a1_1310_0420 = 0.7535d0)
  parameter (a1_1310_0400 = 0.21d0)
c
c
      real*8 a1_1310_1200 , a1_1310_1220 , a1_1310_2000
c
  parameter (a1_1310_1200 = 2.264d0)
  parameter (a1_1310_1220 = 0.3609d0)
  parameter (a1_1310_2000 = 0.074d0)
c
c
      real*8 a1_1310_0111,a1_1310_1000,a1_1310_0220,a1_1310_0200 
c
  parameter (a1_1310_0111 = 1.64d-4)
  parameter (a1_1310_1000 = 1.433d-3)
  parameter (a1_1310_0220 = 2.4d-4)
  parameter (a1_1310_0200 = 9.719d-3)
c
c
      real*8 a1_2110_000 , a1_2110_0420 , a1_2110_0400
c
  parameter (a1_2110_000  = 3.265d-4)
  parameter (a1_2110_0420 = 6.48d-3)
  parameter (a1_2110_0400 = 1.54d-3)
c
c
      real*8 a1_2110_1200,a1_2110_1220,a1_2110_2000,a1_2110_0111
c
  parameter (a1_2110_1200 = 0.0738d0)
  parameter (a1_2110_1220 = 1.957d0)
  parameter (a1_2110_2000 = 3.289d0)
  parameter (a1_2110_0111 = 2.68d-4)
c
c
      real*8 a1_2110_1000 , a1_2110_0220 , a1_2110_0200
c
  parameter (a1_2110_1000 = 3.445d-2)
  parameter (a1_2110_0220 = 1.976d-2)
  parameter (a1_2110_0200 = 6.539d-3)
c
c
      real*8 a1_0403_0402, a1_0403_0202
c
  parameter (a1_0403_0402 = 728.7d0)
  parameter (a1_0403_0202 = 28.9d0)
c
      real*8 a1_1203_1202, a1_1203_1002, a1_1203_0202
c
c
  parameter (a1_1203_1202 = 716.3d0)
  parameter (a1_1203_1002 = 27.25d0)
  parameter (a1_1203_0202 = 43.68d0)
c
c
      real*8 a1_2003_2002, a1_2003_1002, a1_2003_0202
c
  parameter (a1_2003_2002 = 995.7d0)
  parameter (a1_2003_1002 = 83.27d0)
  parameter (a1_2003_0202 = 1.15d0)
c
c
      real*8 a1_0402_0401, a1_0402_0201
c
  parameter (a1_0402_0401 = 700.0d0)
  parameter (a1_0402_0201 = 28.97d0)
c
c
      real*8 a1_1202_1201, a1_1202_1001, a1_1202_0201
c
  parameter (a1_1202_1201 = 687.0d0)  
  parameter (a1_1202_1001 =  26.07d0)
  parameter (a1_1202_0201 =  39.57d0)
c
c
      real*8 a1_2002_2001, a1_2002_1001
c
  parameter (a1_2002_2001 = 686.1d0)
  parameter (a1_2002_1001 = 54.67d0)