(file) Return to hcf2r.f CVS log (file) (dir) Up to [HallC] / pol_hms_single

Diff for /pol_hms_single/hcf2r.f between version 1.2 and 1.3

version 1.2, 2003/12/19 13:37:28 version 1.3, 2005/04/01 16:51:48
Line 2 
Line 2 
       subroutine hcf2r(qq,ww,f2,r)       subroutine hcf2r(qq,ww,f2,r)
  
 c      include 'model_liang.cmn' c      include 'model_liang.cmn'
       implicit real*8 (a-h,o-z)        implicit none
       integer i,j,j1,i1,tmp       integer i,j,j1,i1,tmp
       real*8 qq,ww,w1,w2,ww1,ww2  !,w11,w22       real*8 qq,ww,w1,w2,ww1,ww2  !,w11,w22
       real*8 r,f2,nu,mp,mp2,r1,f21,r2,f22       real*8 r,f2,nu,mp,mp2,r1,f21,r2,f22
Line 22 
Line 22 
  
  
 c      write(*,*) qq,ww c      write(*,*) qq,ww
       if (ww.lt.1.15) then        if (ww.lt.1.15d0) then
          w1 = 0.d0          w1 = 0.d0
          w2 = 0.d0          w2 = 0.d0
       else if (ww.gt.3.9d0) then       else if (ww.gt.3.9d0) then
