Changes in gen_rand_conf, add prng.F - now DiL works
authorBartlomiej Zaborowski <bartek.zaborowski@chem.univ.gda.pl>
Wed, 7 Nov 2012 13:45:10 +0000 (08:45 -0500)
committerBartlomiej Zaborowski <bartek.zaborowski@chem.univ.gda.pl>
Wed, 7 Nov 2012 13:45:10 +0000 (08:45 -0500)
PARAM/torsion_631Gdp_old.parm [new file with mode: 0644]
PARAM/torsion_double_631Gdp_old.parm [new file with mode: 0644]
source/unres/src_CSA/CMakeLists.txt
source/unres/src_CSA/Makefile
source/unres/src_CSA/Makefile_4P
source/unres/src_CSA/cinfo.f [new file with mode: 0644]
source/unres/src_CSA/prng.F [new file with mode: 0644]
source/unres/src_CSA_DiL/CMakeLists.txt
source/unres/src_CSA_DiL/gen_rand_conf.F
source/unres/src_CSA_DiL/prng.F [new file with mode: 0644]
source/unres/src_CSA_DiL/together.F

diff --git a/PARAM/torsion_631Gdp_old.parm b/PARAM/torsion_631Gdp_old.parm
new file mode 100644 (file)
index 0000000..faa8083
--- /dev/null
@@ -0,0 +1,180 @@
+3  *** Parameters derived by integrating MP2/6-31G** local energy surfaces ***
+1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 2
+ 0 0 ** Gly-D-Pro-reg
+ 0 0 ** Gly-D-Ala-reg
+ 6 0 ************ G-G
+     1                    4.03053E-03       0.0
+     2                   -4.90330E-03       0.0
+     3                    2.70952E-01       0.0
+     4                   -7.20927E-02       0.0
+     5                   -3.66258E-03       0.0
+     6                   -6.34294E-03       0.0
+ 6 0 ************ G-A
+     1                    1.27465E-02      -1.01763E-01
+     2                   -8.41985E-03       1.18716E-02
+     3                    2.50582E-01       3.30235E-02
+     4                   -1.02631E-01       1.46337E-01
+     5                    6.63161E-03       8.86275E-03
+     6                   -6.81049E-03       8.71641E-03
+10 0 ************ G-P
+     1                   -1.14953E+00       1.29720E+00
+     2                   -1.08859E-03      -1.20321E-02
+     3                    4.04229E-01       2.92750E-01
+     4                   -2.41994E-01       5.38295E-02
+     5                    4.31804E-02      -6.56306E-02
+     6                    4.29080E-02       1.38685E-01
+     7                   -6.31286E-02      -2.65338E-02
+     8                    1.46743E-02      -3.26658E-03
+     9                   -1.33460E-02       3.01470E-02
+     0                   -1.11258E-02      -3.15723E-02
+ 0 0 ****** Ala-D-Pro
+ 0 0 ****** Ala-D-Ala
+ 6 0 ************ A-G
+     1                    5.14830E-03       9.48478E-02
+     2                   -8.85552E-02       1.20661E-01
+     3                    2.65903E-01      -8.99192E-02
+     4                   -6.16552E-02      -6.07952E-02
+     5                    2.19209E-02      -3.09487E-02
+     6                   -1.58136E-02      -1.74252E-02
+ 6 0 ************ A-A
+     1                    1.05084E+00      -1.83406E-03
+     2                    2.56198E-01       3.75825E-02
+     3                    2.97224E-01      -6.65886E-03
+     4                   -9.30144E-02       2.78499E-02
+     5                    2.69127E-02      -2.79813E-03
+     6                   -5.69348E-03      -5.76679E-03
+10 0 ************ A-P
+     1                    1.97348E+00       3.61809E+00
+     2                    3.25172E-01       6.93324E-01
+     3                    7.44811E-01       2.28853E-01
+     4                   -1.11515E-02       4.21766E-01
+     5                   -1.09160E-01       4.09945E-02
+     6                    1.08618E-01       6.71048E-03
+     7                    1.21785E-02       9.37813E-02
+     8                   -2.77282E-02       2.84981E-02
+     9                    2.15879E-02       1.85528E-02
+     0                   -8.56980E-03       3.92482E-02
+ 0 0 ** Pro-D-Por
+ 0 0 ** Pro-D-Ala
+10 0 ************ P-G
+     1                    2.12359E-02       4.05506E-02
+     2                   -1.19424E-01       2.91087E-01
+     3                    3.81770E-01      -1.11839E-01
+     4                   -9.75339E-02      -7.60089E-02
+     5                   -5.00214E-02       3.42037E-02
+     6                    8.12296E-02       3.98713E-02
+     7                   -1.21613E-02      -4.57682E-02
+     8                    6.95870E-03       2.28245E-02
+     9                    2.02598E-02      -2.00100E-02
+     0                   -1.48685E-02      -1.09627E-02
+10 0 ************ P-A
+     1                    1.40539E+00      -7.42879E-01
+     2                    9.21075E-02      -1.17205E-01
+     3                    4.49479E-01       3.82888E-02
+     4                    8.11975E-02      -3.35594E-02
+     5                   -9.88964E-03      -4.63759E-02
+     6                    8.58396E-02       4.22362E-02
+     7                    2.73049E-02      -5.43392E-02
+     8                   -4.18620E-03      -1.04889E-02
+     9                    1.06804E-02      -7.89066E-03
+     0                   -5.01879E-03      -1.20569E-02
+10 0 ************ P-P
+     1                    4.23219E+00       1.42318E+00
+     2                    1.24967E+00      -2.32575E-01
+     3                    5.38690E-01      -4.25507E-01
+     4                    4.01096E-01      -6.42186E-01
+     5                    2.52852E-01      -2.32257E-01
+     6                   -3.06406E-02      -2.48617E-01
+     7                   -1.68086E-02      -1.55069E-01
+     8                   -9.68421E-02      -1.80718E-01
+     9                   -2.73541E-02      -1.36915E-01
+     0                   -6.08829E-02      -6.71426E-02
+ 0 0 ** Gly-D-Pro-reg
+ 0 0 ** Gly-D-Ala-reg
+ 6 0 ************ G-G
+     1                    4.03053E-03       0.0
+     2                   -4.90330E-03       0.0
+     3                    2.70952E-01       0.0
+     4                   -7.20927E-02       0.0
+     5                   -3.66258E-03       0.0
+     6                   -6.34294E-03       0.0
+ 6 0 ************ G-A
+     1                    1.27465E-02      -1.01763E-01
+     2                   -8.41985E-03       1.18716E-02
+     3                    2.50582E-01       3.30235E-02
+     4                   -1.02631E-01       1.46337E-01
+     5                    6.63161E-03       8.86275E-03
+     6                   -6.81049E-03       8.71641E-03
+10 0 ************ G-P
+     1                   -1.14953E+00       1.29720E+00
+     2                   -1.08859E-03      -1.20321E-02
+     3                    4.04229E-01       2.92750E-01
+     4                   -2.41994E-01       5.38295E-02
+     5                    4.31804E-02      -6.56306E-02
+     6                    4.29080E-02       1.38685E-01
+     7                   -6.31286E-02      -2.65338E-02
+     8                    1.46743E-02      -3.26658E-03
+     9                   -1.33460E-02       3.01470E-02
+     0                   -1.11258E-02      -3.15723E-02
+ 0 0 ****** Ala-D-Pro
+ 0 0 ****** Ala-D-Ala
+ 6 0 ************ A-G
+     1                    5.14830E-03       9.48478E-02
+     2                   -8.85552E-02       1.20661E-01
+     3                    2.65903E-01      -8.99192E-02
+     4                   -6.16552E-02      -6.07952E-02
+     5                    2.19209E-02      -3.09487E-02
+     6                   -1.58136E-02      -1.74252E-02
+ 6 0 ************ A-A
+     1                    1.05084E+00      -1.83406E-03
+     2                    2.56198E-01       3.75825E-02
+     3                    2.97224E-01      -6.65886E-03
+     4                   -9.30144E-02       2.78499E-02
+     5                    2.69127E-02      -2.79813E-03
+     6                   -5.69348E-03      -5.76679E-03
+10 0 ************ A-P
+     1                    1.97348E+00       3.61809E+00
+     2                    3.25172E-01       6.93324E-01
+     3                    7.44811E-01       2.28853E-01
+     4                   -1.11515E-02       4.21766E-01
+     5                   -1.09160E-01       4.09945E-02
+     6                    1.08618E-01       6.71048E-03
+     7                    1.21785E-02       9.37813E-02
+     8                   -2.77282E-02       2.84981E-02
+     9                    2.15879E-02       1.85528E-02
+     0                   -8.56980E-03       3.92482E-02
+ 0 0 ** Pro-D-Por
+ 0 0 ** Pro-D-Ala
+10 0 ************ P-G
+     1                    2.12359E-02       4.05506E-02
+     2                   -1.19424E-01       2.91087E-01
+     3                    3.81770E-01      -1.11839E-01
+     4                   -9.75339E-02      -7.60089E-02
+     5                   -5.00214E-02       3.42037E-02
+     6                    8.12296E-02       3.98713E-02
+     7                   -1.21613E-02      -4.57682E-02
+     8                    6.95870E-03       2.28245E-02
+     9                    2.02598E-02      -2.00100E-02
+     0                   -1.48685E-02      -1.09627E-02
+10 0 ************ P-A
+     1                    1.40539E+00      -7.42879E-01
+     2                    9.21075E-02      -1.17205E-01
+     3                    4.49479E-01       3.82888E-02
+     4                    8.11975E-02      -3.35594E-02
+     5                   -9.88964E-03      -4.63759E-02
+     6                    8.58396E-02       4.22362E-02
+     7                    2.73049E-02      -5.43392E-02
+     8                   -4.18620E-03      -1.04889E-02
+     9                    1.06804E-02      -7.89066E-03
+     0                   -5.01879E-03      -1.20569E-02
+10 0 ************ P-P
+     1                    4.23219E+00       1.42318E+00
+     2                    1.24967E+00      -2.32575E-01
+     3                    5.38690E-01      -4.25507E-01
+     4                    4.01096E-01      -6.42186E-01
+     5                    2.52852E-01      -2.32257E-01
+     6                   -3.06406E-02      -2.48617E-01
+     7                   -1.68086E-02      -1.55069E-01
+     8                   -9.68421E-02      -1.80718E-01
+     9                   -2.73541E-02      -1.36915E-01
+     0                   -6.08829E-02      -6.71426E-02
diff --git a/PARAM/torsion_double_631Gdp_old.parm b/PARAM/torsion_double_631Gdp_old.parm
new file mode 100644 (file)
index 0000000..9eb0747
--- /dev/null
@@ -0,0 +1,5317 @@
+Gpp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Gpa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GpG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GpA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GpP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Gap
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Gaa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GaG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GaA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GaP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GGp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GGa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GGG_abinitio.fitout
+6 6
+-4.46036E-03
+-5.52737E-05
+-2.18852E-02
+4.33960E-03
+1.85025E-03
+4.95855E-04
+-3.24597E-18
+-2.02014E-17
+4.26966E-17
+1.49283E-17
+2.82040E-17
+1.43925E-16
+-7.18411E-03
+7.92072E-03
+-2.17806E-02
+1.67057E-03
+3.84725E-03
+-7.17231E-05
+-7.20379E-18
+-4.20384E-17
+3.03891E-17
+-2.44087E-18
+9.05009E-17
+1.14122E-19
+-7.23813E-03
+3.28293E-03
+-1.39801E-17
+-1.89800E-17
+-1.30490E-03
+1.38772E-02
+-7.10794E-18
+3.31613E-17
+1.99808E-04
+-2.17604E-03
+3.79714E-18
+5.99313E-17
+1.26587E-02
+-2.05259E-02
+1.93061E-17
+2.66731E-17
+-3.66073E-04
+3.55213E-03
+-1.11906E-17
+-3.51248E-18
+-4.59544E-03
+-9.54163E-03
+-1.52654E-17
+3.48214E-18
+4.86862E-03
+-9.46975E-05
+-1.03323E-17
+-5.01265E-17
+-5.04727E-03
+2.44410E-03
+5.33373E-17
+-1.92922E-17
+2.77079E-03
+2.66047E-02
+2.43917E-17
+-3.41779E-17
+3.19743E-03
+-5.00879E-03
+-1.09699E-17
+3.58392E-17
+4.68052E-03
+6.60091E-04
+3.91014E-17
+3.76848E-17
+-3.19326E-04
+2.01276E-04
+-2.20174E-17
+4.87189E-17
+-2.79146E-02
+-6.79853E-02
+1.69022E-17
+-5.18855E-17
+2.15742E-02
+-4.87803E-02
+2.70102E-17
+-5.03522E-18
+1.86440E-04
+-1.17222E-04
+4.25151E-18
+-1.13210E-17
+GGA_abinitio.fitout
+6 6
+-8.46686E-03
+2.62291E-03
+-3.30608E-02
+7.96177E-03
+7.15867E-03
+3.66237E-03
+-7.72008E-04
+-2.02845E-04
+1.62280E-03
+8.54876E-04
+-8.49190E-05
+-5.78018E-04
+-8.84141E-03
+1.13979E-02
+-2.03566E-02
+3.55919E-03
+-3.02160E-04
+-1.26198E-03
+3.52817E-02
+8.28932E-03
+4.25960E-03
+-6.41946E-03
+3.60893E-04
+3.94028E-04
+-1.35465E-02
+4.17716E-03
+8.28258E-02
+1.61355E-02
+8.31401E-03
+1.44081E-02
+1.30831E-02
+-1.86039E-04
+4.43630E-04
+-2.39715E-03
+-3.58254E-03
+-6.00048E-03
+9.14174E-03
+-8.39490E-03
+-6.77698E-04
+1.68410E-03
+-1.67024E-03
+-2.02726E-03
+5.39827E-04
+-8.64739E-04
+-7.91196E-03
+-1.68047E-02
+6.92503E-02
+-8.67821E-02
+3.44170E-03
+8.48263E-04
+3.84754E-03
+2.18871E-03
+-4.09461E-03
+-4.34824E-04
+9.37754E-04
+-2.35904E-03
+1.07912E-02
+3.20561E-02
+-9.95653E-04
+-1.99915E-02
+9.98327E-03
+-8.10801E-03
+-7.96601E-02
+-5.92419E-02
+-4.49290E-03
+1.00540E-03
+-1.71946E-03
+1.14316E-03
+7.92241E-05
+6.05998E-04
+3.78783E-04
+2.15709E-03
+-1.93128E-02
+-5.56518E-02
+-1.58545E-03
+8.20476E-03
+1.36469E-02
+-5.13838E-02
+2.85169E-03
+3.67851E-02
+1.00136E-03
+-2.29434E-03
+-4.83317E-03
+-3.92162E-03
+GGP_abinitio.fitout
+6 6
+6.50508E-02
+-5.89761E-03
+-8.27980E-02
+3.10265E-02
+-3.38962E-03
+6.73743E-04
+-4.16459E-04
+1.00721E-04
+-4.62636E-05
+-3.49098E-03
+-9.41654E-05
+2.71000E-04
+-2.72775E-02
+-8.83340E-04
+-4.37583E-03
+8.50974E-04
+-2.05161E-03
+-2.06744E-03
+3.15235E-02
+-2.21357E-03
+-3.64561E-03
+5.73169E-05
+2.49261E-03
+-7.51165E-03
+-5.84943E-02
+3.97047E-02
+6.60434E-02
+4.52294E-02
+3.58092E-03
+4.17625E-03
+3.31956E-02
+-2.90994E-02
+-4.75490E-04
+-3.86681E-03
+9.36791E-04
+-4.89921E-03
+1.88598E-02
+-4.42295E-03
+1.27901E-02
+2.50468E-03
+-1.10166E-03
+-6.43327E-04
+-7.04067E-03
+6.64965E-04
+1.12632E-02
+1.15798E-01
+-6.89545E-03
+1.25577E-01
+-1.23962E-02
+-6.95734E-03
+3.95449E-03
+-1.02346E-03
+-3.67717E-03
+-3.65381E-03
+-2.61415E-03
+3.03406E-03
+-3.43625E-03
+1.68217E-02
+1.47508E-02
+-1.73684E-01
+-1.75592E-02
+-2.20142E-02
+2.19531E-02
+-1.87179E-02
+-3.30625E-04
+-2.47359E-03
+-7.45116E-04
+-5.56124E-03
+2.54898E-03
+6.05867E-03
+-6.14297E-04
+1.24640E-03
+-1.17324E-02
+-1.19361E-02
+-1.28656E-02
+8.59711E-03
+-5.18731E-03
+7.58531E-03
+-5.10999E-02
+-4.52666E-02
+3.56032E-03
+3.85053E-03
+-3.65556E-03
+4.00514E-03
+GAp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GAa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GAG_abinitio.fitout
+6 6
+-6.97988E-03
+2.11326E-04
+-3.29929E-02
+2.08293E-03
+-4.84598E-04
+2.83986E-03
+9.70254E-04
+-7.61772E-04
+8.32483E-03
+-1.25922E-02
+-3.45241E-04
+-1.52497E-03
+-5.32058E-03
+5.03716E-03
+-2.67407E-02
+-2.11211E-05
+-1.71645E-03
+1.11761E-03
+-2.54093E-03
+-1.64341E-02
+-4.08188E-03
+1.68737E-03
+3.46628E-03
+1.20169E-03
+-1.55672E-03
+-3.27573E-04
+1.18295E-03
+-2.74125E-03
+6.89273E-03
+5.44083E-03
+-1.06816E-02
+8.72845E-03
+-1.39272E-03
+-7.16988E-04
+-5.68857E-04
+4.75969E-04
+-7.24911E-03
+-1.12901E-02
+-1.80540E-02
+4.99824E-05
+8.68770E-04
+-6.51696E-05
+-2.98370E-03
+2.10040E-03
+-2.04330E-04
+-7.57116E-03
+2.10429E-02
+-7.44830E-03
+-9.90521E-04
+1.59562E-03
+2.77189E-03
+-4.99094E-03
+-3.89075E-03
+-2.77476E-03
+-3.39128E-04
+-1.78779E-03
+-1.03301E-02
+3.30321E-02
+-4.51590E-02
+4.76619E-02
+-8.44719E-03
+4.88342E-03
+-1.77256E-03
+-4.07904E-03
+1.13761E-03
+2.28837E-03
+-1.92714E-03
+1.78535E-03
+-4.62403E-05
+6.19471E-04
+-1.29644E-03
+1.96271E-03
+-4.78698E-02
+-8.10150E-02
+-1.16577E-02
+2.11953E-02
+1.00520E-02
+-2.99724E-02
+-4.03191E-03
+2.14660E-03
+1.38202E-04
+-5.32398E-04
+-1.11994E-04
+-1.14554E-03
+GAA_abinitio.fitout
+6 6
+-2.48643E-02
+4.81526E-03
+-4.51777E-02
+2.76581E-02
+-3.04518E-04
+4.36235E-03
+2.17373E-02
+-2.41799E-03
+5.37657E-03
+-2.98047E-02
+-2.52519E-03
+-5.47296E-03
+-1.58140E-02
+-6.16724E-03
+-2.25130E-02
+3.11349E-03
+-3.17153E-03
+-1.31443E-03
+1.86030E-02
+-6.75956E-03
+-7.49626E-03
+-4.44882E-03
+-1.60960E-03
+9.74007E-04
+-1.32605E-02
+-1.96368E-02
+4.43606E-02
+-1.22538E-02
+1.29574E-02
+9.39390E-04
+1.40595E-02
+-4.89717E-04
+-1.80880E-03
+2.35373E-03
+4.11257E-03
+-2.89455E-03
+5.20837E-03
+-7.82527E-03
+-8.10209E-03
+2.09822E-03
+-5.59362E-04
+6.21779E-04
+1.19200E-03
+3.27965E-03
+6.34701E-02
+-9.55843E-02
+7.89666E-02
+-8.69198E-02
+3.85345E-03
+-1.60660E-03
+6.48502E-03
+-2.57321E-04
+-2.77198E-03
+-2.22505E-03
+8.32341E-04
+4.61294E-04
+4.76449E-02
+7.77390E-03
+-4.89485E-02
+9.50176E-04
+-5.87483E-03
+4.95926E-02
+-1.88864E-02
+-2.20620E-02
+3.61988E-03
+3.72923E-04
+-9.01274E-04
+2.99111E-03
+3.87878E-04
+1.49061E-03
+-7.52049E-04
+2.88234E-03
+-2.78368E-02
+-5.76616E-02
+-4.44516E-02
+3.52391E-02
+1.62484E-02
+-1.33273E-02
+-2.03980E-02
+5.05344E-03
+-3.74544E-03
+1.78525E-03
+-2.40457E-03
+-3.94525E-03
+GAP_abinitio.fitout
+6 6
+2.16431E-02
+1.01017E-02
+1.08876E-01
+7.34409E-02
+-1.09853E-02
+9.79247E-04
+4.98205E-02
+-6.29565E-03
+-1.37799E-01
+-4.10759E-02
+-2.46017E-03
+4.80559E-04
+1.58473E-03
+1.75081E-02
+-1.54409E-03
+2.42574E-03
+4.26959E-03
+-1.02854E-03
+1.03385E-02
+-1.57001E-02
+2.02765E-03
+1.87894E-04
+2.12208E-04
+3.96201E-03
+-2.27239E-02
+7.46691E-03
+4.17848E-02
+2.41565E-02
+7.74483E-03
+1.63369E-03
+3.93859E-02
+-2.02702E-02
+-2.30982E-03
+-1.11121E-03
+3.03768E-03
+-1.78059E-03
+-1.92422E-03
+1.47293E-03
+1.85671E-02
+-7.15950E-03
+1.11092E-03
+4.77112E-03
+-3.82005E-03
+6.27417E-04
+1.72187E-01
+8.48193E-02
+5.13335E-02
+2.56405E-02
+-1.86126E-02
+-5.59752E-03
+1.12581E-02
+-9.64361E-04
+6.46262E-04
+-2.31221E-03
+-2.40846E-03
+1.63797E-03
+4.54519E-02
+-5.91857E-02
+-1.00474E-02
+-1.40591E-01
+1.81836E-02
+-2.04773E-02
+6.49684E-02
+3.16978E-02
+-9.31686E-03
+7.77812E-03
+3.71733E-03
+-8.83382E-03
+1.26015E-03
+7.59054E-04
+4.31544E-04
+-8.23154E-04
+1.76513E-02
+-3.27815E-03
+-1.17345E-01
+5.62252E-02
+-1.81312E-02
+4.04796E-02
+-5.55516E-02
+-1.14203E-02
+9.65567E-04
+1.29682E-03
+-5.17527E-03
+1.59823E-03
+GPp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GPa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GPG_abinitio.fitout
+6 6
+-7.59230E-03
+-1.94552E-04
+-5.12944E-03
+3.20316E-03
+8.28212E-04
+-6.12307E-03
+-4.85436E-03
+-3.86687E-03
+5.36871E-03
+-3.81511E-03
+1.48924E-03
+-4.03381E-03
+-1.52037E-03
+5.72898E-03
+-2.07652E-02
+6.52919E-03
+-1.05532E-03
+-1.52037E-03
+-4.69660E-03
+-6.37003E-03
+6.17507E-03
+8.39488E-03
+-1.76904E-03
+8.05108E-03
+-7.98866E-04
+3.32411E-03
+5.85637E-03
+-2.98884E-03
+-1.20298E-02
+-3.13323E-02
+-1.89038E-02
+-1.18356E-02
+-3.02395E-03
+2.43800E-03
+1.38867E-03
+3.67447E-04
+-2.66576E-02
+5.35667E-02
+3.22822E-02
+-4.20099E-02
+-2.76279E-03
+2.73121E-03
+-2.92878E-03
+-3.97641E-03
+5.10833E-03
+-2.13602E-03
+1.80687E-03
+-2.80887E-03
+1.75818E-02
+4.47119E-03
+1.94978E-03
+3.01372E-02
+2.73577E-03
+2.07557E-04
+1.02968E-02
+1.06706E-02
+-1.07292E-02
+-4.42466E-03
+1.03471E-02
+3.26791E-02
+-3.88377E-03
+1.21678E-03
+-7.07759E-04
+1.76439E-03
+-9.08375E-03
+1.74263E-03
+2.04763E-03
+-1.17705E-02
+2.75141E-03
+-6.02603E-03
+-3.58871E-03
+-2.79755E-03
+2.32403E-02
+-5.25103E-02
+1.58724E-02
+-4.78943E-02
+-6.49652E-04
+-1.30365E-02
+-1.24654E-02
+-1.82200E-02
+-1.17657E-04
+-1.48990E-03
+-1.15237E-03
+-1.71433E-04
+GPA_abinitio.fitout
+6 6
+1.76833E-02
+9.87411E-04
+-1.16391E-02
+1.75059E-02
+-6.24857E-03
+-5.22125E-03
+5.86116E-03
+-2.32421E-03
+-1.58225E-02
+7.43433E-03
+4.93394E-03
+-1.89631E-02
+-3.25799E-02
+-2.32415E-02
+-3.10454E-02
+-8.91640E-03
+-7.65479E-03
+-1.27982E-03
+-9.57172E-03
+8.77571E-03
+6.27899E-03
+9.87208E-03
+1.16529E-02
+8.66130E-03
+1.92294E-02
+7.29907E-02
+5.17669E-02
+-2.27035E-02
+6.60383E-03
+2.08489E-02
+-4.78442E-03
+-1.14125E-02
+-4.55136E-03
+7.67130E-03
+1.21326E-02
+1.30754E-02
+-2.73899E-02
+6.76535E-02
+3.83851E-02
+-5.38365E-02
+5.65271E-03
+-5.37810E-04
+6.97088E-03
+7.72153E-03
+2.74100E-02
+-4.70346E-02
+-2.35609E-02
+-5.97889E-02
+3.49906E-02
+3.87800E-02
+2.36671E-02
+3.21217E-02
+3.57922E-03
+3.74125E-03
+1.26977E-02
+1.32543E-02
+-1.83147E-02
+-2.62166E-02
+-5.76715E-03
+-1.83617E-02
+-6.49200E-03
+4.85116E-02
+2.89452E-02
+1.35099E-02
+8.75516E-03
+2.02179E-02
+4.64394E-03
+7.27739E-03
+4.54111E-03
+-3.11138E-03
+2.01092E-03
+2.90504E-03
+2.15560E-02
+-7.36577E-02
+1.91254E-02
+-5.61390E-02
+8.99074E-03
+1.35028E-02
+4.72925E-04
+-3.22296E-03
+-9.02731E-03
+-1.73413E-02
+-9.12097E-03
+1.71334E-02
+GPP_abinitio.fitout
+6 6
+3.55782E-01
+-2.50590E-02
+-2.02781E-02
+1.88976E-01
+-6.17588E-02
+-2.91385E-02
+2.05426E-01
+-5.87770E-03
+-3.30234E-01
+7.61622E-02
+5.16617E-02
+-1.35195E-01
+-3.32316E-03
+4.73134E-03
+9.83600E-03
+1.38053E-02
+2.51654E-03
+6.42227E-03
+-1.41234E-02
+1.49387E-02
+-1.78714E-04
+-5.87439E-03
+1.08191E-03
+-1.97401E-04
+1.05038E-01
+2.91828E-01
+3.05333E-01
+-1.00742E-01
+-1.01012E-02
+-1.49328E-02
+3.57843E-02
+-4.10710E-02
+-2.83700E-02
+-2.30223E-03
+-2.43599E-03
+2.05701E-02
+-1.78170E-02
+-4.82013E-03
+6.73937E-03
+-3.38595E-03
+3.28851E-03
+-8.08060E-03
+-7.65673E-03
+5.09805E-03
+1.57284E-01
+-2.00083E-01
+-2.14765E-01
+-1.45042E-01
+3.67444E-02
+1.99603E-02
+2.72203E-02
+3.10215E-02
+-2.08581E-03
+-1.66514E-03
+-4.79677E-03
+-3.68899E-03
+1.95962E-02
+-5.81681E-03
+-5.65331E-02
+5.75209E-02
+6.98914E-02
+1.28428E-01
+1.38299E-01
+-6.32083E-02
+-1.13264E-02
+-1.45981E-02
+-6.42446E-03
+-1.94911E-02
+-8.84309E-03
+-2.46011E-03
+-2.28773E-03
+-5.77668E-03
+3.89634E-03
+9.03634E-03
+9.92089E-03
+-6.47315E-03
+4.81635E-02
+-5.01397E-02
+3.40332E-02
+-2.73013E-02
+-6.50338E-02
+9.57053E-04
+3.41943E-03
+5.80621E-02
+App
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Apa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+ApG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+ApA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+ApP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Aap
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Aaa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AaG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AaA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AaP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AGp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AGa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AGG_abinitio.fitout
+6 6
+-7.39715E-03
+3.65556E-03
+-1.91078E-02
+2.38136E-03
+-5.37961E-03
+3.61998E-03
+-2.18464E-02
+-4.41896E-03
+5.08921E-03
+3.21624E-03
+5.89060E-03
+1.19802E-03
+-2.19281E-02
+1.49645E-02
+-3.56136E-02
+9.88455E-03
+1.22315E-02
+-3.00099E-03
+2.21176E-03
+3.15110E-03
+-6.61995E-03
+-3.20593E-03
+1.54183E-03
+-8.59009E-04
+-5.99348E-03
+4.65389E-05
+-5.95310E-02
+1.06111E-02
+1.57623E-02
+1.59041E-02
+1.00177E-04
+1.29357E-01
+4.43290E-03
+2.08145E-02
+-4.15715E-03
+1.26761E-03
+1.04441E-02
+-2.89679E-02
+6.04392E-02
+-1.59174E-01
+-5.14944E-02
+1.87340E-02
+3.85570E-02
+-3.73199E-02
+7.51987E-04
+-4.77918E-03
+4.26491E-03
+4.13725E-03
+9.71818E-03
+9.73285E-04
+3.70014E-02
+-2.94368E-04
+-4.04141E-03
+3.29676E-03
+4.40040E-02
+-7.37677E-03
+-2.53272E-03
+1.62907E-02
+1.14612E-03
+-1.48734E-02
+-3.33955E-03
+-2.38826E-03
+-8.32102E-03
+1.27157E-03
+4.85104E-03
+4.73026E-04
+2.63678E-02
+-1.07751E-02
+-2.39741E-03
+-1.84784E-02
+5.03233E-03
+1.06810E-02
+-2.90718E-02
+-5.45636E-02
+4.72574E-03
+3.35612E-02
+1.56482E-02
+-2.86180E-02
+7.62139E-03
+-1.87081E-02
+-1.26071E-03
+2.55160E-03
+3.08408E-03
+6.64625E-04
+AGA_abinitio.fitout
+6 6
+-2.10939E-02
+6.03383E-02
+-5.52248E-02
+1.30827E-02
+-1.40094E-02
+1.47428E-02
+-4.01256E-02
+-2.37278E-02
+-2.87652E-02
+3.68246E-02
+2.93335E-02
+4.16440E-04
+1.17436E-02
+1.16773E-01
+-5.62224E-02
+2.90645E-02
+9.23893E-03
+9.04317E-04
+8.25772E-02
+-8.68232E-03
+5.76473E-02
+-5.38481E-02
+4.24119E-03
+-2.74517E-03
+-7.20633E-01
+-1.75249E-01
+2.57130E-03
+1.87130E-02
+-3.59686E-02
+6.77989E-02
+2.08938E-02
+1.32158E-01
+-1.08429E-02
+7.44245E-02
+-1.44520E-03
+1.01815E-01
+-3.94152E-02
+-3.71396E-02
+5.52084E-02
+-1.15046E-01
+1.03393E-01
+3.15823E-02
+4.11136E-03
+2.05022E-02
+3.50590E-02
+-3.83513E-02
+2.49373E-02
+-8.93949E-02
+1.20974E-02
+-4.74375E-03
+1.72335E-02
+1.47956E-02
+-7.11183E-03
+-1.84173E-03
+2.06808E-02
+-7.58410E-03
+1.52949E-02
+-5.81264E-03
+-5.96499E-04
+-2.82883E-02
+-1.02857E-02
+2.25503E-02
+-2.57388E-02
+-2.87509E-02
+-8.97673E-03
+-8.26320E-03
+-5.34274E-03
+-1.16850E-02
+2.07672E-02
+-3.82975E-03
+-1.01681E-02
+6.05167E-03
+-5.72931E-02
+-4.97570E-02
+1.03005E-03
+4.70162E-02
+-1.06617E-03
+-2.36271E-02
+1.88343E-03
+-3.47428E-03
+1.31956E-03
+-2.15396E-02
+9.27436E-03
+-7.35804E-03
+AGP_abinitio.fitout
+6 6
+1.31640E-01
+2.58075E-01
+-8.58162E-02
+1.38005E-02
+-1.71596E-03
+5.98072E-03
+3.33712E-01
+2.93859E-02
+-3.00524E-02
+7.07210E-03
+7.75509E-03
+-8.65250E-03
+-8.74605E-02
+4.48412E-03
+-4.66310E-02
+3.06797E-02
+-2.19807E-02
+-1.58779E-02
+1.05553E-01
+6.43617E-02
+-4.08730E-02
+-2.03012E-03
+3.47495E-02
+-7.22030E-02
+-7.17920E-01
+-5.04379E-01
+-5.06488E-01
+5.31765E-01
+-1.50917E-01
+4.25707E-01
+5.91572E-02
+7.35844E-03
+-6.86600E-02
+8.34443E-03
+2.77907E-02
+2.06332E-02
+-1.13316E-01
+-6.99060E-02
+2.08245E-01
+-8.67383E-02
+-2.09219E-02
+2.28999E-02
+1.01289E-01
+-7.87242E-02
+-9.80531E-04
+1.26618E-01
+-4.27857E-02
+1.02030E-01
+-5.09267E-02
+-7.30891E-03
+-1.41543E-01
+-1.16047E-02
+6.00033E-02
+-7.06215E-02
+1.06774E-01
+4.95735E-02
+6.20532E-03
+-1.58945E-02
+1.13562E-02
+-1.70360E-01
+-2.32910E-02
+-6.32569E-03
+1.29374E-02
+-2.33028E-02
+7.41635E-02
+6.90855E-02
+3.93931E-02
+-4.32746E-02
+6.13297E-02
+3.73997E-02
+-4.94826E-03
+-3.71951E-03
+-2.28195E-03
+-6.09570E-03
+-8.25585E-03
+2.80583E-02
+5.84276E-03
+-3.18833E-02
+-3.95278E-02
+-3.23312E-02
+5.37754E-03
+3.70124E-03
+-5.30540E-05
+2.20955E-02
+AAp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AAa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AAG_abinitio.fitout
+6 6
+-1.70504E-02
+-1.25725E-02
+-2.16572E-02
+9.17676E-03
+-7.82222E-03
+4.46235E-04
+-2.89066E-02
+-3.18047E-03
+1.97174E-03
+8.13914E-05
+2.44927E-03
+-7.40994E-04
+-1.96244E-02
+4.33134E-02
+-3.37035E-02
+7.64787E-03
+-2.13924E-02
+9.52598E-03
+-3.70523E-02
+-7.40036E-02
+2.40140E-02
+-5.10030E-03
+1.59236E-02
+7.43181E-03
+-1.31910E-02
+-7.98080E-03
+-3.09633E-02
+-7.51476E-03
+1.27383E-01
+-5.50401E-02
+5.28821E-03
+1.30486E-01
+7.11103E-03
+3.19300E-03
+-6.44252E-03
+-1.02063E-02
+6.87598E-02
+-8.99360E-02
+-6.87923E-02
+-1.23161E-01
+2.60656E-03
+-1.74383E-02
+3.18957E-02
+1.71278E-02
+6.27281E-03
+-3.34598E-03
+9.42044E-03
+-7.47981E-03
+-8.29916E-03
+2.89364E-02
+8.84389E-03
+-1.37194E-03
+3.56810E-02
+9.28219E-03
+3.85979E-02
+-2.62633E-02
+-1.93497E-02
+2.38321E-02
+-2.48270E-02
+4.03903E-02
+-3.67611E-04
+-4.36283E-03
+-9.47872E-03
+3.38944E-04
+3.18600E-03
+-4.88922E-03
+5.87363E-03
+-1.14359E-03
+-1.15826E-03
+-3.43521E-03
+7.57689E-03
+3.94852E-03
+-3.77955E-02
+-5.78792E-02
+1.21638E-02
+2.92502E-02
+1.65337E-02
+-2.18462E-02
+-9.40339E-03
+5.38053E-03
+3.03565E-03
+-1.07824E-04
+4.67967E-03
+-2.39506E-03
+AAA_abinitio.fitout
+6 6
+-2.05481E-01
+-8.26191E-02
+-6.95965E-02
+2.35665E-02
+-2.42723E-02
+1.32714E-03
+-8.89860E-02
+3.10180E-02
+-1.78645E-02
+-1.17453E-02
+6.66812E-03
+2.24833E-03
+-2.12827E-01
+-8.67822E-02
+-6.42470E-02
+2.51073E-02
+-1.87957E-02
+-2.67943E-03
+2.43460E-02
+-6.09457E-02
+3.04366E-02
+1.66922E-04
+-4.91981E-04
+4.11041E-03
+-5.57253E-01
+3.53999E-02
+-3.54366E-02
+-2.08256E-01
+-1.00701E-01
+5.79432E-03
+8.35129E-02
+7.78467E-02
+-1.23620E-01
+5.25424E-02
+-7.34795E-02
+3.79700E-02
+3.56303E-02
+-1.04002E-01
+1.77074E-02
+-6.40671E-02
+1.61121E-02
+5.81681E-02
+2.14007E-02
+6.03380E-02
+3.14908E-02
+-8.44041E-02
+-1.35241E-02
+-5.63823E-02
+-1.98125E-02
+-1.07643E-02
+3.85091E-02
+2.61695E-03
+2.67479E-02
+2.91781E-02
+3.32585E-02
+3.74364E-03
+3.95839E-02
+3.44728E-02
+-3.61017E-02
+2.51821E-02
+-1.16023E-02
+9.06885E-03
+-3.49017E-02
+-2.81423E-02
+-7.23186E-03
+-2.86460E-02
+1.58622E-02
+-1.63989E-03
+-1.90603E-02
+-6.01905E-03
+1.24698E-02
+2.47730E-03
+-3.53859E-02
+-1.77117E-02
+-2.28246E-02
+4.82700E-02
+1.54564E-02
+-5.67317E-03
+-7.01419E-03
+4.19916E-03
+7.22849E-03
+-3.30823E-02
+-8.92923E-03
+-1.67293E-02
+AAP_abinitio.fitout
+6 6
+-2.60341E-01
+-1.76464E-01
+-2.58402E-02
+5.08365E-02
+9.37785E-03
+1.34683E-02
+9.21716E-02
+4.11687E-02
+-8.52362E-02
+-1.16991E-02
+-2.83520E-02
+-2.85315E-02
+-1.15678E-01
+1.39526E-01
+1.43150E-03
+-1.24315E-02
+3.67833E-02
+-1.29502E-02
+-1.49832E-01
+-1.70965E-01
+-5.62158E-03
+-2.47031E-02
+-2.08181E-02
+1.00187E-02
+-5.09741E-01
+-3.04707E-01
+-1.77654E-01
+6.17878E-02
+-3.70183E-01
+3.16793E-01
+2.12524E-02
+2.52440E-02
+-6.09105E-02
+-8.32189E-02
+-3.45513E-02
+9.52014E-03
+-1.82716E-01
+1.57791E-02
+5.23627E-02
+-4.44659E-03
+-6.04463E-02
+7.35341E-02
+6.70677E-02
+-5.01091E-02
+1.58910E-01
+7.43892E-02
+-4.84020E-03
+2.64253E-02
+-1.74718E-01
+2.79791E-02
+-1.48718E-01
+-1.55775E-02
+2.66950E-02
+-1.89974E-02
+9.19990E-02
+-1.34598E-02
+2.58915E-02
+-8.32852E-02
+-1.21375E-02
+-1.35189E-01
+-2.34513E-03
+2.60660E-02
+2.74800E-02
+2.42119E-02
+-1.05766E-02
+1.07458E-01
+-9.49010E-02
+4.41061E-02
+2.37077E-04
+1.65542E-02
+5.17467E-03
+-2.63125E-02
+4.40376E-03
+-5.38761E-03
+-9.77473E-02
+6.71620E-02
+2.48552E-03
+-2.73268E-03
+-6.93468E-02
+-2.85697E-02
+-4.28702E-04
+1.93884E-02
+3.32264E-03
+6.53017E-03
+APp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+APa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+APG_abinitio.fitout
+6 6
+2.11374E-02
+2.96387E-02
+7.90197E-03
+2.06850E-02
+2.34088E-02
+6.25406E-03
+-6.08740E-03
+-1.55375E-03
+-5.36785E-03
+-3.11711E-02
+-9.42745E-03
+-9.88549E-03
+-2.33313E-03
+1.05580E-02
+-6.05641E-02
+2.53598E-02
+-3.71104E-03
+-5.01697E-03
+-2.89535E-02
+-2.77652E-02
+3.61965E-02
+2.00866E-02
+-3.16730E-03
+1.42707E-02
+-2.45496E-03
+-6.39974E-03
+-5.31496E-03
+1.49525E-02
+-5.33448E-02
+-6.12883E-02
+1.64008E-02
+5.17439E-02
+-1.32435E-02
+8.86931E-04
+-1.93406E-02
+1.29406E-02
+3.93005E-02
+-4.91735E-02
+8.90276E-02
+-1.29607E-01
+-4.02511E-02
+2.03071E-03
+-3.23640E-03
+1.57962E-02
+2.20466E-03
+4.80277E-03
+-7.51544E-03
+1.37773E-02
+2.79080E-02
+6.38630E-02
+-3.29278E-02
+6.21700E-03
+2.75330E-02
+-6.78977E-02
+7.03002E-02
+-2.20345E-02
+-1.90882E-02
+1.59821E-03
+1.08734E-02
+6.23326E-02
+-1.95481E-02
+6.49369E-03
+1.63273E-03
+1.99442E-02
+1.58751E-03
+-2.23369E-02
+2.25455E-02
+-1.50345E-02
+2.31505E-02
+2.36684E-02
+-2.15136E-02
+-2.13738E-02
+4.16797E-02
+-1.06406E-01
+4.31283E-02
+-5.98608E-02
+-2.21497E-02
+-4.07299E-02
+-4.56594E-04
+-1.37761E-03
+-1.48161E-02
+3.23533E-03
+-1.19902E-02
+1.20310E-02
+APA_abinitio.fitout
+6 6
+8.58406E-02
+7.11000E-02
+-2.46047E-03
+6.50547E-02
+6.13304E-02
+4.29823E-03
+-6.05849E-02
+-2.44367E-03
+-3.77551E-02
+-7.23133E-02
+-3.86123E-03
+-1.07154E-02
+-1.03709E-01
+-7.05574E-02
+-9.49350E-02
+-1.53363E-02
+-1.18499E-02
+-7.61708E-03
+-2.91869E-02
+4.65680E-02
+3.62772E-02
+3.39042E-02
+3.16662E-02
+1.47463E-02
+1.52687E-01
+3.95754E-02
+-1.32083E-02
+-1.89615E-01
+-1.83152E-02
+-3.72732E-03
+-1.41424E-02
+-4.35396E-02
+3.71739E-02
+2.14553E-02
+6.33433E-02
+1.58020E-02
+4.23314E-02
+-7.04643E-02
+1.00041E-01
+-1.66457E-01
+1.68853E-02
+-4.71492E-02
+4.53786E-02
+1.36936E-02
+1.65425E-02
+-1.09111E-01
+-4.96983E-02
+-7.14614E-02
+7.25503E-02
+1.07859E-01
+-6.38044E-02
+-5.08115E-02
+2.52785E-02
+-7.96546E-02
+8.77169E-02
+-1.86505E-02
+-4.12172E-02
+-4.92979E-02
+2.69538E-02
+-2.66457E-02
+6.73594E-02
+8.20067E-02
+1.23491E-02
+-7.82045E-02
+5.59781E-03
+1.69580E-02
+-3.36624E-03
+-2.40832E-02
+5.55433E-02
+1.45973E-02
+-1.31928E-02
+-1.57636E-02
+3.78321E-02
+-1.46667E-01
+5.68449E-02
+-5.81680E-02
+1.44393E-02
+9.01460E-03
+-2.10025E-02
+-1.45872E-02
+4.07156E-03
+3.95902E-02
+4.00640E-02
+4.22339E-02
+APP_abinitio.fitout
+6 6
+7.67218E-01
+2.33457E-01
+-1.46363E-01
+2.05224E-01
+1.06638E-01
+-6.77313E-02
+-6.30041E-01
+-3.81159E-01
+-5.15514E-01
+-2.88904E-01
+-2.90076E-02
+-2.70035E-02
+3.70005E-03
+1.95656E-02
+1.25674E-02
+3.59916E-02
+1.37024E-02
+1.01801E-02
+-3.80359E-02
+5.47337E-02
+1.01848E-02
+-2.57535E-02
+2.66764E-03
+-4.98453E-04
+7.52962E-01
+2.85588E-02
+-2.30346E-03
+-7.23259E-01
+8.10079E-02
+-1.05684E-01
+3.68628E-02
+-1.91424E-03
+2.12266E-01
+-6.32729E-02
+-1.02146E-01
+-2.18836E-01
+-6.31142E-03
+-2.91560E-02
+5.34587E-02
+-7.50095E-03
+7.50583E-02
+-8.13426E-02
+-2.32451E-02
+1.11562E-02
+1.13419E-01
+-3.85175E-01
+-4.29197E-01
+-1.16652E-01
+7.97811E-02
+9.57693E-02
+-8.52841E-02
+-2.21538E-02
+-2.39487E-03
+-8.70199E-03
+1.31160E-02
+-1.57490E-02
+4.10501E-03
+1.15439E-02
+-1.21521E-01
+1.06008E-01
+2.17808E-01
+-4.09860E-02
+-4.67063E-02
+-2.20038E-01
+-1.73742E-02
+-5.15139E-02
+1.45133E-02
+1.43836E-02
+1.12055E-02
+4.37722E-02
+-6.08277E-02
+-2.48175E-02
+1.34523E-02
+-5.61869E-03
+2.54898E-02
+-1.15690E-02
+8.28139E-02
+-9.28803E-02
+-7.82454E-02
+8.53303E-02
+-8.14192E-03
+4.37089E-02
+4.66433E-02
+-7.54954E-03
+Ppp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Ppa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PpG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PpA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PpP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Pap
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Paa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0
+PaG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PaA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PaP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PGp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PGa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PGG_abinitio.fitout
+6 6
+-1.17152E-02
+8.08314E-04
+-1.84929E-02
+7.36520E-03
+6.78035E-04
+1.47704E-03
+-1.70250E-02
+-1.55551E-03
+2.76943E-03
+1.08343E-02
+4.12013E-03
+-3.24136E-03
+-3.01008E-02
+2.32977E-02
+-5.15750E-02
+2.27035E-02
+6.65240E-03
+-2.18380E-03
+5.05778E-04
+1.64632E-03
+-3.06580E-03
+-1.65941E-03
+1.19256E-03
+1.70662E-04
+-3.20902E-02
+1.66855E-03
+-5.58858E-02
+3.53936E-03
+3.31006E-02
+7.91879E-02
+4.61083E-02
+1.38623E-01
+3.86473E-03
+1.37982E-02
+-9.42935E-03
+-2.21267E-02
+3.39327E-02
+-1.03076E-01
+5.55447E-02
+-1.74210E-01
+-5.80200E-02
+4.68277E-02
+1.06477E-01
+-8.69627E-02
+6.78160E-03
+-3.06871E-03
+-7.58285E-04
+-2.00872E-04
+2.83759E-02
+-1.86440E-04
+4.72643E-02
+-3.65989E-03
+-3.13701E-02
+9.28177E-03
+7.65182E-02
+-1.78706E-02
+-8.74290E-03
+7.16994E-03
+3.33572E-03
+-3.65890E-03
+-7.11420E-03
+2.06385E-03
+-8.99078E-03
+2.57823E-03
+1.33229E-02
+-6.60596E-03
+2.05639E-02
+-1.16590E-02
+-7.63912E-03
+-1.97280E-02
+1.62325E-02
+3.47846E-02
+-4.83151E-02
+-6.75034E-02
+3.99377E-03
+9.70541E-03
+5.68454E-03
+-7.85137E-03
+3.89044E-03
+2.39758E-03
+-2.02212E-03
+4.78289E-03
+6.07226E-03
+-9.28848E-03
+PGA_abinitio.fitout
+6 6
+-4.21674E-02
+-1.10869E-02
+-1.91467E-01
+3.91775E-02
+1.94958E-02
+-3.54878E-02
+-3.77851E-02
+1.27993E-02
+6.84937E-03
+5.87173E-02
+-2.19318E-02
+-1.22346E-02
+3.61826E-02
+1.58457E-01
+-8.77671E-02
+5.07088E-02
+1.99345E-02
+-5.69296E-03
+1.80979E-01
+-1.48099E-02
+4.34955E-02
+-1.06913E-01
+-4.23758E-03
+-1.12771E-02
+-9.14668E-01
+-1.27576E-01
+4.31803E-01
+8.48152E-02
+-9.29375E-02
+8.11535E-02
+7.75147E-02
+-1.54831E-02
+1.80980E-02
+3.56840E-01
+1.06163E-03
+1.84371E-01
+-3.63221E-02
+-1.54518E-01
+8.01580E-02
+-1.37715E-01
+3.75783E-02
+6.25028E-02
+-9.71596E-02
+-7.76123E-02
+-7.34366E-03
+-2.53251E-02
+-8.84402E-02
+-1.33588E-01
+2.93982E-02
+-2.04265E-02
+3.61310E-02
+2.56799E-02
+-7.53693E-02
+3.97849E-03
+9.38605E-03
+-1.03295E-02
+1.08952E-02
+-1.28963E-02
+-9.99238E-03
+-1.83259E-02
+-5.49337E-02
+-3.29320E-02
+1.76373E-02
+1.68899E-02
+-1.38793E-02
+-6.36390E-03
+1.10934E-02
+-2.37010E-02
+7.99957E-03
+-8.11723E-03
+-1.10377E-02
+1.49324E-02
+-5.20932E-02
+-5.07061E-02
+9.24310E-03
+3.59607E-02
+8.39780E-03
+2.03269E-02
+3.52271E-02
+2.07302E-02
+2.99236E-02
+1.25733E-02
+2.46004E-02
+2.57522E-02
+PGP_abinitio.fitout
+6 6
+2.03597E-01
+8.31483E-02
+-1.71381E-01
+8.54226E-02
+4.25626E-02
+-5.85819E-02
+2.22794E-01
+-2.02475E-01
+3.25632E-02
+7.33475E-02
+-7.06541E-02
+4.51681E-03
+-1.67080E-01
+2.09618E-02
+-3.46297E-02
+8.43392E-02
+-5.34880E-02
+-6.10693E-02
+2.21688E-01
+2.10843E-01
+-3.28216E-02
+-2.72329E-02
+9.86024E-02
+-1.76610E-01
+-1.60565E+00
+-4.44138E-01
+-2.78194E-01
+1.38773E+00
+-3.28544E-01
+5.61299E-01
+2.68172E-01
+-2.49111E-01
+-3.04849E-02
+3.24789E-01
+-7.02413E-02
+-5.66952E-02
+-1.99543E-02
+-1.54812E-01
+3.44999E-01
+-8.80924E-02
+2.92163E-01
+-1.06140E-01
+1.13480E-01
+-6.77349E-02
+-1.00651E-02
+7.06482E-04
+-1.16447E-02
+-9.76706E-03
+-1.99690E-01
+-4.61642E-05
+-1.88271E-01
+2.39097E-02
+8.30479E-02
+-3.04097E-03
+-5.22889E-02
+-3.72105E-02
+1.89411E-02
+5.14450E-03
+1.73122E-04
+2.81198E-02
+-5.64760E-02
+-2.87088E-02
+1.98950E-02
+-1.16766E-02
+1.01974E-01
+6.66616E-02
+-1.87044E-03
+-1.15424E-01
+4.27475E-02
+1.81373E-02
+-1.20504E-01
+-1.84873E-02
+-1.08167E-01
+1.51483E-02
+-3.84445E-02
+-1.61519E-02
+8.86730E-03
+1.79342E-02
+-6.56259E-02
+-7.10179E-03
+2.27041E-02
+-1.46227E-02
+4.06429E-02
+6.08438E-03
+PAp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PAa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PAG_abinitio.fitout
+6 6
+-1.86727E-02
+-4.32929E-03
+-2.85958E-02
+1.50859E-03
+-3.07202E-03
+-1.39642E-03
+-3.69846E-02
+1.04882E-03
+-8.93575E-03
+6.91102E-03
+8.76224E-03
+-9.37333E-03
+-2.83489E-02
+9.17861E-02
+-4.06504E-02
+3.05162E-03
+-2.76511E-02
+1.12820E-02
+-6.41226E-02
+-1.64995E-01
+1.11209E-01
+2.21353E-03
+2.39645E-02
+2.01533E-02
+-2.47345E-02
+-2.50684E-03
+-3.41213E-02
+-7.94485E-03
+1.05852E-01
+8.78292E-02
+-5.17054E-02
+1.69855E-01
+1.04445E-02
+9.78832E-04
+-2.29385E-02
+-1.47634E-02
+1.18421E-01
+-1.88867E-01
+-6.21589E-02
+-1.26185E-01
+-5.33812E-02
+5.50859E-02
+7.14998E-02
+-8.04550E-02
+-2.26018E-03
+4.68748E-03
+-2.27517E-03
+-1.43078E-02
+-1.17454E-03
+1.23677E-02
+3.28704E-02
+-1.18357E-03
+6.66522E-02
+1.62106E-02
+4.96241E-02
+2.33720E-02
+-3.02511E-02
+9.67404E-03
+2.23762E-02
+-7.04732E-03
+-4.78693E-03
+-4.81439E-03
+-1.62583E-02
+-5.39745E-03
+1.03614E-02
+-1.48080E-03
+2.00905E-03
+-1.32309E-03
+2.19363E-02
+-1.24035E-02
+2.52544E-02
+5.18711E-02
+-2.76865E-02
+-5.74752E-02
+2.65025E-02
+3.94510E-02
+2.18987E-02
+-2.19298E-02
+7.84386E-04
+1.00090E-02
+-1.75104E-05
+-7.15526E-04
+4.07597E-05
+-8.03326E-03
+PAA_abinitio.fitout
+6 6
+-1.48004E-01
+3.55924E-02
+-1.63137E-01
+-8.28753E-02
+-1.06713E-02
+-6.45243E-03
+-3.21653E-02
+6.01919E-02
+-5.05636E-02
+5.67676E-02
+7.86295E-02
+-2.22937E-02
+-5.37140E-01
+-1.70759E-01
+-9.59335E-02
+6.43856E-02
+-8.81145E-03
+-1.00916E-02
+-4.40811E-02
+-3.73393E-02
+1.03283E-01
+-3.42948E-03
+-2.17269E-03
+9.17541E-03
+-9.27782E-01
+-1.11939E-02
+3.02174E-01
+-1.36657E-02
+-7.86916E-02
+6.96968E-03
+1.56723E-01
+2.32634E-02
+-1.70987E-01
+3.76170E-01
+9.99686E-02
+2.18486E-01
+1.06825E-01
+-2.06877E-01
+-2.97818E-03
+-4.52041E-02
+3.13196E-02
+7.85678E-02
+-1.01412E-02
+-3.73421E-02
+1.72860E-02
+1.22020E-01
+2.19112E-02
+-1.30464E-01
+-1.06877E-02
+-3.86686E-02
+3.99139E-02
+5.35222E-04
+1.14431E-02
+6.04451E-02
+3.05464E-02
+1.50339E-02
+6.80800E-02
+3.31187E-03
+-6.31445E-02
+-1.72000E-02
+-7.22567E-02
+2.69832E-02
+-2.09692E-02
+-6.82591E-02
+-1.64572E-02
+-1.33289E-02
+1.19150E-02
+-1.13424E-02
+-1.07905E-02
+-1.69497E-02
+4.81347E-02
+3.76561E-02
+-5.18859E-02
+-3.22088E-02
+-1.06826E-02
+3.06852E-02
+3.58685E-04
+1.38110E-02
+-8.67692E-03
+2.35927E-02
+-1.07325E-02
+-3.63660E-03
+1.88840E-02
+2.41626E-02
+PAP_abinitio.fitout
+6 6
+-6.62363E-01
+7.32861E-02
+1.94710E-01
+9.73458E-02
+-6.28769E-03
+-4.85752E-02
+3.97715E-01
+5.49366E-01
+2.30814E-02
+-8.22262E-02
+-8.64146E-02
+-2.43265E-02
+-4.11363E-01
+3.39575E-01
+1.05049E-02
+-6.65599E-02
+1.35397E-01
+-4.58790E-02
+-5.10344E-01
+-4.70674E-01
+1.31064E-02
+-1.19879E-01
+-1.29575E-02
+1.81257E-02
+-1.09088E+00
+-4.71902E-01
+2.36092E-02
+5.38778E-01
+-6.85489E-01
+4.92470E-01
+4.01670E-01
+-2.32399E-01
+-7.00634E-02
+2.83989E-01
+-1.32175E-01
+1.66960E-02
+-2.21384E-01
+-5.96587E-02
+2.72734E-01
+-4.13095E-02
+2.44356E-01
+6.36661E-02
+1.48425E-01
+-6.93542E-02
+-3.99401E-02
+5.24989E-02
+7.15320E-02
+-1.61202E-01
+-3.72458E-01
+1.25361E-02
+-4.14212E-02
+2.01276E-02
+1.48321E-01
+4.87296E-02
+3.08984E-02
+-7.43849E-02
+2.22397E-02
+-6.30931E-02
+-2.58342E-02
+-1.17253E-01
+1.16469E-03
+-4.30153E-02
+1.00187E-01
+-3.68603E-02
+-2.84632E-02
+1.25879E-01
+-7.91799E-02
+-5.81097E-02
+1.66423E-02
+6.62823E-02
+-2.10941E-02
+-4.56423E-02
+-4.47566E-02
+-2.91895E-02
+-9.58987E-02
+5.13975E-03
+-2.11538E-02
+-6.48730E-03
+-5.39683E-02
+1.34026E-02
+5.00738E-02
+4.52056E-04
+3.10756E-02
+1.90776E-02
+PPp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PPa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PPG_abinitio.fitout
+6 6
+5.16416E-02
+8.62095E-02
+-4.51303E-02
+-1.00017E-02
+5.57916E-02
+-2.82344E-03
+1.31516E-02
+-4.56497E-02
+-2.65691E-02
+-5.08100E-02
+5.49723E-03
+-2.54623E-02
+-1.16607E-02
+8.71292E-03
+-1.04028E-01
+4.88786E-02
+4.79630E-03
+-2.09226E-02
+-6.61022E-02
+-5.72385E-02
+8.20553E-02
+2.70668E-02
+-5.44738E-03
+1.03208E-02
+-1.03496E-02
+-2.79784E-04
+-4.19006E-02
+4.74221E-02
+-5.69082E-02
+-4.47927E-02
+8.07646E-03
+1.26023E-01
+-2.52358E-02
+1.21320E-02
+-3.71100E-02
+4.15186E-02
+5.02204E-02
+-1.58502E-01
+1.32449E-01
+-1.77700E-01
+-3.00666E-02
+-1.30149E-02
+1.33461E-02
+9.21573E-02
+-2.56421E-02
+2.27861E-02
+-1.73030E-02
+3.10550E-02
+3.71478E-02
+9.60699E-02
+-4.83687E-02
+-2.61395E-02
+4.01929E-02
+-1.28377E-01
+9.48191E-02
+-8.40461E-02
+-8.84160E-03
+-1.60301E-02
+1.77380E-02
+8.24025E-02
+-2.67488E-02
+2.84035E-02
+2.91455E-03
+1.70812E-02
+1.54656E-02
+-2.40176E-02
+2.76831E-02
+-5.27066E-03
+2.65944E-02
+4.76662E-02
+-4.29954E-02
+-4.04722E-02
+3.02425E-02
+-8.97502E-02
+5.58394E-02
+-5.79944E-02
+-1.14607E-02
+-5.29796E-02
+2.89066E-02
+7.05200E-02
+-3.06923E-02
+1.81037E-02
+-7.19964E-03
+1.24096E-02
+PPA_abinitio.fitout
+6 6
+1.32899E-01
+1.38710E-01
+-3.14724E-02
+1.18637E-03
+7.58328E-02
+-1.42436E-02
+-6.58168E-02
+-8.66061E-02
+-6.79304E-02
+-1.19298E-01
+-1.61300E-02
+-4.37837E-02
+-1.78927E-01
+-1.06051E-01
+-1.70089E-01
+-8.92743E-03
+-1.36536E-02
+-3.11752E-02
+-6.57058E-02
+8.25739E-02
+8.75466E-02
+4.94767E-02
+3.98260E-02
+7.59713E-03
+7.87024E-02
+-1.31469E-01
+-1.19822E-01
+-1.88926E-01
+-5.37974E-02
+-5.92828E-02
+3.66331E-02
+-7.89996E-02
+3.94025E-02
+-1.22489E-01
+-1.04510E-01
+-8.92355E-02
+3.76168E-02
+-2.21729E-01
+1.45467E-01
+-2.07780E-01
+-2.21517E-02
+-6.26420E-02
+5.66877E-02
+-8.08278E-03
+2.41190E-02
+-1.27733E-01
+-1.22023E-01
+-1.08070E-01
+6.95849E-02
+1.10424E-01
+-8.50267E-02
+-1.19373E-01
+3.44974E-02
+-1.88762E-01
+1.16889E-01
+-8.09305E-02
+-2.55243E-02
+-4.70948E-02
+3.42463E-02
+-6.88564E-03
+8.19412E-02
+-5.04839E-02
+-1.18235E-01
+-1.96336E-01
+1.05609E-02
+1.98306E-02
+-3.37354E-03
+-3.73736E-02
+4.17335E-02
+4.33113E-02
+-5.87118E-02
+-6.90520E-02
+2.89217E-02
+-1.39956E-01
+8.07782E-02
+-5.61360E-02
+-1.10854E-02
+-2.35454E-02
+-2.55672E-02
+-1.54173E-02
+5.08292E-02
+2.02335E-03
+-6.94101E-02
+-7.15398E-02
+PPP_abinitio.fitout
+6 6
+9.55141E-01
+2.82825E-01
+-1.25053E-01
+-7.72476E-02
+-1.53768E-02
+-1.23068E-01
+-1.07860E+00
+-8.60254E-01
+-5.96195E-01
+-5.73846E-01
+-1.56914E-01
+-2.25099E-01
+3.39553E-02
+5.27174E-02
+-4.70984E-04
+6.95022E-02
+3.54170E-02
+-6.13267E-04
+-7.18476E-02
+1.27301E-01
+3.65459E-02
+-6.67390E-02
+-1.04061E-02
+-1.72909E-02
+9.06180E-01
+-3.16174E-01
+-4.36987E-01
+-8.48768E-01
+1.45459E-01
+-1.03920E-01
+6.02209E-02
+-2.74110E-02
+3.85731E-01
+-4.98738E-01
+-6.17808E-01
+-3.88753E-01
+1.82417E-02
+-5.85814E-02
+8.79629E-02
+-1.17158E-02
+9.11518E-02
+-5.60357E-02
+-8.69024E-02
+7.33401E-02
+1.80120E-01
+-4.92279E-01
+-5.69457E-01
+-2.06515E-01
+7.60602E-02
+1.37305E-01
+-1.78078E-01
+-3.68494E-02
+3.13809E-02
+-4.97757E-02
+6.29741E-02
+-1.27850E-02
+3.55620E-03
+6.26058E-03
+-1.59718E-01
+1.40642E-01
+2.54613E-01
+-4.16776E-01
+-4.60517E-01
+-2.65052E-01
+4.02722E-03
+-3.47336E-02
+1.38117E-02
+4.25497E-02
+-4.07847E-03
+8.37183E-02
+-1.45515E-01
+-4.97774E-02
+2.16109E-02
+-4.19394E-02
+5.29746E-02
+-7.24450E-03
+1.28857E-03
+5.46137E-03
+-2.15724E-01
+2.27414E-01
+9.50069E-02
+-1.60267E-01
+-1.77017E-01
+-1.06693E-01
+Gpp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Gpa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GpG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GpA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GpP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Gap
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Gaa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GaG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GaA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GaP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GGp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GGa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GGG_abinitio.fitout
+6 6
+-4.46036E-03
+-5.52737E-05
+-2.18852E-02
+4.33960E-03
+1.85025E-03
+4.95855E-04
+-3.24597E-18
+-2.02014E-17
+4.26966E-17
+1.49283E-17
+2.82040E-17
+1.43925E-16
+-7.18411E-03
+7.92072E-03
+-2.17806E-02
+1.67057E-03
+3.84725E-03
+-7.17231E-05
+-7.20379E-18
+-4.20384E-17
+3.03891E-17
+-2.44087E-18
+9.05009E-17
+1.14122E-19
+-7.23813E-03
+3.28293E-03
+-1.39801E-17
+-1.89800E-17
+-1.30490E-03
+1.38772E-02
+-7.10794E-18
+3.31613E-17
+1.99808E-04
+-2.17604E-03
+3.79714E-18
+5.99313E-17
+1.26587E-02
+-2.05259E-02
+1.93061E-17
+2.66731E-17
+-3.66073E-04
+3.55213E-03
+-1.11906E-17
+-3.51248E-18
+-4.59544E-03
+-9.54163E-03
+-1.52654E-17
+3.48214E-18
+4.86862E-03
+-9.46975E-05
+-1.03323E-17
+-5.01265E-17
+-5.04727E-03
+2.44410E-03
+5.33373E-17
+-1.92922E-17
+2.77079E-03
+2.66047E-02
+2.43917E-17
+-3.41779E-17
+3.19743E-03
+-5.00879E-03
+-1.09699E-17
+3.58392E-17
+4.68052E-03
+6.60091E-04
+3.91014E-17
+3.76848E-17
+-3.19326E-04
+2.01276E-04
+-2.20174E-17
+4.87189E-17
+-2.79146E-02
+-6.79853E-02
+1.69022E-17
+-5.18855E-17
+2.15742E-02
+-4.87803E-02
+2.70102E-17
+-5.03522E-18
+1.86440E-04
+-1.17222E-04
+4.25151E-18
+-1.13210E-17
+GGA_abinitio.fitout
+6 6
+-8.46686E-03
+2.62291E-03
+-3.30608E-02
+7.96177E-03
+7.15867E-03
+3.66237E-03
+-7.72008E-04
+-2.02845E-04
+1.62280E-03
+8.54876E-04
+-8.49190E-05
+-5.78018E-04
+-8.84141E-03
+1.13979E-02
+-2.03566E-02
+3.55919E-03
+-3.02160E-04
+-1.26198E-03
+3.52817E-02
+8.28932E-03
+4.25960E-03
+-6.41946E-03
+3.60893E-04
+3.94028E-04
+-1.35465E-02
+4.17716E-03
+8.28258E-02
+1.61355E-02
+8.31401E-03
+1.44081E-02
+1.30831E-02
+-1.86039E-04
+4.43630E-04
+-2.39715E-03
+-3.58254E-03
+-6.00048E-03
+9.14174E-03
+-8.39490E-03
+-6.77698E-04
+1.68410E-03
+-1.67024E-03
+-2.02726E-03
+5.39827E-04
+-8.64739E-04
+-7.91196E-03
+-1.68047E-02
+6.92503E-02
+-8.67821E-02
+3.44170E-03
+8.48263E-04
+3.84754E-03
+2.18871E-03
+-4.09461E-03
+-4.34824E-04
+9.37754E-04
+-2.35904E-03
+1.07912E-02
+3.20561E-02
+-9.95653E-04
+-1.99915E-02
+9.98327E-03
+-8.10801E-03
+-7.96601E-02
+-5.92419E-02
+-4.49290E-03
+1.00540E-03
+-1.71946E-03
+1.14316E-03
+7.92241E-05
+6.05998E-04
+3.78783E-04
+2.15709E-03
+-1.93128E-02
+-5.56518E-02
+-1.58545E-03
+8.20476E-03
+1.36469E-02
+-5.13838E-02
+2.85169E-03
+3.67851E-02
+1.00136E-03
+-2.29434E-03
+-4.83317E-03
+-3.92162E-03
+GGP_abinitio.fitout
+6 6
+6.50508E-02
+-5.89761E-03
+-8.27980E-02
+3.10265E-02
+-3.38962E-03
+6.73743E-04
+-4.16459E-04
+1.00721E-04
+-4.62636E-05
+-3.49098E-03
+-9.41654E-05
+2.71000E-04
+-2.72775E-02
+-8.83340E-04
+-4.37583E-03
+8.50974E-04
+-2.05161E-03
+-2.06744E-03
+3.15235E-02
+-2.21357E-03
+-3.64561E-03
+5.73169E-05
+2.49261E-03
+-7.51165E-03
+-5.84943E-02
+3.97047E-02
+6.60434E-02
+4.52294E-02
+3.58092E-03
+4.17625E-03
+3.31956E-02
+-2.90994E-02
+-4.75490E-04
+-3.86681E-03
+9.36791E-04
+-4.89921E-03
+1.88598E-02
+-4.42295E-03
+1.27901E-02
+2.50468E-03
+-1.10166E-03
+-6.43327E-04
+-7.04067E-03
+6.64965E-04
+1.12632E-02
+1.15798E-01
+-6.89545E-03
+1.25577E-01
+-1.23962E-02
+-6.95734E-03
+3.95449E-03
+-1.02346E-03
+-3.67717E-03
+-3.65381E-03
+-2.61415E-03
+3.03406E-03
+-3.43625E-03
+1.68217E-02
+1.47508E-02
+-1.73684E-01
+-1.75592E-02
+-2.20142E-02
+2.19531E-02
+-1.87179E-02
+-3.30625E-04
+-2.47359E-03
+-7.45116E-04
+-5.56124E-03
+2.54898E-03
+6.05867E-03
+-6.14297E-04
+1.24640E-03
+-1.17324E-02
+-1.19361E-02
+-1.28656E-02
+8.59711E-03
+-5.18731E-03
+7.58531E-03
+-5.10999E-02
+-4.52666E-02
+3.56032E-03
+3.85053E-03
+-3.65556E-03
+4.00514E-03
+GAp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GAa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GAG_abinitio.fitout
+6 6
+-6.97988E-03
+2.11326E-04
+-3.29929E-02
+2.08293E-03
+-4.84598E-04
+2.83986E-03
+9.70254E-04
+-7.61772E-04
+8.32483E-03
+-1.25922E-02
+-3.45241E-04
+-1.52497E-03
+-5.32058E-03
+5.03716E-03
+-2.67407E-02
+-2.11211E-05
+-1.71645E-03
+1.11761E-03
+-2.54093E-03
+-1.64341E-02
+-4.08188E-03
+1.68737E-03
+3.46628E-03
+1.20169E-03
+-1.55672E-03
+-3.27573E-04
+1.18295E-03
+-2.74125E-03
+6.89273E-03
+5.44083E-03
+-1.06816E-02
+8.72845E-03
+-1.39272E-03
+-7.16988E-04
+-5.68857E-04
+4.75969E-04
+-7.24911E-03
+-1.12901E-02
+-1.80540E-02
+4.99824E-05
+8.68770E-04
+-6.51696E-05
+-2.98370E-03
+2.10040E-03
+-2.04330E-04
+-7.57116E-03
+2.10429E-02
+-7.44830E-03
+-9.90521E-04
+1.59562E-03
+2.77189E-03
+-4.99094E-03
+-3.89075E-03
+-2.77476E-03
+-3.39128E-04
+-1.78779E-03
+-1.03301E-02
+3.30321E-02
+-4.51590E-02
+4.76619E-02
+-8.44719E-03
+4.88342E-03
+-1.77256E-03
+-4.07904E-03
+1.13761E-03
+2.28837E-03
+-1.92714E-03
+1.78535E-03
+-4.62403E-05
+6.19471E-04
+-1.29644E-03
+1.96271E-03
+-4.78698E-02
+-8.10150E-02
+-1.16577E-02
+2.11953E-02
+1.00520E-02
+-2.99724E-02
+-4.03191E-03
+2.14660E-03
+1.38202E-04
+-5.32398E-04
+-1.11994E-04
+-1.14554E-03
+GAA_abinitio.fitout
+6 6
+-2.48643E-02
+4.81526E-03
+-4.51777E-02
+2.76581E-02
+-3.04518E-04
+4.36235E-03
+2.17373E-02
+-2.41799E-03
+5.37657E-03
+-2.98047E-02
+-2.52519E-03
+-5.47296E-03
+-1.58140E-02
+-6.16724E-03
+-2.25130E-02
+3.11349E-03
+-3.17153E-03
+-1.31443E-03
+1.86030E-02
+-6.75956E-03
+-7.49626E-03
+-4.44882E-03
+-1.60960E-03
+9.74007E-04
+-1.32605E-02
+-1.96368E-02
+4.43606E-02
+-1.22538E-02
+1.29574E-02
+9.39390E-04
+1.40595E-02
+-4.89717E-04
+-1.80880E-03
+2.35373E-03
+4.11257E-03
+-2.89455E-03
+5.20837E-03
+-7.82527E-03
+-8.10209E-03
+2.09822E-03
+-5.59362E-04
+6.21779E-04
+1.19200E-03
+3.27965E-03
+6.34701E-02
+-9.55843E-02
+7.89666E-02
+-8.69198E-02
+3.85345E-03
+-1.60660E-03
+6.48502E-03
+-2.57321E-04
+-2.77198E-03
+-2.22505E-03
+8.32341E-04
+4.61294E-04
+4.76449E-02
+7.77390E-03
+-4.89485E-02
+9.50176E-04
+-5.87483E-03
+4.95926E-02
+-1.88864E-02
+-2.20620E-02
+3.61988E-03
+3.72923E-04
+-9.01274E-04
+2.99111E-03
+3.87878E-04
+1.49061E-03
+-7.52049E-04
+2.88234E-03
+-2.78368E-02
+-5.76616E-02
+-4.44516E-02
+3.52391E-02
+1.62484E-02
+-1.33273E-02
+-2.03980E-02
+5.05344E-03
+-3.74544E-03
+1.78525E-03
+-2.40457E-03
+-3.94525E-03
+GAP_abinitio.fitout
+6 6
+2.16431E-02
+1.01017E-02
+1.08876E-01
+7.34409E-02
+-1.09853E-02
+9.79247E-04
+4.98205E-02
+-6.29565E-03
+-1.37799E-01
+-4.10759E-02
+-2.46017E-03
+4.80559E-04
+1.58473E-03
+1.75081E-02
+-1.54409E-03
+2.42574E-03
+4.26959E-03
+-1.02854E-03
+1.03385E-02
+-1.57001E-02
+2.02765E-03
+1.87894E-04
+2.12208E-04
+3.96201E-03
+-2.27239E-02
+7.46691E-03
+4.17848E-02
+2.41565E-02
+7.74483E-03
+1.63369E-03
+3.93859E-02
+-2.02702E-02
+-2.30982E-03
+-1.11121E-03
+3.03768E-03
+-1.78059E-03
+-1.92422E-03
+1.47293E-03
+1.85671E-02
+-7.15950E-03
+1.11092E-03
+4.77112E-03
+-3.82005E-03
+6.27417E-04
+1.72187E-01
+8.48193E-02
+5.13335E-02
+2.56405E-02
+-1.86126E-02
+-5.59752E-03
+1.12581E-02
+-9.64361E-04
+6.46262E-04
+-2.31221E-03
+-2.40846E-03
+1.63797E-03
+4.54519E-02
+-5.91857E-02
+-1.00474E-02
+-1.40591E-01
+1.81836E-02
+-2.04773E-02
+6.49684E-02
+3.16978E-02
+-9.31686E-03
+7.77812E-03
+3.71733E-03
+-8.83382E-03
+1.26015E-03
+7.59054E-04
+4.31544E-04
+-8.23154E-04
+1.76513E-02
+-3.27815E-03
+-1.17345E-01
+5.62252E-02
+-1.81312E-02
+4.04796E-02
+-5.55516E-02
+-1.14203E-02
+9.65567E-04
+1.29682E-03
+-5.17527E-03
+1.59823E-03
+GPp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GPa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+GPG_abinitio.fitout
+6 6
+-7.59230E-03
+-1.94552E-04
+-5.12944E-03
+3.20316E-03
+8.28212E-04
+-6.12307E-03
+-4.85436E-03
+-3.86687E-03
+5.36871E-03
+-3.81511E-03
+1.48924E-03
+-4.03381E-03
+-1.52037E-03
+5.72898E-03
+-2.07652E-02
+6.52919E-03
+-1.05532E-03
+-1.52037E-03
+-4.69660E-03
+-6.37003E-03
+6.17507E-03
+8.39488E-03
+-1.76904E-03
+8.05108E-03
+-7.98866E-04
+3.32411E-03
+5.85637E-03
+-2.98884E-03
+-1.20298E-02
+-3.13323E-02
+-1.89038E-02
+-1.18356E-02
+-3.02395E-03
+2.43800E-03
+1.38867E-03
+3.67447E-04
+-2.66576E-02
+5.35667E-02
+3.22822E-02
+-4.20099E-02
+-2.76279E-03
+2.73121E-03
+-2.92878E-03
+-3.97641E-03
+5.10833E-03
+-2.13602E-03
+1.80687E-03
+-2.80887E-03
+1.75818E-02
+4.47119E-03
+1.94978E-03
+3.01372E-02
+2.73577E-03
+2.07557E-04
+1.02968E-02
+1.06706E-02
+-1.07292E-02
+-4.42466E-03
+1.03471E-02
+3.26791E-02
+-3.88377E-03
+1.21678E-03
+-7.07759E-04
+1.76439E-03
+-9.08375E-03
+1.74263E-03
+2.04763E-03
+-1.17705E-02
+2.75141E-03
+-6.02603E-03
+-3.58871E-03
+-2.79755E-03
+2.32403E-02
+-5.25103E-02
+1.58724E-02
+-4.78943E-02
+-6.49652E-04
+-1.30365E-02
+-1.24654E-02
+-1.82200E-02
+-1.17657E-04
+-1.48990E-03
+-1.15237E-03
+-1.71433E-04
+GPA_abinitio.fitout
+6 6
+1.76833E-02
+9.87411E-04
+-1.16391E-02
+1.75059E-02
+-6.24857E-03
+-5.22125E-03
+5.86116E-03
+-2.32421E-03
+-1.58225E-02
+7.43433E-03
+4.93394E-03
+-1.89631E-02
+-3.25799E-02
+-2.32415E-02
+-3.10454E-02
+-8.91640E-03
+-7.65479E-03
+-1.27982E-03
+-9.57172E-03
+8.77571E-03
+6.27899E-03
+9.87208E-03
+1.16529E-02
+8.66130E-03
+1.92294E-02
+7.29907E-02
+5.17669E-02
+-2.27035E-02
+6.60383E-03
+2.08489E-02
+-4.78442E-03
+-1.14125E-02
+-4.55136E-03
+7.67130E-03
+1.21326E-02
+1.30754E-02
+-2.73899E-02
+6.76535E-02
+3.83851E-02
+-5.38365E-02
+5.65271E-03
+-5.37810E-04
+6.97088E-03
+7.72153E-03
+2.74100E-02
+-4.70346E-02
+-2.35609E-02
+-5.97889E-02
+3.49906E-02
+3.87800E-02
+2.36671E-02
+3.21217E-02
+3.57922E-03
+3.74125E-03
+1.26977E-02
+1.32543E-02
+-1.83147E-02
+-2.62166E-02
+-5.76715E-03
+-1.83617E-02
+-6.49200E-03
+4.85116E-02
+2.89452E-02
+1.35099E-02
+8.75516E-03
+2.02179E-02
+4.64394E-03
+7.27739E-03
+4.54111E-03
+-3.11138E-03
+2.01092E-03
+2.90504E-03
+2.15560E-02
+-7.36577E-02
+1.91254E-02
+-5.61390E-02
+8.99074E-03
+1.35028E-02
+4.72925E-04
+-3.22296E-03
+-9.02731E-03
+-1.73413E-02
+-9.12097E-03
+1.71334E-02
+GPP_abinitio.fitout
+6 6
+3.55782E-01
+-2.50590E-02
+-2.02781E-02
+1.88976E-01
+-6.17588E-02
+-2.91385E-02
+2.05426E-01
+-5.87770E-03
+-3.30234E-01
+7.61622E-02
+5.16617E-02
+-1.35195E-01
+-3.32316E-03
+4.73134E-03
+9.83600E-03
+1.38053E-02
+2.51654E-03
+6.42227E-03
+-1.41234E-02
+1.49387E-02
+-1.78714E-04
+-5.87439E-03
+1.08191E-03
+-1.97401E-04
+1.05038E-01
+2.91828E-01
+3.05333E-01
+-1.00742E-01
+-1.01012E-02
+-1.49328E-02
+3.57843E-02
+-4.10710E-02
+-2.83700E-02
+-2.30223E-03
+-2.43599E-03
+2.05701E-02
+-1.78170E-02
+-4.82013E-03
+6.73937E-03
+-3.38595E-03
+3.28851E-03
+-8.08060E-03
+-7.65673E-03
+5.09805E-03
+1.57284E-01
+-2.00083E-01
+-2.14765E-01
+-1.45042E-01
+3.67444E-02
+1.99603E-02
+2.72203E-02
+3.10215E-02
+-2.08581E-03
+-1.66514E-03
+-4.79677E-03
+-3.68899E-03
+1.95962E-02
+-5.81681E-03
+-5.65331E-02
+5.75209E-02
+6.98914E-02
+1.28428E-01
+1.38299E-01
+-6.32083E-02
+-1.13264E-02
+-1.45981E-02
+-6.42446E-03
+-1.94911E-02
+-8.84309E-03
+-2.46011E-03
+-2.28773E-03
+-5.77668E-03
+3.89634E-03
+9.03634E-03
+9.92089E-03
+-6.47315E-03
+4.81635E-02
+-5.01397E-02
+3.40332E-02
+-2.73013E-02
+-6.50338E-02
+9.57053E-04
+3.41943E-03
+5.80621E-02
+App
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Apa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+ApG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+ApA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+ApP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Aap
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Aaa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AaG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AaA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AaP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AGp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AGa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AGG_abinitio.fitout
+6 6
+-7.39715E-03
+3.65556E-03
+-1.91078E-02
+2.38136E-03
+-5.37961E-03
+3.61998E-03
+-2.18464E-02
+-4.41896E-03
+5.08921E-03
+3.21624E-03
+5.89060E-03
+1.19802E-03
+-2.19281E-02
+1.49645E-02
+-3.56136E-02
+9.88455E-03
+1.22315E-02
+-3.00099E-03
+2.21176E-03
+3.15110E-03
+-6.61995E-03
+-3.20593E-03
+1.54183E-03
+-8.59009E-04
+-5.99348E-03
+4.65389E-05
+-5.95310E-02
+1.06111E-02
+1.57623E-02
+1.59041E-02
+1.00177E-04
+1.29357E-01
+4.43290E-03
+2.08145E-02
+-4.15715E-03
+1.26761E-03
+1.04441E-02
+-2.89679E-02
+6.04392E-02
+-1.59174E-01
+-5.14944E-02
+1.87340E-02
+3.85570E-02
+-3.73199E-02
+7.51987E-04
+-4.77918E-03
+4.26491E-03
+4.13725E-03
+9.71818E-03
+9.73285E-04
+3.70014E-02
+-2.94368E-04
+-4.04141E-03
+3.29676E-03
+4.40040E-02
+-7.37677E-03
+-2.53272E-03
+1.62907E-02
+1.14612E-03
+-1.48734E-02
+-3.33955E-03
+-2.38826E-03
+-8.32102E-03
+1.27157E-03
+4.85104E-03
+4.73026E-04
+2.63678E-02
+-1.07751E-02
+-2.39741E-03
+-1.84784E-02
+5.03233E-03
+1.06810E-02
+-2.90718E-02
+-5.45636E-02
+4.72574E-03
+3.35612E-02
+1.56482E-02
+-2.86180E-02
+7.62139E-03
+-1.87081E-02
+-1.26071E-03
+2.55160E-03
+3.08408E-03
+6.64625E-04
+AGA_abinitio.fitout
+6 6
+-2.10939E-02
+6.03383E-02
+-5.52248E-02
+1.30827E-02
+-1.40094E-02
+1.47428E-02
+-4.01256E-02
+-2.37278E-02
+-2.87652E-02
+3.68246E-02
+2.93335E-02
+4.16440E-04
+1.17436E-02
+1.16773E-01
+-5.62224E-02
+2.90645E-02
+9.23893E-03
+9.04317E-04
+8.25772E-02
+-8.68232E-03
+5.76473E-02
+-5.38481E-02
+4.24119E-03
+-2.74517E-03
+-7.20633E-01
+-1.75249E-01
+2.57130E-03
+1.87130E-02
+-3.59686E-02
+6.77989E-02
+2.08938E-02
+1.32158E-01
+-1.08429E-02
+7.44245E-02
+-1.44520E-03
+1.01815E-01
+-3.94152E-02
+-3.71396E-02
+5.52084E-02
+-1.15046E-01
+1.03393E-01
+3.15823E-02
+4.11136E-03
+2.05022E-02
+3.50590E-02
+-3.83513E-02
+2.49373E-02
+-8.93949E-02
+1.20974E-02
+-4.74375E-03
+1.72335E-02
+1.47956E-02
+-7.11183E-03
+-1.84173E-03
+2.06808E-02
+-7.58410E-03
+1.52949E-02
+-5.81264E-03
+-5.96499E-04
+-2.82883E-02
+-1.02857E-02
+2.25503E-02
+-2.57388E-02
+-2.87509E-02
+-8.97673E-03
+-8.26320E-03
+-5.34274E-03
+-1.16850E-02
+2.07672E-02
+-3.82975E-03
+-1.01681E-02
+6.05167E-03
+-5.72931E-02
+-4.97570E-02
+1.03005E-03
+4.70162E-02
+-1.06617E-03
+-2.36271E-02
+1.88343E-03
+-3.47428E-03
+1.31956E-03
+-2.15396E-02
+9.27436E-03
+-7.35804E-03
+AGP_abinitio.fitout
+6 6
+1.31640E-01
+2.58075E-01
+-8.58162E-02
+1.38005E-02
+-1.71596E-03
+5.98072E-03
+3.33712E-01
+2.93859E-02
+-3.00524E-02
+7.07210E-03
+7.75509E-03
+-8.65250E-03
+-8.74605E-02
+4.48412E-03
+-4.66310E-02
+3.06797E-02
+-2.19807E-02
+-1.58779E-02
+1.05553E-01
+6.43617E-02
+-4.08730E-02
+-2.03012E-03
+3.47495E-02
+-7.22030E-02
+-7.17920E-01
+-5.04379E-01
+-5.06488E-01
+5.31765E-01
+-1.50917E-01
+4.25707E-01
+5.91572E-02
+7.35844E-03
+-6.86600E-02
+8.34443E-03
+2.77907E-02
+2.06332E-02
+-1.13316E-01
+-6.99060E-02
+2.08245E-01
+-8.67383E-02
+-2.09219E-02
+2.28999E-02
+1.01289E-01
+-7.87242E-02
+-9.80531E-04
+1.26618E-01
+-4.27857E-02
+1.02030E-01
+-5.09267E-02
+-7.30891E-03
+-1.41543E-01
+-1.16047E-02
+6.00033E-02
+-7.06215E-02
+1.06774E-01
+4.95735E-02
+6.20532E-03
+-1.58945E-02
+1.13562E-02
+-1.70360E-01
+-2.32910E-02
+-6.32569E-03
+1.29374E-02
+-2.33028E-02
+7.41635E-02
+6.90855E-02
+3.93931E-02
+-4.32746E-02
+6.13297E-02
+3.73997E-02
+-4.94826E-03
+-3.71951E-03
+-2.28195E-03
+-6.09570E-03
+-8.25585E-03
+2.80583E-02
+5.84276E-03
+-3.18833E-02
+-3.95278E-02
+-3.23312E-02
+5.37754E-03
+3.70124E-03
+-5.30540E-05
+2.20955E-02
+AAp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AAa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+AAG_abinitio.fitout
+6 6
+-1.70504E-02
+-1.25725E-02
+-2.16572E-02
+9.17676E-03
+-7.82222E-03
+4.46235E-04
+-2.89066E-02
+-3.18047E-03
+1.97174E-03
+8.13914E-05
+2.44927E-03
+-7.40994E-04
+-1.96244E-02
+4.33134E-02
+-3.37035E-02
+7.64787E-03
+-2.13924E-02
+9.52598E-03
+-3.70523E-02
+-7.40036E-02
+2.40140E-02
+-5.10030E-03
+1.59236E-02
+7.43181E-03
+-1.31910E-02
+-7.98080E-03
+-3.09633E-02
+-7.51476E-03
+1.27383E-01
+-5.50401E-02
+5.28821E-03
+1.30486E-01
+7.11103E-03
+3.19300E-03
+-6.44252E-03
+-1.02063E-02
+6.87598E-02
+-8.99360E-02
+-6.87923E-02
+-1.23161E-01
+2.60656E-03
+-1.74383E-02
+3.18957E-02
+1.71278E-02
+6.27281E-03
+-3.34598E-03
+9.42044E-03
+-7.47981E-03
+-8.29916E-03
+2.89364E-02
+8.84389E-03
+-1.37194E-03
+3.56810E-02
+9.28219E-03
+3.85979E-02
+-2.62633E-02
+-1.93497E-02
+2.38321E-02
+-2.48270E-02
+4.03903E-02
+-3.67611E-04
+-4.36283E-03
+-9.47872E-03
+3.38944E-04
+3.18600E-03
+-4.88922E-03
+5.87363E-03
+-1.14359E-03
+-1.15826E-03
+-3.43521E-03
+7.57689E-03
+3.94852E-03
+-3.77955E-02
+-5.78792E-02
+1.21638E-02
+2.92502E-02
+1.65337E-02
+-2.18462E-02
+-9.40339E-03
+5.38053E-03
+3.03565E-03
+-1.07824E-04
+4.67967E-03
+-2.39506E-03
+AAA_abinitio.fitout
+6 6
+-2.05481E-01
+-8.26191E-02
+-6.95965E-02
+2.35665E-02
+-2.42723E-02
+1.32714E-03
+-8.89860E-02
+3.10180E-02
+-1.78645E-02
+-1.17453E-02
+6.66812E-03
+2.24833E-03
+-2.12827E-01
+-8.67822E-02
+-6.42470E-02
+2.51073E-02
+-1.87957E-02
+-2.67943E-03
+2.43460E-02
+-6.09457E-02
+3.04366E-02
+1.66922E-04
+-4.91981E-04
+4.11041E-03
+-5.57253E-01
+3.53999E-02
+-3.54366E-02
+-2.08256E-01
+-1.00701E-01
+5.79432E-03
+8.35129E-02
+7.78467E-02
+-1.23620E-01
+5.25424E-02
+-7.34795E-02
+3.79700E-02
+3.56303E-02
+-1.04002E-01
+1.77074E-02
+-6.40671E-02
+1.61121E-02
+5.81681E-02
+2.14007E-02
+6.03380E-02
+3.14908E-02
+-8.44041E-02
+-1.35241E-02
+-5.63823E-02
+-1.98125E-02
+-1.07643E-02
+3.85091E-02
+2.61695E-03
+2.67479E-02
+2.91781E-02
+3.32585E-02
+3.74364E-03
+3.95839E-02
+3.44728E-02
+-3.61017E-02
+2.51821E-02
+-1.16023E-02
+9.06885E-03
+-3.49017E-02
+-2.81423E-02
+-7.23186E-03
+-2.86460E-02
+1.58622E-02
+-1.63989E-03
+-1.90603E-02
+-6.01905E-03
+1.24698E-02
+2.47730E-03
+-3.53859E-02
+-1.77117E-02
+-2.28246E-02
+4.82700E-02
+1.54564E-02
+-5.67317E-03
+-7.01419E-03
+4.19916E-03
+7.22849E-03
+-3.30823E-02
+-8.92923E-03
+-1.67293E-02
+AAP_abinitio.fitout
+6 6
+-2.60341E-01
+-1.76464E-01
+-2.58402E-02
+5.08365E-02
+9.37785E-03
+1.34683E-02
+9.21716E-02
+4.11687E-02
+-8.52362E-02
+-1.16991E-02
+-2.83520E-02
+-2.85315E-02
+-1.15678E-01
+1.39526E-01
+1.43150E-03
+-1.24315E-02
+3.67833E-02
+-1.29502E-02
+-1.49832E-01
+-1.70965E-01
+-5.62158E-03
+-2.47031E-02
+-2.08181E-02
+1.00187E-02
+-5.09741E-01
+-3.04707E-01
+-1.77654E-01
+6.17878E-02
+-3.70183E-01
+3.16793E-01
+2.12524E-02
+2.52440E-02
+-6.09105E-02
+-8.32189E-02
+-3.45513E-02
+9.52014E-03
+-1.82716E-01
+1.57791E-02
+5.23627E-02
+-4.44659E-03
+-6.04463E-02
+7.35341E-02
+6.70677E-02
+-5.01091E-02
+1.58910E-01
+7.43892E-02
+-4.84020E-03
+2.64253E-02
+-1.74718E-01
+2.79791E-02
+-1.48718E-01
+-1.55775E-02
+2.66950E-02
+-1.89974E-02
+9.19990E-02
+-1.34598E-02
+2.58915E-02
+-8.32852E-02
+-1.21375E-02
+-1.35189E-01
+-2.34513E-03
+2.60660E-02
+2.74800E-02
+2.42119E-02
+-1.05766E-02
+1.07458E-01
+-9.49010E-02
+4.41061E-02
+2.37077E-04
+1.65542E-02
+5.17467E-03
+-2.63125E-02
+4.40376E-03
+-5.38761E-03
+-9.77473E-02
+6.71620E-02
+2.48552E-03
+-2.73268E-03
+-6.93468E-02
+-2.85697E-02
+-4.28702E-04
+1.93884E-02
+3.32264E-03
+6.53017E-03
+APp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+APa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+APG_abinitio.fitout
+6 6
+2.11374E-02
+2.96387E-02
+7.90197E-03
+2.06850E-02
+2.34088E-02
+6.25406E-03
+-6.08740E-03
+-1.55375E-03
+-5.36785E-03
+-3.11711E-02
+-9.42745E-03
+-9.88549E-03
+-2.33313E-03
+1.05580E-02
+-6.05641E-02
+2.53598E-02
+-3.71104E-03
+-5.01697E-03
+-2.89535E-02
+-2.77652E-02
+3.61965E-02
+2.00866E-02
+-3.16730E-03
+1.42707E-02
+-2.45496E-03
+-6.39974E-03
+-5.31496E-03
+1.49525E-02
+-5.33448E-02
+-6.12883E-02
+1.64008E-02
+5.17439E-02
+-1.32435E-02
+8.86931E-04
+-1.93406E-02
+1.29406E-02
+3.93005E-02
+-4.91735E-02
+8.90276E-02
+-1.29607E-01
+-4.02511E-02
+2.03071E-03
+-3.23640E-03
+1.57962E-02
+2.20466E-03
+4.80277E-03
+-7.51544E-03
+1.37773E-02
+2.79080E-02
+6.38630E-02
+-3.29278E-02
+6.21700E-03
+2.75330E-02
+-6.78977E-02
+7.03002E-02
+-2.20345E-02
+-1.90882E-02
+1.59821E-03
+1.08734E-02
+6.23326E-02
+-1.95481E-02
+6.49369E-03
+1.63273E-03
+1.99442E-02
+1.58751E-03
+-2.23369E-02
+2.25455E-02
+-1.50345E-02
+2.31505E-02
+2.36684E-02
+-2.15136E-02
+-2.13738E-02
+4.16797E-02
+-1.06406E-01
+4.31283E-02
+-5.98608E-02
+-2.21497E-02
+-4.07299E-02
+-4.56594E-04
+-1.37761E-03
+-1.48161E-02
+3.23533E-03
+-1.19902E-02
+1.20310E-02
+APA_abinitio.fitout
+6 6
+8.58406E-02
+7.11000E-02
+-2.46047E-03
+6.50547E-02
+6.13304E-02
+4.29823E-03
+-6.05849E-02
+-2.44367E-03
+-3.77551E-02
+-7.23133E-02
+-3.86123E-03
+-1.07154E-02
+-1.03709E-01
+-7.05574E-02
+-9.49350E-02
+-1.53363E-02
+-1.18499E-02
+-7.61708E-03
+-2.91869E-02
+4.65680E-02
+3.62772E-02
+3.39042E-02
+3.16662E-02
+1.47463E-02
+1.52687E-01
+3.95754E-02
+-1.32083E-02
+-1.89615E-01
+-1.83152E-02
+-3.72732E-03
+-1.41424E-02
+-4.35396E-02
+3.71739E-02
+2.14553E-02
+6.33433E-02
+1.58020E-02
+4.23314E-02
+-7.04643E-02
+1.00041E-01
+-1.66457E-01
+1.68853E-02
+-4.71492E-02
+4.53786E-02
+1.36936E-02
+1.65425E-02
+-1.09111E-01
+-4.96983E-02
+-7.14614E-02
+7.25503E-02
+1.07859E-01
+-6.38044E-02
+-5.08115E-02
+2.52785E-02
+-7.96546E-02
+8.77169E-02
+-1.86505E-02
+-4.12172E-02
+-4.92979E-02
+2.69538E-02
+-2.66457E-02
+6.73594E-02
+8.20067E-02
+1.23491E-02
+-7.82045E-02
+5.59781E-03
+1.69580E-02
+-3.36624E-03
+-2.40832E-02
+5.55433E-02
+1.45973E-02
+-1.31928E-02
+-1.57636E-02
+3.78321E-02
+-1.46667E-01
+5.68449E-02
+-5.81680E-02
+1.44393E-02
+9.01460E-03
+-2.10025E-02
+-1.45872E-02
+4.07156E-03
+3.95902E-02
+4.00640E-02
+4.22339E-02
+APP_abinitio.fitout
+6 6
+7.67218E-01
+2.33457E-01
+-1.46363E-01
+2.05224E-01
+1.06638E-01
+-6.77313E-02
+-6.30041E-01
+-3.81159E-01
+-5.15514E-01
+-2.88904E-01
+-2.90076E-02
+-2.70035E-02
+3.70005E-03
+1.95656E-02
+1.25674E-02
+3.59916E-02
+1.37024E-02
+1.01801E-02
+-3.80359E-02
+5.47337E-02
+1.01848E-02
+-2.57535E-02
+2.66764E-03
+-4.98453E-04
+7.52962E-01
+2.85588E-02
+-2.30346E-03
+-7.23259E-01
+8.10079E-02
+-1.05684E-01
+3.68628E-02
+-1.91424E-03
+2.12266E-01
+-6.32729E-02
+-1.02146E-01
+-2.18836E-01
+-6.31142E-03
+-2.91560E-02
+5.34587E-02
+-7.50095E-03
+7.50583E-02
+-8.13426E-02
+-2.32451E-02
+1.11562E-02
+1.13419E-01
+-3.85175E-01
+-4.29197E-01
+-1.16652E-01
+7.97811E-02
+9.57693E-02
+-8.52841E-02
+-2.21538E-02
+-2.39487E-03
+-8.70199E-03
+1.31160E-02
+-1.57490E-02
+4.10501E-03
+1.15439E-02
+-1.21521E-01
+1.06008E-01
+2.17808E-01
+-4.09860E-02
+-4.67063E-02
+-2.20038E-01
+-1.73742E-02
+-5.15139E-02
+1.45133E-02
+1.43836E-02
+1.12055E-02
+4.37722E-02
+-6.08277E-02
+-2.48175E-02
+1.34523E-02
+-5.61869E-03
+2.54898E-02
+-1.15690E-02
+8.28139E-02
+-9.28803E-02
+-7.82454E-02
+8.53303E-02
+-8.14192E-03
+4.37089E-02
+4.66433E-02
+-7.54954E-03
+Ppp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Ppa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PpG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PpA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PpP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Pap
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+Paa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0
+PaG
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PaA
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PaP
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PGp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PGa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PGG_abinitio.fitout
+6 6
+-1.17152E-02
+8.08314E-04
+-1.84929E-02
+7.36520E-03
+6.78035E-04
+1.47704E-03
+-1.70250E-02
+-1.55551E-03
+2.76943E-03
+1.08343E-02
+4.12013E-03
+-3.24136E-03
+-3.01008E-02
+2.32977E-02
+-5.15750E-02
+2.27035E-02
+6.65240E-03
+-2.18380E-03
+5.05778E-04
+1.64632E-03
+-3.06580E-03
+-1.65941E-03
+1.19256E-03
+1.70662E-04
+-3.20902E-02
+1.66855E-03
+-5.58858E-02
+3.53936E-03
+3.31006E-02
+7.91879E-02
+4.61083E-02
+1.38623E-01
+3.86473E-03
+1.37982E-02
+-9.42935E-03
+-2.21267E-02
+3.39327E-02
+-1.03076E-01
+5.55447E-02
+-1.74210E-01
+-5.80200E-02
+4.68277E-02
+1.06477E-01
+-8.69627E-02
+6.78160E-03
+-3.06871E-03
+-7.58285E-04
+-2.00872E-04
+2.83759E-02
+-1.86440E-04
+4.72643E-02
+-3.65989E-03
+-3.13701E-02
+9.28177E-03
+7.65182E-02
+-1.78706E-02
+-8.74290E-03
+7.16994E-03
+3.33572E-03
+-3.65890E-03
+-7.11420E-03
+2.06385E-03
+-8.99078E-03
+2.57823E-03
+1.33229E-02
+-6.60596E-03
+2.05639E-02
+-1.16590E-02
+-7.63912E-03
+-1.97280E-02
+1.62325E-02
+3.47846E-02
+-4.83151E-02
+-6.75034E-02
+3.99377E-03
+9.70541E-03
+5.68454E-03
+-7.85137E-03
+3.89044E-03
+2.39758E-03
+-2.02212E-03
+4.78289E-03
+6.07226E-03
+-9.28848E-03
+PGA_abinitio.fitout
+6 6
+-4.21674E-02
+-1.10869E-02
+-1.91467E-01
+3.91775E-02
+1.94958E-02
+-3.54878E-02
+-3.77851E-02
+1.27993E-02
+6.84937E-03
+5.87173E-02
+-2.19318E-02
+-1.22346E-02
+3.61826E-02
+1.58457E-01
+-8.77671E-02
+5.07088E-02
+1.99345E-02
+-5.69296E-03
+1.80979E-01
+-1.48099E-02
+4.34955E-02
+-1.06913E-01
+-4.23758E-03
+-1.12771E-02
+-9.14668E-01
+-1.27576E-01
+4.31803E-01
+8.48152E-02
+-9.29375E-02
+8.11535E-02
+7.75147E-02
+-1.54831E-02
+1.80980E-02
+3.56840E-01
+1.06163E-03
+1.84371E-01
+-3.63221E-02
+-1.54518E-01
+8.01580E-02
+-1.37715E-01
+3.75783E-02
+6.25028E-02
+-9.71596E-02
+-7.76123E-02
+-7.34366E-03
+-2.53251E-02
+-8.84402E-02
+-1.33588E-01
+2.93982E-02
+-2.04265E-02
+3.61310E-02
+2.56799E-02
+-7.53693E-02
+3.97849E-03
+9.38605E-03
+-1.03295E-02
+1.08952E-02
+-1.28963E-02
+-9.99238E-03
+-1.83259E-02
+-5.49337E-02
+-3.29320E-02
+1.76373E-02
+1.68899E-02
+-1.38793E-02
+-6.36390E-03
+1.10934E-02
+-2.37010E-02
+7.99957E-03
+-8.11723E-03
+-1.10377E-02
+1.49324E-02
+-5.20932E-02
+-5.07061E-02
+9.24310E-03
+3.59607E-02
+8.39780E-03
+2.03269E-02
+3.52271E-02
+2.07302E-02
+2.99236E-02
+1.25733E-02
+2.46004E-02
+2.57522E-02
+PGP_abinitio.fitout
+6 6
+2.03597E-01
+8.31483E-02
+-1.71381E-01
+8.54226E-02
+4.25626E-02
+-5.85819E-02
+2.22794E-01
+-2.02475E-01
+3.25632E-02
+7.33475E-02
+-7.06541E-02
+4.51681E-03
+-1.67080E-01
+2.09618E-02
+-3.46297E-02
+8.43392E-02
+-5.34880E-02
+-6.10693E-02
+2.21688E-01
+2.10843E-01
+-3.28216E-02
+-2.72329E-02
+9.86024E-02
+-1.76610E-01
+-1.60565E+00
+-4.44138E-01
+-2.78194E-01
+1.38773E+00
+-3.28544E-01
+5.61299E-01
+2.68172E-01
+-2.49111E-01
+-3.04849E-02
+3.24789E-01
+-7.02413E-02
+-5.66952E-02
+-1.99543E-02
+-1.54812E-01
+3.44999E-01
+-8.80924E-02
+2.92163E-01
+-1.06140E-01
+1.13480E-01
+-6.77349E-02
+-1.00651E-02
+7.06482E-04
+-1.16447E-02
+-9.76706E-03
+-1.99690E-01
+-4.61642E-05
+-1.88271E-01
+2.39097E-02
+8.30479E-02
+-3.04097E-03
+-5.22889E-02
+-3.72105E-02
+1.89411E-02
+5.14450E-03
+1.73122E-04
+2.81198E-02
+-5.64760E-02
+-2.87088E-02
+1.98950E-02
+-1.16766E-02
+1.01974E-01
+6.66616E-02
+-1.87044E-03
+-1.15424E-01
+4.27475E-02
+1.81373E-02
+-1.20504E-01
+-1.84873E-02
+-1.08167E-01
+1.51483E-02
+-3.84445E-02
+-1.61519E-02
+8.86730E-03
+1.79342E-02
+-6.56259E-02
+-7.10179E-03
+2.27041E-02
+-1.46227E-02
+4.06429E-02
+6.08438E-03
+PAp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PAa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PAG_abinitio.fitout
+6 6
+-1.86727E-02
+-4.32929E-03
+-2.85958E-02
+1.50859E-03
+-3.07202E-03
+-1.39642E-03
+-3.69846E-02
+1.04882E-03
+-8.93575E-03
+6.91102E-03
+8.76224E-03
+-9.37333E-03
+-2.83489E-02
+9.17861E-02
+-4.06504E-02
+3.05162E-03
+-2.76511E-02
+1.12820E-02
+-6.41226E-02
+-1.64995E-01
+1.11209E-01
+2.21353E-03
+2.39645E-02
+2.01533E-02
+-2.47345E-02
+-2.50684E-03
+-3.41213E-02
+-7.94485E-03
+1.05852E-01
+8.78292E-02
+-5.17054E-02
+1.69855E-01
+1.04445E-02
+9.78832E-04
+-2.29385E-02
+-1.47634E-02
+1.18421E-01
+-1.88867E-01
+-6.21589E-02
+-1.26185E-01
+-5.33812E-02
+5.50859E-02
+7.14998E-02
+-8.04550E-02
+-2.26018E-03
+4.68748E-03
+-2.27517E-03
+-1.43078E-02
+-1.17454E-03
+1.23677E-02
+3.28704E-02
+-1.18357E-03
+6.66522E-02
+1.62106E-02
+4.96241E-02
+2.33720E-02
+-3.02511E-02
+9.67404E-03
+2.23762E-02
+-7.04732E-03
+-4.78693E-03
+-4.81439E-03
+-1.62583E-02
+-5.39745E-03
+1.03614E-02
+-1.48080E-03
+2.00905E-03
+-1.32309E-03
+2.19363E-02
+-1.24035E-02
+2.52544E-02
+5.18711E-02
+-2.76865E-02
+-5.74752E-02
+2.65025E-02
+3.94510E-02
+2.18987E-02
+-2.19298E-02
+7.84386E-04
+1.00090E-02
+-1.75104E-05
+-7.15526E-04
+4.07597E-05
+-8.03326E-03
+PAA_abinitio.fitout
+6 6
+-1.48004E-01
+3.55924E-02
+-1.63137E-01
+-8.28753E-02
+-1.06713E-02
+-6.45243E-03
+-3.21653E-02
+6.01919E-02
+-5.05636E-02
+5.67676E-02
+7.86295E-02
+-2.22937E-02
+-5.37140E-01
+-1.70759E-01
+-9.59335E-02
+6.43856E-02
+-8.81145E-03
+-1.00916E-02
+-4.40811E-02
+-3.73393E-02
+1.03283E-01
+-3.42948E-03
+-2.17269E-03
+9.17541E-03
+-9.27782E-01
+-1.11939E-02
+3.02174E-01
+-1.36657E-02
+-7.86916E-02
+6.96968E-03
+1.56723E-01
+2.32634E-02
+-1.70987E-01
+3.76170E-01
+9.99686E-02
+2.18486E-01
+1.06825E-01
+-2.06877E-01
+-2.97818E-03
+-4.52041E-02
+3.13196E-02
+7.85678E-02
+-1.01412E-02
+-3.73421E-02
+1.72860E-02
+1.22020E-01
+2.19112E-02
+-1.30464E-01
+-1.06877E-02
+-3.86686E-02
+3.99139E-02
+5.35222E-04
+1.14431E-02
+6.04451E-02
+3.05464E-02
+1.50339E-02
+6.80800E-02
+3.31187E-03
+-6.31445E-02
+-1.72000E-02
+-7.22567E-02
+2.69832E-02
+-2.09692E-02
+-6.82591E-02
+-1.64572E-02
+-1.33289E-02
+1.19150E-02
+-1.13424E-02
+-1.07905E-02
+-1.69497E-02
+4.81347E-02
+3.76561E-02
+-5.18859E-02
+-3.22088E-02
+-1.06826E-02
+3.06852E-02
+3.58685E-04
+1.38110E-02
+-8.67692E-03
+2.35927E-02
+-1.07325E-02
+-3.63660E-03
+1.88840E-02
+2.41626E-02
+PAP_abinitio.fitout
+6 6
+-6.62363E-01
+7.32861E-02
+1.94710E-01
+9.73458E-02
+-6.28769E-03
+-4.85752E-02
+3.97715E-01
+5.49366E-01
+2.30814E-02
+-8.22262E-02
+-8.64146E-02
+-2.43265E-02
+-4.11363E-01
+3.39575E-01
+1.05049E-02
+-6.65599E-02
+1.35397E-01
+-4.58790E-02
+-5.10344E-01
+-4.70674E-01
+1.31064E-02
+-1.19879E-01
+-1.29575E-02
+1.81257E-02
+-1.09088E+00
+-4.71902E-01
+2.36092E-02
+5.38778E-01
+-6.85489E-01
+4.92470E-01
+4.01670E-01
+-2.32399E-01
+-7.00634E-02
+2.83989E-01
+-1.32175E-01
+1.66960E-02
+-2.21384E-01
+-5.96587E-02
+2.72734E-01
+-4.13095E-02
+2.44356E-01
+6.36661E-02
+1.48425E-01
+-6.93542E-02
+-3.99401E-02
+5.24989E-02
+7.15320E-02
+-1.61202E-01
+-3.72458E-01
+1.25361E-02
+-4.14212E-02
+2.01276E-02
+1.48321E-01
+4.87296E-02
+3.08984E-02
+-7.43849E-02
+2.22397E-02
+-6.30931E-02
+-2.58342E-02
+-1.17253E-01
+1.16469E-03
+-4.30153E-02
+1.00187E-01
+-3.68603E-02
+-2.84632E-02
+1.25879E-01
+-7.91799E-02
+-5.81097E-02
+1.66423E-02
+6.62823E-02
+-2.10941E-02
+-4.56423E-02
+-4.47566E-02
+-2.91895E-02
+-9.58987E-02
+5.13975E-03
+-2.11538E-02
+-6.48730E-03
+-5.39683E-02
+1.34026E-02
+5.00738E-02
+4.52056E-04
+3.10756E-02
+1.90776E-02
+PPp
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PPa
+0 0
+ 0 0
+ 0 0
+ 0 0
+ 0 0 
+ 0 0
+PPG_abinitio.fitout
+6 6
+5.16416E-02
+8.62095E-02
+-4.51303E-02
+-1.00017E-02
+5.57916E-02
+-2.82344E-03
+1.31516E-02
+-4.56497E-02
+-2.65691E-02
+-5.08100E-02
+5.49723E-03
+-2.54623E-02
+-1.16607E-02
+8.71292E-03
+-1.04028E-01
+4.88786E-02
+4.79630E-03
+-2.09226E-02
+-6.61022E-02
+-5.72385E-02
+8.20553E-02
+2.70668E-02
+-5.44738E-03
+1.03208E-02
+-1.03496E-02
+-2.79784E-04
+-4.19006E-02
+4.74221E-02
+-5.69082E-02
+-4.47927E-02
+8.07646E-03
+1.26023E-01
+-2.52358E-02
+1.21320E-02
+-3.71100E-02
+4.15186E-02
+5.02204E-02
+-1.58502E-01
+1.32449E-01
+-1.77700E-01
+-3.00666E-02
+-1.30149E-02
+1.33461E-02
+9.21573E-02
+-2.56421E-02
+2.27861E-02
+-1.73030E-02
+3.10550E-02
+3.71478E-02
+9.60699E-02
+-4.83687E-02
+-2.61395E-02
+4.01929E-02
+-1.28377E-01
+9.48191E-02
+-8.40461E-02
+-8.84160E-03
+-1.60301E-02
+1.77380E-02
+8.24025E-02
+-2.67488E-02
+2.84035E-02
+2.91455E-03
+1.70812E-02
+1.54656E-02
+-2.40176E-02
+2.76831E-02
+-5.27066E-03
+2.65944E-02
+4.76662E-02
+-4.29954E-02
+-4.04722E-02
+3.02425E-02
+-8.97502E-02
+5.58394E-02
+-5.79944E-02
+-1.14607E-02
+-5.29796E-02
+2.89066E-02
+7.05200E-02
+-3.06923E-02
+1.81037E-02
+-7.19964E-03
+1.24096E-02
+PPA_abinitio.fitout
+6 6
+1.32899E-01
+1.38710E-01
+-3.14724E-02
+1.18637E-03
+7.58328E-02
+-1.42436E-02
+-6.58168E-02
+-8.66061E-02
+-6.79304E-02
+-1.19298E-01
+-1.61300E-02
+-4.37837E-02
+-1.78927E-01
+-1.06051E-01
+-1.70089E-01
+-8.92743E-03
+-1.36536E-02
+-3.11752E-02
+-6.57058E-02
+8.25739E-02
+8.75466E-02
+4.94767E-02
+3.98260E-02
+7.59713E-03
+7.87024E-02
+-1.31469E-01
+-1.19822E-01
+-1.88926E-01
+-5.37974E-02
+-5.92828E-02
+3.66331E-02
+-7.89996E-02
+3.94025E-02
+-1.22489E-01
+-1.04510E-01
+-8.92355E-02
+3.76168E-02
+-2.21729E-01
+1.45467E-01
+-2.07780E-01
+-2.21517E-02
+-6.26420E-02
+5.66877E-02
+-8.08278E-03
+2.41190E-02
+-1.27733E-01
+-1.22023E-01
+-1.08070E-01
+6.95849E-02
+1.10424E-01
+-8.50267E-02
+-1.19373E-01
+3.44974E-02
+-1.88762E-01
+1.16889E-01
+-8.09305E-02
+-2.55243E-02
+-4.70948E-02
+3.42463E-02
+-6.88564E-03
+8.19412E-02
+-5.04839E-02
+-1.18235E-01
+-1.96336E-01
+1.05609E-02
+1.98306E-02
+-3.37354E-03
+-3.73736E-02
+4.17335E-02
+4.33113E-02
+-5.87118E-02
+-6.90520E-02
+2.89217E-02
+-1.39956E-01
+8.07782E-02
+-5.61360E-02
+-1.10854E-02
+-2.35454E-02
+-2.55672E-02
+-1.54173E-02
+5.08292E-02
+2.02335E-03
+-6.94101E-02
+-7.15398E-02
+PPP_abinitio.fitout
+6 6
+9.55141E-01
+2.82825E-01
+-1.25053E-01
+-7.72476E-02
+-1.53768E-02
+-1.23068E-01
+-1.07860E+00
+-8.60254E-01
+-5.96195E-01
+-5.73846E-01
+-1.56914E-01
+-2.25099E-01
+3.39553E-02
+5.27174E-02
+-4.70984E-04
+6.95022E-02
+3.54170E-02
+-6.13267E-04
+-7.18476E-02
+1.27301E-01
+3.65459E-02
+-6.67390E-02
+-1.04061E-02
+-1.72909E-02
+9.06180E-01
+-3.16174E-01
+-4.36987E-01
+-8.48768E-01
+1.45459E-01
+-1.03920E-01
+6.02209E-02
+-2.74110E-02
+3.85731E-01
+-4.98738E-01
+-6.17808E-01
+-3.88753E-01
+1.82417E-02
+-5.85814E-02
+8.79629E-02
+-1.17158E-02
+9.11518E-02
+-5.60357E-02
+-8.69024E-02
+7.33401E-02
+1.80120E-01
+-4.92279E-01
+-5.69457E-01
+-2.06515E-01
+7.60602E-02
+1.37305E-01
+-1.78078E-01
+-3.68494E-02
+3.13809E-02
+-4.97757E-02
+6.29741E-02
+-1.27850E-02
+3.55620E-03
+6.26058E-03
+-1.59718E-01
+1.40642E-01
+2.54613E-01
+-4.16776E-01
+-4.60517E-01
+-2.65052E-01
+4.02722E-03
+-3.47336E-02
+1.38117E-02
+4.25497E-02
+-4.07847E-03
+8.37183E-02
+-1.45515E-01
+-4.97774E-02
+2.16109E-02
+-4.19394E-02
+5.29746E-02
+-7.24450E-03
+1.28857E-03
+5.46137E-03
+-2.15724E-01
+2.27414E-01
+9.50069E-02
+-1.60267E-01
+-1.77017E-01
+-1.06693E-01
+
index 8ee5798..72e2266 100644 (file)
@@ -47,7 +47,7 @@ set(UNRES_CSA_SRC0
        parmread.F
        pinorm.f
        printmat.f
-       prng_32.F
+       prng.F
        randgens.f
        ran.f
        readpdb.F
index 74b5dd0..4b6b59c 120000 (symlink)
@@ -1 +1 @@
-Makefile-DFA-NEWPARM.piasek
\ No newline at end of file
+Makefile_4P
\ No newline at end of file
index 915eec2..082a3da 100644 (file)
@@ -1,7 +1,6 @@
 CPPFLAGS = -DPROCOR -DLINUX -DPGI -DISNAN -DMP -DMPI -DUNRES \
            -DSPLITELE -DAMD64 -DLANG0 -DPROCOR \
-           -DCRYST_BOND -DCRYST_THETA -DCRYST_SC 
-# -DPROCOR
+           -DCRYST_BOND -DCRYST_THETA -DCRYST_SC -DPROCOR
 #           -DTSCSC
 #-DTIMING \
 # -DCRYST_BOND -DCRYST_THETA -DCRYST_SC 
@@ -51,8 +50,8 @@ object = unres_csa.o arcos.o cartprint.o chainbuild.o initialize_p.o \
         distfit.o banach.o TMscore_subroutine.o minim_mult.o
 
 unres: ${object} 
-#      cc -o compinfo compinfo.c
-#      ./compinfo | true
+       cc -o compinfo compinfo.c
+       ./compinfo | true
        ${FC} ${FFLAGS} cinfo.f
        ${FC} ${OPT} ${object} cinfo.o ${LIBS}  -o ${BIN}
 
diff --git a/source/unres/src_CSA/cinfo.f b/source/unres/src_CSA/cinfo.f
new file mode 100644 (file)
index 0000000..34dcb6d
--- /dev/null
@@ -0,0 +1,29 @@
+C 0 4196011 1
+      subroutine cinfo
+      include 'COMMON.IOUNITS'
+      write(iout,*)'++++ Compile info ++++'
+      write(iout,*)'Version 0.4196011 build 1'
+      write(iout,*)'compiled Mon Nov  5 08:27:57 2012'
+      write(iout,*)'compiled by bz94@matrix.chem.cornell.edu'
+      write(iout,*)'OS name:    Linux '
+      write(iout,*)'OS release: 2.6.34.9-69.fc13.x86_64 '
+      write(iout,*)'OS version:',
+     & ' #1 SMP Tue May 3 09:23:03 UTC 2011 '
+      write(iout,*)'flags:'
+      write(iout,*)'CPPFLAGS = -DPROCOR -DLINUX -DPGI -DISNAN -DMP ...'
+      write(iout,*)'INSTALL_DIR = /users/software/mpich-1.2.7p1_int...'
+      write(iout,*)'FC= ifort '
+      write(iout,*)'OPT =  -O3 -ip -w '
+      write(iout,*)'FFLAGS = -c ${OPT} -I$(INSTALL_DIR)/include'
+      write(iout,*)'FFLAGS1 = -c -w -g -O0 -d2 -CA -CB -I$(INSTALL_...'
+      write(iout,*)'FFLAGS2 = -c -w -g -O0 -I$(INSTALL_DIR)/include'
+      write(iout,*)'FFLAGS3 = -c -w -O3 -mp'
+      write(iout,*)'FFLAGSE = -c -w -O3 -ipo -ipo_obj  -opt_report ...'
+      write(iout,*)'BIN = ../../../bin/unres/CSA/unres_csa_ifort_mp...'
+      write(iout,*)'LIBS =  -lpthread -L$(INSTALL_DIR)/lib -lmpich'
+      write(iout,*)'ARCH = LINUX'
+      write(iout,*)'PP = /lib/cpp -P'
+      write(iout,*)'object = unres_csa.o arcos.o cartprint.o chainb...'
+      write(iout,*)'++++ End of compile info ++++'
+      return
+      end
diff --git a/source/unres/src_CSA/prng.F b/source/unres/src_CSA/prng.F
new file mode 100644 (file)
index 0000000..73f6766
--- /dev/null
@@ -0,0 +1,525 @@
+      real*8 function prng_next(me)
+      implicit none
+      integer me
+c
+c Calling sequence:
+c      <new random number> = prng_next ( <ordinal of generator desired> )
+c      <vector of random #s> = vprng ( <ordinal>, <vector>, <length> )
+c
+c This code is based on a sequential algorithm provided by Mal Kalos.
+c This version uses a single 64-bit word to store the initial seeds
+c and additive constants.
+c A 64-bit floating point number is returned.
+c
+c The array "iparam" is full-word aligned, being padded by zeros to
+c let each generator be on a subpage boundary.
+c That is, rows 1 and 2 in a given column of the array are for real, 
+c rows 3-16 are bogus.
+c
+c July 12, 1993: double the number of sequences.  We should have been
+c                using two packets per seed, rather than four
+c October 31, 1993: merge the two arrays of seeds and constants,
+c                and switch to 64-bit arithmetic.
+c June 1994: port to RS6K.  Internal state is kept as 2 64-bit integers
+c The ishft function is defined only on 32-bit integers, so we will
+c shift numbers by dividing by 2**11 and then adding on 2**53-1.
+c
+c November 1994: ishift now works on 64-bit numbers (though it gives a
+c warning).  Thus we go back to using it.  John Zollweg also added the
+c vprng() routine to return vectors of real*8 random numbers.
+c
+      real*8 recip53
+      parameter ( recip53 = 2.0D0**(-53) )
+      integer*8 two
+      parameter ( two = 2**11)
+      integer*8 m,ishift
+c      parameter ( m = 34522712143931 )          ! 11**13
+c      parameter ( ishift = 9007199254740991 )   ! 2**53-1
+
+      integer nmax 
+      integer*8 iparam
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      integer*8 next
+
+crc  g77 doesn't support integer*8 constants
+      m = dint(34522712143931.0d0)
+      ishift = dint(9007199254740991.0d0)
+
+c RS6K porting note: ishift now takes 64-bit integers , with a warning
+      if ( 0.le.me .and. me.le.nmax ) then
+         next = iparam(1,me)*m + iparam(2,me)
+         iparam(1,me) = next
+         prng_next = recip53 * ishft( next, -11 )
+      else
+         prng_next=-1.0D0
+      endif
+
+      end
+c
+c   vprng(me, rn, num)       Get a vector of random numbers
+c
+      subroutine vprng(me,rn,num)
+      real*8 recip53, rn(1)
+      parameter ( recip53 = 2.0D0**(-53) )
+      integer*8 m,iparam
+c      parameter ( m = 34522712143931 )          ! 11**13
+      integer nmax, num, me
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      integer*8 next
+
+crc  g77 doesn't support integer*8 constants
+      m = dint(34522712143931.0d0)
+
+      if ( 0.le.me .and. me.le.nmax ) then
+         do 1 i=1,num
+            next = iparam(1,me)*m + iparam(2,me)
+            iparam(1,me) = next
+            rn(i) = recip53 * ishft( next, -11 )
+    1    continue
+      else
+         rn(1)=-1.0D0
+      endif
+      return
+      end
+
+c
+c   prng_chkpnt          Get the current state of a generator
+c
+c Calling sequence:
+c   logical prng_chkpnt, status
+c   status = prng_chkpnt (me, iseed)    where
+c
+c     me is the particular generator whose state is being gotten
+c     seed is an 4-element integer array where the "l"-values will be saved
+c
+      logical function prng_chkpnt (me, iseed)
+      implicit none
+      integer me 
+      integer*8 iseed
+
+      integer nmax 
+      integer*8 iparam
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      if (me .lt. 0 .or. me .gt. nmax) then
+         prng_chkpnt=.false.
+      else
+        prng_chkpnt=.true.
+        iseed=iparam(1,me)
+      endif
+      end
+c
+c   prng_restart          Restart generator from a saved state
+c
+c Calling sequence:
+c   logical prng_restart, status
+c   status = prng_restart (me, iseed)    where
+c
+c     me is the particular generator being restarted
+c     iseed is a 8-byte integer containing the "l"-values
+c
+      logical function prng_restart (me, iseed)
+      implicit none
+      integer me
+      integer*8 iseed
+
+      integer nmax 
+      integer*8 iparam
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      if (me .lt. 0 .or. me .gt. nmax) then
+         prng_restart=.false.
+        return
+      else
+        prng_restart=.true.
+        iparam(1,me)=iseed
+      endif
+      end
+
+      block data prngblk
+      parameter(nmax=1021)
+      integer*8 iparam
+      common/ksrprng/iparam(2,0:nmax)
+      data (iparam(1,i),iparam(2,i),i=   0,  29) /
+     + 11848219, 11848219, 11848237, 11848237, 11848241, 11848241,
+     + 11848247, 11848247, 11848253, 11848253, 11848271, 11848271,
+     + 11848297, 11848297, 11848313, 11848313, 11848339, 11848339,
+     + 11848351, 11848351, 11848357, 11848357, 11848363, 11848363,
+     + 11848367, 11848367, 11848373, 11848373, 11848379, 11848379,
+     + 11848393, 11848393, 11848433, 11848433, 11848451, 11848451,
+     + 11848469, 11848469, 11848477, 11848477, 11848489, 11848489,
+     + 11848493, 11848493, 11848513, 11848513, 11848523, 11848523,
+     + 11848531, 11848531, 11848537, 11848537, 11848553, 11848553,
+     + 11848589, 11848589, 11848591, 11848591, 11848601, 11848601 /
+      data (iparam(1,i),iparam(2,i),i=  30,  59) /
+     + 11848619, 11848619, 11848637, 11848637, 11848663, 11848663,
+     + 11848673, 11848673, 11848679, 11848679, 11848691, 11848691,
+     + 11848699, 11848699, 11848709, 11848709, 11848717, 11848717,
+     + 11848721, 11848721, 11848729, 11848729, 11848741, 11848741,
+     + 11848751, 11848751, 11848757, 11848757, 11848787, 11848787,
+     + 11848801, 11848801, 11848829, 11848829, 11848853, 11848853,
+     + 11848861, 11848861, 11848867, 11848867, 11848873, 11848873,
+     + 11848891, 11848891, 11848909, 11848909, 11848919, 11848919,
+     + 11848931, 11848931, 11848937, 11848937, 11848961, 11848961,
+     + 11848981, 11848981, 11849021, 11849021, 11849039, 11849039 /
+      data (iparam(1,i),iparam(2,i),i=  60,  89) /
+     + 11849053, 11849053, 11849059, 11849059, 11849069, 11849069,
+     + 11849077, 11849077, 11849087, 11849087, 11849093, 11849093,
+     + 11849107, 11849107, 11849111, 11849111, 11849129, 11849129,
+     + 11849137, 11849137, 11849177, 11849177, 11849183, 11849183,
+     + 11849203, 11849203, 11849231, 11849231, 11849237, 11849237,
+     + 11849239, 11849239, 11849249, 11849249, 11849251, 11849251,
+     + 11849269, 11849269, 11849273, 11849273, 11849291, 11849291,
+     + 11849297, 11849297, 11849309, 11849309, 11849339, 11849339,
+     + 11849359, 11849359, 11849363, 11849363, 11849399, 11849399,
+     + 11849401, 11849401, 11849413, 11849413, 11849417, 11849417 /
+      data (iparam(1,i),iparam(2,i),i=  90, 119) /
+     + 11849437, 11849437, 11849443, 11849443, 11849473, 11849473,
+     + 11849491, 11849491, 11849503, 11849503, 11849507, 11849507,
+     + 11849557, 11849557, 11849567, 11849567, 11849569, 11849569,
+     + 11849573, 11849573, 11849587, 11849587, 11849599, 11849599,
+     + 11849633, 11849633, 11849641, 11849641, 11849653, 11849653,
+     + 11849659, 11849659, 11849671, 11849671, 11849683, 11849683,
+     + 11849689, 11849689, 11849693, 11849693, 11849699, 11849699,
+     + 11849701, 11849701, 11849707, 11849707, 11849713, 11849713,
+     + 11849723, 11849723, 11849741, 11849741, 11849743, 11849743,
+     + 11849759, 11849759, 11849767, 11849767, 11849771, 11849771 /
+      data (iparam(1,i),iparam(2,i),i= 120, 149) /
+     + 11849791, 11849791, 11849801, 11849801, 11849809, 11849809,
+     + 11849813, 11849813, 11849869, 11849869, 11849881, 11849881,
+     + 11849891, 11849891, 11849909, 11849909, 11849923, 11849923,
+     + 11849933, 11849933, 11849947, 11849947, 11849987, 11849987,
+     + 11850001, 11850001, 11850011, 11850011, 11850019, 11850019,
+     + 11850023, 11850023, 11850031, 11850031, 11850049, 11850049,
+     + 11850061, 11850061, 11850073, 11850073, 11850077, 11850077,
+     + 11850103, 11850103, 11850109, 11850109, 11850121, 11850121,
+     + 11850127, 11850127, 11850133, 11850133, 11850149, 11850149,
+     + 11850161, 11850161, 11850169, 11850169, 11850191, 11850191 /
+      data (iparam(1,i),iparam(2,i),i= 150, 179) /
+     + 11850233, 11850233, 11850247, 11850247, 11850259, 11850259,
+     + 11850269, 11850269, 11850283, 11850283, 11850301, 11850301,
+     + 11850341, 11850341, 11850347, 11850347, 11850367, 11850367,
+     + 11850373, 11850373, 11850379, 11850379, 11850389, 11850389,
+     + 11850407, 11850407, 11850427, 11850427, 11850437, 11850437,
+     + 11850469, 11850469, 11850481, 11850481, 11850511, 11850511,
+     + 11850529, 11850529, 11850541, 11850541, 11850557, 11850557,
+     + 11850607, 11850607, 11850611, 11850611, 11850667, 11850667,
+     + 11850677, 11850677, 11850679, 11850679, 11850701, 11850701,
+     + 11850731, 11850731, 11850739, 11850739, 11850749, 11850749 /
+      data (iparam(1,i),iparam(2,i),i= 180, 209) /
+     + 11850791, 11850791, 11850803, 11850803, 11850829, 11850829,
+     + 11850833, 11850833, 11850859, 11850859, 11850877, 11850877,
+     + 11850899, 11850899, 11850907, 11850907, 11850913, 11850913,
+     + 11850919, 11850919, 11850931, 11850931, 11850941, 11850941,
+     + 11850947, 11850947, 11850953, 11850953, 11850961, 11850961,
+     + 11850983, 11850983, 11850991, 11850991, 11850997, 11850997,
+     + 11851031, 11851031, 11851033, 11851033, 11851051, 11851051,
+     + 11851061, 11851061, 11851067, 11851067, 11851093, 11851093,
+     + 11851109, 11851109, 11851123, 11851123, 11851127, 11851127,
+     + 11851139, 11851139, 11851157, 11851157, 11851163, 11851163 /
+      data (iparam(1,i),iparam(2,i),i= 210, 239) /
+     + 11851181, 11851181, 11851201, 11851201, 11851219, 11851219,
+     + 11851291, 11851291, 11851303, 11851303, 11851309, 11851309,
+     + 11851313, 11851313, 11851319, 11851319, 11851349, 11851349,
+     + 11851351, 11851351, 11851361, 11851361, 11851373, 11851373,
+     + 11851403, 11851403, 11851409, 11851409, 11851423, 11851423,
+     + 11851447, 11851447, 11851451, 11851451, 11851481, 11851481,
+     + 11851493, 11851493, 11851519, 11851519, 11851523, 11851523,
+     + 11851529, 11851529, 11851547, 11851547, 11851549, 11851549,
+     + 11851559, 11851559, 11851577, 11851577, 11851589, 11851589,
+     + 11851591, 11851591, 11851597, 11851597, 11851603, 11851603 /
+      data (iparam(1,i),iparam(2,i),i= 240, 269) /
+     + 11851607, 11851607, 11851613, 11851613, 11851621, 11851621,
+     + 11851627, 11851627, 11851639, 11851639, 11851673, 11851673,
+     + 11851681, 11851681, 11851727, 11851727, 11851753, 11851753,
+     + 11851759, 11851759, 11851787, 11851787, 11851793, 11851793,
+     + 11851799, 11851799, 11851813, 11851813, 11851841, 11851841,
+     + 11851859, 11851859, 11851867, 11851867, 11851891, 11851891,
+     + 11851909, 11851909, 11851919, 11851919, 11851927, 11851927,
+     + 11851933, 11851933, 11851949, 11851949, 11851967, 11851967,
+     + 11851997, 11851997, 11852017, 11852017, 11852051, 11852051,
+     + 11852053, 11852053, 11852059, 11852059, 11852083, 11852083 /
+      data (iparam(1,i),iparam(2,i),i= 270, 299) /
+     + 11852089, 11852089, 11852129, 11852129, 11852147, 11852147,
+     + 11852149, 11852149, 11852161, 11852161, 11852171, 11852171,
+     + 11852177, 11852177, 11852209, 11852209, 11852221, 11852221,
+     + 11852237, 11852237, 11852251, 11852251, 11852263, 11852263,
+     + 11852273, 11852273, 11852279, 11852279, 11852287, 11852287,
+     + 11852293, 11852293, 11852297, 11852297, 11852303, 11852303,
+     + 11852311, 11852311, 11852327, 11852327, 11852339, 11852339,
+     + 11852341, 11852341, 11852359, 11852359, 11852369, 11852369,
+     + 11852437, 11852437, 11852453, 11852453, 11852459, 11852459,
+     + 11852473, 11852473, 11852513, 11852513, 11852531, 11852531 /
+      data (iparam(1,i),iparam(2,i),i= 300, 329) /
+     + 11852537, 11852537, 11852539, 11852539, 11852557, 11852557,
+     + 11852573, 11852573, 11852579, 11852579, 11852591, 11852591,
+     + 11852609, 11852609, 11852611, 11852611, 11852623, 11852623,
+     + 11852641, 11852641, 11852647, 11852647, 11852657, 11852657,
+     + 11852663, 11852663, 11852717, 11852717, 11852719, 11852719,
+     + 11852741, 11852741, 11852759, 11852759, 11852767, 11852767,
+     + 11852773, 11852773, 11852803, 11852803, 11852807, 11852807,
+     + 11852809, 11852809, 11852831, 11852831, 11852833, 11852833,
+     + 11852837, 11852837, 11852857, 11852857, 11852873, 11852873,
+     + 11852879, 11852879, 11852891, 11852891, 11852917, 11852917 /
+      data (iparam(1,i),iparam(2,i),i= 330, 359) /
+     + 11852921, 11852921, 11852957, 11852957, 11852959, 11852959,
+     + 11852969, 11852969, 11852983, 11852983, 11852989, 11852989,
+     + 11853001, 11853001, 11853013, 11853013, 11853019, 11853019,
+     + 11853031, 11853031, 11853089, 11853089, 11853133, 11853133,
+     + 11853157, 11853157, 11853161, 11853161, 11853181, 11853181,
+     + 11853203, 11853203, 11853217, 11853217, 11853221, 11853221,
+     + 11853227, 11853227, 11853241, 11853241, 11853307, 11853307,
+     + 11853319, 11853319, 11853323, 11853323, 11853329, 11853329,
+     + 11853367, 11853367, 11853383, 11853383, 11853419, 11853419,
+     + 11853421, 11853421, 11853427, 11853427, 11853449, 11853449 /
+      data (iparam(1,i),iparam(2,i),i= 360, 389) /
+     + 11853451, 11853451, 11853463, 11853463, 11853529, 11853529,
+     + 11853557, 11853557, 11853571, 11853571, 11853601, 11853601,
+     + 11853613, 11853613, 11853617, 11853617, 11853629, 11853629,
+     + 11853649, 11853649, 11853659, 11853659, 11853679, 11853679,
+     + 11853689, 11853689, 11853719, 11853719, 11853731, 11853731,
+     + 11853757, 11853757, 11853761, 11853761, 11853773, 11853773,
+     + 11853791, 11853791, 11853817, 11853817, 11853839, 11853839,
+     + 11853847, 11853847, 11853857, 11853857, 11853869, 11853869,
+     + 11853883, 11853883, 11853887, 11853887, 11853889, 11853889,
+     + 11853893, 11853893, 11853899, 11853899, 11853911, 11853911 /
+      data (iparam(1,i),iparam(2,i),i= 390, 419) /
+     + 11853931, 11853931, 11853943, 11853943, 11853979, 11853979,
+     + 11853991, 11853991, 11854001, 11854001, 11854009, 11854009,
+     + 11854019, 11854019, 11854057, 11854057, 11854061, 11854061,
+     + 11854147, 11854147, 11854159, 11854159, 11854163, 11854163,
+     + 11854169, 11854169, 11854211, 11854211, 11854247, 11854247,
+     + 11854261, 11854261, 11854267, 11854267, 11854279, 11854279,
+     + 11854303, 11854303, 11854327, 11854327, 11854331, 11854331,
+     + 11854333, 11854333, 11854363, 11854363, 11854379, 11854379,
+     + 11854399, 11854399, 11854411, 11854411, 11854429, 11854429,
+     + 11854433, 11854433, 11854439, 11854439, 11854441, 11854441 /
+      data (iparam(1,i),iparam(2,i),i= 420, 449) /
+     + 11854463, 11854463, 11854477, 11854477, 11854489, 11854489,
+     + 11854517, 11854517, 11854519, 11854519, 11854523, 11854523,
+     + 11854529, 11854529, 11854567, 11854567, 11854571, 11854571,
+     + 11854573, 11854573, 11854603, 11854603, 11854607, 11854607,
+     + 11854681, 11854681, 11854691, 11854691, 11854709, 11854709,
+     + 11854723, 11854723, 11854757, 11854757, 11854783, 11854783,
+     + 11854793, 11854793, 11854813, 11854813, 11854847, 11854847,
+     + 11854853, 11854853, 11854873, 11854873, 11854877, 11854877,
+     + 11854883, 11854883, 11854891, 11854891, 11854897, 11854897,
+     + 11854901, 11854901, 11854919, 11854919, 11854937, 11854937 /
+      data (iparam(1,i),iparam(2,i),i= 450, 479) /
+     + 11854961, 11854961, 11854963, 11854963, 11854979, 11854979,
+     + 11855003, 11855003, 11855017, 11855017, 11855023, 11855023,
+     + 11855029, 11855029, 11855033, 11855033, 11855111, 11855111,
+     + 11855141, 11855141, 11855147, 11855147, 11855149, 11855149,
+     + 11855159, 11855159, 11855177, 11855177, 11855203, 11855203,
+     + 11855213, 11855213, 11855219, 11855219, 11855231, 11855231,
+     + 11855267, 11855267, 11855269, 11855269, 11855303, 11855303,
+     + 11855309, 11855309, 11855321, 11855321, 11855329, 11855329,
+     + 11855339, 11855339, 11855351, 11855351, 11855353, 11855353,
+     + 11855357, 11855357, 11855359, 11855359, 11855381, 11855381 /
+      data (iparam(1,i),iparam(2,i),i= 480, 509) /
+     + 11855383, 11855383, 11855387, 11855387, 11855399, 11855399,
+     + 11855407, 11855407, 11855413, 11855413, 11855489, 11855489,
+     + 11855491, 11855491, 11855507, 11855507, 11855521, 11855521,
+     + 11855531, 11855531, 11855549, 11855549, 11855551, 11855551,
+     + 11855567, 11855567, 11855581, 11855581, 11855587, 11855587,
+     + 11855593, 11855593, 11855633, 11855633, 11855653, 11855653,
+     + 11855663, 11855663, 11855687, 11855687, 11855689, 11855689,
+     + 11855699, 11855699, 11855713, 11855713, 11855731, 11855731,
+     + 11855737, 11855737, 11855743, 11855743, 11855747, 11855747,
+     + 11855759, 11855759, 11855773, 11855773, 11855801, 11855801 /
+      data (iparam(1,i),iparam(2,i),i= 510, 539) /
+     + 11855807, 11855807, 11855813, 11855813, 11855827, 11855827,
+     + 11855839, 11855839, 11855869, 11855869, 11855881, 11855881,
+     + 11855903, 11855903, 11855911, 11855911, 11855933, 11855933,
+     + 11855959, 11855959, 11855989, 11855989, 11855993, 11855993,
+     + 11855999, 11855999, 11856001, 11856001, 11856023, 11856023,
+     + 11856049, 11856049, 11856071, 11856071, 11856101, 11856101,
+     + 11856107, 11856107, 11856113, 11856113, 11856139, 11856139,
+     + 11856151, 11856151, 11856161, 11856161, 11856179, 11856179,
+     + 11856193, 11856193, 11856199, 11856199, 11856223, 11856223,
+     + 11856239, 11856239, 11856263, 11856263, 11856269, 11856269 /
+      data (iparam(1,i),iparam(2,i),i= 540, 569) /
+     + 11856281, 11856281, 11856287, 11856287, 11856307, 11856307,
+     + 11856311, 11856311, 11856329, 11856329, 11856343, 11856343,
+     + 11856359, 11856359, 11856371, 11856371, 11856373, 11856373,
+     + 11856409, 11856409, 11856419, 11856419, 11856461, 11856461,
+     + 11856469, 11856469, 11856473, 11856473, 11856479, 11856479,
+     + 11856511, 11856511, 11856517, 11856517, 11856541, 11856541,
+     + 11856547, 11856547, 11856553, 11856553, 11856583, 11856583,
+     + 11856629, 11856629, 11856641, 11856641, 11856653, 11856653,
+     + 11856659, 11856659, 11856673, 11856673, 11856697, 11856697,
+     + 11856709, 11856709, 11856727, 11856727, 11856731, 11856731 /
+      data (iparam(1,i),iparam(2,i),i= 570, 599) /
+     + 11856763, 11856763, 11856809, 11856809, 11856811, 11856811,
+     + 11856821, 11856821, 11856841, 11856841, 11856857, 11856857,
+     + 11856877, 11856877, 11856883, 11856883, 11856899, 11856899,
+     + 11856919, 11856919, 11856947, 11856947, 11856953, 11856953,
+     + 11856979, 11856979, 11857003, 11857003, 11857033, 11857033,
+     + 11857037, 11857037, 11857039, 11857039, 11857049, 11857049,
+     + 11857061, 11857061, 11857067, 11857067, 11857073, 11857073,
+     + 11857081, 11857081, 11857091, 11857091, 11857093, 11857093,
+     + 11857099, 11857099, 11857123, 11857123, 11857127, 11857127,
+     + 11857147, 11857147, 11857151, 11857151, 11857193, 11857193 /
+      data (iparam(1,i),iparam(2,i),i= 600, 629) /
+     + 11857217, 11857217, 11857229, 11857229, 11857243, 11857243,
+     + 11857249, 11857249, 11857267, 11857267, 11857277, 11857277,
+     + 11857291, 11857291, 11857303, 11857303, 11857309, 11857309,
+     + 11857327, 11857327, 11857331, 11857331, 11857333, 11857333,
+     + 11857361, 11857361, 11857367, 11857367, 11857369, 11857369,
+     + 11857393, 11857393, 11857399, 11857399, 11857409, 11857409,
+     + 11857421, 11857421, 11857423, 11857423, 11857451, 11857451,
+     + 11857453, 11857453, 11857457, 11857457, 11857477, 11857477,
+     + 11857481, 11857481, 11857493, 11857493, 11857499, 11857499,
+     + 11857519, 11857519, 11857523, 11857523, 11857529, 11857529 /
+      data (iparam(1,i),iparam(2,i),i= 630, 659) /
+     + 11857543, 11857543, 11857561, 11857561, 11857589, 11857589,
+     + 11857591, 11857591, 11857613, 11857613, 11857621, 11857621,
+     + 11857661, 11857661, 11857667, 11857667, 11857693, 11857693,
+     + 11857697, 11857697, 11857709, 11857709, 11857711, 11857711,
+     + 11857751, 11857751, 11857753, 11857753, 11857759, 11857759,
+     + 11857763, 11857763, 11857777, 11857777, 11857787, 11857787,
+     + 11857793, 11857793, 11857801, 11857801, 11857817, 11857817,
+     + 11857819, 11857819, 11857831, 11857831, 11857837, 11857837,
+     + 11857873, 11857873, 11857877, 11857877, 11857883, 11857883,
+     + 11857889, 11857889, 11857907, 11857907, 11857913, 11857913 /
+      data (iparam(1,i),iparam(2,i),i= 660, 689) /
+     + 11857931, 11857931, 11857969, 11857969, 11857991, 11857991,
+     + 11857999, 11857999, 11858009, 11858009, 11858017, 11858017,
+     + 11858023, 11858023, 11858029, 11858029, 11858039, 11858039,
+     + 11858051, 11858051, 11858057, 11858057, 11858059, 11858059,
+     + 11858101, 11858101, 11858111, 11858111, 11858131, 11858131,
+     + 11858149, 11858149, 11858159, 11858159, 11858177, 11858177,
+     + 11858191, 11858191, 11858201, 11858201, 11858227, 11858227,
+     + 11858243, 11858243, 11858267, 11858267, 11858269, 11858269,
+     + 11858279, 11858279, 11858281, 11858281, 11858291, 11858291,
+     + 11858311, 11858311, 11858323, 11858323, 11858359, 11858359 /
+      data (iparam(1,i),iparam(2,i),i= 690, 719) /
+     + 11858377, 11858377, 11858381, 11858381, 11858387, 11858387,
+     + 11858423, 11858423, 11858443, 11858443, 11858447, 11858447,
+     + 11858479, 11858479, 11858533, 11858533, 11858543, 11858543,
+     + 11858551, 11858551, 11858557, 11858557, 11858569, 11858569,
+     + 11858573, 11858573, 11858579, 11858579, 11858597, 11858597,
+     + 11858599, 11858599, 11858629, 11858629, 11858657, 11858657,
+     + 11858659, 11858659, 11858683, 11858683, 11858701, 11858701,
+     + 11858719, 11858719, 11858723, 11858723, 11858729, 11858729,
+     + 11858747, 11858747, 11858779, 11858779, 11858783, 11858783,
+     + 11858801, 11858801, 11858807, 11858807, 11858813, 11858813 /
+      data (iparam(1,i),iparam(2,i),i= 720, 749) /
+     + 11858839, 11858839, 11858851, 11858851, 11858893, 11858893,
+     + 11858897, 11858897, 11858921, 11858921, 11858947, 11858947,
+     + 11858953, 11858953, 11858969, 11858969, 11858971, 11858971,
+     + 11858989, 11858989, 11859017, 11859017, 11859031, 11859031,
+     + 11859049, 11859049, 11859061, 11859061, 11859073, 11859073,
+     + 11859077, 11859077, 11859079, 11859079, 11859083, 11859083,
+     + 11859101, 11859101, 11859109, 11859109, 11859137, 11859137,
+     + 11859139, 11859139, 11859151, 11859151, 11859157, 11859157,
+     + 11859163, 11859163, 11859167, 11859167, 11859179, 11859179,
+     + 11859187, 11859187, 11859229, 11859229, 11859233, 11859233 /
+      data (iparam(1,i),iparam(2,i),i= 750, 779) /
+     + 11859241, 11859241, 11859247, 11859247, 11859269, 11859269,
+     + 11859293, 11859293, 11859307, 11859307, 11859311, 11859311,
+     + 11859349, 11859349, 11859359, 11859359, 11859371, 11859371,
+     + 11859377, 11859377, 11859383, 11859383, 11859427, 11859427,
+     + 11859433, 11859433, 11859451, 11859451, 11859457, 11859457,
+     + 11859461, 11859461, 11859473, 11859473, 11859481, 11859481,
+     + 11859487, 11859487, 11859493, 11859493, 11859503, 11859503,
+     + 11859509, 11859509, 11859539, 11859539, 11859541, 11859541,
+     + 11859563, 11859563, 11859569, 11859569, 11859571, 11859571,
+     + 11859583, 11859583, 11859599, 11859599, 11859611, 11859611 /
+      data (iparam(1,i),iparam(2,i),i= 780, 809) /
+     + 11859643, 11859643, 11859707, 11859707, 11859713, 11859713,
+     + 11859719, 11859719, 11859739, 11859739, 11859751, 11859751,
+     + 11859791, 11859791, 11859817, 11859817, 11859821, 11859821,
+     + 11859833, 11859833, 11859847, 11859847, 11859853, 11859853,
+     + 11859877, 11859877, 11859889, 11859889, 11859893, 11859893,
+     + 11859901, 11859901, 11859907, 11859907, 11859917, 11859917,
+     + 11859923, 11859923, 11859929, 11859929, 11859961, 11859961,
+     + 11859979, 11859979, 11859989, 11859989, 11859997, 11859997,
+     + 11860021, 11860021, 11860031, 11860031, 11860039, 11860039,
+     + 11860049, 11860049, 11860081, 11860081, 11860087, 11860087 /
+      data (iparam(1,i),iparam(2,i),i= 810, 839) /
+     + 11860097, 11860097, 11860103, 11860103, 11860109, 11860109,
+     + 11860117, 11860117, 11860133, 11860133, 11860151, 11860151,
+     + 11860171, 11860171, 11860207, 11860207, 11860223, 11860223,
+     + 11860231, 11860231, 11860243, 11860243, 11860267, 11860267,
+     + 11860301, 11860301, 11860307, 11860307, 11860327, 11860327,
+     + 11860379, 11860379, 11860397, 11860397, 11860411, 11860411,
+     + 11860469, 11860469, 11860477, 11860477, 11860483, 11860483,
+     + 11860487, 11860487, 11860489, 11860489, 11860493, 11860493,
+     + 11860517, 11860517, 11860547, 11860547, 11860567, 11860567,
+     + 11860573, 11860573, 11860613, 11860613, 11860619, 11860619 /
+      data (iparam(1,i),iparam(2,i),i= 840, 869) /
+     + 11860627, 11860627, 11860637, 11860637, 11860643, 11860643,
+     + 11860649, 11860649, 11860661, 11860661, 11860669, 11860669,
+     + 11860687, 11860687, 11860691, 11860691, 11860697, 11860697,
+     + 11860699, 11860699, 11860703, 11860703, 11860727, 11860727,
+     + 11860741, 11860741, 11860753, 11860753, 11860777, 11860777,
+     + 11860787, 11860787, 11860789, 11860789, 11860811, 11860811,
+     + 11860837, 11860837, 11860859, 11860859, 11860867, 11860867,
+     + 11860889, 11860889, 11860897, 11860897, 11860963, 11860963,
+     + 11860969, 11860969, 11860973, 11860973, 11860993, 11860993,
+     + 11861011, 11861011, 11861033, 11861033, 11861071, 11861071 /
+      data (iparam(1,i),iparam(2,i),i= 870, 899) /
+     + 11861081, 11861081, 11861089, 11861089, 11861093, 11861093,
+     + 11861099, 11861099, 11861107, 11861107, 11861131, 11861131,
+     + 11861141, 11861141, 11861159, 11861159, 11861167, 11861167,
+     + 11861191, 11861191, 11861197, 11861197, 11861207, 11861207,
+     + 11861219, 11861219, 11861221, 11861221, 11861231, 11861231,
+     + 11861237, 11861237, 11861273, 11861273, 11861293, 11861293,
+     + 11861299, 11861299, 11861303, 11861303, 11861327, 11861327,
+     + 11861351, 11861351, 11861357, 11861357, 11861363, 11861363,
+     + 11861371, 11861371, 11861401, 11861401, 11861407, 11861407,
+     + 11861411, 11861411, 11861413, 11861413, 11861429, 11861429 /
+      data (iparam(1,i),iparam(2,i),i= 900, 929) /
+     + 11861441, 11861441, 11861467, 11861467, 11861527, 11861527,
+     + 11861539, 11861539, 11861543, 11861543, 11861557, 11861557,
+     + 11861569, 11861569, 11861573, 11861573, 11861579, 11861579,
+     + 11861581, 11861581, 11861599, 11861599, 11861611, 11861611,
+     + 11861617, 11861617, 11861627, 11861627, 11861639, 11861639,
+     + 11861651, 11861651, 11861659, 11861659, 11861671, 11861671,
+     + 11861683, 11861683, 11861687, 11861687, 11861693, 11861693,
+     + 11861701, 11861701, 11861711, 11861711, 11861713, 11861713,
+     + 11861749, 11861749, 11861791, 11861791, 11861803, 11861803,
+     + 11861819, 11861819, 11861827, 11861827, 11861849, 11861849 /
+      data (iparam(1,i),iparam(2,i),i= 930, 959) /
+     + 11861873, 11861873, 11861879, 11861879, 11861887, 11861887,
+     + 11861911, 11861911, 11861917, 11861917, 11861921, 11861921,
+     + 11861923, 11861923, 11861953, 11861953, 11861959, 11861959,
+     + 11861987, 11861987, 11862007, 11862007, 11862013, 11862013,
+     + 11862029, 11862029, 11862031, 11862031, 11862049, 11862049,
+     + 11862077, 11862077, 11862083, 11862083, 11862157, 11862157,
+     + 11862167, 11862167, 11862199, 11862199, 11862203, 11862203,
+     + 11862217, 11862217, 11862223, 11862223, 11862229, 11862229,
+     + 11862233, 11862233, 11862239, 11862239, 11862241, 11862241,
+     + 11862259, 11862259, 11862269, 11862269, 11862271, 11862271 /
+      data (iparam(1,i),iparam(2,i),i= 960, 989) /
+     + 11862293, 11862293, 11862307, 11862307, 11862313, 11862313,
+     + 11862317, 11862317, 11862343, 11862343, 11862353, 11862353,
+     + 11862373, 11862373, 11862391, 11862391, 11862439, 11862439,
+     + 11862469, 11862469, 11862493, 11862493, 11862527, 11862527,
+     + 11862547, 11862547, 11862563, 11862563, 11862569, 11862569,
+     + 11862577, 11862577, 11862581, 11862581, 11862611, 11862611,
+     + 11862623, 11862623, 11862661, 11862661, 11862673, 11862673,
+     + 11862679, 11862679, 11862701, 11862701, 11862703, 11862703,
+     + 11862713, 11862713, 11862761, 11862761, 11862791, 11862791,
+     + 11862803, 11862803, 11862839, 11862839, 11862841, 11862841 /
+      data (iparam(1,i),iparam(2,i),i= 990,1019) /
+     + 11862857, 11862857, 11862869, 11862869, 11862881, 11862881,
+     + 11862911, 11862911, 11862919, 11862919, 11862959, 11862959,
+     + 11862979, 11862979, 11862989, 11862989, 11862997, 11862997,
+     + 11863021, 11863021, 11863031, 11863031, 11863037, 11863037,
+     + 11863039, 11863039, 11863057, 11863057, 11863067, 11863067,
+     + 11863073, 11863073, 11863099, 11863099, 11863109, 11863109,
+     + 11863121, 11863121, 11863123, 11863123, 11863133, 11863133,
+     + 11863151, 11863151, 11863153, 11863153, 11863171, 11863171,
+     + 11863183, 11863183, 11863207, 11863207, 11863213, 11863213,
+     + 11863237, 11863237, 11863249, 11863249, 11863253, 11863253 /
+      data (iparam(1,i),iparam(2,i),i=1020,1021) /
+     + 11863259, 11863259, 11863279, 11863279 /
+      end
index 3a628ac..ac396e2 100644 (file)
@@ -47,7 +47,7 @@ set(UNRES_CSA_DiL_SRC0
        parmread.F
        pinorm.f
        printmat.f
-       prng_32.F
+       prng.F
        randgens.f
        ran.f
        readpdb.F
@@ -86,7 +86,7 @@ set(UNRES_CSA_DiL_PP_SRC
        MP.F
        newconf.F
        parmread.F
-       prng_32.F
+       prng.F
        readpdb.F
        readrtns_csa.F
        rmsd.F
index 78d4cca..f212780 100644 (file)
@@ -12,14 +12,18 @@ C Generate random conformation or chain cut and regrowth.
       include 'COMMON.CONTROL'
       logical overlap,back,fail
 cd    print *,' CG Processor',me,' maxgen=',maxgen
+c      write(iout,*) "czy kurwa wogole wchodze"
       maxsi=100
 cd    write (iout,*) 'Gen_Rand_conf: nstart=',nstart
       if (nstart.lt.5) then
         it1=iabs(itype(2))
         phi(4)=gen_phi(4,iabs(itype(2)),abs(itype(3)))
-c       write(iout,*)'phi(4)=',rad2deg*phi(4)
-        if (nstart.lt.3) theta(3)=gen_theta(itype(2),pi,phi(4))
-c       write(iout,*)'theta(3)=',rad2deg*theta(3) 
+c        write(iout,*)'phi(4)=',rad2deg*phi(4)
+        ichir1=isign(1,itype(1))
+        ichir2=isign(1,itype(3))
+        if (nstart.lt.3) theta(3)=gen_theta(itype(2),ichir1,ichir2,
+     &                                     pi,phi(4))
+       write(iout,*)'theta(3)=',rad2deg*theta(3) 
         if (it1.ne.10) then
           nsi=0
           fail=.true.
@@ -38,10 +42,10 @@ c       write(iout,*)'theta(3)=',rad2deg*theta(3)
       endif
 
       maxnit=0
-
+      iprint=10 
       nit=0
       niter=0
-      back=.false.
+      back=.true.
       do while (i.le.nres .and. niter.lt.maxgen)
         if (i.lt.nstart) then
           if(iprint.gt.1) then
@@ -54,16 +58,20 @@ c       write(iout,*)'theta(3)=',rad2deg*theta(3)
           endif
           return1
         endif
-       it1=abs(itype(i-1))
-       it2=abs(itype(i-2))
-       it=abs(itype(i))
+       it1=itype(i-1)
+       it2=itype(i-2)
+       it=itype(i)
+        ichir3=isign(1,itype(i))
+        ichir2=isign(1,itype(i-1))
+        ichir0=isign(1,itype(i-3))
+        ichir1=isign(1,itype(i-2))
 c       print *,'Gen_Rand_Conf: i=',i,' it=',it,' it1=',it1,' it2=',it2,
 c     &    ' nit=',nit,' niter=',niter,' maxgen=',maxgen
        phi(i+1)=gen_phi(i+1,it1,it)
        if (back) then
           phi(i)=gen_phi(i+1,it2,it1)
           print *,'phi(',i,')=',phi(i)
-         theta(i-1)=gen_theta(it2,phi(i-1),phi(i))
+         theta(i-1)=gen_theta(it2,ichir0,ichir2,phi(i-1),phi(i))
          if (it2.ne.10) then
             nsi=0
             fail=.true.
@@ -75,7 +83,7 @@ c     &    ' nit=',nit,' niter=',niter,' maxgen=',maxgen
           endif
          call locate_next_res(i-1)
        endif
-       theta(i)=gen_theta(it1,phi(i),phi(i+1))
+       theta(i)=gen_theta(it1,ichir1,ichir3,phi(i),phi(i+1))
         if (it1.ne.10) then 
         nsi=0
         fail=.true.
@@ -212,14 +220,14 @@ C 8/13/98 Generate phi using pre-defined boundaries
       return
       end
 c---------------------------------------------------------------------------
-      double precision function gen_theta(it,gama,gama1)
+      double precision function gen_theta(it,ichir1,ichir2,gama,gama1)
       implicit real*8 (a-h,o-z)
       include 'DIMENSIONS'
       include 'COMMON.LOCAL'
       include 'COMMON.GEO'
       double precision y(2),z(2)
       double precision theta_max,theta_min
-c     print *,'gen_theta: it=',it
+      print *,'gen_theta: it=',it
       theta_min=0.05D0*pi
       theta_max=0.95D0*pi
       if (dabs(gama).gt.dwapi) then
@@ -238,8 +246,8 @@ c     print *,'gen_theta: it=',it
       endif  
       thet_pred_mean=a0thet(it)
       do k=1,2
-        thet_pred_mean=thet_pred_mean+athet(k,it,1,1)*y(k)
-     &                +bthet(k,it,1,1)*z(k)
+        thet_pred_mean=thet_pred_mean+athet(k,it,ichir1,ichir2)
+     &               *y(k)+bthet(k,it,ichir1,ichir2)*z(k)
       enddo
       sig=polthet(3,it)
       do j=2,0,-1
@@ -247,14 +255,16 @@ c     print *,'gen_theta: it=',it
       enddo
       sig=0.5D0/(sig*sig+sigc0(it))
       ak=dexp(gthet(1,it)-
-     &0.5D0*((gthet(2,it)-thet_pred_mean)/gthet(3,it))**2)
-c     print '(i5,5(1pe14.4))',it,(gthet(j,it),j=1,3)
-c     print '(5(1pe14.4))',thet_pred_mean,theta0(it),sig,sig0(it),ak
-      theta_temp=binorm(thet_pred_mean,theta0(it),sig,sig0(it),ak) 
+     &0.5D0*((gthet(2,it)-thet_pred_mean)
+     &/gthet(3,it))**2)
+      print '(i5,5(1pe14.4))',it,(gthet(j,it),j=1,3)
+      print '(5(1pe14.4))',thet_pred_mean,theta0(it),sig,sig0(it),ak
+      theta_temp=binorm(thet_pred_mean,theta0(it),sig
+     &,sig0(it),ak) 
       if (theta_temp.lt.theta_min) theta_temp=theta_min
       if (theta_temp.gt.theta_max) theta_temp=theta_max
       gen_theta=theta_temp
-c     print '(a)','Exiting GENTHETA.'
+      print '(a)','Exiting GENTHETA.'
       return
       end
 c-------------------------------------------------------------------------
@@ -287,7 +297,7 @@ cd     &   'Error in GenSide: it=',it,' theta=',the
         return
       endif
       tant=dtan(the-pipol)
-      nlobit=nlob(it)
+      nlobit=nlob(iabs(it))
       if (lprint) then
 #ifdef MPI
         print '(a,i4,a)','CG Processor:',me,' Enter Gen_Side.'
@@ -324,7 +334,7 @@ cd     &   'Error in GenSide: it=',it,' theta=',the
           enddo
         enddo
         W1i=a(1,1)-W1i
-        W1(i)=dexp(bsc(i,it)-0.5D0*W1i*zz1*zz1)
+        W1(i)=dexp(bsc(i,iabs(it))-0.5D0*W1i*zz1*zz1)
 c        if (lprint) write(*,'(a,3(1pe15.5)/)')
 c     &          'detAp, W1, anormi',detApi,W1i,anormi
        do k=2,3
@@ -576,17 +586,22 @@ C Calculate a random integer number from the range (n1,n2).
 c--------------------------------------------------------------------------
       double precision function binorm(x1,x2,sigma1,sigma2,ak)
       implicit real*8 (a-h,o-z)
-c     print '(a)','Enter BINORM.'
+      print *,'Enter BINORM.',x1,x2,sigma1,sigma2,ak
       alowb=dmin1(x1-3.0D0*sigma1,x2-3.0D0*sigma2)
       aupb=dmax1(x1+3.0D0*sigma1,x2+3.0D0*sigma2)
       seg=sigma1/(sigma1+ak*sigma2)
       alen=ran_number(0.0D0,1.0D0)
       if (alen.lt.seg) then
+c        print *,'przed anorm',x1,sigma1,alowb,aupb
+c        print *, 'anorm',anorm_distr(x1,sigma1,alowb,aupb)
         binorm=anorm_distr(x1,sigma1,alowb,aupb)
+      
       else
+c        print *,'przed anorm',x2,sigma2,alowb,aupb
+c        print *, 'anorm',anorm_distr(x2,sigma2,alowb,aupb)
         binorm=anorm_distr(x2,sigma2,alowb,aupb)
       endif
-c     print '(a)','Exiting BINORM.'
+      print '(a)','Exiting BINORM.'
       return
       end
 c-----------------------------------------------------------------------
@@ -613,6 +628,7 @@ c
         if(iset.eq.0) then
 1               v1=2.0d0*ran_number(0.0d0,1.0d0)-1.0d0
                 v2=2.0d0*ran_number(0.0d0,1.0d0)-1.0d0
+c                print *,"anorm: iset",iset," v1",v1," v2",v2," rsq",rsq
                 rsq=v1**2+v2**2
                 if(rsq.ge.1.d0.or.rsq.eq.0.0d0) goto 1
                 fac=sqrt(-2.0d0*log(rsq)/rsq)
diff --git a/source/unres/src_CSA_DiL/prng.F b/source/unres/src_CSA_DiL/prng.F
new file mode 100644 (file)
index 0000000..73f6766
--- /dev/null
@@ -0,0 +1,525 @@
+      real*8 function prng_next(me)
+      implicit none
+      integer me
+c
+c Calling sequence:
+c      <new random number> = prng_next ( <ordinal of generator desired> )
+c      <vector of random #s> = vprng ( <ordinal>, <vector>, <length> )
+c
+c This code is based on a sequential algorithm provided by Mal Kalos.
+c This version uses a single 64-bit word to store the initial seeds
+c and additive constants.
+c A 64-bit floating point number is returned.
+c
+c The array "iparam" is full-word aligned, being padded by zeros to
+c let each generator be on a subpage boundary.
+c That is, rows 1 and 2 in a given column of the array are for real, 
+c rows 3-16 are bogus.
+c
+c July 12, 1993: double the number of sequences.  We should have been
+c                using two packets per seed, rather than four
+c October 31, 1993: merge the two arrays of seeds and constants,
+c                and switch to 64-bit arithmetic.
+c June 1994: port to RS6K.  Internal state is kept as 2 64-bit integers
+c The ishft function is defined only on 32-bit integers, so we will
+c shift numbers by dividing by 2**11 and then adding on 2**53-1.
+c
+c November 1994: ishift now works on 64-bit numbers (though it gives a
+c warning).  Thus we go back to using it.  John Zollweg also added the
+c vprng() routine to return vectors of real*8 random numbers.
+c
+      real*8 recip53
+      parameter ( recip53 = 2.0D0**(-53) )
+      integer*8 two
+      parameter ( two = 2**11)
+      integer*8 m,ishift
+c      parameter ( m = 34522712143931 )          ! 11**13
+c      parameter ( ishift = 9007199254740991 )   ! 2**53-1
+
+      integer nmax 
+      integer*8 iparam
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      integer*8 next
+
+crc  g77 doesn't support integer*8 constants
+      m = dint(34522712143931.0d0)
+      ishift = dint(9007199254740991.0d0)
+
+c RS6K porting note: ishift now takes 64-bit integers , with a warning
+      if ( 0.le.me .and. me.le.nmax ) then
+         next = iparam(1,me)*m + iparam(2,me)
+         iparam(1,me) = next
+         prng_next = recip53 * ishft( next, -11 )
+      else
+         prng_next=-1.0D0
+      endif
+
+      end
+c
+c   vprng(me, rn, num)       Get a vector of random numbers
+c
+      subroutine vprng(me,rn,num)
+      real*8 recip53, rn(1)
+      parameter ( recip53 = 2.0D0**(-53) )
+      integer*8 m,iparam
+c      parameter ( m = 34522712143931 )          ! 11**13
+      integer nmax, num, me
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      integer*8 next
+
+crc  g77 doesn't support integer*8 constants
+      m = dint(34522712143931.0d0)
+
+      if ( 0.le.me .and. me.le.nmax ) then
+         do 1 i=1,num
+            next = iparam(1,me)*m + iparam(2,me)
+            iparam(1,me) = next
+            rn(i) = recip53 * ishft( next, -11 )
+    1    continue
+      else
+         rn(1)=-1.0D0
+      endif
+      return
+      end
+
+c
+c   prng_chkpnt          Get the current state of a generator
+c
+c Calling sequence:
+c   logical prng_chkpnt, status
+c   status = prng_chkpnt (me, iseed)    where
+c
+c     me is the particular generator whose state is being gotten
+c     seed is an 4-element integer array where the "l"-values will be saved
+c
+      logical function prng_chkpnt (me, iseed)
+      implicit none
+      integer me 
+      integer*8 iseed
+
+      integer nmax 
+      integer*8 iparam
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      if (me .lt. 0 .or. me .gt. nmax) then
+         prng_chkpnt=.false.
+      else
+        prng_chkpnt=.true.
+        iseed=iparam(1,me)
+      endif
+      end
+c
+c   prng_restart          Restart generator from a saved state
+c
+c Calling sequence:
+c   logical prng_restart, status
+c   status = prng_restart (me, iseed)    where
+c
+c     me is the particular generator being restarted
+c     iseed is a 8-byte integer containing the "l"-values
+c
+      logical function prng_restart (me, iseed)
+      implicit none
+      integer me
+      integer*8 iseed
+
+      integer nmax 
+      integer*8 iparam
+      parameter(nmax=1021)
+      common/ksrprng/iparam(2,0:nmax)
+
+      if (me .lt. 0 .or. me .gt. nmax) then
+         prng_restart=.false.
+        return
+      else
+        prng_restart=.true.
+        iparam(1,me)=iseed
+      endif
+      end
+
+      block data prngblk
+      parameter(nmax=1021)
+      integer*8 iparam
+      common/ksrprng/iparam(2,0:nmax)
+      data (iparam(1,i),iparam(2,i),i=   0,  29) /
+     + 11848219, 11848219, 11848237, 11848237, 11848241, 11848241,
+     + 11848247, 11848247, 11848253, 11848253, 11848271, 11848271,
+     + 11848297, 11848297, 11848313, 11848313, 11848339, 11848339,
+     + 11848351, 11848351, 11848357, 11848357, 11848363, 11848363,
+     + 11848367, 11848367, 11848373, 11848373, 11848379, 11848379,
+     + 11848393, 11848393, 11848433, 11848433, 11848451, 11848451,
+     + 11848469, 11848469, 11848477, 11848477, 11848489, 11848489,
+     + 11848493, 11848493, 11848513, 11848513, 11848523, 11848523,
+     + 11848531, 11848531, 11848537, 11848537, 11848553, 11848553,
+     + 11848589, 11848589, 11848591, 11848591, 11848601, 11848601 /
+      data (iparam(1,i),iparam(2,i),i=  30,  59) /
+     + 11848619, 11848619, 11848637, 11848637, 11848663, 11848663,
+     + 11848673, 11848673, 11848679, 11848679, 11848691, 11848691,
+     + 11848699, 11848699, 11848709, 11848709, 11848717, 11848717,
+     + 11848721, 11848721, 11848729, 11848729, 11848741, 11848741,
+     + 11848751, 11848751, 11848757, 11848757, 11848787, 11848787,
+     + 11848801, 11848801, 11848829, 11848829, 11848853, 11848853,
+     + 11848861, 11848861, 11848867, 11848867, 11848873, 11848873,
+     + 11848891, 11848891, 11848909, 11848909, 11848919, 11848919,
+     + 11848931, 11848931, 11848937, 11848937, 11848961, 11848961,
+     + 11848981, 11848981, 11849021, 11849021, 11849039, 11849039 /
+      data (iparam(1,i),iparam(2,i),i=  60,  89) /
+     + 11849053, 11849053, 11849059, 11849059, 11849069, 11849069,
+     + 11849077, 11849077, 11849087, 11849087, 11849093, 11849093,
+     + 11849107, 11849107, 11849111, 11849111, 11849129, 11849129,
+     + 11849137, 11849137, 11849177, 11849177, 11849183, 11849183,
+     + 11849203, 11849203, 11849231, 11849231, 11849237, 11849237,
+     + 11849239, 11849239, 11849249, 11849249, 11849251, 11849251,
+     + 11849269, 11849269, 11849273, 11849273, 11849291, 11849291,
+     + 11849297, 11849297, 11849309, 11849309, 11849339, 11849339,
+     + 11849359, 11849359, 11849363, 11849363, 11849399, 11849399,
+     + 11849401, 11849401, 11849413, 11849413, 11849417, 11849417 /
+      data (iparam(1,i),iparam(2,i),i=  90, 119) /
+     + 11849437, 11849437, 11849443, 11849443, 11849473, 11849473,
+     + 11849491, 11849491, 11849503, 11849503, 11849507, 11849507,
+     + 11849557, 11849557, 11849567, 11849567, 11849569, 11849569,
+     + 11849573, 11849573, 11849587, 11849587, 11849599, 11849599,
+     + 11849633, 11849633, 11849641, 11849641, 11849653, 11849653,
+     + 11849659, 11849659, 11849671, 11849671, 11849683, 11849683,
+     + 11849689, 11849689, 11849693, 11849693, 11849699, 11849699,
+     + 11849701, 11849701, 11849707, 11849707, 11849713, 11849713,
+     + 11849723, 11849723, 11849741, 11849741, 11849743, 11849743,
+     + 11849759, 11849759, 11849767, 11849767, 11849771, 11849771 /
+      data (iparam(1,i),iparam(2,i),i= 120, 149) /
+     + 11849791, 11849791, 11849801, 11849801, 11849809, 11849809,
+     + 11849813, 11849813, 11849869, 11849869, 11849881, 11849881,
+     + 11849891, 11849891, 11849909, 11849909, 11849923, 11849923,
+     + 11849933, 11849933, 11849947, 11849947, 11849987, 11849987,
+     + 11850001, 11850001, 11850011, 11850011, 11850019, 11850019,
+     + 11850023, 11850023, 11850031, 11850031, 11850049, 11850049,
+     + 11850061, 11850061, 11850073, 11850073, 11850077, 11850077,
+     + 11850103, 11850103, 11850109, 11850109, 11850121, 11850121,
+     + 11850127, 11850127, 11850133, 11850133, 11850149, 11850149,
+     + 11850161, 11850161, 11850169, 11850169, 11850191, 11850191 /
+      data (iparam(1,i),iparam(2,i),i= 150, 179) /
+     + 11850233, 11850233, 11850247, 11850247, 11850259, 11850259,
+     + 11850269, 11850269, 11850283, 11850283, 11850301, 11850301,
+     + 11850341, 11850341, 11850347, 11850347, 11850367, 11850367,
+     + 11850373, 11850373, 11850379, 11850379, 11850389, 11850389,
+     + 11850407, 11850407, 11850427, 11850427, 11850437, 11850437,
+     + 11850469, 11850469, 11850481, 11850481, 11850511, 11850511,
+     + 11850529, 11850529, 11850541, 11850541, 11850557, 11850557,
+     + 11850607, 11850607, 11850611, 11850611, 11850667, 11850667,
+     + 11850677, 11850677, 11850679, 11850679, 11850701, 11850701,
+     + 11850731, 11850731, 11850739, 11850739, 11850749, 11850749 /
+      data (iparam(1,i),iparam(2,i),i= 180, 209) /
+     + 11850791, 11850791, 11850803, 11850803, 11850829, 11850829,
+     + 11850833, 11850833, 11850859, 11850859, 11850877, 11850877,
+     + 11850899, 11850899, 11850907, 11850907, 11850913, 11850913,
+     + 11850919, 11850919, 11850931, 11850931, 11850941, 11850941,
+     + 11850947, 11850947, 11850953, 11850953, 11850961, 11850961,
+     + 11850983, 11850983, 11850991, 11850991, 11850997, 11850997,
+     + 11851031, 11851031, 11851033, 11851033, 11851051, 11851051,
+     + 11851061, 11851061, 11851067, 11851067, 11851093, 11851093,
+     + 11851109, 11851109, 11851123, 11851123, 11851127, 11851127,
+     + 11851139, 11851139, 11851157, 11851157, 11851163, 11851163 /
+      data (iparam(1,i),iparam(2,i),i= 210, 239) /
+     + 11851181, 11851181, 11851201, 11851201, 11851219, 11851219,
+     + 11851291, 11851291, 11851303, 11851303, 11851309, 11851309,
+     + 11851313, 11851313, 11851319, 11851319, 11851349, 11851349,
+     + 11851351, 11851351, 11851361, 11851361, 11851373, 11851373,
+     + 11851403, 11851403, 11851409, 11851409, 11851423, 11851423,
+     + 11851447, 11851447, 11851451, 11851451, 11851481, 11851481,
+     + 11851493, 11851493, 11851519, 11851519, 11851523, 11851523,
+     + 11851529, 11851529, 11851547, 11851547, 11851549, 11851549,
+     + 11851559, 11851559, 11851577, 11851577, 11851589, 11851589,
+     + 11851591, 11851591, 11851597, 11851597, 11851603, 11851603 /
+      data (iparam(1,i),iparam(2,i),i= 240, 269) /
+     + 11851607, 11851607, 11851613, 11851613, 11851621, 11851621,
+     + 11851627, 11851627, 11851639, 11851639, 11851673, 11851673,
+     + 11851681, 11851681, 11851727, 11851727, 11851753, 11851753,
+     + 11851759, 11851759, 11851787, 11851787, 11851793, 11851793,
+     + 11851799, 11851799, 11851813, 11851813, 11851841, 11851841,
+     + 11851859, 11851859, 11851867, 11851867, 11851891, 11851891,
+     + 11851909, 11851909, 11851919, 11851919, 11851927, 11851927,
+     + 11851933, 11851933, 11851949, 11851949, 11851967, 11851967,
+     + 11851997, 11851997, 11852017, 11852017, 11852051, 11852051,
+     + 11852053, 11852053, 11852059, 11852059, 11852083, 11852083 /
+      data (iparam(1,i),iparam(2,i),i= 270, 299) /
+     + 11852089, 11852089, 11852129, 11852129, 11852147, 11852147,
+     + 11852149, 11852149, 11852161, 11852161, 11852171, 11852171,
+     + 11852177, 11852177, 11852209, 11852209, 11852221, 11852221,
+     + 11852237, 11852237, 11852251, 11852251, 11852263, 11852263,
+     + 11852273, 11852273, 11852279, 11852279, 11852287, 11852287,
+     + 11852293, 11852293, 11852297, 11852297, 11852303, 11852303,
+     + 11852311, 11852311, 11852327, 11852327, 11852339, 11852339,
+     + 11852341, 11852341, 11852359, 11852359, 11852369, 11852369,
+     + 11852437, 11852437, 11852453, 11852453, 11852459, 11852459,
+     + 11852473, 11852473, 11852513, 11852513, 11852531, 11852531 /
+      data (iparam(1,i),iparam(2,i),i= 300, 329) /
+     + 11852537, 11852537, 11852539, 11852539, 11852557, 11852557,
+     + 11852573, 11852573, 11852579, 11852579, 11852591, 11852591,
+     + 11852609, 11852609, 11852611, 11852611, 11852623, 11852623,
+     + 11852641, 11852641, 11852647, 11852647, 11852657, 11852657,
+     + 11852663, 11852663, 11852717, 11852717, 11852719, 11852719,
+     + 11852741, 11852741, 11852759, 11852759, 11852767, 11852767,
+     + 11852773, 11852773, 11852803, 11852803, 11852807, 11852807,
+     + 11852809, 11852809, 11852831, 11852831, 11852833, 11852833,
+     + 11852837, 11852837, 11852857, 11852857, 11852873, 11852873,
+     + 11852879, 11852879, 11852891, 11852891, 11852917, 11852917 /
+      data (iparam(1,i),iparam(2,i),i= 330, 359) /
+     + 11852921, 11852921, 11852957, 11852957, 11852959, 11852959,
+     + 11852969, 11852969, 11852983, 11852983, 11852989, 11852989,
+     + 11853001, 11853001, 11853013, 11853013, 11853019, 11853019,
+     + 11853031, 11853031, 11853089, 11853089, 11853133, 11853133,
+     + 11853157, 11853157, 11853161, 11853161, 11853181, 11853181,
+     + 11853203, 11853203, 11853217, 11853217, 11853221, 11853221,
+     + 11853227, 11853227, 11853241, 11853241, 11853307, 11853307,
+     + 11853319, 11853319, 11853323, 11853323, 11853329, 11853329,
+     + 11853367, 11853367, 11853383, 11853383, 11853419, 11853419,
+     + 11853421, 11853421, 11853427, 11853427, 11853449, 11853449 /
+      data (iparam(1,i),iparam(2,i),i= 360, 389) /
+     + 11853451, 11853451, 11853463, 11853463, 11853529, 11853529,
+     + 11853557, 11853557, 11853571, 11853571, 11853601, 11853601,
+     + 11853613, 11853613, 11853617, 11853617, 11853629, 11853629,
+     + 11853649, 11853649, 11853659, 11853659, 11853679, 11853679,
+     + 11853689, 11853689, 11853719, 11853719, 11853731, 11853731,
+     + 11853757, 11853757, 11853761, 11853761, 11853773, 11853773,
+     + 11853791, 11853791, 11853817, 11853817, 11853839, 11853839,
+     + 11853847, 11853847, 11853857, 11853857, 11853869, 11853869,
+     + 11853883, 11853883, 11853887, 11853887, 11853889, 11853889,
+     + 11853893, 11853893, 11853899, 11853899, 11853911, 11853911 /
+      data (iparam(1,i),iparam(2,i),i= 390, 419) /
+     + 11853931, 11853931, 11853943, 11853943, 11853979, 11853979,
+     + 11853991, 11853991, 11854001, 11854001, 11854009, 11854009,
+     + 11854019, 11854019, 11854057, 11854057, 11854061, 11854061,
+     + 11854147, 11854147, 11854159, 11854159, 11854163, 11854163,
+     + 11854169, 11854169, 11854211, 11854211, 11854247, 11854247,
+     + 11854261, 11854261, 11854267, 11854267, 11854279, 11854279,
+     + 11854303, 11854303, 11854327, 11854327, 11854331, 11854331,
+     + 11854333, 11854333, 11854363, 11854363, 11854379, 11854379,
+     + 11854399, 11854399, 11854411, 11854411, 11854429, 11854429,
+     + 11854433, 11854433, 11854439, 11854439, 11854441, 11854441 /
+      data (iparam(1,i),iparam(2,i),i= 420, 449) /
+     + 11854463, 11854463, 11854477, 11854477, 11854489, 11854489,
+     + 11854517, 11854517, 11854519, 11854519, 11854523, 11854523,
+     + 11854529, 11854529, 11854567, 11854567, 11854571, 11854571,
+     + 11854573, 11854573, 11854603, 11854603, 11854607, 11854607,
+     + 11854681, 11854681, 11854691, 11854691, 11854709, 11854709,
+     + 11854723, 11854723, 11854757, 11854757, 11854783, 11854783,
+     + 11854793, 11854793, 11854813, 11854813, 11854847, 11854847,
+     + 11854853, 11854853, 11854873, 11854873, 11854877, 11854877,
+     + 11854883, 11854883, 11854891, 11854891, 11854897, 11854897,
+     + 11854901, 11854901, 11854919, 11854919, 11854937, 11854937 /
+      data (iparam(1,i),iparam(2,i),i= 450, 479) /
+     + 11854961, 11854961, 11854963, 11854963, 11854979, 11854979,
+     + 11855003, 11855003, 11855017, 11855017, 11855023, 11855023,
+     + 11855029, 11855029, 11855033, 11855033, 11855111, 11855111,
+     + 11855141, 11855141, 11855147, 11855147, 11855149, 11855149,
+     + 11855159, 11855159, 11855177, 11855177, 11855203, 11855203,
+     + 11855213, 11855213, 11855219, 11855219, 11855231, 11855231,
+     + 11855267, 11855267, 11855269, 11855269, 11855303, 11855303,
+     + 11855309, 11855309, 11855321, 11855321, 11855329, 11855329,
+     + 11855339, 11855339, 11855351, 11855351, 11855353, 11855353,
+     + 11855357, 11855357, 11855359, 11855359, 11855381, 11855381 /
+      data (iparam(1,i),iparam(2,i),i= 480, 509) /
+     + 11855383, 11855383, 11855387, 11855387, 11855399, 11855399,
+     + 11855407, 11855407, 11855413, 11855413, 11855489, 11855489,
+     + 11855491, 11855491, 11855507, 11855507, 11855521, 11855521,
+     + 11855531, 11855531, 11855549, 11855549, 11855551, 11855551,
+     + 11855567, 11855567, 11855581, 11855581, 11855587, 11855587,
+     + 11855593, 11855593, 11855633, 11855633, 11855653, 11855653,
+     + 11855663, 11855663, 11855687, 11855687, 11855689, 11855689,
+     + 11855699, 11855699, 11855713, 11855713, 11855731, 11855731,
+     + 11855737, 11855737, 11855743, 11855743, 11855747, 11855747,
+     + 11855759, 11855759, 11855773, 11855773, 11855801, 11855801 /
+      data (iparam(1,i),iparam(2,i),i= 510, 539) /
+     + 11855807, 11855807, 11855813, 11855813, 11855827, 11855827,
+     + 11855839, 11855839, 11855869, 11855869, 11855881, 11855881,
+     + 11855903, 11855903, 11855911, 11855911, 11855933, 11855933,
+     + 11855959, 11855959, 11855989, 11855989, 11855993, 11855993,
+     + 11855999, 11855999, 11856001, 11856001, 11856023, 11856023,
+     + 11856049, 11856049, 11856071, 11856071, 11856101, 11856101,
+     + 11856107, 11856107, 11856113, 11856113, 11856139, 11856139,
+     + 11856151, 11856151, 11856161, 11856161, 11856179, 11856179,
+     + 11856193, 11856193, 11856199, 11856199, 11856223, 11856223,
+     + 11856239, 11856239, 11856263, 11856263, 11856269, 11856269 /
+      data (iparam(1,i),iparam(2,i),i= 540, 569) /
+     + 11856281, 11856281, 11856287, 11856287, 11856307, 11856307,
+     + 11856311, 11856311, 11856329, 11856329, 11856343, 11856343,
+     + 11856359, 11856359, 11856371, 11856371, 11856373, 11856373,
+     + 11856409, 11856409, 11856419, 11856419, 11856461, 11856461,
+     + 11856469, 11856469, 11856473, 11856473, 11856479, 11856479,
+     + 11856511, 11856511, 11856517, 11856517, 11856541, 11856541,
+     + 11856547, 11856547, 11856553, 11856553, 11856583, 11856583,
+     + 11856629, 11856629, 11856641, 11856641, 11856653, 11856653,
+     + 11856659, 11856659, 11856673, 11856673, 11856697, 11856697,
+     + 11856709, 11856709, 11856727, 11856727, 11856731, 11856731 /
+      data (iparam(1,i),iparam(2,i),i= 570, 599) /
+     + 11856763, 11856763, 11856809, 11856809, 11856811, 11856811,
+     + 11856821, 11856821, 11856841, 11856841, 11856857, 11856857,
+     + 11856877, 11856877, 11856883, 11856883, 11856899, 11856899,
+     + 11856919, 11856919, 11856947, 11856947, 11856953, 11856953,
+     + 11856979, 11856979, 11857003, 11857003, 11857033, 11857033,
+     + 11857037, 11857037, 11857039, 11857039, 11857049, 11857049,
+     + 11857061, 11857061, 11857067, 11857067, 11857073, 11857073,
+     + 11857081, 11857081, 11857091, 11857091, 11857093, 11857093,
+     + 11857099, 11857099, 11857123, 11857123, 11857127, 11857127,
+     + 11857147, 11857147, 11857151, 11857151, 11857193, 11857193 /
+      data (iparam(1,i),iparam(2,i),i= 600, 629) /
+     + 11857217, 11857217, 11857229, 11857229, 11857243, 11857243,
+     + 11857249, 11857249, 11857267, 11857267, 11857277, 11857277,
+     + 11857291, 11857291, 11857303, 11857303, 11857309, 11857309,
+     + 11857327, 11857327, 11857331, 11857331, 11857333, 11857333,
+     + 11857361, 11857361, 11857367, 11857367, 11857369, 11857369,
+     + 11857393, 11857393, 11857399, 11857399, 11857409, 11857409,
+     + 11857421, 11857421, 11857423, 11857423, 11857451, 11857451,
+     + 11857453, 11857453, 11857457, 11857457, 11857477, 11857477,
+     + 11857481, 11857481, 11857493, 11857493, 11857499, 11857499,
+     + 11857519, 11857519, 11857523, 11857523, 11857529, 11857529 /
+      data (iparam(1,i),iparam(2,i),i= 630, 659) /
+     + 11857543, 11857543, 11857561, 11857561, 11857589, 11857589,
+     + 11857591, 11857591, 11857613, 11857613, 11857621, 11857621,
+     + 11857661, 11857661, 11857667, 11857667, 11857693, 11857693,
+     + 11857697, 11857697, 11857709, 11857709, 11857711, 11857711,
+     + 11857751, 11857751, 11857753, 11857753, 11857759, 11857759,
+     + 11857763, 11857763, 11857777, 11857777, 11857787, 11857787,
+     + 11857793, 11857793, 11857801, 11857801, 11857817, 11857817,
+     + 11857819, 11857819, 11857831, 11857831, 11857837, 11857837,
+     + 11857873, 11857873, 11857877, 11857877, 11857883, 11857883,
+     + 11857889, 11857889, 11857907, 11857907, 11857913, 11857913 /
+      data (iparam(1,i),iparam(2,i),i= 660, 689) /
+     + 11857931, 11857931, 11857969, 11857969, 11857991, 11857991,
+     + 11857999, 11857999, 11858009, 11858009, 11858017, 11858017,
+     + 11858023, 11858023, 11858029, 11858029, 11858039, 11858039,
+     + 11858051, 11858051, 11858057, 11858057, 11858059, 11858059,
+     + 11858101, 11858101, 11858111, 11858111, 11858131, 11858131,
+     + 11858149, 11858149, 11858159, 11858159, 11858177, 11858177,
+     + 11858191, 11858191, 11858201, 11858201, 11858227, 11858227,
+     + 11858243, 11858243, 11858267, 11858267, 11858269, 11858269,
+     + 11858279, 11858279, 11858281, 11858281, 11858291, 11858291,
+     + 11858311, 11858311, 11858323, 11858323, 11858359, 11858359 /
+      data (iparam(1,i),iparam(2,i),i= 690, 719) /
+     + 11858377, 11858377, 11858381, 11858381, 11858387, 11858387,
+     + 11858423, 11858423, 11858443, 11858443, 11858447, 11858447,
+     + 11858479, 11858479, 11858533, 11858533, 11858543, 11858543,
+     + 11858551, 11858551, 11858557, 11858557, 11858569, 11858569,
+     + 11858573, 11858573, 11858579, 11858579, 11858597, 11858597,
+     + 11858599, 11858599, 11858629, 11858629, 11858657, 11858657,
+     + 11858659, 11858659, 11858683, 11858683, 11858701, 11858701,
+     + 11858719, 11858719, 11858723, 11858723, 11858729, 11858729,
+     + 11858747, 11858747, 11858779, 11858779, 11858783, 11858783,
+     + 11858801, 11858801, 11858807, 11858807, 11858813, 11858813 /
+      data (iparam(1,i),iparam(2,i),i= 720, 749) /
+     + 11858839, 11858839, 11858851, 11858851, 11858893, 11858893,
+     + 11858897, 11858897, 11858921, 11858921, 11858947, 11858947,
+     + 11858953, 11858953, 11858969, 11858969, 11858971, 11858971,
+     + 11858989, 11858989, 11859017, 11859017, 11859031, 11859031,
+     + 11859049, 11859049, 11859061, 11859061, 11859073, 11859073,
+     + 11859077, 11859077, 11859079, 11859079, 11859083, 11859083,
+     + 11859101, 11859101, 11859109, 11859109, 11859137, 11859137,
+     + 11859139, 11859139, 11859151, 11859151, 11859157, 11859157,
+     + 11859163, 11859163, 11859167, 11859167, 11859179, 11859179,
+     + 11859187, 11859187, 11859229, 11859229, 11859233, 11859233 /
+      data (iparam(1,i),iparam(2,i),i= 750, 779) /
+     + 11859241, 11859241, 11859247, 11859247, 11859269, 11859269,
+     + 11859293, 11859293, 11859307, 11859307, 11859311, 11859311,
+     + 11859349, 11859349, 11859359, 11859359, 11859371, 11859371,
+     + 11859377, 11859377, 11859383, 11859383, 11859427, 11859427,
+     + 11859433, 11859433, 11859451, 11859451, 11859457, 11859457,
+     + 11859461, 11859461, 11859473, 11859473, 11859481, 11859481,
+     + 11859487, 11859487, 11859493, 11859493, 11859503, 11859503,
+     + 11859509, 11859509, 11859539, 11859539, 11859541, 11859541,
+     + 11859563, 11859563, 11859569, 11859569, 11859571, 11859571,
+     + 11859583, 11859583, 11859599, 11859599, 11859611, 11859611 /
+      data (iparam(1,i),iparam(2,i),i= 780, 809) /
+     + 11859643, 11859643, 11859707, 11859707, 11859713, 11859713,
+     + 11859719, 11859719, 11859739, 11859739, 11859751, 11859751,
+     + 11859791, 11859791, 11859817, 11859817, 11859821, 11859821,
+     + 11859833, 11859833, 11859847, 11859847, 11859853, 11859853,
+     + 11859877, 11859877, 11859889, 11859889, 11859893, 11859893,
+     + 11859901, 11859901, 11859907, 11859907, 11859917, 11859917,
+     + 11859923, 11859923, 11859929, 11859929, 11859961, 11859961,
+     + 11859979, 11859979, 11859989, 11859989, 11859997, 11859997,
+     + 11860021, 11860021, 11860031, 11860031, 11860039, 11860039,
+     + 11860049, 11860049, 11860081, 11860081, 11860087, 11860087 /
+      data (iparam(1,i),iparam(2,i),i= 810, 839) /
+     + 11860097, 11860097, 11860103, 11860103, 11860109, 11860109,
+     + 11860117, 11860117, 11860133, 11860133, 11860151, 11860151,
+     + 11860171, 11860171, 11860207, 11860207, 11860223, 11860223,
+     + 11860231, 11860231, 11860243, 11860243, 11860267, 11860267,
+     + 11860301, 11860301, 11860307, 11860307, 11860327, 11860327,
+     + 11860379, 11860379, 11860397, 11860397, 11860411, 11860411,
+     + 11860469, 11860469, 11860477, 11860477, 11860483, 11860483,
+     + 11860487, 11860487, 11860489, 11860489, 11860493, 11860493,
+     + 11860517, 11860517, 11860547, 11860547, 11860567, 11860567,
+     + 11860573, 11860573, 11860613, 11860613, 11860619, 11860619 /
+      data (iparam(1,i),iparam(2,i),i= 840, 869) /
+     + 11860627, 11860627, 11860637, 11860637, 11860643, 11860643,
+     + 11860649, 11860649, 11860661, 11860661, 11860669, 11860669,
+     + 11860687, 11860687, 11860691, 11860691, 11860697, 11860697,
+     + 11860699, 11860699, 11860703, 11860703, 11860727, 11860727,
+     + 11860741, 11860741, 11860753, 11860753, 11860777, 11860777,
+     + 11860787, 11860787, 11860789, 11860789, 11860811, 11860811,
+     + 11860837, 11860837, 11860859, 11860859, 11860867, 11860867,
+     + 11860889, 11860889, 11860897, 11860897, 11860963, 11860963,
+     + 11860969, 11860969, 11860973, 11860973, 11860993, 11860993,
+     + 11861011, 11861011, 11861033, 11861033, 11861071, 11861071 /
+      data (iparam(1,i),iparam(2,i),i= 870, 899) /
+     + 11861081, 11861081, 11861089, 11861089, 11861093, 11861093,
+     + 11861099, 11861099, 11861107, 11861107, 11861131, 11861131,
+     + 11861141, 11861141, 11861159, 11861159, 11861167, 11861167,
+     + 11861191, 11861191, 11861197, 11861197, 11861207, 11861207,
+     + 11861219, 11861219, 11861221, 11861221, 11861231, 11861231,
+     + 11861237, 11861237, 11861273, 11861273, 11861293, 11861293,
+     + 11861299, 11861299, 11861303, 11861303, 11861327, 11861327,
+     + 11861351, 11861351, 11861357, 11861357, 11861363, 11861363,
+     + 11861371, 11861371, 11861401, 11861401, 11861407, 11861407,
+     + 11861411, 11861411, 11861413, 11861413, 11861429, 11861429 /
+      data (iparam(1,i),iparam(2,i),i= 900, 929) /
+     + 11861441, 11861441, 11861467, 11861467, 11861527, 11861527,
+     + 11861539, 11861539, 11861543, 11861543, 11861557, 11861557,
+     + 11861569, 11861569, 11861573, 11861573, 11861579, 11861579,
+     + 11861581, 11861581, 11861599, 11861599, 11861611, 11861611,
+     + 11861617, 11861617, 11861627, 11861627, 11861639, 11861639,
+     + 11861651, 11861651, 11861659, 11861659, 11861671, 11861671,
+     + 11861683, 11861683, 11861687, 11861687, 11861693, 11861693,
+     + 11861701, 11861701, 11861711, 11861711, 11861713, 11861713,
+     + 11861749, 11861749, 11861791, 11861791, 11861803, 11861803,
+     + 11861819, 11861819, 11861827, 11861827, 11861849, 11861849 /
+      data (iparam(1,i),iparam(2,i),i= 930, 959) /
+     + 11861873, 11861873, 11861879, 11861879, 11861887, 11861887,
+     + 11861911, 11861911, 11861917, 11861917, 11861921, 11861921,
+     + 11861923, 11861923, 11861953, 11861953, 11861959, 11861959,
+     + 11861987, 11861987, 11862007, 11862007, 11862013, 11862013,
+     + 11862029, 11862029, 11862031, 11862031, 11862049, 11862049,
+     + 11862077, 11862077, 11862083, 11862083, 11862157, 11862157,
+     + 11862167, 11862167, 11862199, 11862199, 11862203, 11862203,
+     + 11862217, 11862217, 11862223, 11862223, 11862229, 11862229,
+     + 11862233, 11862233, 11862239, 11862239, 11862241, 11862241,
+     + 11862259, 11862259, 11862269, 11862269, 11862271, 11862271 /
+      data (iparam(1,i),iparam(2,i),i= 960, 989) /
+     + 11862293, 11862293, 11862307, 11862307, 11862313, 11862313,
+     + 11862317, 11862317, 11862343, 11862343, 11862353, 11862353,
+     + 11862373, 11862373, 11862391, 11862391, 11862439, 11862439,
+     + 11862469, 11862469, 11862493, 11862493, 11862527, 11862527,
+     + 11862547, 11862547, 11862563, 11862563, 11862569, 11862569,
+     + 11862577, 11862577, 11862581, 11862581, 11862611, 11862611,
+     + 11862623, 11862623, 11862661, 11862661, 11862673, 11862673,
+     + 11862679, 11862679, 11862701, 11862701, 11862703, 11862703,
+     + 11862713, 11862713, 11862761, 11862761, 11862791, 11862791,
+     + 11862803, 11862803, 11862839, 11862839, 11862841, 11862841 /
+      data (iparam(1,i),iparam(2,i),i= 990,1019) /
+     + 11862857, 11862857, 11862869, 11862869, 11862881, 11862881,
+     + 11862911, 11862911, 11862919, 11862919, 11862959, 11862959,
+     + 11862979, 11862979, 11862989, 11862989, 11862997, 11862997,
+     + 11863021, 11863021, 11863031, 11863031, 11863037, 11863037,
+     + 11863039, 11863039, 11863057, 11863057, 11863067, 11863067,
+     + 11863073, 11863073, 11863099, 11863099, 11863109, 11863109,
+     + 11863121, 11863121, 11863123, 11863123, 11863133, 11863133,
+     + 11863151, 11863151, 11863153, 11863153, 11863171, 11863171,
+     + 11863183, 11863183, 11863207, 11863207, 11863213, 11863213,
+     + 11863237, 11863237, 11863249, 11863249, 11863253, 11863253 /
+      data (iparam(1,i),iparam(2,i),i=1020,1021) /
+     + 11863259, 11863259, 11863279, 11863279 /
+      end
index 8bc9d7a..e4c042a 100644 (file)
@@ -756,6 +756,7 @@ cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
         call mpi_abort(mpi_comm_world,ierror,ierrcode)
       endif
       do n=1,nconf
+        print*,"n=",n
 c  pull out external and internal variables for next start
         call putx(xin,n,rad)
 !        write (iout,*) 'XIN from FEEDIN N=',n
@@ -813,6 +814,7 @@ c  retrieve latest results
          if(iw_pdb.gt.0) 
      &        call write_csa_pdb(xout,eout,nft,ind(1),iw_pdb)
         endif
+        print *,"koniec petli n=",n
       enddo
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 c  no more input