Line 36 
Line 36 
  
  
       data rp/       data rp/
      >   0.33129,   1.32645,  -0.79900,  -3.11559,    2.5681,  -1.82534,       > 0.33129d0, 1.32645d0,-0.79900d0,  -3.11559d0,    2.5681d0,  -1.82534d0,
      >   0.28813,   0.90973,  -0.54869,  -3.76651,    3.4846,  -2.02225,       > 0.28813d0, 0.90973d0,-0.54869d0,  -3.76651d0,    3.4846d0,  -2.02225d0,
      >   0.56382,   0.43268,   0.08276,  23.37011,  -25.4063,  13.71927,       > 0.56382d0, 0.43268d0, 0.08276d0,  23.37011d0,  -25.4063d0,  13.71927d0,
      >   0.67816,  -0.09742,   1.91168,  51.09415,  -61.0272,  17.29566,       > 0.67816d0,-0.09742d0, 1.91168d0,  51.09415d0,  -61.0272d0,  17.29566d0,
      >   0.64090,   0.51581,   0.28952, -69.99436,  -98.0096, 133.50420,       > 0.64090d0, 0.51581d0, 0.28952d0, -69.99436d0,  -98.0096d0, 133.50420d0,
      >   0.46245,  -0.02240,   0.31779,  -3.26504,   77.9481,  14.49329,       > 0.46245d0,-0.02240d0, 0.31779d0,  -3.26504d0,   77.9481d0,  14.49329d0,
      >   0.39436,  -0.02232,   0.34328,  24.69035,   50.0479,  13.91268,       > 0.39436d0,-0.02232d0, 0.34328d0,  24.69035d0,   50.0479d0,  13.91268d0,
      >   0.39805,   0.03919,   0.07716,   3.75484,   10.3324,   8.08673,       > 0.39805d0, 0.03919d0, 0.07716d0,   3.75484d0,   10.3324d0,   8.08673d0,
      >   0.38997,   0.00900,   0.15343,  16.36940,    4.9903,   7.44212,       > 0.38997d0, 0.00900d0, 0.15343d0,  16.36940d0,    4.9903d0,   7.44212d0,
      >   0.29674,   0.02783,   0.06118,  27.36369,  -17.4484,   4.32732,       > 0.29674d0, 0.02783d0, 0.06118d0,  27.36369d0,  -17.4484d0,   4.32732d0,
      >   0.35439,  -0.01126,   0.23150,  53.46022,  -18.2759,   7.18313,       > 0.35439d0,-0.01126d0, 0.23150d0,  53.46022d0,  -18.2759d0,   7.18313d0,
      >   0.31987,   0.01455,   0.23862,  81.55396,  -74.7982,   5.96615,       > 0.31987d0, 0.01455d0, 0.23862d0,  81.55396d0,  -74.7982d0,   5.96615d0,
      >   0.31163,  -0.00699,   0.12995,  51.23706,   18.1217,   4.82674,       > 0.31163d0,-0.00699d0, 0.12995d0,  51.23706d0,   18.1217d0,   4.82674d0,
      >   0.32176,   0.12831,   0.00514,  68.39391,  -73.3401,   3.86676,       > 0.32176d0, 0.12831d0, 0.00514d0,  68.39391d0,  -73.3401d0,   3.86676d0,
      >   0.32514,   0.20568,   0.02398,  31.46726,  -34.1291,   3.42527,       > 0.32514d0, 0.20568d0, 0.02398d0,  31.46726d0,  -34.1291d0,   3.42527d0,
      >   0.48734,   0.32263,   0.07471,  42.71339,  -49.2712,   7.19561,       > 0.48734d0, 0.32263d0, 0.07471d0,  42.71339d0,  -49.2712d0,   7.19561d0,
      >   0.38366,   0.16378,   0.00110,  26.58418,  -29.0223,   3.80760,       > 0.38366d0, 0.16378d0, 0.00110d0,  26.58418d0,  -29.0223d0,   3.80760d0,
      >   0.52817,  -0.37320,  10.22177,  44.87227,  -53.3864,  15.03388,       > 0.52817d0,-0.37320d0,10.22177d0,  44.87227d0,  -53.3864d0,  15.03388d0,
      >   0.34971,   0.24432,   0.03942,  15.46841,  -17.0274,   2.87509,       > 0.34971d0, 0.24432d0, 0.03942d0,  15.46841d0,  -17.0274d0,   2.87509d0,
      >   0.32309,   0.05971,   0.23399,  81.18824,  -85.6957,   4.59427,       > 0.32309d0, 0.05971d0, 0.23399d0,  81.18824d0,  -85.6957d0,   4.59427d0,
      >   0.48408,   0.14606,  -0.12127,  78.59016,  -92.2077,   2.15488,       > 0.48408d0, 0.14606d0,-0.12127d0,  78.59016d0,  -92.2077d0,   2.15488d0,
      >   0.44745,   0.08725,  -0.07341,  96.91505, -113.3328,   1.61255,       > 0.44745d0, 0.08725d0,-0.07341d0,  96.91505d0, -113.3328d0,   1.61255d0,
      >   0.17621,   0.00267,   0.16446, 290.88310, -284.5832,   4.78794,       > 0.17621d0, 0.00267d0, 0.16446d0, 290.88310d0, -284.5832d0,   4.78794d0,
      >   0.19235,   0.03305,   0.26250,  86.32166,  -86.1215,   4.45001,       > 0.19235d0, 0.03305d0, 0.26250d0,  86.32166d0,  -86.1215d0,   4.45001d0,
      >   0.25966,   0.02309,   0.06332,  54.29687,  -49.3171,   3.33319,       > 0.25966d0, 0.02309d0, 0.06332d0,  54.29687d0,  -49.3171d0,   3.33319d0,
      >   0.35243,   0.02105,  -0.18868, 110.69520, -226.4427,   4.95000,       > 0.35243d0, 0.02105d0,-0.18868d0, 110.69520d0, -226.4427d0,   4.95000d0,
      >   0.36534,  -0.21314,   2.81021, -59.05300,  -37.7424,  31.42299,       > 0.36534d0,-0.21314d0, 2.81021d0, -59.05300d0,  -37.7424d0,  31.42299d0,
      >   0.28897,  -5.27934,  62.75448,  15.15152,    5.7736,  29.82406,       > 0.28897d0,-5.27934d0,62.75448d0,  15.15152d0,    5.7736d0,  29.82406d0,
      >   0.31192,  -0.03546,   0.09286,  19.06862,  -33.4685,   1.56198,       > 0.31192d0,-0.03546d0, 0.09286d0,  19.06862d0,  -33.4685d0,   1.56198d0,
      >   0.35713,  -0.20097,   2.05940, 436.08360, -517.4285,  14.75464,       > 0.35713d0,-0.20097d0, 2.05940d0, 436.08360d0, -517.4285d0,  14.75464d0,
      >   0.32423,   0.05549,   0.04984,  74.70284,  -85.5771,   3.61358,       > 0.32423d0, 0.05549d0, 0.04984d0,  74.70284d0,  -85.5771d0,   3.61358d0,
      >   0.26321,   0.02532,   0.04165, 291.81850, -330.3480,   3.42817,       > 0.26321d0, 0.02532d0, 0.04165d0, 291.81850d0, -330.3480d0,   3.42817d0,
      >   0.19235,  -0.00521,   0.30864, 262.23130, -271.3294,   4.52254,       > 0.19235d0,-0.00521d0, 0.30864d0, 262.23130d0, -271.3294d0,   4.52254d0,
      >   0.17858,  -0.00018,   0.28614, 195.28650, -195.6441,   3.92928,       > 0.17858d0,-0.00018d0, 0.28614d0, 195.28650d0, -195.6441d0,   3.92928d0,
      >   0.14180,   0.00801,   0.21545, 131.05210, -121.0251,   3.20372,       > 0.14180d0, 0.00801d0, 0.21545d0, 131.05210d0, -121.0251d0,   3.20372d0,
      >   0.19810,   0.00240,   0.29005, 122.32390, -121.0733,   3.60232,       > 0.19810d0, 0.00240d0, 0.29005d0, 122.32390d0, -121.0733d0,   3.60232d0,
      >   0.24832,  -0.00999,   0.46397,  81.22332,  -85.5181,   3.40233,       > 0.24832d0,-0.00999d0, 0.46397d0,  81.22332d0,  -85.5181d0,   3.40233d0,
      >   0.31740,   0.12344,   0.07026,  53.01398,  -64.9467,   3.08821,       > 0.31740d0, 0.12344d0, 0.07026d0,  53.01398d0,  -64.9467d0,   3.08821d0,
      >   0.32739,  -0.06064,   0.79821, 104.36790, -124.3242,   4.74774,       > 0.32739d0,-0.06064d0, 0.79821d0, 104.36790d0, -124.3242d0,   4.74774d0,
      >   0.34597,   0.11624,  -0.04088,  50.29695,  -64.8436,   2.35262,       > 0.34597d0, 0.11624d0,-0.04088d0,  50.29695d0,  -64.8436d0,   2.35262d0,
      >   0.33840,   0.11035,   0.03128,  39.27810,  -51.3317,   2.38903,       > 0.33840d0, 0.11035d0, 0.03128d0,  39.27810d0,  -51.3317d0,   2.38903d0,
      >   0.38341,   0.12239,  -0.12886,  59.12659,  -80.0135,   1.94849,       > 0.38341d0, 0.12239d0,-0.12886d0,  59.12659d0,  -80.0135d0,   1.94849d0,
      >   0.34759,  -0.06730,   0.99427,  87.77917, -116.0158,   4.31371,       > 0.34759d0,-0.06730d0, 0.99427d0,  87.77917d0, -116.0158d0,   4.31371d0,
      >   0.35035,   0.13825,   0.05372,  37.00152,  -51.0832,   2.35058,       > 0.35035d0, 0.13825d0, 0.05372d0,  37.00152d0,  -51.0832d0,   2.35058d0,
      >   0.31877,   0.09654,   0.12551,  70.29255,  -93.9217,   2.72640,       > 0.31877d0, 0.09654d0, 0.12551d0,  70.29255d0,  -93.9217d0,   2.72640d0,
      >   0.32985,   0.25111,  -0.24644,  11.54706,  -17.2718,   1.00506/       > 0.32985d0, 0.25111d0,-0.24644d0,  11.54706d0,  -17.2718d0,   1.00506d0/
  
       data f2p/       data f2p/
      >   0.00270,  -0.02563,   0.06438,   0.13288,  -0.18442,  59.57852,       >   0.00270d0,  -0.02563d0,  0.06438d0,   0.13288d0,  -0.18442d0,  59.57852d0,
      > 126.68970,-156.52300,  259.6862,  45.71608,  51.94866,-337.30260,       > 126.68970d0,-156.52300d0,259.6862d0,   45.71608d0,  51.94866d0,-337.30260d0,
      >  -0.01266,   0.53825,  -1.14310,  -1.10755,   2.67523, -96.68052,       >  -0.01266d0,   0.53825d0, -1.14310d0,  -1.10755d0,   2.67523d0, -96.68052d0,
      > -65.71029, 144.62170, -350.7377,  61.76398, 116.55460, 172.52990,       > -65.71029d0, 144.62170d0,-350.7377d0,  61.76398d0, 116.55460d0, 172.52990d0,
      >  -0.00820,   0.30865,  -0.10750,  -0.79825,   1.01081, -37.89112,       >  -0.00820d0,   0.30865d0, -0.10750d0,  -0.79825d0,   1.01081d0, -37.89112d0,
      > -99.42540, 120.79090, -161.3490, -27.73087, -37.25240, 223.31690,       > -99.42540d0, 120.79090d0,-161.3490d0, -27.73087d0, -37.25240d0, 223.31690d0,
      >  -0.01626,   0.68263,  -2.08330,   0.93416,   7.66115, -20.34949,       >  -0.01626d0,   0.68263d0, -2.08330d0,   0.93416d0,   7.66115d0, -20.34949d0,
      > -89.79055, 100.32260,  -89.9365, -81.08472,  33.33122, 139.55320,       > -89.79055d0, 100.32260d0, -89.9365d0, -81.08472d0,  33.33122d0, 139.55320d0,
      >   0.10085,  -0.39180,   4.39375,   2.96299,  -9.43805,   5.53277,       >   0.10085d0,  -0.39180d0,  4.39375d0,   2.96299d0,  -9.43805d0,   5.53277d0,
      > -49.30489,  43.73740,   12.3122,-123.36040, 117.44570,  -2.55460,       > -49.30489d0,  43.73740d0,  12.3122d0,-123.36040d0, 117.44570d0,  -2.55460d0,
      >   0.16801,   1.66346,  -3.58292,  -1.43295,   6.54949,   9.59234,       >   0.16801d0,   1.66346d0, -3.58292d0,  -1.43295d0,   6.54949d0,   9.59234d0,
      > -15.43190,   7.56539,   32.3811, -79.52365,  79.85619, -30.67446,       > -15.43190d0,   7.56539d0,  32.3811d0, -79.52365d0,  79.85619d0, -30.67446d0,
      >   0.24036,   0.66839,  -1.02483,   2.73454,  -3.12791,   4.67220,       >   0.24036d0,   0.66839d0, -1.02483d0,   2.73454d0,  -3.12791d0,   4.67220d0,
      >  -1.94161,  -1.02171,   14.8272, -21.21147,  20.23464, -11.97409,       >  -1.94161d0,  -1.02171d0,  14.8272d0, -21.21147d0,  20.23464d0, -11.97409d0,
      >   0.14395,  -0.24303,   0.47600,   5.20955,   4.07727,  -6.50083,       >   0.14395d0,  -0.24303d0,  0.47600d0,   5.20955d0,   4.07727d0,  -6.50083d0,
      > -59.81225,  63.16642,  -30.5782, -70.20249,  50.86131,  55.80522,       > -59.81225d0,  63.16642d0, -30.5782d0, -70.20249d0,  50.86131d0,  55.80522d0,
      >   0.19595,   0.49078,   0.09058,   2.08639,   0.61890,   3.84467,       >   0.19595d0,   0.49078d0,  0.09058d0,   2.08639d0,   0.61890d0,   3.84467d0,
      >  30.61177, -30.33699,   17.4414,  32.96692, -20.06578, -29.88399,       >  30.61177d0, -30.33699d0,  17.4414d0,  32.96692d0, -20.06578d0, -29.88399d0,
      >   0.16363,   0.21846,   0.20264,   2.53441,  -0.76028,   3.90782,       >   0.16363d0,   0.21846d0,  0.20264d0,   2.53441d0,  -0.76028d0,   3.90782d0,
      >  36.58707, -35.68809,   18.2261,  40.25367, -24.84555, -33.50865,       >  36.58707d0, -35.68809d0,  18.2261d0,  40.25367d0, -24.84555d0, -33.50865d0,
      >   0.11976,  -0.45730,   1.08627,   2.63121,  -0.58440,  -5.72117,       >   0.11976d0,  -0.45730d0,  1.08627d0,   2.63121d0,  -0.58440d0,  -5.72117d0,
      > -53.80934,  57.73899,  -25.7978, -52.77444,  37.38223,  48.12931,       > -53.80934d0,  57.73899d0, -25.7978d0, -52.77444d0,  37.38223d0,  48.12931d0,
      >   0.03129,   0.08877,   0.31201,   0.52214,   0.56099,   9.45694,       >   0.03129d0,   0.08877d0,  0.31201d0,   0.52214d0,   0.56099d0,   9.45694d0,
      >  80.42830, -77.88046,   41.7567,  69.77966, -24.70381, -89.21936,       >  80.42830d0, -77.88046d0,  41.7567d0,  69.77966d0, -24.70381d0, -89.21936d0,
      >   0.04149,   0.05422,   0.14645,   0.46049,   0.70058,   7.52857,       >   0.04149d0,   0.05422d0,  0.14645d0,   0.46049d0,   0.70058d0,   7.52857d0,
      >  73.36507, -69.99093,   33.8845,  68.44820, -34.26609, -69.93432,       >  73.36507d0, -69.99093d0,  33.8845d0,  68.44820d0, -34.26609d0, -69.93432d0,
      >   0.03743,   0.05538,   0.06989,   0.42875,   0.74432,   9.91687,       >   0.03743d0,   0.05538d0,  0.06989d0,   0.42875d0,   0.74432d0,   9.91687d0,
      >  99.10915, -95.21824,   44.3231,  88.44396, -44.56128, -92.21440,       >  99.10915d0, -95.21824d0,  44.3231d0,  88.44396d0, -44.56128d0, -92.21440d0,
      >   0.04214,   0.07361,  -0.05129,   0.40464,   1.03509,   8.78332,       >   0.04214d0,   0.07361d0, -0.05129d0,   0.40464d0,   1.03509d0,   8.78332d0,
      >  90.67270, -86.53600,   39.0777,  79.55918, -42.40906, -79.76051,       >  90.67270d0, -86.53600d0,  39.0777d0,  79.55918d0, -42.40906d0, -79.76051d0,
      >   0.02977,   0.24989,  -0.36480,   0.43282,   1.66379,  10.58702,       >   0.02977d0,   0.24989d0, -0.36480d0,   0.43282d0,   1.66379d0,  10.58702d0,
      >  99.35179, -95.66632,   45.5028,  74.12631, -26.76438, -97.86623,       >  99.35179d0, -95.66632d0,  45.5028d0,  74.12631d0, -26.76438d0, -97.86623d0,
      >   0.04183,   0.17483,  -0.23191,   0.45452,   1.41034,   8.17936,       >   0.04183d0,   0.17483d0, -0.23191d0,   0.45452d0,   1.41034d0,   8.17936d0,
      >  74.28743, -69.96535,   34.5675,  50.35936, -11.35428, -76.36961,       >  74.28743d0, -69.96535d0,  34.5675d0,  50.35936d0, -11.35428d0, -76.36961d0,
      >   0.05173,   0.04448,   0.13424,   0.38647,   0.54959,   6.16662,       >   0.05173d0,   0.04448d0,  0.13424d0,   0.38647d0,   0.54959d0,   6.16662d0,
      >  51.08549, -46.36961,   25.3451,  28.23540,   3.70411, -58.01252,       >  51.08549d0, -46.36961d0,  25.3451d0,  28.23540d0,   3.70411d0, -58.01252d0,
      >   0.08047,   0.07636,  -0.05934,   0.34274,   0.77572,   4.27257,       >   0.08047d0,   0.07636d0, -0.05934d0,   0.34274d0,   0.77572d0,   4.27257d0,
      >  38.90837, -35.54622,   17.7386,  24.76801,  -9.17321, -33.32258,       >  38.90837d0, -35.54622d0,  17.7386d0,  24.76801d0,  -9.17321d0, -33.32258d0,
      >   0.08916,   0.08720,  -0.04772,   0.36316,   0.77081,   4.10228,       >   0.08916d0,   0.08720d0, -0.04772d0,   0.36316d0,   0.77081d0,   4.10228d0,
      >  40.26231, -37.09061,   17.1693,  26.59951, -14.07698, -29.96852,       >  40.26231d0, -37.09061d0,  17.1693d0,  26.59951d0, -14.07698d0, -29.96852d0,
      >   0.02688,   0.20865,   0.24902,   0.19521,   0.37927,  -0.57453,       >   0.02688d0,   0.20865d0,  0.24902d0,   0.19521d0,   0.37927d0,  -0.57453d0,
      >  -9.00651,  14.38374,   -2.4497,  -6.51968,   9.73428,   3.89960,       >  -9.00651d0,  14.38374d0,  -2.4497d0,  -6.51968d0,   9.73428d0,   3.89960d0,
      >   0.04139,   0.15301,   0.32248,   0.36027,   0.17856,   3.32167,       >   0.04139d0,   0.15301d0,  0.32248d0,   0.36027d0,   0.17856d0,   3.32167d0,
      >  28.88752, -23.90850,   13.4160,  12.82948,   3.90580, -29.60664,       >  28.88752d0, -23.90850d0,  13.4160d0,  12.82948d0,   3.90580d0, -29.60664d0,
      >   0.05169,  -0.19792,   0.70477,   0.98245,   0.19970,   1.07366,       >   0.05169d0,  -0.19792d0,  0.70477d0,   0.98245d0,   0.19970d0,   1.07366d0,
      >  -0.96632,  13.46862,    3.5916, -10.17772,  41.48280, -29.53078,       >  -0.96632d0,  13.46862d0,   3.5916d0, -10.17772d0,  41.48280d0, -29.53078d0,
      >   0.04810,  -0.14599,   0.44319,   0.86999,   0.59196,   2.37207,       >   0.04810d0,  -0.14599d0,  0.44319d0,   0.86999d0,   0.59196d0,   2.37207d0,
      >  16.49698,  -3.58537,    9.1814,   2.29473,  31.25257, -39.15542,       >  16.49698d0,  -3.58537d0,   9.1814d0,   2.29473d0,  31.25257d0, -39.15542d0,
      >   0.04731,  -0.14069,   0.33899,   0.74002,   0.62827,  -0.83045,       >   0.04731d0,  -0.14069d0,  0.33899d0,   0.74002d0,   0.62827d0,  -0.83045d0,
      > -17.17104,  32.13834,   -3.8414, -13.50958,  37.48258, -11.74448,       > -17.17104d0,  32.13834d0,  -3.8414d0, -13.50958d0,  37.48258d0, -11.74448d0,
      >   0.05221,  -0.11177,   0.27687,   0.68277,   0.62925,   2.10930,       >   0.05221d0,  -0.11177d0,  0.27687d0,   0.68277d0,   0.62925d0,   2.10930d0,
      >  13.86561,  -1.66573,    7.9977,   0.35896,  27.99104, -32.29613,       >  13.86561d0,  -1.66573d0,   7.9977d0,   0.35896d0,  27.99104d0, -32.29613d0,
      >   0.06442,  -0.07257,   0.11952,   0.50396,   0.53317,   4.29581,       >   0.06442d0,  -0.07257d0,  0.11952d0,   0.50396d0,   0.53317d0,   4.29581d0,
      >  44.07500, -35.61928,   17.2312,  19.12265,   0.86142, -37.08731,       >  44.07500d0, -35.61928d0,  17.2312d0,  19.12265d0,   0.86142d0, -37.08731d0,
      >   0.07351,  -0.10911,   0.12156,   0.46602,   0.48331,   5.17804,       >   0.07351d0,  -0.10911d0,  0.12156d0,   0.46602d0,   0.48331d0,   5.17804d0,
      >  58.31429, -50.40657,   21.0203,  27.73335, -10.71966, -39.46624,       >  58.31429d0, -50.40657d0,  21.0203d0,  27.73335d0, -10.71966d0, -39.46624d0,
      >   0.06943,  -0.07336,   0.15002,   0.46675,   0.51645,   3.31844,       >   0.06943d0,  -0.07336d0,  0.15002d0,   0.46675d0,   0.51645d0,   3.31844d0,
      >  37.30048, -28.68561,   13.4842,  17.16954,  -2.77701, -26.91567,       >  37.30048d0, -28.68561d0,  13.4842d0,  17.16954d0,  -2.77701d0, -26.91567d0,
      >   0.05361,  -0.02633,   0.19376,   0.44691,   0.44814,   2.88220,       >   0.05361d0,  -0.02633d0,  0.19376d0,   0.44691d0,   0.44814d0,   2.88220d0,
      >  29.95785, -19.93872,   11.4738,  10.83239,   7.47675, -28.12408,       >  29.95785d0, -19.93872d0,  11.4738d0,  10.83239d0,   7.47675d0, -28.12408d0,
      >   0.03795,  -0.04353,   0.30496,   0.53973,   0.44511,   1.82131,       >   0.03795d0,  -0.04353d0,  0.30496d0,   0.53973d0,   0.44511d0,   1.82131d0,
      >  14.62905,  -0.35599,    6.9752,   1.48548,  24.99437, -29.88984,       >  14.62905d0,  -0.35599d0,   6.9752d0,   1.48548d0,  24.99437d0, -29.88984d0,
      >   0.04790,   0.05632,   0.20660,   0.39304,   0.45673,   4.19462,       >   0.04790d0,   0.05632d0,  0.20660d0,   0.39304d0,   0.45673d0,   4.19462d0,
      >  47.41709, -39.83233,   16.6487,  17.18184,  -2.52659, -32.65720,       >  47.41709d0, -39.83233d0,  16.6487d0,  17.18184d0,  -2.52659d0, -32.65720d0,
      >   0.05140,   0.11201,   0.19756,   0.26992,   0.30311,   3.20459,       >   0.05140d0,   0.11201d0,  0.19756d0,   0.26992d0,   0.30311d0,   3.20459d0,
      >  37.47848, -31.60155,   12.7851,  13.97016,  -8.18266, -18.63125,       >  37.47848d0, -31.60155d0,  12.7851d0,  13.97016d0,  -8.18266d0, -18.63125d0,
      >   0.04703,   0.17203,   0.21807,   0.23841,   0.26405,   1.17976,       >   0.04703d0,   0.17203d0,  0.21807d0,   0.23841d0,   0.26405d0,   1.17976d0,
      >  10.84169,  -4.87255,    4.5850,   2.60271,  -1.08836,  -2.73421,       >  10.84169d0,  -4.87255d0,   4.5850d0,   2.60271d0,  -1.08836d0,  -2.73421d0,
      >   0.09651,  -0.07097,  -0.01140,   0.41291,   0.67944,   8.63177,       >   0.09651d0,  -0.07097d0, -0.01140d0,   0.41291d0,   0.67944d0,   8.63177d0,
      > 116.48230,-115.22170,   34.6775,  50.71781, -62.71606, -31.72991,       > 116.48230d0,-115.22170d0,  34.6775d0,  50.71781d0, -62.71606d0, -31.72991d0,
      >   0.05407,  -0.01853,   0.18484,   0.39172,   0.39883,   4.54709,       >   0.05407d0,  -0.01853d0,  0.18484d0,   0.39172d0,   0.39883d0,   4.54709d0,
      >  57.66357, -47.37784,   18.0009,  22.01861, -16.19442, -25.28908,       >  57.66357d0, -47.37784d0,  18.0009d0,  22.01861d0, -16.19442d0, -25.28908d0,
      >   0.05646,  -0.02226,   0.19051,   0.38210,   0.35937,   2.35732,       >   0.05646d0,  -0.02226d0,  0.19051d0,   0.38210d0,   0.35937d0,   2.35732d0,
      >  31.71869, -21.30330,    9.5296,  14.48209, -14.30969,  -7.90519,       >  31.71869d0, -21.30330d0,   9.5296d0,  14.48209d0, -14.30969d0,  -7.90519d0,
      >   0.04477,   0.00480,   0.15076,   0.29189,   0.31529,   3.12081,       >   0.04477d0,   0.00480d0,  0.15076d0,   0.29189d0,   0.31529d0,   3.12081d0,
      >  41.81401, -31.97399,   12.4581,  15.71344, -15.68024, -13.09181,       >  41.81401d0, -31.97399d0,  12.4581d0,  15.71344d0, -15.68024d0, -13.09181d0,
      >   0.03143,   0.06733,   0.18032,   0.27882,   0.27721,   1.09240,       >   0.03143d0,   0.06733d0,  0.18032d0,   0.27882d0,   0.27721d0,   1.09240d0,
      >  15.14939,  -7.18577,    4.4895,   7.48708, -11.64642,   1.23596,       >  15.14939d0,  -7.18577d0,   4.4895d0,   7.48708d0, -11.64642d0,   1.23596d0,
      >   0.09641,  -0.07037,  -0.04597,   0.26390,   0.48362,   3.68185,       >   0.09641d0,  -0.07037d0, -0.04597d0,   0.26390d0,   0.48362d0,   3.68185d0,
      >  61.90438, -66.11752,   15.0943,  29.88299, -76.76054,  30.73642,       >  61.90438d0, -66.11752d0,  15.0943d0,  29.88299d0, -76.76054d0,  30.73642d0,
      >   0.05348,  -0.06567,   0.26431,   0.51672,  -0.39508,  -0.17103,       >   0.05348d0,  -0.06567d0,  0.26431d0,   0.51672d0,  -0.39508d0,  -0.17103d0,
      >   0.38781,   9.77783,   -0.2076,   3.74785, -13.89592,  15.93324,       >   0.38781d0,   9.77783d0,  -0.2076d0,   3.74785d0, -13.89592d0,  15.93324d0,
      >   0.01978,   0.02381,   0.16282,   0.27720,   0.28227,   1.10671,       >   0.01978d0,   0.02381d0,  0.16282d0,   0.27720d0,   0.28227d0,   1.10671d0,
      >  19.57988,  -8.77431,    4.6923,  10.90089, -21.92500,   4.82789,       >  19.57988d0,  -8.77431d0,   4.6923d0,  10.90089d0, -21.92500d0,   4.82789d0,
      >   0.02816,   0.08035,   0.15177,   0.19596,   0.20729,   0.16046,       >   0.02816d0,   0.08035d0,  0.15177d0,   0.19596d0,   0.20729d0,   0.16046d0,
      >   4.85393,   1.60025,    0.9924,   6.34038, -23.91168,  19.81668,       >   4.85393d0,   1.60025d0,   0.9924d0,   6.34038d0, -23.91168d0,  19.81668d0,
      >   0.00810,   0.06904,   0.19159,   0.22518,   0.16841,   0.22909,       >   0.00810d0,   0.06904d0,  0.19159d0,   0.22518d0,   0.16841d0,   0.22909d0,
      >   4.46808,  -1.67027,    1.1692,   4.17188, -22.18001,  13.56817,       >   4.46808d0,  -1.67027d0,   1.1692d0,   4.17188d0, -22.18001d0,  13.56817d0,
      >  -0.00459,   0.08266,   0.20067,   0.22163,   0.15112,   0.67785,       >  -0.00459d0,   0.08266d0,  0.20067d0,   0.22163d0,   0.15112d0,   0.67785d0,
      >  10.03290, -23.97605,    2.8133,   1.49803, -32.05834,  14.49983,       >  10.03290d0, -23.97605d0,   2.8133d0,   1.49803d0, -32.05834d0,  14.49983d0,
      >  -0.00869,   0.07988,   0.20642,   0.23439,   0.16011,   1.58356,       >  -0.00869d0,   0.07988d0,  0.20642d0,   0.23439d0,   0.16011d0,   1.58356d0,
      >  23.98258, -56.34460,    6.2639,  -0.04471, -43.27986,  13.87953/       >  23.98258d0, -56.34460d0,   6.2639d0,  -0.04471d0, -43.27986d0,  13.87953d0/
         endif         endif
  
 c      if (ww .ge. 3.9) then c      if (ww .ge. 3.9) then
Line 229 
Line 229 
  
 c      if (ww.le.1.2 .or. ww.ge.3.85 c      if (ww.le.1.2 .or. ww.ge.3.85
 c     >.or. abs(ww-ww1).le.0.001) then c     >.or. abs(ww-ww1).le.0.001) then
       if (ww.le.1.2 .or. abs(ww-ww1).le.0.001) then        if (ww.le.1.2d0 .or. abs(ww-ww1).le.0.001d0) then
           r = r1           r = r1
           f2 = f21           f2 = f21
       else       else
Line 268 
Line 268 
       rlog = fac/tmp2       rlog = fac/tmp2
  
       r2 = rp(i)*rlog       r2 = rp(i)*rlog
      &+(rp(i+1)/qq+rp(i+2)/(q4+.3**2))       &+(rp(i+1)/qq+rp(i+2)/(q4+.3d0**2))
      &*(1.+rp(i+3)*xx+rp(i+4)*(xx**2))       &*(1.d0+rp(i+3)*xx+rp(i+4)*(xx**2))
      &*(xx**rp(i+5))      &*(xx**rp(i+5))
  
       q = log(qq)       q = log(qq)
       qth = .2d0+3.2d0*xx       qth = .2d0+3.2d0*xx
       dq = q-qth       dq = q-qth
       y = 1.-xx        y = 1.d0-xx
       f2th = f2p(i1)*y**1+f2p(i1+1)*y**2+f2p(i1+2)*y**3+       f2th = f2p(i1)*y**1+f2p(i1+1)*y**2+f2p(i1+2)*y**3+
      >       f2p(i1+3)*y**4+f2p(i1+4)*y**5      >       f2p(i1+3)*y**4+f2p(i1+4)*y**5
       quad = (f2p(i1+5)+f2p(i1+6)*xx+f2p(i1+7)*xx**2)*dq**2       quad = (f2p(i1+5)+f2p(i1+6)*xx+f2p(i1+7)*xx**2)*dq**2


Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3

Analyzer/Replay: Mark Jones, Documents: Stephen Wood
Powered by
ViewCVS 0.9.2-cvsgraph-1.4.0