--- /dev/null
+2 0 0 0 0 0 0 2 0 0 2 2 2 2 -1 -1 2 1 1 0 0 0 0 0 -1 -1 -1
+1 1 0
+1 1 0
+1 1 0
+1 1 0
+1 1 0
+18 3 0
+1 1 0
+23 3 1
+24 3 2
+39 3 1
+40 3 2
+35.5 3 -1
+65 3 2
+ 0.58768 4.70723 -0.56523 0.00000 -0.37777 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.10164 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.11870 0.10011 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Cys_Wat
+ 0.48757 4.36438 -0.10008 0.00000 -0.62321 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Met_Wat
+ 0.03940 5.02694 0.01422 0.00000 -0.68367 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Phe_Wat
+ 0.03466 4.57927 0.03194 0.00000 0.48274 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Ile_Wat
+ 0.03466 4.57927 0.03194 0.00000 0.48274 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Leu_Wat
+ 0.46647 5.05492 -0.83472 0.00000 -0.84517 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Val_Wat
+ 0.03471 5.23211 0.14279 0.00000 -0.46298 0.00000 0.00000 0.00000 0.00000 0.00000 0.29795 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Trp_Wat
+ 0.03469 5.80599 -0.39290 0.00000 -0.67892 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Tyr_Wat
+ 0.33886 4.80831 -0.70632 0.00000 -0.17959 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Ala_Wat
+ 3.09494 2.74193 0.72755 0.00000 0.82658 0.00000 0.00000 0.00000 0.00000 0.00000 0.29795 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00324 0.11062 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Gly_Wat
+ 0.42486 4.30589 -0.63453 0.00000 -0.84053 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.11356 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.22062 0.09877 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Thr_Wat
+ 2.89301 4.13052 -0.29291 0.00000 0.75500 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.10017 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.10179 0.10003 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Ser_Wat
+ 0.05200 3.32602 0.44346 0.00000 -0.08444 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.10152 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.13018 0.09971 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Gln_Wat
+ 2.90457 4.19597 -0.27898 0.00000 0.73064 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.10013 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.10140 0.10002 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Asn_Wat
+ 0.35034 4.85725 -0.58884 0.00000 -0.12883 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00331 0.10419 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Glu_Wat
+ 0.26653 3.95618 -0.50716 0.00000 -0.70774 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 -0.45179 0.08149 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Asp_Wat
+ 0.17487 3.34571 0.42472 0.00000 -0.06630 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.10049 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.11562 0.09976 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 His_Wat
+ 0.03652 3.73385 0.07213 0.00000 -0.27415 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.07695 0.10099 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Arg_Wat
+ 0.07463 4.77335 -0.21276 0.00000 -0.74350 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 -0.18116 0.09247 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Lys_Wat
+ 0.03471 6.34515 -0.70002 0.00000 -0.80060 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Pro_Wat
+ 0.48757 4.36438 -0.10008 0.00000 -0.62321 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 SeMet_Wa
+ 0.03940 5.02694 0.01422 0.00000 -0.68367 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 DapBz_Wa
+ 0.33886 4.80831 -0.70632 0.00000 -0.17959 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Aib_Wat
+ 0.33886 4.80831 -0.70632 0.00000 -0.17959 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Abu_Wat
+ 0.35034 4.85725 -0.58884 0.00000 -0.12883 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00331 0.10419 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Glu_Wat
+ 0.35034 4.85725 -0.58884 0.00000 -0.12883 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00331 0.10419 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Glu_Wat
+ 0.26653 3.95618 -0.50716 0.00000 -0.70774 0.00000 0.00000 0.00000 0.00000 0.00000 2.97950 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 -0.45179 0.08149 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 1.00000 1.00000 0.00000 Asp_Wat
+5.69891E-02 3.17002E+00 4.95421E-01 0.00000E+00 8.09856E-01 0.00000E+00 5.95645E-01 5.61069E+00 2.58276E+00 5.23464E+00 4.42980E+00 0.00000E+00 4.14020E-01 0.00000E+00 1 0 0 0 0 4.12992E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -9.04131E+00 3.62811E-02 9.69329E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.35052E+00 0.00000E+00
+2.79444E+00 2.91727E+00 8.33501E-01 0.00000E+00 9.99999E-01 0.00000E+00 9.81661E-01 6.93235E+00 5.97322E-01 4.87889E+00 4.43795E+00 0.00000E+00 4.32100E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.45596E+00 2.83670E+00 7.43295E-01 0.00000E+00 9.99920E-01 0.00000E+00 5.98655E-01 6.55996E+00 3.46739E-02 4.74838E+00 4.56099E+00 0.00000E+00 3.76099E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.48735E-02 4.36914E+00 2.58957E-01 0.00000E+00 -9.67829E-02 0.00000E+00 1.34649E+00 8.42957E-01 3.46574E-02 2.60246E-01 3.82376E+00 0.00000E+00 2.49457E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.52118E-02 4.63748E+00 3.64272E-01 0.00000E+00 1.59087E-01 0.00000E+00 9.53869E-01 4.77396E-01 3.46872E-02 2.68034E-01 4.11193E+00 0.00000E+00 4.00815E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.49660E-02 4.34679E+00 2.09480E-01 0.00000E+00 6.64390E-03 0.00000E+00 1.36813E+00 7.14527E-01 3.46574E-02 6.55911E-01 4.11720E+00 0.00000E+00 1.61904E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46667E-02 2.97236E+00 4.04152E-01 0.00000E+00 -9.99966E-01 0.00000E+00 5.18590E-01 1.20064E+01 3.46665E-02 4.46463E+00 4.08858E+00 0.00000E+00 6.95524E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+4.22623E-02 3.55300E+00 7.52451E-01 0.00000E+00 9.99989E-01 0.00000E+00 1.46027E+00 8.56620E-01 1.45909E-01 2.11520E+00 4.15491E+00 0.00000E+00 4.92352E-01 0.00000E+00 1 0 0 0 0 3.46593E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 0.00000E+00 3.46574E-02 -3.18395E-02 1.03355E-01 3.46574E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E-01 0.00000E+00 1.00000E+00 0.00000E+00
+2.79183E-01 3.38227E+00 4.58897E-01 0.00000E+00 9.10411E-01 0.00000E+00 8.34463E+00 4.99499E+00 6.90949E+00 8.14074E-02 2.51032E+00 0.00000E+00 1.48920E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.84139E+00 2.72458E+00 7.07073E-01 0.00000E+00 9.99751E-01 0.00000E+00 1.55244E+00 5.23148E+00 3.53037E+00 5.29844E+00 4.18614E+00 0.00000E+00 4.88630E-01 0.00000E+00 1 0 0 0 0 3.53127E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.32639E-02 2.00298E+00 4.09554E-01 1.42549E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.15490E+00 0.00000E+00
+7.57785E+00 2.69464E+00 6.59173E-01 0.00000E+00 9.96763E-01 0.00000E+00 2.25237E+00 5.92011E+00 2.68553E+00 5.20529E+00 4.39604E+00 0.00000E+00 4.44385E-02 0.00000E+00 1 0 0 0 0 4.02319E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -3.67928E+00 1.79969E-01 1.50196E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.08211E+00 0.00000E+00
+2.33132E-01 2.96945E+00 4.53728E-01 0.00000E+00 7.24649E-01 0.00000E+00 8.22898E+00 5.43525E+00 4.88545E+00 4.96245E+00 3.78540E+00 0.00000E+00 2.45870E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -7.67410E-01 1.90266E-01 4.59049E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00010E+00 0.00000E+00
+9.48215E+00 2.89500E+00 7.84420E-01 0.00000E+00 9.98547E-01 0.00000E+00 3.43103E+00 6.23121E+00 3.00319E+00 5.18071E+00 4.41926E+00 0.00000E+00 3.36059E-01 0.00000E+00 1 0 0 0 0 1.46750E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 3.95183E+00 4.00136E-01 1.90755E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.47250E-01 0.00000E+00
+1.89742E-01 1.95673E+00 9.00000E-01 0.00000E+00 9.98063E-01 0.00000E+00 7.37857E-01 6.04538E+00 6.70453E-01 5.20022E+00 4.33545E+00 0.00000E+00 4.02223E-01 0.00000E+00 1 0 0 0 0 4.72482E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -1.50741E+00 1.59389E-01 1.88996E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.66766E-01 0.00000E+00
+9.92793E+00 2.64623E+00 2.65809E-01 0.00000E+00 7.31160E-01 0.00000E+00 1.71506E+00 5.95064E+00 3.71304E-02 5.37684E+00 4.25598E+00 0.00000E+00 -7.82987E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.61327E+00
+6.55245E+00 2.72566E+00 8.01881E-01 0.00000E+00 1.00000E+00 0.00000E+00 1.11117E+00 8.74351E+00 1.93383E+00 5.21161E+00 4.43586E+00 0.00000E+00 1.34392E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.82464E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 5.49894E+00
+3.46816E-02 3.48560E+00 1.78424E-01 0.00000E+00 -4.61362E-01 0.00000E+00 7.66510E-01 5.49102E-01 1.78264E-01 3.46574E-02 3.29851E+00 0.00000E+00 4.50032E-01 0.00000E+00 1 0 0 0 0 5.53438E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 0.00000E+00 3.46574E-02 -1.04360E+01 6.05739E-01 3.46574E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 0.00000E+00 9.99998E-01 0.00000E+00
+3.47440E-02 4.93389E+00 6.59651E-01 0.00000E+00 9.09560E-01 0.00000E+00 1.45067E+00 5.85525E+00 5.87716E-01 5.05488E+00 2.87499E+00 0.00000E+00 8.80085E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 0.00000E+00 1.00907E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.58498E-01
+3.47157E-02 3.88251E+00 -1.92325E-02 0.00000E+00 1.43402E-02 0.00000E+00 1.52175E+00 6.22560E+00 7.96165E-01 5.37775E+00 4.03779E+00 0.00000E+00 6.32819E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.83448E+00
+3.47578E-02 4.01099E+00 1.98246E-01 0.00000E+00 -8.62837E-01 0.00000E+00 1.49163E+00 9.06598E-01 2.76008E-01 3.93459E-01 3.70271E+00 0.00000E+00 2.87385E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.79444E+00 2.91727E+00 8.33501E-01 0.00000E+00 9.99999E-01 0.00000E+00 9.81661E-01 6.93235E+00 5.97322E-01 4.87889E+00 4.43795E+00 0.00000E+00 4.32100E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.45596E+00 2.83670E+00 7.43295E-01 0.00000E+00 9.99920E-01 0.00000E+00 5.98655E-01 6.55996E+00 3.46739E-02 4.74838E+00 4.56099E+00 0.00000E+00 3.76099E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.79183E-01 3.38227E+00 4.58897E-01 0.00000E+00 9.10411E-01 0.00000E+00 8.34463E+00 4.99499E+00 6.90949E+00 8.14074E-02 2.51032E+00 0.00000E+00 1.48920E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.79183E-01 3.38227E+00 4.58897E-01 0.00000E+00 9.10411E-01 0.00000E+00 8.34463E+00 4.99499E+00 6.90949E+00 8.14074E-02 2.51032E+00 0.00000E+00 1.48920E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.84139E+00 2.72458E+00 7.07073E-01 0.00000E+00 9.99751E-01 0.00000E+00 1.55244E+00 5.23148E+00 3.53037E+00 5.29844E+00 4.18614E+00 0.00000E+00 4.88630E-01 0.00000E+00 1 0 0 0 0 3.53127E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.32639E-02 2.00298E+00 4.09554E-01 1.42549E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.15490E+00 0.00000E+00
+3.84139E+00 2.72458E+00 7.07073E-01 0.00000E+00 9.99751E-01 0.00000E+00 1.55244E+00 5.23148E+00 3.53037E+00 5.29844E+00 4.18614E+00 0.00000E+00 4.88630E-01 0.00000E+00 1 0 0 0 0 3.53127E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.32639E-02 2.00298E+00 4.09554E-01 1.42549E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.15490E+00 0.00000E+00
+3.84139E+00 2.72458E+00 7.07073E-01 0.00000E+00 9.99751E-01 0.00000E+00 1.55244E+00 5.23148E+00 3.53037E+00 5.29844E+00 4.18614E+00 0.00000E+00 4.88630E-01 0.00000E+00 1 0 0 0 0 3.53127E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.32639E-02 2.00298E+00 4.09554E-01 1.42549E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.15490E+00 0.00000E+00
+5.40033E-01 4.41744E+00 6.40857E-01 0.00000E+00 9.57375E-01 0.00000E+00 1.41585E+00 5.84146E+00 2.22249E+00 5.28349E+00 3.26342E+00 0.00000E+00 -5.41440E-01 0.00000E+00 1 0 0 0 0 5.60413E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -2.40876E+00 3.96712E-02 2.97530E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.87733E+00 0.00000E+00
+5.12100E+00 2.96715E+00 7.87559E-01 0.00000E+00 2.06425E-01 0.00000E+00 3.01066E+00 7.14712E+00 8.28676E-01 4.84834E+00 4.32378E+00 0.00000E+00 5.87310E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.48244E-02 5.14144E+00 9.01893E-02 0.00000E+00 -7.75830E-02 0.00000E+00 5.96794E-01 6.73438E+00 3.16339E+00 4.88579E+00 4.64210E+00 0.00000E+00 5.22546E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46604E-02 5.80191E+00 3.28443E-01 0.00000E+00 3.22840E-01 0.00000E+00 7.85743E-01 4.19477E+00 1.94123E+00 4.02330E-01 4.51534E+00 0.00000E+00 -3.32144E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46578E-02 5.80193E+00 3.29777E-01 0.00000E+00 2.93787E-02 0.00000E+00 1.55415E+00 9.03230E-01 3.66901E-02 2.24085E-01 4.12118E+00 0.00000E+00 1.86773E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46590E-02 5.53520E+00 9.56586E-02 0.00000E+00 -8.07031E-01 0.00000E+00 1.02853E-01 1.92892E-01 2.99391E-01 2.61753E-01 4.91832E+00 0.00000E+00 5.31986E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46585E-02 6.44209E+00 5.92565E-01 0.00000E+00 7.38064E-01 0.00000E+00 5.00150E+00 5.00182E+00 4.99940E+00 4.99850E+00 3.06919E+00 0.00000E+00 4.89695E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+6.03185E+00 2.59077E+00 8.10272E-01 0.00000E+00 1.00000E+00 0.00000E+00 3.47984E+00 6.83831E+00 8.29494E-02 5.19389E+00 4.09247E+00 0.00000E+00 6.48867E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -3.41371E+00 3.46971E-02 1.69559E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 2.79208E+00 0.00000E+00
+2.07431E-01 4.25937E+00 3.76895E-01 0.00000E+00 5.43405E-01 0.00000E+00 8.37959E+00 5.12688E+00 6.82556E+00 3.47687E-02 2.55561E+00 0.00000E+00 6.33778E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.02661E+01 3.01002E+00 3.54124E-01 0.00000E+00 3.33249E-01 0.00000E+00 8.06931E+00 5.46078E+00 4.08705E+00 4.78382E+00 4.56983E+00 0.00000E+00 2.81656E-01 0.00000E+00 1 0 0 0 0 3.52831E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -3.82122E+00 3.46574E-02 3.50807E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.99964E-01 0.00000E+00
+2.88710E+00 2.91588E+00 1.06001E-01 0.00000E+00 1.90717E-01 0.00000E+00 4.94190E+00 6.67311E+00 4.69960E+00 5.17698E+00 4.49051E+00 0.00000E+00 7.36208E-02 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -3.62781E+00 1.27188E-01 3.53663E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00364E+01 0.00000E+00
+3.50268E-02 3.43342E+00 2.69955E-01 0.00000E+00 4.40375E-01 0.00000E+00 3.50964E+00 6.12026E+00 4.18020E+00 5.24424E+00 4.33443E+00 0.00000E+00 3.66855E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -2.64593E+00 3.65549E-02 9.43631E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00243E+01 0.00000E+00
+3.66584E+00 3.04726E+00 3.54686E-01 0.00000E+00 3.76939E-01 0.00000E+00 3.61242E+00 5.62115E+00 4.28055E+00 5.26053E+00 4.63323E+00 0.00000E+00 2.52640E-01 0.00000E+00 1 0 0 0 0 3.15535E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -7.79305E+00 6.48937E-01 5.34975E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.35220E+00 0.00000E+00
+4.00006E+00 3.02283E+00 3.36720E-01 0.00000E+00 2.61431E-01 0.00000E+00 4.02255E+00 5.60722E+00 4.16026E+00 5.20224E+00 4.72663E+00 0.00000E+00 2.66122E-01 0.00000E+00 1 0 0 0 0 1.38418E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -7.73594E+00 1.88595E+00 3.46997E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.57392E+00 0.00000E+00
+4.42356E-02 5.03283E+00 4.68392E-01 0.00000E+00 8.84686E-01 0.00000E+00 3.92084E+00 4.03821E+00 6.66349E+00 5.50903E+00 5.39108E+00 0.00000E+00 6.92130E-02 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.38418E+00
+8.07755E+00 3.45780E+00 7.18077E-01 0.00000E+00 7.50121E-01 0.00000E+00 1.38508E+01 1.11095E+01 1.00815E+01 4.69231E+00 4.45843E+00 0.00000E+00 2.41181E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.82464E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.36851E+00
+5.16850E-01 4.62992E+00 7.30337E-01 0.00000E+00 9.99922E-01 0.00000E+00 2.42528E+00 5.10463E+00 1.36519E+00 5.27152E+00 4.41810E+00 0.00000E+00 2.84141E-01 0.00000E+00 1 0 0 0 0 1.55118E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -6.71818E+00 1.79649E+01 1.63310E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.68335E+00 0.00000E+00
+3.47644E-02 5.68514E+00 7.31836E-01 0.00000E+00 9.96188E-01 0.00000E+00 2.96272E-01 7.16537E+00 4.09320E+00 5.10387E+00 7.45972E+00 0.00000E+00 -5.87795E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.73051E+00
+3.39499E+00 1.98986E+00 8.76452E-01 0.00000E+00 -9.89326E-01 0.00000E+00 8.97340E+00 1.06200E+01 1.30437E+00 4.77487E+00 3.86883E+00 0.00000E+00 5.75799E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 2.03908E+00
+5.16940E-01 4.02533E+00 6.11267E-02 0.00000E+00 -9.99091E-01 0.00000E+00 3.53497E+00 2.84895E+00 4.14480E-02 1.74708E-01 3.22577E+00 0.00000E+00 3.49203E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+5.12100E+00 2.96715E+00 7.87559E-01 0.00000E+00 2.06425E-01 0.00000E+00 3.01066E+00 7.14712E+00 8.28676E-01 4.84834E+00 4.32378E+00 0.00000E+00 5.87310E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.48244E-02 5.14144E+00 9.01893E-02 0.00000E+00 -7.75830E-02 0.00000E+00 5.96794E-01 6.73438E+00 3.16339E+00 4.88579E+00 4.64210E+00 0.00000E+00 5.22546E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.07431E-01 4.25937E+00 3.76895E-01 0.00000E+00 5.43405E-01 0.00000E+00 8.37959E+00 5.12688E+00 6.82556E+00 3.47687E-02 2.55561E+00 0.00000E+00 6.33778E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.07431E-01 4.25937E+00 3.76895E-01 0.00000E+00 5.43405E-01 0.00000E+00 8.37959E+00 5.12688E+00 6.82556E+00 3.47687E-02 2.55561E+00 0.00000E+00 6.33778E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.02661E+01 3.01002E+00 3.54124E-01 0.00000E+00 3.33249E-01 0.00000E+00 8.06931E+00 5.46078E+00 4.08705E+00 4.78382E+00 4.56983E+00 0.00000E+00 2.81656E-01 0.00000E+00 1 0 0 0 0 3.52831E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -3.82122E+00 3.46574E-02 3.50807E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.99964E-01 0.00000E+00
+1.02661E+01 3.01002E+00 3.54124E-01 0.00000E+00 3.33249E-01 0.00000E+00 8.06931E+00 5.46078E+00 4.08705E+00 4.78382E+00 4.56983E+00 0.00000E+00 2.81656E-01 0.00000E+00 1 0 0 0 0 3.52831E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -3.82122E+00 3.46574E-02 3.50807E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.99964E-01 0.00000E+00
+1.02661E+01 3.01002E+00 3.54124E-01 0.00000E+00 3.33249E-01 0.00000E+00 8.06931E+00 5.46078E+00 4.08705E+00 4.78382E+00 4.56983E+00 0.00000E+00 2.81656E-01 0.00000E+00 1 0 0 0 0 3.52831E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -3.82122E+00 3.46574E-02 3.50807E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.99964E-01 0.00000E+00
+2.28417E+00 3.11967E+00 7.53444E-01 0.00000E+00 9.96355E-01 0.00000E+00 2.17928E+00 5.59535E+00 3.51811E+00 5.22672E+00 4.58445E+00 0.00000E+00 7.68478E-02 0.00000E+00 1 0 0 0 0 2.75910E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -7.74441E+00 9.50957E-02 1.88460E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.95370E+00 0.00000E+00
+5.61517E-02 4.42618E+00 5.27067E-01 0.00000E+00 -9.35872E-01 0.00000E+00 3.94007E+00 5.62577E+00 5.69252E+00 5.01986E+00 4.39277E+00 0.00000E+00 6.44335E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.62482E+00 2.89552E+00 7.19111E-01 0.00000E+00 9.99654E-01 0.00000E+00 6.22191E-01 6.71472E+00 7.59725E-01 4.71116E+00 4.92977E+00 0.00000E+00 2.83076E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.58977E-02 4.71829E+00 2.02388E-01 0.00000E+00 -2.35139E-01 0.00000E+00 7.11334E-01 4.36250E-01 6.12703E-02 7.84648E-02 4.03229E+00 0.00000E+00 8.91540E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.49641E-02 5.18602E+00 2.64948E-01 0.00000E+00 1.08992E-01 0.00000E+00 6.58099E-01 3.00470E-01 6.08084E-02 2.57582E-01 4.15735E+00 0.00000E+00 3.35986E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+5.67675E-01 3.87476E+00 3.04352E-01 0.00000E+00 -5.44993E-01 0.00000E+00 1.19380E+00 6.77261E-01 3.46574E-02 3.73156E-01 4.49257E+00 0.00000E+00 3.42866E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+4.51330E-02 4.55712E+00 5.85804E-01 0.00000E+00 6.42133E-01 0.00000E+00 2.97864E+00 5.31736E+00 5.45513E+00 4.92960E+00 4.77883E+00 0.00000E+00 4.17574E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.29390E+00 3.57431E+00 7.03840E-01 0.00000E+00 1.75902E-01 0.00000E+00 8.02108E+00 7.46898E+00 7.07329E+00 5.15783E+00 4.28321E+00 0.00000E+00 5.38351E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -1.79887E+01 5.32217E+01 3.46684E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00352E+01 0.00000E+00
+4.32837E-02 4.32014E+00 4.68194E-01 0.00000E+00 9.27389E-01 0.00000E+00 6.98569E+00 4.18288E+00 5.94434E+00 4.44280E-01 3.11342E+00 0.00000E+00 2.42507E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+6.09339E+00 3.10691E+00 7.44248E-01 0.00000E+00 9.99993E-01 0.00000E+00 1.15828E+00 5.50927E+00 2.87480E+00 5.26883E+00 5.10064E+00 0.00000E+00 -2.41577E-02 0.00000E+00 1 0 0 0 0 3.79731E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -6.35055E+00 3.87025E-02 2.19172E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.20089E+00 0.00000E+00
+3.06949E+00 3.15302E+00 6.82317E-01 0.00000E+00 9.99435E-01 0.00000E+00 2.60888E+00 5.53742E+00 3.91430E+00 5.21037E+00 4.88989E+00 0.00000E+00 -2.36374E-01 0.00000E+00 1 0 0 0 0 2.00875E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -5.63044E+00 1.03176E-01 1.60774E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.96490E+00 0.00000E+00
+4.56308E-02 3.98470E+00 4.43004E-01 0.00000E+00 8.59186E-01 0.00000E+00 1.10114E+01 8.19441E+00 8.79990E+00 4.89299E+00 3.54455E+00 0.00000E+00 1.96321E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 9.99942E-01 -3.32776E+00 2.37226E-01 4.01340E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.83870E-01 0.00000E+00
+3.79299E+00 3.04104E+00 8.27103E-01 0.00000E+00 9.94705E-01 0.00000E+00 2.04203E+00 5.82234E+00 2.85011E+00 5.21114E+00 4.62378E+00 0.00000E+00 2.16705E-01 0.00000E+00 1 0 0 0 0 3.47082E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -5.14175E+00 9.65693E-02 2.86577E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.88727E+00 0.00000E+00
+4.56223E+00 3.20950E+00 7.63502E-01 0.00000E+00 9.99714E-01 0.00000E+00 1.89230E+00 5.77012E+00 2.86322E+00 5.17611E+00 5.02402E+00 0.00000E+00 -7.59549E-02 0.00000E+00 1 0 0 0 0 3.46868E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -3.98524E+00 9.95547E-02 2.01432E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.94608E+00 0.00000E+00
+3.30017E+00 2.92509E+00 7.27062E-01 0.00000E+00 9.99966E-01 0.00000E+00 8.26725E-01 6.17445E+00 5.73256E-01 5.35144E+00 4.57875E+00 0.00000E+00 3.25776E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.37409E+00
+5.42960E+00 3.06271E+00 7.23838E-01 0.00000E+00 9.99995E-01 0.00000E+00 1.61783E+00 5.98429E+00 2.27894E+00 5.33877E+00 4.78185E+00 0.00000E+00 -9.02454E-02 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 8.68574E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.18496E+00
+4.32522E-02 3.66063E+00 3.07488E-01 0.00000E+00 -9.45380E-01 0.00000E+00 3.41055E+00 2.23532E+00 3.08684E+00 2.18157E-01 2.92952E+00 0.00000E+00 4.69712E-01 0.00000E+00 1 0 0 0 0 1.85456E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 0.00000E+00 3.46574E-02 -1.74782E+01 1.10687E+00 3.46574E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 0.00000E+00 1.00000E+00 0.00000E+00
+6.92086E-01 2.93371E+00 6.49846E-01 0.00000E+00 9.99924E-01 0.00000E+00 2.09511E+00 5.91998E+00 2.49017E+00 5.41138E+00 3.97807E+00 0.00000E+00 6.85924E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55288E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 4.92867E-02
+3.54380E+00 2.87461E+00 7.92366E-01 0.00000E+00 8.09536E-01 0.00000E+00 2.81889E+00 6.56658E+00 1.68978E+00 5.34487E+00 4.06904E+00 0.00000E+00 4.72316E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 -8.95252E-06
+3.51189E-02 4.15569E+00 2.58552E-01 0.00000E+00 -7.36162E-01 0.00000E+00 1.15411E+00 5.45723E-01 2.59335E-01 4.43090E-01 3.73618E+00 0.00000E+00 2.34780E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+5.61517E-02 4.42618E+00 5.27067E-01 0.00000E+00 -9.35872E-01 0.00000E+00 3.94007E+00 5.62577E+00 5.69252E+00 5.01986E+00 4.39277E+00 0.00000E+00 6.44335E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.62482E+00 2.89552E+00 7.19111E-01 0.00000E+00 9.99654E-01 0.00000E+00 6.22191E-01 6.71472E+00 7.59725E-01 4.71116E+00 4.92977E+00 0.00000E+00 2.83076E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+4.32837E-02 4.32014E+00 4.68194E-01 0.00000E+00 9.27389E-01 0.00000E+00 6.98569E+00 4.18288E+00 5.94434E+00 4.44280E-01 3.11342E+00 0.00000E+00 2.42507E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+4.32837E-02 4.32014E+00 4.68194E-01 0.00000E+00 9.27389E-01 0.00000E+00 6.98569E+00 4.18288E+00 5.94434E+00 4.44280E-01 3.11342E+00 0.00000E+00 2.42507E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+6.09339E+00 3.10691E+00 7.44248E-01 0.00000E+00 9.99993E-01 0.00000E+00 1.15828E+00 5.50927E+00 2.87480E+00 5.26883E+00 5.10064E+00 0.00000E+00 -2.41577E-02 0.00000E+00 1 0 0 0 0 3.79731E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -6.35055E+00 3.87025E-02 2.19172E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.20089E+00 0.00000E+00
+6.09339E+00 3.10691E+00 7.44248E-01 0.00000E+00 9.99993E-01 0.00000E+00 1.15828E+00 5.50927E+00 2.87480E+00 5.26883E+00 5.10064E+00 0.00000E+00 -2.41577E-02 0.00000E+00 1 0 0 0 0 3.79731E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -6.35055E+00 3.87025E-02 2.19172E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.20089E+00 0.00000E+00
+6.09339E+00 3.10691E+00 7.44248E-01 0.00000E+00 9.99993E-01 0.00000E+00 1.15828E+00 5.50927E+00 2.87480E+00 5.26883E+00 5.10064E+00 0.00000E+00 -2.41577E-02 0.00000E+00 1 0 0 0 0 3.79731E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -6.35055E+00 3.87025E-02 2.19172E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.20089E+00 0.00000E+00
+9.97525E-01 4.18439E+00 6.77393E-01 0.00000E+00 9.97401E-01 0.00000E+00 7.31614E-01 5.99156E+00 5.12616E-01 5.16280E+00 4.78825E+00 0.00000E+00 9.56238E-02 0.00000E+00 1 0 0 0 0 4.65837E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 3.56508E-02 -4.26339E+00 3.60916E-02 1.77410E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.16288E+00 0.00000E+00
+3.46727E-02 5.69741E+00 4.04784E-01 0.00000E+00 -9.99866E-01 0.00000E+00 1.16515E-01 5.94188E+00 3.69978E+00 4.62284E+00 8.10860E+00 0.00000E+00 5.24562E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.47274E-02 5.72829E+00 4.22155E-01 0.00000E+00 7.66092E-01 0.00000E+00 1.96847E+00 1.27875E+00 3.46574E-02 7.89190E-01 3.22939E+00 0.00000E+00 6.50326E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.20559E+00 4.06572E+00 4.77602E-01 0.00000E+00 5.40577E-01 0.00000E+00 2.38951E+00 1.87918E+00 8.19911E-02 5.90656E-02 3.60279E+00 0.00000E+00 8.77579E-03 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.96601E-01 5.10122E+00 2.68703E-01 0.00000E+00 -4.52724E-01 0.00000E+00 3.99839E+00 2.99985E+00 4.64968E-02 5.05432E-01 2.69873E+00 0.00000E+00 6.80297E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.12752E+00 3.88423E+00 4.97048E-01 0.00000E+00 4.07795E-01 0.00000E+00 5.56897E+00 4.56273E+00 2.80514E+00 3.47279E-02 2.95142E+00 0.00000E+00 8.60369E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+4.43247E-02 5.96363E+00 6.22133E-01 0.00000E+00 6.59525E-01 0.00000E+00 4.99795E+00 4.99646E+00 5.00199E+00 5.00205E+00 3.00902E+00 0.00000E+00 4.87953E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.60283E+00 3.69545E+00 4.46478E-01 0.00000E+00 -4.18371E-01 0.00000E+00 2.10068E+01 1.58383E+01 1.51080E+01 3.43730E+00 3.75758E+00 0.00000E+00 5.41827E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -1.30256E+01 1.84492E+00 7.79864E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.33937E-01 0.00000E+00
+1.68072E-01 4.16175E+00 5.68603E-01 0.00000E+00 9.48088E-01 0.00000E+00 8.97167E+00 5.79892E+00 6.95293E+00 3.46913E-02 2.57100E+00 0.00000E+00 -8.29498E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+8.58958E+00 2.90846E+00 4.75166E-01 0.00000E+00 4.68647E-01 0.00000E+00 5.27341E+00 3.05368E+00 2.01220E+00 4.96228E+00 4.67037E+00 0.00000E+00 3.26370E-01 0.00000E+00 1 0 0 0 0 1.65824E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -4.64188E+00 3.46574E-02 3.49844E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00002E+00 0.00000E+00
+5.09309E+00 3.14518E+00 7.42449E-02 0.00000E+00 1.80707E-01 0.00000E+00 4.80775E+00 6.52559E+00 4.52505E+00 5.16510E+00 4.84967E+00 0.00000E+00 5.71291E-02 0.00000E+00 1 0 0 0 0 5.03834E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -4.06644E+00 4.08528E-01 3.60491E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00425E+01 0.00000E+00
+6.76474E+00 2.75846E+00 2.46733E-01 0.00000E+00 -1.08755E-01 0.00000E+00 2.14001E+00 4.84231E+00 3.44439E-01 5.01561E+00 4.57000E+00 0.00000E+00 2.36119E-01 0.00000E+00 1 0 0 0 0 3.77635E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 9.99460E-01 -2.97361E+00 1.67156E+00 7.49271E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.36961E+00 0.00000E+00
+2.44973E+00 3.89948E+00 6.40772E-01 0.00000E+00 8.05096E-01 0.00000E+00 3.20174E+00 6.31066E+00 1.82636E+00 5.30654E+00 2.21824E+00 0.00000E+00 -5.53286E-02 0.00000E+00 1 0 0 0 0 3.57955E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -3.35319E+00 1.15949E-01 1.86536E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.18263E+00 0.00000E+00
+3.87252E+00 2.92265E+00 8.36180E-01 0.00000E+00 9.99800E-01 0.00000E+00 3.59703E+00 5.59916E+00 4.47886E+00 5.21866E+00 4.76985E+00 0.00000E+00 3.57409E-01 0.00000E+00 1 0 0 0 0 7.33078E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -2.18070E+00 3.95470E-02 1.84526E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.03630E+00 0.00000E+00
+4.65424E-02 4.62683E+00 5.19170E-01 0.00000E+00 2.67693E-01 0.00000E+00 2.95628E+00 4.56700E+00 8.15507E+00 5.29944E+00 4.25301E+00 0.00000E+00 3.03135E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 3.10373E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 3.27171E+00
+1.28156E-01 3.60914E+00 1.12677E-02 0.00000E+00 -2.96058E-01 0.00000E+00 2.00930E+00 2.14180E+00 1.15580E+01 1.92910E+04 1.48840E+00 0.00000E+00 1.91371E-01 0.00000E+00 1 0 0 0 0 2.65936E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 2.88450E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 4.13140E+01
+1.34610E+00 3.11693E+00 5.92165E-01 0.00000E+00 9.99955E-01 0.00000E+00 1.23508E+00 4.69920E+00 1.44058E+00 5.28377E+00 4.60458E+00 0.00000E+00 -5.18828E-01 0.00000E+00 1 0 0 0 0 4.44947E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -8.01020E+00 1.81209E+01 3.24712E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.87027E+00 0.00000E+00
+2.07983E+00 7.53143E-01 6.80598E-01 0.00000E+00 9.99918E-01 0.00000E+00 2.26601E+00 6.07154E+00 6.53911E-01 5.32884E+00 4.96206E+00 0.00000E+00 5.34943E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55288E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.60319E+00
+3.52541E-01 1.81925E+00 9.00000E-01 0.00000E+00 -8.45490E-01 0.00000E+00 4.30142E+00 7.33512E+00 1.11363E+00 5.24191E+00 4.32761E+00 0.00000E+00 5.17484E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 6.45011E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.76514E+00
+1.66728E+00 2.45263E+00 5.26602E-01 0.00000E+00 -9.99369E-01 0.00000E+00 3.75332E+00 2.65059E+00 3.46574E-02 3.46911E-02 2.78295E+00 0.00000E+00 3.36168E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46727E-02 5.69741E+00 4.04784E-01 0.00000E+00 -9.99866E-01 0.00000E+00 1.16515E-01 5.94188E+00 3.69978E+00 4.62284E+00 8.10860E+00 0.00000E+00 5.24562E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.47274E-02 5.72829E+00 4.22155E-01 0.00000E+00 7.66092E-01 0.00000E+00 1.96847E+00 1.27875E+00 3.46574E-02 7.89190E-01 3.22939E+00 0.00000E+00 6.50326E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.68072E-01 4.16175E+00 5.68603E-01 0.00000E+00 9.48088E-01 0.00000E+00 8.97167E+00 5.79892E+00 6.95293E+00 3.46913E-02 2.57100E+00 0.00000E+00 -8.29498E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+1.68072E-01 4.16175E+00 5.68603E-01 0.00000E+00 9.48088E-01 0.00000E+00 8.97167E+00 5.79892E+00 6.95293E+00 3.46913E-02 2.57100E+00 0.00000E+00 -8.29498E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+8.58958E+00 2.90846E+00 4.75166E-01 0.00000E+00 4.68647E-01 0.00000E+00 5.27341E+00 3.05368E+00 2.01220E+00 4.96228E+00 4.67037E+00 0.00000E+00 3.26370E-01 0.00000E+00 1 0 0 0 0 1.65824E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -4.64188E+00 3.46574E-02 3.49844E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00002E+00 0.00000E+00
+8.58958E+00 2.90846E+00 4.75166E-01 0.00000E+00 4.68647E-01 0.00000E+00 5.27341E+00 3.05368E+00 2.01220E+00 4.96228E+00 4.67037E+00 0.00000E+00 3.26370E-01 0.00000E+00 1 0 0 0 0 1.65824E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -4.64188E+00 3.46574E-02 3.49844E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00002E+00 0.00000E+00
+8.58958E+00 2.90846E+00 4.75166E-01 0.00000E+00 4.68647E-01 0.00000E+00 5.27341E+00 3.05368E+00 2.01220E+00 4.96228E+00 4.67037E+00 0.00000E+00 3.26370E-01 0.00000E+00 1 0 0 0 0 1.65824E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 -4.64188E+00 3.46574E-02 3.49844E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00002E+00 0.00000E+00
+8.62098E-01 3.73125E+00 7.00692E-01 0.00000E+00 9.88889E-01 0.00000E+00 2.88466E-01 5.33382E+00 3.04897E+00 5.12037E+00 6.22961E+00 0.00000E+00 -9.31924E-01 0.00000E+00 1 0 0 0 0 3.86053E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 2.82330E+00 9.04100E-02 1.80415E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.16114E+00 0.00000E+00
+3.46578E-02 4.46387E+00 5.12970E-01 0.00000E+00 -8.91112E-01 0.00000E+00 6.41980E+00 7.07251E+00 6.09954E+00 4.90064E+00 3.89398E+00 0.00000E+00 5.20810E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.72698E+00 2.58793E+00 9.00000E-01 0.00000E+00 1.00000E+00 0.00000E+00 1.43984E+01 5.14046E+00 6.67706E-01 1.14234E+00 3.01613E+00 0.00000E+00 2.44037E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.47035E-02 5.09875E+00 5.04219E-01 0.00000E+00 8.91590E-01 0.00000E+00 6.59671E-01 3.60321E-01 3.49066E-02 4.60251E-02 4.05916E+00 0.00000E+00 -2.08514E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46599E-02 5.09831E+00 4.71043E-01 0.00000E+00 4.96533E-01 0.00000E+00 6.71879E-01 3.31190E-01 7.55743E-02 8.68816E-02 4.24439E+00 0.00000E+00 2.25158E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.74845E-02 5.08895E+00 1.55634E-01 0.00000E+00 -7.28288E-01 0.00000E+00 5.35658E-01 6.28217E-02 3.46574E-02 4.21658E-02 4.76316E+00 0.00000E+00 2.43749E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46579E-02 2.78072E+00 9.00000E-01 0.00000E+00 -9.99516E-01 0.00000E+00 2.45365E+00 6.02447E+00 3.53157E-02 5.50721E+00 4.25319E+00 0.00000E+00 5.55443E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.47629E-02 5.15825E+00 6.75281E-01 0.00000E+00 1.00000E+00 0.00000E+00 5.29977E-01 5.94692E+00 2.12632E+00 5.32595E+00 5.02401E+00 0.00000E+00 6.98022E-01 0.00000E+00 1 0 0 0 0 3.84865E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -7.68570E+00 3.92230E-02 8.86179E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.99600E-01 0.00000E+00
+3.46589E-02 5.14054E+00 3.09344E-01 0.00000E+00 7.55465E-01 0.00000E+00 8.64675E+00 4.56216E+00 7.60156E+00 2.16120E-01 3.18280E+00 0.00000E+00 7.71626E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46626E-02 4.59525E+00 6.70482E-01 0.00000E+00 9.98615E-01 0.00000E+00 2.27779E+00 6.15674E+00 2.47068E+00 5.34326E+00 1.31224E+00 0.00000E+00 -1.25169E-01 0.00000E+00 1 0 0 0 0 3.46578E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.37751E-02 -4.92279E-01 1.39322E+00 1.53857E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.42315E+00 0.00000E+00
+4.82225E+00 3.40841E+00 6.65188E-01 0.00000E+00 9.96584E-01 0.00000E+00 1.57438E+00 5.83803E+00 2.52051E+00 5.14874E+00 5.28915E+00 0.00000E+00 -3.29340E-01 0.00000E+00 1 0 0 0 0 3.71778E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 3.10474E+00 8.71437E-02 1.46444E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.11272E+00 0.00000E+00
+3.48212E-02 4.52622E+00 5.71513E-01 0.00000E+00 9.99715E-01 0.00000E+00 1.14151E+00 9.48303E-01 3.46617E-02 9.92175E+00 4.07600E+00 0.00000E+00 4.69311E-01 0.00000E+00 1 0 0 0 0 3.15315E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00907E-01 7.49595E+00 5.05569E+01 2.25084E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+01 0.00000E+00
+2.20138E+00 2.96657E+00 5.24561E-01 0.00000E+00 9.94153E-01 0.00000E+00 2.25014E+00 6.17356E+00 1.60783E+00 5.22631E+00 3.87712E+00 0.00000E+00 6.15718E-01 0.00000E+00 1 0 0 0 0 1.23959E-01 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 9.56493E-01 6.79183E-02 1.49726E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.43421E+00 0.00000E+00
+1.60652E+00 3.72149E+00 7.24008E-01 0.00000E+00 9.99546E-01 0.00000E+00 1.48146E-01 6.01322E+00 3.57132E-02 5.00817E+00 6.48242E+00 0.00000E+00 -9.96397E-01 0.00000E+00 1 0 0 0 0 3.46797E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -1.57206E+00 3.50256E-02 1.87839E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 9.74397E-01 0.00000E+00
+3.47326E-02 4.55046E+00 7.27910E-01 0.00000E+00 9.95261E-01 0.00000E+00 8.19747E-01 5.98971E+00 3.47676E-02 5.37011E+00 4.41816E+00 0.00000E+00 4.43236E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 5.55281E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 1.57257E+00
+3.27904E+00 3.26999E+00 7.54352E-01 0.00000E+00 9.99576E-01 0.00000E+00 1.64589E+00 6.15598E+00 2.96046E-01 5.34135E+00 4.44316E+00 0.00000E+00 2.52188E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.82464E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 -4.74670E-03
+3.42380E+00 1.22887E+00 9.00000E-01 0.00000E+00 9.99672E-01 0.00000E+00 4.86590E+00 6.32326E+00 1.12727E-01 5.18911E+00 3.56132E+00 0.00000E+00 5.75126E-01 0.00000E+00 1 0 0 0 0 6.94340E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.55619E-02 -1.70312E+00 1.80914E+01 4.35390E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.88453E+00 0.00000E+00
+3.46710E-02 4.95669E+00 6.04046E-01 0.00000E+00 -6.76790E-01 0.00000E+00 1.70302E+00 5.91142E+00 4.60400E-01 6.09915E+00 3.37648E+00 0.00000E+00 6.51286E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.10924E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 3.41926E-01
+3.47258E-02 5.24509E+00 6.56205E-01 0.00000E+00 9.67802E-01 0.00000E+00 3.00032E+00 8.31601E+00 2.80541E+00 4.46214E+00 1.36151E+00 0.00000E+00 6.82412E-01 0.00000E+00 1 0 0 0 0 3.46574E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00908E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 8.54380E+00
+1.20595E-01 4.37444E+00 2.52255E-01 0.00000E+00 -9.99908E-01 0.00000E+00 8.30188E-01 3.60469E-02 1.53706E-01 1.22822E-01 4.44863E+00 0.00000E+00 3.82358E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46578E-02 4.46387E+00 5.12970E-01 0.00000E+00 -8.91112E-01 0.00000E+00 6.41980E+00 7.07251E+00 6.09954E+00 4.90064E+00 3.89398E+00 0.00000E+00 5.20810E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+2.72698E+00 2.58793E+00 9.00000E-01 0.00000E+00 1.00000E+00 0.00000E+00 1.43984E+01 5.14046E+00 6.67706E-01 1.14234E+00 3.01613E+00 0.00000E+00 2.44037E-01 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46589E-02 5.14054E+00 3.09344E-01 0.00000E+00 7.55465E-01 0.00000E+00 8.64675E+00 4.56216E+00 7.60156E+00 2.16120E-01 3.18280E+00 0.00000E+00 7.71626E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46589E-02 5.14054E+00 3.09344E-01 0.00000E+00 7.55465E-01 0.00000E+00 8.64675E+00 4.56216E+00 7.60156E+00 2.16120E-01 3.18280E+00 0.00000E+00 7.71626E-02 0.00000E+00 1 0 0 0 0 0.00000E+00 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.00000E+00 0.00000E+00
+3.46626E-02 4.59525E+00 6.70482E-01 0.00000E+00 9.98615E-01 0.00000E+00 2.27779E+00 6.15674E+00 2.47068E+00 5.34326E+00 1.31224E+00 0.00000E+00 -1.25169E-01 0.00000E+00 1 0 0 0 0 3.46578E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.37751E-02 -4.92279E-01 1.39322E+00 1.53857E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.42315E+00 0.00000E+00
+3.46626E-02 4.59525E+00 6.70482E-01 0.00000E+00 9.98615E-01 0.00000E+00 2.27779E+00 6.15674E+00 2.47068E+00 5.34326E+00 1.31224E+00 0.00000E+00 -1.25169E-01 0.00000E+00 1 0 0 0 0 3.46578E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.37751E-02 -4.92279E-01 1.39322E+00 1.53857E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.42315E+00 0.00000E+00
+3.46626E-02 4.59525E+00 6.70482E-01 0.00000E+00 9.98615E-01 0.00000E+00 2.27779E+00 6.15674E+00 2.47068E+00 5.34326E+00 1.31224E+00 0.00000E+00 -1.25169E-01 0.00000E+00 1 0 0 0 0 3.46578E-02 0 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 0 1.00000E+00 7.37751E-02 -4.92279E-01 1.39322E+00 1.53857E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00 1.00000E+00 1.42315E+00 0.00000E+00
--- /dev/null
+5.87676E-01 4.70723E+00 -5.65227E-01 0.00000E+00 -3.77774E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 1.01639E-01 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 1.18696E-01 1.00112E-01 Cys_Wat
+4.87567E-01 4.36438E+00 -1.00079E-01 0.00000E+00 -6.23208E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Met_Wat
+3.94000E-02 5.02694E+00 1.42243E-02 0.00000E+00 -6.83672E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Phe_Wat
+3.46595E-02 4.57927E+00 3.19434E-02 0.00000E+00 4.82737E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Ile_Wat
+3.46595E-02 4.57927E+00 3.19434E-02 0.00000E+00 4.82737E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Leu_Wat
+4.66467E-01 5.05492E+00 -8.34722E-01 0.00000E+00 -8.45171E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Val_Wat
+3.47087E-02 5.23211E+00 1.42793E-01 0.00000E+00 -4.62977E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E-01 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Trp_Wat
+3.46916E-02 5.80599E+00 -3.92895E-01 0.00000E+00 -6.78919E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Tyr_Wat
+3.38861E-01 4.80831E+00 -7.06324E-01 0.00000E+00 -1.79590E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Ala_Wat
+3.09494E+00 2.74193E+00 7.27546E-01 0.00000E+00 8.26579E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E-01 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 3.24110E-03 1.10622E-01 Gly_Wat
+4.24862E-01 4.30589E+00 -6.34529E-01 0.00000E+00 -8.40527E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 1.13556E-01 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 2.20617E-01 9.87730E-02 Thr_Wat
+2.89301E+00 4.13052E+00 -2.92911E-01 0.00000E+00 7.55002E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 1.00167E-01 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 1.01790E-01 1.00034E-01 Ser_Wat
+5.20032E-02 3.32602E+00 4.43461E-01 0.00000E+00 -8.44420E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 1.01520E-01 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 1.30177E-01 9.97089E-02 Gln_Wat
+2.90457E+00 4.19597E+00 -2.78983E-01 0.00000E+00 7.30642E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 1.00128E-01 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 1.01402E-01 1.00021E-01 Asn_Wat
+3.50339E-01 4.85725E+00 -5.88844E-01 0.00000E+00 -1.28830E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 3.31407E-03 1.04186E-01 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Glu_Wat
+2.66525E-01 3.95618E+00 -5.07160E-01 0.00000E+00 -7.07736E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 -4.51793E-01 8.14863E-02 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Asp_Wat
+1.74874E-01 3.34571E+00 4.24724E-01 0.00000E+00 -6.62958E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 1.00493E-01 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 1.15621E-01 9.97586E-02 His_Wat
+3.65235E-02 3.73385E+00 7.21323E-02 0.00000E+00 -2.74153E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 7.69543E-02 1.00987E-01 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Arg_Wat
+7.46277E-02 4.77335E+00 -2.12759E-01 0.00000E+00 -7.43504E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 -1.81159E-01 9.24723E-02 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Lys_Wat
+3.47102E-02 6.34515E+00 -7.00023E-01 0.00000E+00 -8.00596E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Pro_Wat
+4.87567E-01 4.36438E+00 -1.00079E-01 0.00000E+00 -6.23208E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 SeMet_Wat
+3.94000E-02 5.02694E+00 1.42243E-02 0.00000E+00 -6.83672E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 DapBz_Wat
+3.38861E-01 4.80831E+00 -7.06324E-01 0.00000E+00 -1.79590E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Aib_Wat
+3.38861E-01 4.80831E+00 -7.06324E-01 0.00000E+00 -1.79590E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Abu_Wat
+3.50339E-01 4.85725E+00 -5.88844E-01 0.00000E+00 -1.28830E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 3.31407E-03 1.04186E-01 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Glu_Wat
+3.50339E-01 4.85725E+00 -5.88844E-01 0.00000E+00 -1.28830E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 3.31407E-03 1.04186E-01 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Glu_Wat
+2.66525E-01 3.95618E+00 -5.07160E-01 0.00000E+00 -7.07736E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.97950E+00 0.00000E+00 0.00000E+00 0.00000E+00 1 0.00 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00 0.00 0.00 0.00 0.00 1.0 1.0 1.0 -4.51793E-01 8.14863E-02 0.00 0.00 0.00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.0 1.0 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 Asp_Wat
mnum=molnum(i)
iti=iabs(itype(i,mnum))
if (mnum.eq.5) iti=itype(i,mnum)
+! if (itype(i,mnum).eq.ntyp1_molec(mnum)) then
+! do j=1,3
+! incr(j)=d_t(j,i)
+! enddo
+! endif
if (itype(i,1).eq.10 .or. itype(i,mnum).eq.ntyp1_molec(mnum)&
.or.mnum.ge.3) then
do j=1,3
mnum=molnum(i)
if (itype(i,mnum).ne.ntyp1_molec(mnum)&
.and.itype(i+1,mnum).ne.ntyp1_molec(mnum)) then
- if (mnum.eq.5) Ip(mnum)=Isc(itype(i,mnum),mnum)
+ if (mnum.eq.5) Ip(mnum)=0.0
! write (iout,*) "i",i
! write (iout,*) "i",i," mag1",mag1," mag2",mag2
do j=1,3
do i=innt,inct-1
mnum=molnum(i)
- if (mnum.eq.5) mp(mnum)=msc(itype(i,mnum),mnum)
- if (mnum.eq.5) mp(mnum)=msc(itype(i,mnum),mnum)
+ if (mnum.eq.5) mp(mnum)=0.0d0
+! if (mnum.eq.5) mp(mnum)=msc(itype(i,mnum),mnum)
!c write (iout,*) i,(d_t(j,i),j=1,3),(d_t(j,i+1),j=1,3)
do j=1,3
v(j)=0.5d0*(d_t(j,i)+d_t(j,i+1))
!c to the velocities of the first Calpha.
do i=innt,inct
mnum=molnum(i)
+ if (mnum.eq.5) then
+ iti=itype(i,mnum)
+ else
iti=iabs(itype(i,mnum))
+ endif
if (itype(i,1).eq.10.or.mnum.ge.3.or. itype(i,mnum).eq.ntyp1_molec(mnum)) then
!c write (iout,*) i,iti,(d_t(j,i),j=1,3)
do j=1,3
do j=1,3
incr(j)=d_t(j,i+1)-d_t(j,i)
enddo
- if (mnum.eq.5) Ip(mnum)=Isc(itype(i,mnum),mnum)
+ if (mnum.eq.5) Ip(mnum)=0.0d0
!c write (iout,*) i,(incr(j),j=1,3)
!c write (iout,*) "Kinetic rotp:",i,(incr(j),j=1,3)
KEr_p=KEr_p+Ip(mnum)*(incr(1)*incr(1)+incr(2)*incr(2)&
!c The rotational part of the side chain virtual bond
do i=innt,inct
mnum=molnum(i)
- iti=iabs(itype(i,mnum))
+! iti=iabs(itype(i,mnum))
+ if (mnum.eq.5) then
+ iti=itype(i,mnum)
+ else
+ iti=iabs(itype(i,mnum))
+ endif
+
! if (iti.ne.10.and.mnum.lt.3) then
if (itype(i,1).ne.10.and.mnum.lt.3.and. itype(i,mnum).ne.ntyp1_molec(mnum)) then
do j=1,3
write (tytul,'("time",f8.2)') totT
if(mdpdb) then
+ write(iout,*) "before hairpin"
call hairpin(.true.,nharp,iharp)
+ write(iout,*) "before secondary"
call secondary2(.true.)
+ write(iout,*) "before pdbout"
call pdbout(potE,tytul,ipdb)
! call enerprint(potEcomp)
else
! call ginv_mult(fric_work, d_af_work)
! call ginv_mult(stochforcvec, d_as_work)
#ifdef FIVEDIAG
+ write(iout,*) "forces before fivediaginv"
+ do i=1,dimen*3
+ write(iout,*) "fricwork",i,fric_work(i)
+ enddo
call fivediaginv_mult(dimen,fric_work, d_af_work)
call fivediaginv_mult(dimen,stochforcvec, d_as_work)
if (large) then
! include 'COMMON.INTERACT'
! include 'COMMON.IOUNITS'
! include 'COMMON.NAMES'
- real(kind=8),dimension(6*nres) :: stochforcvec,d_as_work1 !(MAXRES6) maxres6=6*maxres
+ real(kind=8),dimension(:),allocatable :: stochforcvec,d_as_work1 !(MAXRES6) maxres6=6*maxres
real(kind=8) :: cos60 = 0.5d0, sin60 = 0.86602540378443864676d0
integer :: i,j,ind,inres,mnum
+ if (.not.allocated(stochforcvec)) allocate(stochforcvec(6*nres))
+ if (.not.allocated(d_as_work1)) allocate(d_as_work1(6*nres))
! Revised 3/31/05 AL: correlation between random contributions to
! position and velocity increments included.
! The correlation coefficients are calculated at low-friction limit.
! include 'COMMON.TIME1'
real(kind=8) :: xv,sigv,lowb,highb ,Ek1
#ifdef FIVEDIAG
- integer ichain,n,innt,inct,ibeg,ierr
+ integer ichain,n,innt,inct,ibeg,ierr,innt_org
real(kind=8) ,allocatable, dimension(:):: work
integer,allocatable,dimension(:) :: iwork
! double precision Ghalf(mmaxres2_chain),Geigen(maxres2_chain),&
EK=0.0d0
Ek3=0.0d0
#ifdef DEBUG
- write(iout,*), nchain
+ write(iout,*), "nchain",nchain
#endif
do ichain=1,nchain
ind=0
! ghalf=0.0d0
n=dimen_chain(ichain)
innt=iposd_chain(ichain)
- if (molnum(innt).eq.5) go to 137
+! innt_org=
+ innt_org=chain_border(1,ichain)
+ if ((molnum(innt_org).eq.5).or.(molnum(innt_org).eq.4)) go to 137
if(.not.allocated(ghalf)) print *,"COCO"
if(.not.allocated(Ghalf)) allocate(Ghalf(1300*(1300+1)/2))
ghalf=0.0d0
enddo
#endif
137 continue
- write(iout,*) "HERE,",n
+ write(iout,*) "HERE,",n,innt
+ innt_org=chain_border(1,ichain)
xv=0.0d0
ii=0
do i=1,n
do k=1,3
ii=ii+1
- if (molnum(innt).eq.5) geigen(i)=1.0/msc(itype(innt+i-1,5),5)
-
+ mnum=molnum(innt_org)
+ if (molnum(innt_org).ge.4) geigen(i)=3.0/msc(itype(innt_org+i-1,mnum),mnum)
+! if (molnum(innt).eq.5) write(iout,*) "typ",i,innt-1+i,itype(innt+i-1,5)
sigv=dsqrt((Rb*t_bath)/geigen(i))
lowb=-5*sigv
highb=5*sigv
d_t_work_new(ii)=anorm_distr(xv,sigv,lowb,highb)
EK=EK+0.5d0*geigen(i)*d_t_work_new(ii)**2
-! write (iout,*) "i",i," ii",ii," geigen",geigen(i), &
-! " d_t_work_new",d_t_work_new(ii)
+ write (iout,*) "i",i," ii",ii," geigen",geigen(i), &
+ " d_t_work_new",d_t_work_new(ii),innt_org+i-1
enddo
enddo
- if (molnum(innt).eq.5) then
-
+ if (molnum(innt_org).ge.4) then
+ mnum=molnum(innt_org)
do k=1,3
do i=1,n
ind=(i-1)*3+k
d_t_work(ind)=0.0d0
- masinv=1.0/msc(itype(innt+i-1,5),5)
+ masinv=1.0d0/msc(itype(innt_org+i-1,mnum),mnum)
d_t_work(ind)=d_t_work(ind)&
+masinv*d_t_work_new((i-1)*3+k)
enddo
d_t(:,0)=d_t(:,1)
d_t(:,1)=0.0d0
endif
+ if (large) then
+ write (iout,*)
+ write (iout,*) "Random vel after 1st transf the Calpha,SC space"
+ write (iout,'(3hORG,1h(,i5,1h),3f10.5)') 0,(d_t(j,0),j=1,3)
+ do i=1,nres
+ mnum=molnum(i)
+ write (iout,'(a3,1h(,i5,1h),3f10.5,3x,3f10.5)')&
+ restyp(itype(i,mnum),mnum),i,(d_t(j,i),j=1,3),(d_t(j,i+nres),j=1,3)
+ enddo
+ endif
+
!c d_a(:,0)=d_a(:,1)
!c d_a(:,1)=0.0d0
!c write (iout,*) "Shifting accelerations"
do ichain=2,nchain
+ write(iout,*) "nchain",ichain,chain_border1(1,ichain),molnum(chain_border1(1,ichain))
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
!c write (iout,*) "ichain",chain_border1(1,ichain)-1,
!c & chain_border1(1,ichain)
d_t(:,chain_border1(1,ichain)-1)=d_t(:,chain_border1(1,ichain))
enddo
!c write (iout,*) "Adding accelerations"
do ichain=2,nchain
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
!c write (iout,*) "chain",ichain,chain_border1(1,ichain)-1,
!c & chain_border(2,ichain-1)
d_t(:,chain_border1(1,ichain)-1)=&
do ichain=2,nchain
write (iout,*) "chain",ichain,chain_border1(1,ichain)-1,&
chain_border(2,ichain-1)
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
+
d_t(:,chain_border1(1,ichain)-1)=&
d_t(:,chain_border1(1,ichain)-1)+d_t(:,chain_border(2,ichain-1))
d_t(:,chain_border(2,ichain-1))=0.0d0
enddo
+ if (large) then
+ write (iout,*)
+ write (iout,*) "Random vel after 2nd transf the Calpha,SC space"
+ write (iout,'(3hORG,1h(,i5,1h),3f10.5)') 0,(d_t(j,0),j=1,3)
+ do i=1,nres
+ mnum=molnum(i)
+ write (iout,'(a3,1h(,i5,1h),3f10.5,3x,3f10.5)')&
+ restyp(itype(i,mnum),mnum),i,(d_t(j,i),j=1,3),(d_t(j,i+nres),j=1,3)
+ enddo
+ endif
+
#else
ibeg=0
!c do j=1,3
M_PEP=0.0d0
do i=nnt,nct-1
mnum=molnum(i)
- if (mnum.ge.5) mp(mnum)=msc(itype(i,mnum),mnum)
+! if (mnum.ge.5) mp(mnum)=msc(itype(i,mnum),mnum)
write(iout,*) "WTF",itype(i,mnum),i,mnum,mp(mnum)
! if (itype(i,mnum).eq.ntyp1_molec(mnum)) cycle
M_PEP=M_PEP+mp(mnum)
M_SC=0.0d0
do i=nnt,nct
mnum=molnum(i)
- if (mnum.ge.5) cycle
+! if (mnum.ge.5) cycle
iti=iabs(itype(i,mnum))
M_SC=M_SC+msc(iabs(iti),mnum)
inres=i+nres
+ if (mnum.ge.4) inres=i
do j=1,3
cm(j)=cm(j)+msc(iabs(iti),mnum)*c(j,inres)
enddo
! write(iout,*) "Center of mass:",cm
do i=nnt,nct-1
mnum=molnum(i)
- if (mnum.ge.5) mp(mnum)=msc(itype(i,mnum),mnum)
+! if (mnum.ge.5) mp(mnum)=msc(itype(i,mnum),mnum)
do j=1,3
pr(j)=c(j,i)+0.5d0*dc(j,i)-cm(j)
enddo
do i=nnt,nct
mnum=molnum(i)
iti=iabs(itype(i,mnum))
- if (mnum.ge.5) cycle
+! if (mnum.ge.5) cycle
inres=i+nres
+ if (mnum.ge.4) inres=i
do j=1,3
pr(j)=c(j,inres)-cm(j)
enddo
summas=0.0d0
do i=nnt,nct
mnum=molnum(i)
- if (mnum.ge.5) mp(mnum)=msc(itype(i,mnum),mnum)
+! if (mnum.ge.4) mp(mnum)=msc(itype(i,mnum),mnum)
if (i.lt.nct) then
summas=summas+mp(mnum)
do j=1,3
! print *,i,j,vv(j),d_t(j,i)
enddo
endif
- if (mnum.ne.4) then
+! if (mnum.ne.4) then
amas=msc(iabs(itype(i,mnum)),mnum)
- else
- amas=0.0d0
- endif
+! else
+! amas=0.0d0
+! endif
summas=summas+amas
if (itype(i,mnum).ne.10 .and. itype(i,mnum).ne.ntyp1_molec(mnum)&
.and.(mnum.lt.4)) then
!c innt=chain_border(1,1)
!c inct=chain_border(2,1)
do i=innt,inct
+ mnum=molnum(i)
vvec(ind+1)=v_work(j,i)
ind=ind+1
! if (iabs(itype(i)).ne.10) then
! include 'COMMON.IOUNITS'
integer :: IERROR
integer :: i,j,ind,ind1,m,ichain,innt,inct
- logical :: lprn = .false.
+ logical :: lprn = .true.
real(kind=8) :: dtdi !el ,gamvec(2*nres)
!el real(kind=8),dimension(2*nres,2*nres) :: ginvfric,fcopy
! real(kind=8),allocatable,dimension(:,:) :: fcopy
enddo
!c DU1fric part
do ichain=1,nchain
- mnum=molnum(i)
-
ind=iposd_chain(ichain)
innt=chain_border(1,ichain)
inct=chain_border(2,ichain)
do i=innt,inct
+ mnum=molnum(i)
if (iabs(itype(i,1)).ne.10.and.mnum.le.2) then
ind=ind+2
else
enddo
!c DU2fric part
do ichain=1,nchain
- mnum=molnum(i)
ind=iposd_chain(ichain)
innt=chain_border(1,ichain)
inct=chain_border(2,ichain)
do i=innt,inct-1
+ mnum=molnum(i)
if (iabs(itype(i,1)).ne.10.and.mnum.le.2) then
DU2fric(ind)=gamvec(i-nnt+1)/4
DU2fric(ind+1)=0.0d0
! the line below might be wrong
#ifdef FIVEDIAG
real(kind=8) :: rs(2*nres),xsolv(2*nres)
-#ifdef CHECK5DSOL
+!#ifdef CHECK5DSOL
real(kind=8) :: rscheck(2*nres),rsold(2*nres)
-#endif
+!#endif
#endif
logical :: lprn = .false.
!el common /cipiszcze/ itime
endif
#ifdef CHECK5DSOL
! Check the solution
- call fivediagmult(n,DMorig(innt),DU1orig(innt),DU2orig(innt),
- & xsolv,rscheck)
+ call fivediagmult(n,DMorig(innt),DU1orig(innt),DU2orig(innt),&
+ xsolv,rscheck)
do i=1,n
- write(iout,*) "i",i,"rsold",rsold(i),"rscheck",rscheck(i),
- & "ratio",rscheck(i)/rsold(i)
+ write(iout,*) "i",i,"rsold",rsold(i),"rscheck",rscheck(i),&
+ "ratio",rscheck(i)/rsold(i)
enddo
! end check
#endif
#define CHUJ
#ifdef CHUJ
do ichain=2,nchain
+!TEST 27.06.2023 godz 16.00
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
!c write (iout,*) "ichain",chain_border1(1,ichain)-1,
!c & chain_border1(1,ichain)
d_a(:,chain_border1(1,ichain)-1)=d_a(:,chain_border1(1,ichain))
enddo
!c write (iout,*) "Adding accelerations"
do ichain=2,nchain
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
!c write (iout,*) "chain",ichain,chain_border1(1,ichain)-1,
!c & chain_border(2,ichain-1)
d_a(:,chain_border1(1,ichain)-1)=&
mnum=molnum(innt)
inct=chain_border(2,ichain)
if (mnum.eq.5) mp(mnum)=0.0
+ if (mnum.eq.5) ip(mnum)=0.0
! if (mnum.eq.5) mp(mnum)=msc(itype(innt,mnum),mnum)
DM(ind)=mp(mnum)/4+ip(mnum)/4
if (iabs(itype(innt,1)).eq.10.or.molnum(innt).gt.2) then
! if (mnum.eq.5) mp(mnum)=msc(itype(i,mnum),mnum)
! if (mnum.eq.5) ip(mnum)=isc(itype(i,mnum),mnum)
if (mnum.eq.5) mp(mnum)=0.0
+ if (mnum.eq.5) ip(mnum)=0.0
! if (mnum.eq.5) mp(mnum)=msc(itype(i,mnum),mnum)
DM(ind)=2*ip(mnum)/4+mp(mnum)/2
if (iabs(itype(i,1)).eq.10.or.molnum(i).gt.2) then
enddo
DMorig=DM
DU1orig=DU1
- DU2orig=DU2
+ DU2orig=DU2
+ gmatout=.true.
if (gmatout) then
write (iout,*)"The upper part of the five-diagonal inertia matrix"
endif
subroutine fivediaginv_mult(ndim,forces,d_a_vec)
use energy_data, only:nchain,chain_border,nct,nnt,molnum,&
chain_border1,itype
+ use geometry_data, only: nside
integer ndim
- double precision forces(3*ndim),accel(3,0:maxres2),rs(ndim), &
- xsolv(ndim),d_a_vec(6*nres)
+
+ real(kind=8),dimension(:),allocatable ::forces,d_a_vec
+ real(kind=8),dimension(:),allocatable :: xsolv,rs
+ real(kind=8),dimension(:,:),allocatable :: accel
integer i,j,ind,ichain,n,iposc,innt,inct,inct_prev,mnum
accel=0.0d0
+!#define DEBUG
+ if (.not.allocated(forces)) allocate(forces(6*nres))
+ if (.not.allocated(d_a_vec)) allocate(d_a_vec(6*nres))
+ if (.not.allocated(xsolv)) allocate(xsolv(3*ndim))
+ if (.not.allocated(rs)) allocate(rs(3*ndim))
+ if (.not.allocated(accel)) allocate(accel(3,0:2*nres))
+
+#ifdef DEBUG
+ do i=1,6*nres
+ write(iout,*) "received forces",i,forces(i)
+ enddo
+#endif
do j=1,3
!Compute accelerations in Calpha and SC
do ichain=1,nchain
innt=chain_border(1,ichain)
inct=chain_border(2,ichain)
do i=iposc,iposc+n-1
+! write(iout,*) "index",3*(i-1)+j,forces(3*(i-1)+j)
rs(i-iposc+1)=forces(3*(i-1)+j)
enddo
#ifdef DEBUG
- write (iout,*) "j",j," chain",ichain
+ write (iout,*) "j",j," chain",ichain,"n",n
+ write (iout,*) "innt",innt,inct,iposc
write (iout,*) "rs"
- write (iout,'(f10.5)') (rs(i),i=1,n)
+ do i=1,n
+ write (iout,'(i5,f10.5)') i,rs(i)
+ enddo
#endif
call FDISYS (n,DM(iposc),DU1(iposc),DU2(iposc),rs,xsolv)
#ifdef DEBUG
#ifdef DEBUG
write (iout,*) "accel in CA-SC basis"
do i=1,nres
- write (iout,'(i5,3f10.5,5x,3f10.5)') i,(accel(j,i),j=1,3),
- & (accel(j,i+nres),j=1,3)
+ write (iout,'(i5,3f10.5,5x,3f10.5)') i,(accel(j,i),j=1,3),&
+ (accel(j,i+nres),j=1,3)
enddo
write (iout,*) "nnt",nnt
#endif
accel(:,1)=0.0d0
endif
do ichain=2,nchain
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
accel(:,chain_border1(1,ichain)-1)= &
accel(:,chain_border1(1,ichain))
accel(:,chain_border1(1,ichain))=0.0d0
enddo
do ichain=2,nchain
+ if (molnum(chain_border1(1,ichain)+1).eq.5) cycle
accel(:,chain_border1(1,ichain)-1)= &
accel(:,chain_border1(1,ichain)-1) &
+accel(:,chain_border(2,ichain-1))
#endif
return
end subroutine
-
+#undef DEBUG
#else
!-----------------------------------------------------------------------------
! include 'COMMON.FFIELD'
! include 'COMMON.NAMES'
integer :: ncont
- integer,dimension(2,100*nres) :: icont !(2,maxcont) (maxcont=12*maxres)
+ integer,dimension(:,:),allocatable :: icont !(2,maxcont) (maxcont=12*maxres)
integer :: nharp
integer,dimension(4,nres) :: iharp !(4,nres/3)(4,maxres/3)
logical :: lprint,not_done
real(kind=8) :: rcomp=6.0d0
!el local variables
integer :: i,j,kkk,k,i1,i2,it1,it2,j1,ii1,jj1
-! allocate(icont(2,100*nres))
-
+ if (.not.allocated(icont)) then
+ allocate(icont(2,100*nres_molec(1)+1))
+ endif
ncont=0
kkk=0
! print *,'nnt=',nnt,' nct=',nct
real(kind=8) :: ael6_i,ael3_i
real(kind=8),dimension(2,2) :: app_,bpp_,rpp_
integer :: ncont
- integer,dimension(2,100*nres) :: icont !(2,100*nres)(2,maxcont) (maxcont=12*maxres)
- real(kind=8),dimension(100*nres) :: econt !(maxcont)
+ integer,dimension(:,:),allocatable :: icont !(2,100*nres)(2,maxcont) (maxcont=12*maxres)
+ real(kind=8),dimension(:),allocatable :: econt !(maxcont)
!el local variables
integer :: i,j,k,iteli,itelj,i1,i2,it1,it2,ic1,ic2
real(kind=8) :: elcutoff,elecutoff_14,rri,ees,evdw
data elpp_3 / 0.0503d0, 0.0000d0, 0.0000d0, 0.0692d0/
!el allocate(econt(100*nres)) !(maxcont)
-
+ if (.not.allocated(icont)) then
+ allocate(icont(2,100*nres_molec(1)+1))
+ endif
+ if (.not.allocated(econt)) then
+ allocate(econt(100*nres_molec(1)+1))
+ endif
elcutoff = -0.3d0
elecutoff_14 = -0.5d0
if (lprint) write (iout,'(a)') &
! include 'COMMON.CONTROL'
integer :: ncont,i,j,i1,j1,nbeta,nstrand,ii1,jj1,ij,nhelix,&
iii1,jjj1
- integer,dimension(2,100*nres) :: icont !(2,maxcont) (maxcont=12*maxres)
+ integer,dimension(:,:),allocatable :: icont !(2,maxcont) (maxcont=12*maxres)
integer,dimension(nres,0:4) :: isec !(maxres,4)
integer,dimension(nres) :: nsec !(maxres)
logical :: lprint,not_done !,freeres
!el external freeres
!el allocate(icont(2,100*nres),isec(nres,4),nsec(nres))
-
+ if (.not.allocated(icont)) then
+ allocate(icont(2,100*nres+1))
+ endif
if(.not.dccart) call chainbuild_cart
if(.not.allocated(hfrag)) allocate(hfrag(2,nres/3)) !(2,maxres/3)
!d call write_pdb(99,'sec structure',0d0)
! include 'COMMON.IOUNITS'
! include 'COMMON.DISTFIT'
- integer :: ncont,icont(2,nres*nres/2),isec(nres,3)
+ integer :: ncont,isec(nres,3)
logical :: lprint,not_done
- real(kind=4) :: dcont(nres*nres/2),d
+ real(kind=4) :: d
real(kind=4) :: rcomp = 7.0
real(kind=4) :: rbeta = 5.2
real(kind=4) :: ralfa = 5.2
real(kind=8),dimension(3) :: xpi,xpj
integer :: i,k,j,i1,j1,nbeta,nstrand,ii1,jj1,ij,iii1,jjj1,&
nhelix
+ integer, dimension(:,:),allocatable :: icont
+ real(kind=4),dimension(:),allocatable :: dcont
+ if (.not.allocated(icont)) then
+ allocate(icont(2,100*nres_molec(1)+1))
+ endif
+ if (.not.allocated(dcont)) then
+ allocate(dcont(100*nres_molec(1)+1))
+ endif
call chainbuild_cart
!d call write_pdb(99,'sec structure',0d0)
ncont=0
(xpi(2)-xpj(2))*(xpi(2)-xpj(2)) + &
(xpi(3)-xpj(3))*(xpi(3)-xpj(3))
if ( d.lt.rcomp*rcomp) then
+ if (ncont.gt.(100*nres_molec(1)+1)) ncont=100*nres_molec(1)+1
ncont=ncont+1
icont(1,ncont)=i
icont(2,ncont)=j
! print *,i,j,gg_lipi(3),gg_lipj(3),sss_ele_cut
! write (iout,*) "gg",(gg(k),k=1,3)
do k=1,3
- gradpepcatx(k,i)=gradpepcatx(k,i)-gg(k) &
+ gradpepcatx(k,i)=gradpepcatx(k,i)-gg(k)*sss_ele_cut &
+(eom12*(dc_norm(k,j)-om12*dc_norm(k,nres+i)) &
- +eom1*(erij(k)-om1*dc_norm(k,nres+i)))*dsci_inv
+ +eom1*(erij(k)-om1*dc_norm(k,nres+i)))*dsci_inv*sss_ele_cut
! gradpepcatx(k,j)=gradpepcatx(k,j)+gg(k) &
! +(eom12*(dc_norm(k,nres+i)-om12*dc_norm(k,j)) &
! Calculate the components of the gradient in DC and X
!
do l=1,3
- gradpepcat(l,i)=gradpepcat(l,i)-gg(l)
- gradpepcat(l,j)=gradpepcat(l,j)+gg(l)
+ gradpepcat(l,i)=gradpepcat(l,i)-gg(l)*sss_ele_cut
+ gradpepcat(l,j)=gradpepcat(l,j)+gg(l)*sss_ele_cut
enddo
end subroutine sc_grad_cat
! eom2=0.0d0
! eom12=evdwij*eps1_om12
! end diagnostics
+! write (iout,*) "gg",(gg(k),k=1,3)
do k=1,3
dcosom1(k) = rij * (dc_norm(k,i) - om1 * erij(k))
dcosom2(k) = rij * (dc_norm(k,nres+j) - om2 * erij(k))
gg(k) = gg(k) + eom1 * dcosom1(k) + eom2 * dcosom2(k)
- gvdwc_pepbase(k,i)= gvdwc_pepbase(k,i) +0.5*(- gg(k)) &
+ gradpepcat(k,i)= gradpepcat(k,i) +sss_ele_cut*(0.5*(- gg(k)) &
+ (-eom12*(dc_norm(k,nres+j)-om12*dc_norm(k,i)))&
*dsci_inv*2.0 &
- - (eom1*(erij(k)-om1*dc_norm(k,i)))*dsci_inv*2.0
- gvdwc_pepbase(k,i+1)= gvdwc_pepbase(k,i+1) +0.5*(- gg(k)) &
+ - (eom1*(erij(k)-om1*dc_norm(k,i)))*dsci_inv*2.0)
+ gradpepcat(k,i+1)= gradpepcat(k,i+1) +sss_ele_cut*(0.5*(- gg(k)) &
- (-eom12*(dc_norm(k,nres+j)-om12*dc_norm(k,i))) &
*dsci_inv*2.0 &
- + (eom1*(erij(k)-om1*dc_norm(k,i)))*dsci_inv*2.0
- gradpepcat(k,j)=gradpepcat(k,j)+gg(k)
+ + (eom1*(erij(k)-om1*dc_norm(k,i)))*dsci_inv*2.0)
+ gradpepcat(k,j)=gradpepcat(k,j)+gg(k)*sss_ele_cut
enddo
end subroutine sc_grad_cat_pep
real(kind=8) :: xj_safe,yj_safe,zj_safe,xj_temp,yj_temp,zj_temp,&
dist_temp, dist_init,ssgradlipi,ssgradlipj, &
sslipi,sslipj,faclip,alpha_sco
- integer :: ii
+ integer :: ii,ki
real(kind=8) :: fracinbuf
real (kind=8) :: escpho
real (kind=8),dimension(4):: ener
enddo
! go to 17
! do i=1,nres_molec(1)-1 ! loop over all peptide groups needs parralelization
- do i=ibond_start,ibond_end
+! do i=ibond_start,ibond_end
+ do ki=g_listcatscnorm_start,g_listcatscnorm_end
+ i=newcontlistcatscnormi(ki)
+ j=newcontlistcatscnormj(ki)
! print *,"I am in EVDW",i
itypi=iabs(itype(i,1))
dyi=dc_norm(2,nres+i)
dzi=dc_norm(3,nres+i)
dsci_inv=vbld_inv(i+nres)
- do j=itmp+1,itmp+nres_molec(5)
+! do j=itmp+1,itmp+nres_molec(5)
! Calculate SC interaction energy.
itypj=iabs(itype(j,5))
zj=boxshift(zj-zi,boxzsize)
! write(iout,*) "xj,yj,zj", xj,yj,zj,boxxsize
+ dxj=0.0
+ dyj=0.0
+ dzj=0.0
! dxj = dc_norm( 1, nres+j )
! dyj = dc_norm( 2, nres+j )
! dzj = dc_norm( 3, nres+j )
! rij holds 1/(distance of Calpha atoms)
rrij = 1.0D0 / ( xj*xj + yj*yj + zj*zj)
rij = dsqrt(rrij)
+ sss_ele_cut=sscale_ele(1.0d0/(rij))
+ sss_ele_grad=sscagrad_ele(1.0d0/(rij))
+! print *,sss_ele_cut,sss_ele_grad,&
+! 1.0d0/(rij),r_cut_ele,rlamb_ele
+ if (sss_ele_cut.le.0.0) cycle
CALL sc_angular
! this should be in elgrad_init but om's are calculated by sc_angular
! which in turn is used by older potentials
! END IF
!#else
evdw = evdw &
- + evdwij
+ + evdwij*sss_ele_cut
!#endif
c1 = c1 * eps1 * eps2rt**2 * eps3rt**2
fac = -expon * (c1 + evdwij) * rij_shift
sigder = fac * sigder
! Calculate distance derivative
- gg(1) = fac
- gg(2) = fac
- gg(3) = fac
+ gg(1) = fac*sss_ele_cut+evdwij*sss_ele_grad
+ gg(2) = fac*sss_ele_cut+evdwij*sss_ele_grad
+ gg(3) = fac*sss_ele_cut+evdwij*sss_ele_grad
! print *,"GG(1),distance grad",gg(1)
fac = chis1 * sqom1 + chis2 * sqom2 &
- 2.0d0 * chis12 * om1 * om2 * om12
dtop = b1cav * ((Lambf / (2.0d0 * eagle)) + (b2cav * Lambf))
dbot = 12.0d0 * b4cav * bat * Lambf
- dFdR = ((dtop * bot - top * dbot) / botsq) * sparrow
-
+ dFdR = ((dtop * bot - top * dbot) / botsq) * sparrow*sss_ele_cut+&
+ Fcav*sss_ele_grad
+ Fcav=Fcav*sss_ele_cut
dtop = b1cav * ((Chif / (2.0d0 * eagle)) + (b2cav * Chif))
dbot = 12.0d0 * b4cav * bat * Chif
eagle = Lambf * pom
isel = iabs(Qi) + 1 ! ion is always charged so iabs(Qj)
! print *,i,itype(i,1),isel
IF (isel.eq.0) THEN
-!c! No charges - do nothing
eheadtail = 0.0d0
-
ELSE IF (isel.eq.1) THEN
-!c! Nonpolar-charge interactions
if ((itype(i,1).eq.27).or.(itype(i,1).eq.26).or.(itype(i,1).eq.25)) then
Qi=Qi*2
Qij=Qij*2
endif
-
CALL enq_cat(epol)
eheadtail = epol
-! eheadtail = 0.0d0
-
ELSE IF (isel.eq.3) THEN
-!c! Dipole-charge interactions
if ((itype(i,1).eq.27).or.(itype(i,1).eq.26).or.(itype(i,1).eq.25)) then
Qi=Qi*2
Qij=Qij*2
endif
-! write(iout,*) "KURWA0",d1
-
CALL edq_cat(ecl, elj, epol)
eheadtail = ECL + elj + epol
-! eheadtail = 0.0d0
-
ELSE IF ((isel.eq.2)) THEN
-
-!c! Same charge-charge interaction ( +/+ or -/- )
if ((itype(i,1).eq.27).or.(itype(i,1).eq.26).or.(itype(i,1).eq.25)) then
Qi=Qi*2
Qij=Qij*2
endif
-
CALL eqq_cat(Ecl,Egb,Epol,Fisocav,Elj)
eheadtail = ECL + Egb + Epol + Fisocav + Elj
-! eheadtail = 0.0d0
-
-! ELSE IF ((isel.eq.2.and. &
-! iabs(Qi).eq.1).and. &
-! nstate(itypi,itypj).ne.1) THEN
-!c! Different charge-charge interaction ( +/- or -/+ )
-! if ((itype(i,1).eq.27).or.(itype(i,1).eq.26).or.(itype(i,1).eq.25)) then
-! Qi=Qi*2
-! Qij=Qij*2
-! endif
-! if ((itype(j,1).eq.27).or.(itype(j,1).eq.26).or.(itype(j,1).eq.25)) then
-! Qj=Qj*2
-! Qij=Qij*2
-! endif
-!
-! CALL energy_quad(istate,eheadtail,Ecl,Egb,Epol,Fisocav,Elj,Equad)
END IF ! this endif ends the "catch the gly-gly" at the beggining of Fcav
- else
- write(iout,*) "not yet implemented",j,itype(j,5)
+ else ! here is water and other molecules
+ isel = iabs(Qi)+2
+! isel=2
+! if (isel.eq.4) isel=2
+ if (isel.eq.2) then
+ eheadtail = 0.0d0
+ else if (isel.eq.3) then
+ if ((itype(i,1).eq.27).or.(itype(i,1).eq.26).or.(itype(i,1).eq.25)) then
+ Qi=Qi*2
+ Qij=Qij*2
+ endif
+ call eqd_cat(ecl,elj,epol)
+ eheadtail = ECL + elj + epol
+ else if (isel.eq.4) then
+ call edd_cat(ecl)
+ eheadtail = ECL
+ endif
+! write(iout,*) "not yet implemented",j,itype(j,5)
endif
!! endif ! turn off electrostatic
evdw = evdw + Fcav + eheadtail
!c!-------------------------------------------------------------------
!c! NAPISY KONCOWE
END DO ! j
- END DO ! i
+! END DO ! i
!c write (iout,*) "Number of loop steps in EGB:",ind
!c energy_dec=.false.
! print *,"EVDW KURW",evdw,nres
!!! return
17 continue
! go to 23
- do i=ibond_start,ibond_end
+! do i=ibond_start,ibond_end
+
+ do ki=g_listcatpnorm_start,g_listcatpnorm_end
+ i=newcontlistcatpnormi(ki)
+ j=newcontlistcatpnormj(ki)
! print *,"I am in EVDW",i
itypi=10 ! the peptide group parameters are for glicine
dyi=dc_norm(2,i)
dzi=dc_norm(3,i)
dsci_inv=vbld_inv(i+1)/2.0
- do j=itmp+1,itmp+nres_molec(5)
+! do j=itmp+1,itmp+nres_molec(5)
! Calculate SC interaction energy.
itypj=iabs(itype(j,5))
dCAVdOM1 = 0.0d0
dCAVdOM2 = 0.0d0
dCAVdOM12 = 0.0d0
- dscj_inv = vbld_inv(j+nres)
+ dscj_inv = 0.0d0 ! vbld_inv(j+nres)
! print *,i,j,dscj_inv,dsci_inv
! rij holds 1/(distance of Calpha atoms)
rrij = 1.0D0 / ( xj*xj + yj*yj + zj*zj)
rij = dsqrt(rrij)
+ sss_ele_cut=sscale_ele(1.0d0/(rij))
+ sss_ele_grad=sscagrad_ele(1.0d0/(rij))
+! print *,sss_ele_cut,sss_ele_grad,&
+! 1.0d0/(rij),r_cut_ele,rlamb_ele
+ if (sss_ele_cut.le.0.0) cycle
CALL sc_angular
! this should be in elgrad_init but om's are calculated by sc_angular
! which in turn is used by older potentials
! om = omega, sqom = om^2
+ om2=0.0d0
+ om12=0.0d0
sqom1 = om1 * om1
sqom2 = om2 * om2
sqom12 = om12 * om12
! END IF
!#else
evdw = evdw &
- + evdwij
+ + evdwij*sss_ele_cut
!#endif
c1 = c1 * eps1 * eps2rt**2 * eps3rt**2
fac = -expon * (c1 + evdwij) * rij_shift
sigder = fac * sigder
! Calculate distance derivative
- gg(1) = fac
- gg(2) = fac
- gg(3) = fac
+ gg(1) = fac*sss_ele_cut+evdwij*sss_ele_grad
+ gg(2) = fac*sss_ele_cut+evdwij*sss_ele_grad
+ gg(3) = fac*sss_ele_cut+evdwij*sss_ele_grad
fac = chis1 * sqom1 + chis2 * sqom2 &
- 2.0d0 * chis12 * om1 * om2 * om12
dtop = b1cav * ((Lambf / (2.0d0 * eagle)) + (b2cav * Lambf))
dbot = 12.0d0 * b4cav * bat * Lambf
- dFdR = ((dtop * bot - top * dbot) / botsq) * sparrow
-
+ dFdR = ((dtop * bot - top * dbot) / botsq) * sparrow*sss_ele_cut+&
+ Fcav*sss_ele_grad
+ Fcav=Fcav*sss_ele_cut
dtop = b1cav * ((Chif / (2.0d0 * eagle)) + (b2cav * Chif))
dbot = 12.0d0 * b4cav * bat * Chif
eagle = Lambf * pom
dFdOM1 = -(chis1 * om1 - chis12 * om2 * om12) / (eagle)
+
dFdOM2 = -(chis2 * om2 - chis12 * om1 * om12) / (eagle)
dFdOM12 = chis12 * (chis1 * om1 * om12 - om2) &
* (chis2 * om2 * om12 - om1) / (eagle * pom)
dFdL = ((dtop * bot - top * dbot) / botsq)
dCAVdOM1 = dFdL * ( dFdOM1 )
- dCAVdOM2 = dFdL * ( dFdOM2 )
- dCAVdOM12 = dFdL * ( dFdOM12 )
+! dCAVdOM2 = dFdL * ( dFdOM2 )
+! dCAVdOM12 = dFdL * ( dFdOM12 )
+ dCAVdOM2=0.0d0
+ dCAVdOM12=0.0d0
DO k= 1, 3
ertail(k) = Rtail_distance(k)/Rtail
! eheadtail = 0.0d0
else
!HERE WATER and other types of molecules solvents will be added
- write(iout,*) "not yet implemented"
+! write(iout,*) "not yet implemented"
+ CALL edd_cat_pep(ecl)
+ eheadtail=ecl
! CALL edd_cat_pep
+! eheadtail=0.0d0
endif
evdw = evdw + Fcav + eheadtail
! if (evdw.gt.1.0d6) then
!c!-------------------------------------------------------------------
!c! NAPISY KONCOWE
END DO ! j
- END DO ! i
+! END DO ! i
!c write (iout,*) "Number of loop steps in EGB:",ind
!c energy_dec=.false.
! print *,"EVDW KURW",evdw,nres
! do i=1,nres_molec(1)-1 ! loop over all peptide groups needs parralelization
do i=ibond_start,ibond_end
! cycle
+
if ((itype(i,1).eq.ntyp1).or.(itype(i+1,1).eq.ntyp1)) cycle ! leave dummy atoms
xi=0.5d0*(c(1,i)+c(1,i+1))
yi=0.5d0*(c(2,i)+c(2,i+1))
real(kind=8) :: xj_safe,yj_safe,zj_safe,xj_temp,yj_temp,zj_temp,&
dist_temp, dist_init,ssgradlipi,ssgradlipj, &
sslipi,sslipj,faclip,alpha_sco
- integer :: ii
+ integer :: ii,icont
real(kind=8) :: fracinbuf
real (kind=8) :: escpho
real (kind=8),dimension(4):: ener
sss_ele_cut=1.0d0
countss=0
! print *,"EVDW KURW",evdw,nres
- do i=iatsc_s,iatsc_e
+! do i=iatsc_s,iatsc_e
! print *,"I am in EVDW",i
+ do icont=g_listscsc_start,g_listscsc_end
+ i=newcontlisti(icont)
+ j=newcontlistj(icont)
+
itypi=iabs(itype(i,1))
! if (i.ne.47) cycle
if (itypi.eq.ntyp1) cycle
!
! Calculate SC interaction energy.
!
- do iint=1,nint_gr(i)
- do j=istart(i,iint),iend(i,iint)
+! do iint=1,nint_gr(i)
+! do j=istart(i,iint),iend(i,iint)
! print *,"JA PIER",i,j,iint,istart(i,iint),iend(i,iint)
IF (dyn_ss_mask(i).and.dyn_ss_mask(j)) THEN
call dyn_ssbond_ene(i,j,evdwij,countss)
! rij holds 1/(distance of Calpha atoms)
rrij = 1.0D0 / ( xj*xj + yj*yj + zj*zj)
rij = dsqrt(rrij)
+ sss_ele_cut=sscale_ele(1.0d0/(rij))
+ sss_ele_grad=sscagrad_ele(1.0d0/(rij))
+! print *,sss_ele_cut,sss_ele_grad,&
+! 1.0d0/(rij),r_cut_ele,rlamb_ele
+ if (sss_ele_cut.le.0.0) cycle
+
!----------------------------
CALL sc_angular
! this should be in elgrad_init but om's are calculated by sc_angular
! END IF
!#else
evdw = evdw &
- + evdwij
+ + evdwij*sss_ele_cut
!#endif
c1 = c1 * eps1 * eps2rt**2 * eps3rt**2
sigder = fac * sigder
! fac = rij * fac
! Calculate distance derivative
- gg(1) = fac
- gg(2) = fac
- gg(3) = fac
+ gg(1) = fac*sss_ele_cut+evdwij*sss_ele_grad
+ gg(2) = fac*sss_ele_cut+evdwij*sss_ele_grad
+ gg(3) = fac*sss_ele_cut+evdwij*sss_ele_grad
! if (b2.gt.0.0) then
fac = chis1 * sqom1 + chis2 * sqom2 &
- 2.0d0 * chis12 * om1 * om2 * om12
dtop = b1cav * ((Lambf / (2.0d0 * eagle)) + (b2cav * Lambf))
dbot = 12.0d0 * b4cav * bat * Lambf
- dFdR = ((dtop * bot - top * dbot) / botsq) * sparrow
-
+ dFdR = ((dtop * bot - top * dbot) / botsq) * sparrow*sss_ele_cut&
+ +Fcav*sss_ele_grad
+ Fcav=Fcav*sss
dtop = b1cav * ((Chif / (2.0d0 * eagle)) + (b2cav * Chif))
dbot = 12.0d0 * b4cav * bat * Chif
eagle = Lambf * pom
END IF
!c!-------------------------------------------------------------------
!c! NAPISY KONCOWE
- END DO ! j
- END DO ! iint
+ ! END DO ! j
+ !END DO ! iint
END DO ! i
!c write (iout,*) "Number of loop steps in EGB:",ind
!c energy_dec=.false.
!c! Coulomb electrostatic interaction
Ecl = (332.0d0 * Qij) / Rhead
!c! derivative of Ecl is Gcl...
- dGCLdR = (-332.0d0 * Qij ) / Rhead_sq
+ dGCLdR = (-332.0d0 * Qij ) / Rhead_sq*sss_ele_cut+ECL*sss_ele_grad
dGCLdOM1 = 0.0d0
dGCLdOM2 = 0.0d0
dGCLdOM12 = 0.0d0
+ ECL=ECL*sss_ele_grad
ee0 = dexp(-( Rhead_sq ) / (4.0d0 * a12sq))
Fgb = sqrt( ( Rhead_sq ) + a12sq * ee0)
debkap=debaykap(itypi,itypj)
-(332.0d0 * Qij *&
(dexp(-debkap*Fgb)*debkap/eps_out))/ Fgb
dFGBdR = ( Rhead * ( 2.0d0 - (0.5d0 * ee0) ) )/ ( 2.0d0 * Fgb )
- dGGBdR = dGGBdFGB * dFGBdR
+ dGGBdR = dGGBdFGB * dFGBdR*sss_ele_cut+Egb*sss_ele_grad
+ Egb=Egb*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Fisocav - isotropic cavity creation term
!c! or "how much energy it costs to put charged head in water"
* ( 2.0d0 - 0.5d0 * ee1) ) / ( 2.0d0 * fgb1 )
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2))&
* ( 2.0d0 - 0.5d0 * ee2) ) / ( 2.0d0 * fgb2 )
- dPOLdR1 = dPOLdFGB1 * dFGBdR1
+ dPOLdR1 = dPOLdFGB1 * dFGBdR1*sss_ele_cut
!c! dPOLdR1 = 0.0d0
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
Elj = 4.0d0 * eps_head * pom * (pom-1.0d0)
!c! derivative of Elj is Glj
dGLJdR = 4.0d0 * eps_head*(((-12.0d0*pis**12.0d0)/(Rhead**13.0d0))&
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))*sss_ele_cut+&
+ (ELJ+epol)*sss_ele_grad
+ epol=epol*sss_ele_cut
+ Elj=Elj*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Return the results
!c! These things do the dRdX derivatives, that is
!c! Coulomb electrostatic interaction
Ecl = (332.0d0 * Qij) / Rhead
!c! derivative of Ecl is Gcl...
- dGCLdR = (-332.0d0 * Qij ) / Rhead_sq
+ dGCLdR = (-332.0d0 * Qij ) / Rhead_sq*sss_ele_cut+ECL*sss_ele_grad
+ ECL=ECL*sss_ele_cut
dGCLdOM1 = 0.0d0
dGCLdOM2 = 0.0d0
dGCLdOM12 = 0.0d0
-(332.0d0 * Qij *&
(dexp(-debkap*Fgb)*debkap/eps_out))/ Fgb
dFGBdR = ( Rhead * ( 2.0d0 - (0.5d0 * ee0) ) )/ ( 2.0d0 * Fgb )
- dGGBdR = dGGBdFGB * dFGBdR
+ dGGBdR = dGGBdFGB * dFGBdR*sss_ele_cut+Egb*sss_ele_grad
+ Egb=Egb*sss_ele_grad
!c!-------------------------------------------------------------------
!c! Fisocav - isotropic cavity creation term
!c! or "how much energy it costs to put charged head in water"
!c! Derivative of Fisocav is GCV...
dtop = al1 * ((1.0d0 / (2.0d0 * dsqrt(pom))) + al2)
dbot = 12.0d0 * al4 * pom ** 11.0d0
- dGCVdR = ((dtop * bot - top * dbot) / botsq) * csig
+ dGCVdR = ((dtop * bot - top * dbot) / botsq) * csig*sss_ele_cut&
+ +FisoCav*sss_ele_grad
+ FisoCav=FisoCav*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Epol
!c! Polarization energy - charged heads polarize hydrophobic "neck"
* ( 2.0d0 - 0.5d0 * ee1) ) / ( 2.0d0 * fgb1 )
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2))&
* ( 2.0d0 - 0.5d0 * ee2) ) / ( 2.0d0 * fgb2 )
- dPOLdR1 = dPOLdFGB1 * dFGBdR1
+ dPOLdR1 = dPOLdFGB1 * dFGBdR1!*sss_ele_cut+epol*sss_ele_grad
!c! dPOLdR1 = 0.0d0
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2!*sss_ele_cut+epol*sss_ele_grad
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
dPOLdOM2 = dPOLdFGB1 * dFGBdOM2
+! epol=epol*sss_ele_cut
!c! dPOLdOM2 = 0.0d0
!c!-------------------------------------------------------------------
!c! Elj
Elj = 4.0d0 * eps_head * pom * (pom-1.0d0)
!c! derivative of Elj is Glj
dGLJdR = 4.0d0 * eps_head*(((-12.0d0*pis**12.0d0)/(Rhead**13.0d0))&
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))*sss_ele_cut&
+ +(Elj+epol)*sss_ele_grad
+ Elj=Elj*sss_ele_cut
+ epol=epol*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Return the results
!c! These things do the dRdX derivatives, that is
dcosom2(k) = rij * (dc_norm(k,nres+j) - om2 * erij(k))
gg(k) = gg(k) + eom1 * dcosom1(k) + eom2 * dcosom2(k)
!c! this acts on hydrophobic center of interaction
- gvdwx(k,i)= gvdwx(k,i) - gg(k) &
+ gvdwx(k,i)= gvdwx(k,i) - gg(k)*sss_ele_cut &
+ (eom12*(dc_norm(k,nres+j)-om12*dc_norm(k,nres+i))&
- + eom1*(erij(k)-om1*dc_norm(k,nres+i)))*dsci_inv
- gvdwx(k,j)= gvdwx(k,j) + gg(k) &
+ + eom1*(erij(k)-om1*dc_norm(k,nres+i)))*dsci_inv*sss_ele_cut
+ gvdwx(k,j)= gvdwx(k,j) + gg(k)*sss_ele_cut &
+ (eom12*(dc_norm(k,nres+i)-om12*dc_norm(k,nres+j))&
- + eom2*(erij(k)-om2*dc_norm(k,nres+j)))*dscj_inv
+ + eom2*(erij(k)-om2*dc_norm(k,nres+j)))*dscj_inv*sss_ele_cut
!c! this acts on Calpha
- gvdwc(k,i)=gvdwc(k,i)-gg(k)
- gvdwc(k,j)=gvdwc(k,j)+gg(k)
+ gvdwc(k,i)=gvdwc(k,i)-gg(k)*sss_ele_cut
+ gvdwc(k,j)=gvdwc(k,j)+gg(k)*sss_ele_cut
END DO
!c! sc_grad is done, now we will compute
eheadtail = 0.0d0
!c! Ecl = 0.0d0
!c! write (*,*) "Ecl = ", Ecl
!c! derivative of Ecl is Gcl...
- dGCLdR = (-332.0d0 * Qij ) / (Rhead_sq * eps_in)
+ dGCLdR = (-332.0d0 * Qij ) / (Rhead_sq * eps_in)*sss_ele_cut+ECL*sss_ele_grad
+ ECL=ecl*sss_ele_cut
!c! dGCLdR = 0.0d0
dGCLdOM1 = 0.0d0
dGCLdOM2 = 0.0d0
dGGBdFGB = -(-332.0d0 * Qij * eps_inout_fac) / (Fgb * Fgb)
dFGBdR = ( Rhead * ( 2.0d0 - (0.5d0 * ee0) ) )&
/ ( 2.0d0 * Fgb )
- dGGBdR = dGGBdFGB * dFGBdR
+ dGGBdR = dGGBdFGB * dFGBdR*sss_ele_cut+Egb*sss_ele_grad
+ Egb=Egb*sss_ele_cut
!c! dGGBdR = 0.0d0
!c!-------------------------------------------------------------------
!c! Fisocav - isotropic cavity creation term
FisoCav = top / bot
dtop = al1 * ((1.0d0 / (2.0d0 * dsqrt(pom))) + al2)
dbot = 12.0d0 * al4 * pom ** 11.0d0
- dGCVdR = ((dtop * bot - top * dbot) / botsq) * csig
+ dGCVdR = ((dtop * bot - top * dbot) / botsq) * csig*sss_ele_cut+FisoCav*sss_ele_grad
+ FisoCav=FisoCav*sss_ele_cut
+
!c! dGCVdR = 0.0d0
!c!-------------------------------------------------------------------
!c! Polarization energy
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2)) &
* ( 2.0d0 - 0.5d0 * ee2) ) &
/ ( 2.0d0 * fgb2 )
- dPOLdR1 = dPOLdFGB1 * dFGBdR1
+ dPOLdR1 = dPOLdFGB1 * dFGBdR1*sss_ele_cut
!c! dPOLdR1 = 0.0d0
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
!c! derivative of Elj is Glj
dGLJdR = 4.0d0 * eps_head &
* (((-12.0d0*pis**12.0d0)/(Rhead**13.0d0)) &
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))*sss_ele_cut+&
+ (epol+Elj)*sss_ele_grad
+ Elj=Elj*sss_ele_cut
+ epol=epol*sss_ele_cut
!c! dGLJdR = 0.0d0
!c!-------------------------------------------------------------------
!c! Equad
Equad = fac * Beta1
!c! Equad = 0.0d0
!c! derivative of Equad...
- dQUADdR = ((2.5d0 * Wqd * Beta1) / (Fgb**6.0d0)) * dFGBdR
+ dQUADdR = ((2.5d0 * Wqd * Beta1) / (Fgb**6.0d0)) * dFGBdR*sss_ele_cut&
+ + Equad*sss_ele_grad
+ Equad=Equad*sss_ele_cut
!c! dQUADdR = 0.0d0
dQUADdOM1 = fac* (-75.0d0*om1 + 315.0d0*om1*sqom2 - 45.0d0*om2*om12)
!c! dQUADdOM1 = 0.0d0
DO k = 1, 3
dcosom1(k) = rij * (dc_norm(k,nres+i) - om1 * erij(k))
dcosom2(k) = rij * (dc_norm(k,nres+j) - om2 * erij(k))
- tuna(k) = eom1 * dcosom1(k) + eom2 * dcosom2(k)
+ tuna(k) = eom1 * dcosom1(k) + eom2 * dcosom2(k)*sss_ele_cut
END DO
!c! Radial stuff
DO k = 1, 3
- dQUADdR * pom&
- tuna(k) &
+ (eom12*(dc_norm(k,nres+j)-om12*dc_norm(k,nres+i))&
- + eom1*(erij(k)-om1*dc_norm(k,nres+i)))*dsci_inv
+ + eom1*(erij(k)-om1*dc_norm(k,nres+i)))*dsci_inv*sss_ele_cut
pom = erhead(k)+facd2*(erhead(k)-erdxj*dC_norm(k,j+nres))
!c! this acts on hydrophobic center of interaction
+ dQUADdR * pom &
+ tuna(k) &
+ (eom12*(dc_norm(k,nres+i)-om12*dc_norm(k,nres+j)) &
- + eom2*(erij(k)-om2*dc_norm(k,nres+j)))*dscj_inv
+ + eom2*(erij(k)-om2*dc_norm(k,nres+j)))*dscj_inv*sss_ele_cut
!c! this acts on Calpha
gheadtail(k,3,1) = gheadtail(k,3,1) &
dFGBdOM2 = (((R1 * R1 * chi1 * om2) / (MomoFac1 * MomoFac1)) &
* (2.0d0 - 0.5d0 * ee1) ) &
/ (2.0d0 * fgb1)
- dPOLdR1 = dPOLdFGB1 * dFGBdR1
+ dPOLdR1 = dPOLdFGB1 * dFGBdR1*sss_ele_cut+epol*sss_ele_grad
+ epol=epol*sss_ele_cut
!c! dPOLdR1 = 0.0d0
dPOLdOM1 = 0.0d0
dPOLdOM2 = dPOLdFGB1 * dFGBdOM2
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2)) &
* (2.0d0 - 0.5d0 * ee2) ) &
/ (2.0d0 * fgb2)
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut+epol*sss_ele_grad
+ epol=epol*sss_ele_cut
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2)) &
* (2.0d0 - 0.5d0 * ee2) ) &
/ (2.0d0 * fgb2)
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut+epol*sss_ele_grad
+ epol=epol*sss_ele_cut
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
hawk = w2 * Qi * Qi * (1.0d0 - sqom2)
Ecl = sparrow / Rhead**2.0d0 &
- hawk / Rhead**4.0d0
- dGCLdR = - 2.0d0 * sparrow / Rhead**3.0d0 &
- + 4.0d0 * hawk / Rhead**5.0d0
+ dGCLdR = (- 2.0d0 * sparrow / Rhead**3.0d0 &
+ + 4.0d0 * hawk / Rhead**5.0d0)*sss_ele_cut+Ecl*sss_ele_grad
+ Ecl=Ecl*sss_ele_cut
!c! dF/dom1
dGCLdOM1 = (w1 * Qi) / (Rhead**2.0d0)
!c! dF/dom2
dFGBdOM2 = (((R1 * R1 * chi1 * om2) / (MomoFac1 * MomoFac1)) &
* (2.0d0 - 0.5d0 * ee1) ) &
/ (2.0d0 * fgb1)
- dPOLdR1 = dPOLdFGB1 * dFGBdR1
+ dPOLdR1 = dPOLdFGB1 * dFGBdR1*sss_ele_cut+epol*sss_ele_grad
!c! dPOLdR1 = 0.0d0
dPOLdOM1 = 0.0d0
dPOLdOM2 = dPOLdFGB1 * dFGBdOM2
!c! derivative of Elj is Glj
dGLJdR = 4.0d0 * eps_head &
* (((-12.0d0*pis**12.0d0)/(Rhead**13.0d0)) &
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))*sss_ele_cut+elj*sss_ele_grad
+ Elj=Elj*sss_ele_cut
DO k = 1, 3
erhead(k) = Rhead_distance(k)/Rhead
erhead_tail(k,1) = ((ctail(k,2)-chead(k,1))/R1)
END DO
RETURN
END SUBROUTINE eqd
+
+ SUBROUTINE eqd_cat(Ecl,Elj,Epol)
+ use calc_data
+ use comm_momo
+ double precision facd4, federmaus,ecl,elj,epol
+ alphapol1 = alphapolcat(itypi,itypj)
+ w1 = wqdipcat(1,itypi,itypj)
+ w2 = wqdipcat(2,itypi,itypj)
+ pis = sig0headcat(itypi,itypj)
+ eps_head = epsheadcat(itypi,itypj)
+! eps_head=0.0d0
+! w2=0.0d0
+! alphapol1=0.0d0
+!c!-------------------------------------------------------------------
+!c! R1 - distance between head of ith side chain and tail of jth sidechain
+ R1 = 0.0d0
+ DO k = 1, 3
+!c! Calculate head-to-tail distances
+ R1=R1+(ctail(k,2)-chead(k,1))**2
+ END DO
+!c! Pitagoras
+ R1 = dsqrt(R1)
+
+!c! R1 = dsqrt((Rtail**2)+((dtail(1,itypi,itypj)
+!c! & +dhead(1,1,itypi,itypj))**2))
+!c! R2 = dsqrt((Rtail**2)+((dtail(2,itypi,itypj)
+!c! & +dhead(2,1,itypi,itypj))**2))
+
+!c!-------------------------------------------------------------------
+!c! ecl
+ sparrow = w1 * Qi * om1
+ hawk = w2 * Qi * Qi * (1.0d0 - sqom2)
+ Ecl = sparrow / Rhead**2.0d0 &
+ - hawk / Rhead**4.0d0
+ dGCLdR =sss_ele_cut*(-2.0d0 * sparrow / Rhead**3.0d0 &
+ + 4.0d0 * hawk / Rhead**5.0d0)+sss_ele_grad*ECL
+ ECL=ECL*sss_ele_cut
+!c! dF/dom1
+ dGCLdOM1 = (w1 * Qi) / (Rhead**2.0d0)
+!c! dF/dom2
+ dGCLdOM2 = 0.0d0 !
+
+!(2.0d0 * w2 * Qi * Qi * om2) / (Rhead ** 4.0d0)
+
+!c--------------------------------------------------------------------
+!c Polarization energy
+!c Epol
+ MomoFac1 = (1.0d0 - chi1 * sqom2)
+ RR1 = R1 * R1 / MomoFac1
+ ee1 = exp(-( RR1 / (4.0d0 * a12sq) ))
+ fgb1 = sqrt( RR1 + a12sq * ee1)
+ epol = 332.0d0 * eps_inout_fac * (( alphapol1 / fgb1 )**4.0d0)
+!c! epol = 0.0d0
+!c!------------------------------------------------------------------
+!c! derivative of Epol is Gpol...
+ dPOLdFGB1 = -(1328.0d0 * eps_inout_fac * alphapol1 ** 4.0d0) &
+ / (fgb1 ** 5.0d0)
+ dFGBdR1 = ( (R1 / MomoFac1) &
+ * ( 2.0d0 - (0.5d0 * ee1) ) ) &
+ / ( 2.0d0 * fgb1 )
+ dFGBdOM2 = 0.0d0 ! as om2 is 0
+! (((R1 * R1 * chi1 * om2) / (MomoFac1 * MomoFac1)) &
+! * (2.0d0 - 0.5d0 * ee1) ) &
+! / (2.0d0 * fgb1)
+ dPOLdR1 = dPOLdFGB1 * dFGBdR1*sss_ele_cut+epol*sss_ele_grad
+!c! dPOLdR1 = 0.0d0
+ dPOLdOM1 = 0.0d0
+! dPOLdOM2 = dPOLdFGB1 * dFGBdOM2
+ dPOLdOM2 = 0.0d0
+ epol=epol*sss_ele_cut
+!c!-------------------------------------------------------------------
+!c! Elj
+ pom = (pis / Rhead)**6.0d0
+ Elj = 4.0d0 * eps_head * pom * (pom-1.0d0)
+!c! derivative of Elj is Glj
+ dGLJdR = 4.0d0 * eps_head*sss_ele_cut &
+ * (((-12.0d0*pis**12.0d0)/(Rhead**13.0d0)) &
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))+Elj*sss_ele_grad
+ Elj=Elj*sss_ele_cut
+ DO k = 1, 3
+ erhead(k) = Rhead_distance(k)/Rhead
+ erhead_tail(k,1) = ((ctail(k,2)-chead(k,1))/R1)
+ END DO
+
+ erdxi = scalar( erhead(1), dC_norm(1,i+nres) )
+ bat = scalar( erhead_tail(1,1), dC_norm(1,i+nres) )
+ facd1 = d1 * vbld_inv(i+nres)
+
+ DO k = 1, 3
+ hawk = (erhead_tail(k,1) + &
+ facd1 * (erhead_tail(k,1) - bat * dC_norm(k,i+nres)))
+
+ pom = erhead(k)+facd1*(erhead(k)-erdxi*dC_norm(k,i+nres))
+ gradpepcatx(k,i) = gradpepcatx(k,i) &
+ - dGCLdR * pom&
+ - dPOLdR1 * hawk &
+ - dGLJdR * pom
+
+! pom = erhead(k)+facd2*(erhead(k)-erdxj*dC_norm(k,j+nres))
+! gradpepcatx(k,j) = gradpepcatx(k,j) &
+! + dGCLdR * pom &
+! + dPOLdR1 * (erhead_tail(k,1) &
+! -facd4 * (erhead_tail(k,1) - federmaus * dC_norm(k,j+nres))) &
+! + dGLJdR * pom
+
+
+ gradpepcat(k,i) = gradpepcat(k,i) &
+ - dGCLdR * erhead(k) &
+ - dPOLdR1 * erhead_tail(k,1) &
+ - dGLJdR * erhead(k)
+
+ gradpepcat(k,j) = gradpepcat(k,j) &
+ + dGCLdR * erhead(k) &
+ + dPOLdR1 * erhead_tail(k,1) &
+ + dGLJdR * erhead(k)
+
+ END DO
+ RETURN
+ END SUBROUTINE eqd_cat
+
SUBROUTINE edq(Ecl,Elj,Epol)
! IMPLICIT NONE
use comm_momo
!c!-------------------------------------------------------------------
!c! derivative of ecl is Gcl
!c! dF/dr part
- dGCLdR = - 2.0d0 * sparrow / Rhead**3.0d0 &
- + 4.0d0 * hawk / Rhead**5.0d0
+ dGCLdR =sss_ele_cut*(- 2.0d0 * sparrow / Rhead**3.0d0 &
+ + 4.0d0 * hawk / Rhead**5.0d0)+Ecl*sss_ele_grad
!c! dF/dom1
dGCLdOM1 = (w1 * Qj) / (Rhead**2.0d0)
!c! dF/dom2
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2)) &
* (2.0d0 - 0.5d0 * ee2) ) &
/ (2.0d0 * fgb2)
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut
+ epol=epol*sss_ele_cut
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
!c! derivative of Elj is Glj
dGLJdR = 4.0d0 * eps_head &
* (((-12.0d0*pis**12.0d0)/(Rhead**13.0d0)) &
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))*sss_ele_cut+Elj*sss_ele_grad
+ elj=elj*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Return the results
!c! (see comments in Eqq)
!c!-------------------------------------------------------------------
!c! derivative of ecl is Gcl
!c! dF/dr part
- dGCLdR = - 2.0d0 * sparrow / Rhead**3.0d0 &
- + 4.0d0 * hawk / Rhead**5.0d0
+ dGCLdR =( - 2.0d0 * sparrow / Rhead**3.0d0 &
+ + 4.0d0 * hawk / Rhead**5.0d0)*sss_ele_cut+ECL*sss_ele_grad
!c! dF/dom1
dGCLdOM1 = (w1 * Qj) / (Rhead**2.0d0)
!c! dF/dom2
dGCLdOM2 = (2.0d0 * w2 * Qj * Qj * om2) / (Rhead ** 4.0d0)
+ ECL=ECL*sss_ele_cut
!c--------------------------------------------------------------------
!c--------------------------------------------------------------------
!c Polarization energy
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2)) &
* (2.0d0 - 0.5d0 * ee2) ) &
/ (2.0d0 * fgb2)
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut+epol*sss_ele_grad
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
dPOLdOM2 = 0.0d0
+ epol=epol*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Elj
pom = (pis / Rhead)**6.0d0
!c! derivative of Elj is Glj
dGLJdR = 4.0d0 * eps_head &
* (((-12.0d0*pis**12.0d0)/(Rhead**13.0d0)) &
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))*sss_ele_cut+&
+ Elj*sss_ele_grad
+ Elj=Elj*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Return the results
!c!-------------------------------------------------------------------
!c! derivative of ecl is Gcl
!c! dF/dr part
- dGCLdR = - 2.0d0 * sparrow / Rhead**3.0d0 &
- + 4.0d0 * hawk / Rhead**5.0d0
+ dGCLdR = (- 2.0d0 * sparrow / Rhead**3.0d0 &
+ + 4.0d0 * hawk / Rhead**5.0d0)*sss_ele_cut+&
+ ECL*sss_ele_grad
+ ECL=ECL*sss_ele_cut
!c! dF/dom1
dGCLdOM1 = (w1 * Qj) / (Rhead**2.0d0)
!c! dF/dom2
dFGBdOM1 = (((R2 * R2 * chi2 * om1) / (MomoFac2 * MomoFac2)) &
* (2.0d0 - 0.5d0 * ee2) ) &
/ (2.0d0 * fgb2)
- dPOLdR2 = dPOLdFGB2 * dFGBdR2
+ dPOLdR2 = dPOLdFGB2 * dFGBdR2*sss_ele_cut+epol*sss_ele_grad
+ epol=epol*sss_ele_grad
!c! dPOLdR2 = 0.0d0
dPOLdOM1 = dPOLdFGB2 * dFGBdOM1
!c! dPOLdOM1 = 0.0d0
pom = (pis / Rhead)**6.0d0
Elj = 4.0d0 * eps_head * pom * (pom-1.0d0)
!c! derivative of Elj is Glj
- dGLJdR = 4.0d0 * eps_head &
+ dGLJdR = 4.0d0 * eps_head*sss_ele_cut &
* (((-12.0d0*pis**12.0d0)/(Rhead**13.0d0)) &
- + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))
+ + (( 6.0d0*pis**6.0d0) /(Rhead**7.0d0)))+Elj*sss_ele_grad
+ Elj=Elj*sss_ele_cut
!c!-------------------------------------------------------------------
!c! Return the results
c1 = (-3.0d0 * w1 * fac) / (Rhead ** 4.0d0)
c2 = (-6.0d0 * w2) / (Rhead ** 7.0d0) &
* (4.0d0 + fac * fac - 3.0d0 * (sqom1 + sqom2))
- dGCLdR = c1 - c2
+ dGCLdR = (c1 - c2)*sss_ele_cut+ECL*sss_ele_grad
+ ECL=ECL*sss_ele_cut
!c! dECL/dom1
c1 = (-3.0d0 * w1 * om2 ) / (Rhead**3.0d0)
c2 = (-6.0d0 * w2) / (Rhead**6.0d0) &
END DO
RETURN
END SUBROUTINE edd
+ SUBROUTINE edd_cat(ECL)
+! IMPLICIT NONE
+ use comm_momo
+ use calc_data
+
+ double precision ecl
+!c! csig = sigiso(itypi,itypj)
+ w1 = wqdipcat(1,itypi,itypj)
+ w2 = wqdipcat(2,itypi,itypj)
+! w2=0.0d0
+!c!-------------------------------------------------------------------
+!c! ECL
+! print *,"om1",om1,om2,om12
+ fac = - 3.0d0 * om1 !after integer and simplify
+ c1 = (w1 / (Rhead**3.0d0)) * fac
+ c2 = (w2 / Rhead ** 6.0d0) &
+ * (4.0d0 + 6.0d0*sqom1 ) !after integration and simplification
+ ECL = c1 - c2
+!c! dervative of ECL is GCL...
+!c! dECL/dr
+ c1 = (-3.0d0 * w1 * fac) / (Rhead ** 4.0d0)
+ c2 = (-6.0d0 * w2) / (Rhead ** 7.0d0) &
+ * (4.0d0 + 6.0d0*sqom1)
+ dGCLdR = (c1 - c2)*sss_ele_cut+ECL*sss_ele_grad
+!c! dECL/dom1
+ c1 = (-3.0d0 * w1) / (Rhead**3.0d0)
+ c2 = (12.0d0 * w2*om1) / (Rhead**6.0d0)
+ dGCLdOM1 = c1 - c2
+!c! dECL/dom2
+! c1 = (-3.0d0 * w1 * om1 ) / (Rhead**3.0d0)
+ c1=0.0 ! this is because om2 is 0
+! c2 = (-6.0d0 * w2) / (Rhead**6.0d0) &
+! * ( om1 * om12 - 3.0d0 * sqom1 * om2 + om2 )
+ c2=0.0 !om is 0
+ dGCLdOM2 = c1 - c2
+!c! dECL/dom12
+! c1 = w1 / (Rhead ** 3.0d0)
+ c1=0.0d0 ! this is because om12 is 0
+! c2 = ( 2.0d0 * w2 * fac ) / Rhead ** 6.0d0
+ c2=0.0d0 !om12 is 0
+ dGCLdOM12 = c1 - c2
+!c!-------------------------------------------------------------------
+!c! Return the results
+!c! (see comments in Eqq)
+ DO k= 1, 3
+ erhead(k) = Rhead_distance(k)/Rhead
+ END DO
+ erdxi = scalar( erhead(1), dC_norm(1,i+nres) )
+ erdxj = scalar( erhead(1), dC_norm(1,j+nres) )
+ facd1 = d1 * vbld_inv(i+nres)
+ facd2 = d2 * vbld_inv(j+nres)
+ DO k = 1, 3
+
+ pom = erhead(k)+facd1*(erhead(k)-erdxi*dC_norm(k,i+nres))
+ gradpepcatx(k,i) = gradpepcatx(k,i) - dGCLdR * pom
+! pom = erhead(k)+facd2*(erhead(k)-erdxj*dC_norm(k,j+nres))
+! gradpepcatx(k,j) = gradpepcatx(k,j) + dGCLdR * pom
+
+ gradpepcat(k,i) = gradpepcat(k,i) - dGCLdR * erhead(k)
+ gradpepcat(k,j) = gradpepcat(k,j) + dGCLdR * erhead(k)
+ END DO
+ RETURN
+ END SUBROUTINE edd_cat
+ SUBROUTINE edd_cat_pep(ECL)
+! IMPLICIT NONE
+ use comm_momo
+ use calc_data
+
+ double precision ecl
+!c! csig = sigiso(itypi,itypj)
+ w1 = wqdipcat(1,itypi,itypj)
+ w2 = wqdipcat(2,itypi,itypj)
+!c!-------------------------------------------------------------------
+!c! ECL
+ fac = (om12 - 3.0d0 * om1 * om2)
+ c1 = (w1 / (Rhead**3.0d0)) * fac
+ c2 = (w2 / Rhead ** 6.0d0) &
+ * (4.0d0 + fac * fac -3.0d0 * (sqom1 + sqom2))
+ ECL = c1 - c2
+!c! dECL/dr
+ c1 = (-3.0d0 * w1 * fac) / (Rhead ** 4.0d0)
+ c2 = (-6.0d0 * w2) / (Rhead ** 7.0d0) &
+ * (4.0d0 + fac * fac - 3.0d0 * (sqom1 + sqom2))
+ dGCLdR = (c1 - c2)*sss_ele_cut+ECL*sss_ele_grad
+ ECL=ECL*sss_ele_cut
+!c! dECL/dom1
+ c1 = (-3.0d0 * w1 * om2 ) / (Rhead**3.0d0)
+ c2 = (-6.0d0 * w2) / (Rhead**6.0d0) &
+ * ( om2 * om12 - 3.0d0 * om1 * sqom2 + om1 )
+ dGCLdOM1 = c1 - c2
+!c! dECL/dom2
+ c1 = (-3.0d0 * w1 * om1 ) / (Rhead**3.0d0)
+ c2 = (-6.0d0 * w2) / (Rhead**6.0d0) &
+ * ( om1 * om12 - 3.0d0 * sqom1 * om2 + om2 )
+ dGCLdOM2 = c1 - c2
+ dGCLdOM2=0.0d0 ! this is because om2=0
+!c! dECL/dom12
+ c1 = w1 / (Rhead ** 3.0d0)
+ c2 = ( 2.0d0 * w2 * fac ) / Rhead ** 6.0d0
+ dGCLdOM12 = c1 - c2
+ dGCLdOM12=0.0d0 !this is because om12=0.0
+!c!-------------------------------------------------------------------
+!c! Return the results
+!c! (see comments in Eqq)
+ DO k= 1, 3
+ erhead(k) = Rhead_distance(k)/Rhead
+ END DO
+ erdxi = scalar( erhead(1), dC_norm(1,i) )
+ erdxj = scalar( erhead(1), dC_norm(1,j+nres) )
+ facd1 = d1 * vbld_inv(i)
+ facd2 = d2 * vbld_inv(j+nres)
+ DO k = 1, 3
+
+ pom = facd1*(erhead(k)-erdxi*dC_norm(k,i))
+ gradpepcat(k,i) = gradpepcat(k,i) + dGCLdR * pom
+ gradpepcat(k,i+1) = gradpepcat(k,i+1) - dGCLdR * pom
+! pom = erhead(k)+facd2*(erhead(k)-erdxj*dC_norm(k,j+nres))
+! gradpepcatx(k,j) = gradpepcatx(k,j) + dGCLdR * pom
+
+ gradpepcat(k,i) = gradpepcat(k,i) - dGCLdR * erhead(k)*0.5d0
+ gradpepcat(k,i+1) = gradpepcat(k,i+1)- dGCLdR * erhead(k)*0.5d0
+ gradpepcat(k,j) = gradpepcat(k,j) + dGCLdR * erhead(k)
+ END DO
+ RETURN
+ END SUBROUTINE edd_cat_pep
+
SUBROUTINE elgrad_init(eheadtail,Egb,Ecl,Elj,Equad,Epol)
! IMPLICIT NONE
use comm_momo
if (itype(j,5).le.5) then
ilist_catscnorm=ilist_catscnorm+1
! this can be substituted by cantor and anti-cantor
+! write(iout,*) "have contact",i,j,ilist_catscnorm
contlistcatscnormi(ilist_catscnorm)=i
contlistcatscnormj(ilist_catscnorm)=j
+! write(iout,*) "have contact2",i,j,ilist_catscnorm,&
+! contlistcatscnormi(ilist_catscnorm),contlistcatscnormj(ilist_catscnorm)
+
else
ilist_catsctran=ilist_catsctran+1
! this can be substituted by cantor and anti-cantor
ilist_catscnorm,ilist_catpnorm,ilist_catscang
do i=1,ilist_catsctran
- write (iout,*) i,contlistcatsctrani(i),contlistcatsctranj(i)
+ write (iout,*) i,contlistcatsctrani(i),contlistcatsctranj(i),&
+ itype(j,contlistcatsctranj(i))
enddo
do i=1,ilist_catptran
write (iout,*) i,contlistcatptrani(i),contlistcatsctranj(i)
enddo
do i=1,ilist_catscnorm
- write (iout,*) i,contlistcatscnormi(i),contlistcatsctranj(i)
+ write (iout,*) i,contlistcatscnormi(i),contlistcatscnormj(i)
enddo
do i=1,ilist_catpnorm
- write (iout,*) i,contlistcatpnormi(i),contlistcatsctranj(i)
+ write (iout,*) i,contlistcatpnormi(i),contlistcatscnormj(i)
enddo
do i=1,ilist_catscang
write (iout,*) i,contlistcatscangi(i),contlistcatscangi(i)
i_ilist_catpnorm(0:nprocs),i_ilist_catsctran(0:nprocs),&
i_ilist_catptran(0:nprocs),i_ilist_catscang(0:nprocs),&
i_ilist_catscangf(0:nprocs),i_ilist_catscangt(0:nprocs)
- write(iout,*),"START make_catcat"
+! write(iout,*),"START make_catcat"
ilist_catpnorm=0
ilist_catscnorm=0
ilist_catptran=0
write (iout,*) i,newcontlistcatcatnormi(i),newcontlistcatcatnormj(i)
enddo
#endif
- write(iout,*),"END make_catcat"
+! write(iout,*),"END make_catcat"
return
end subroutine make_cat_cat_list
boxx(1)=boxxsize
boxx(2)=boxysize
boxx(3)=boxzsize
+ write(iout,*) "start ecattran",g_listcatsctran_start,g_listcatsctran_end
do k=g_listcatsctran_start,g_listcatsctran_end
i=newcontlistcatsctrani(k)
j=newcontlistcatsctranj(k)
!c & itype(chain_border(1,i)+k).eq.itype(chain_border(1,j)+k))
do k=0,chain_length(i)-1
!c k=k+1
- mnum=molnum(k)
+ mnum=molnum(k+1)
if (itype(chain_border(1,i)+k,mnum).ne.&
itype(chain_border(1,j)+k,mnum)) exit
enddo
read(iion,*) msc(i,5),restok(i,5),ichargecat(i)
print *,msc(i,5),restok(i,5)
enddo
- ip(5)=0.2
+! ip(5)=0.2
! isc(5)=0.2
read (iion,*) ncatprotparm
allocate(catprm(ncatprotparm,4))
!HERE THE MASS of MARTINI
write(*,*) "before MARTINI PARAM"
do i=1,ntyp_molec(4)
- msc(i,4)=0.0d0
- mp(4)=72.0d0
+ msc(i,4)=72.0d0
+ mp(4)=0.0d0
isc(i,4)=0.d0
enddo
ip(4)=0.0
+ msc(ntyp_molec(4)+1,4)=0.1d0
!relative dielectric constant = 15 for implicit screening
k_coulomb_lip=332.0d0/15.0d0
!kbond = 1250 kJ/(mol*nm*2)
! Ions by Aga
- allocate(alphapolcat(ntyp,ntyp),epsheadcat(ntyp,ntyp),sig0headcat(ntyp,ntyp))
- allocate(alphapolcat2(ntyp,ntyp))
- allocate(sigiso1cat(ntyp,ntyp),rborn1cat(ntyp,ntyp),rborn2cat(ntyp,ntyp),sigmap1cat(ntyp,ntyp))
- allocate(sigmap2cat(ntyp,ntyp),sigiso2cat(ntyp,ntyp))
- allocate(chis1cat(ntyp,ntyp),chis2cat(ntyp,ntyp),wquadcat(ntyp,ntyp),chipp1cat(ntyp,ntyp),chipp2cat(ntyp,ntyp))
- allocate(epsintabcat(ntyp,ntyp))
- allocate(dtailcat(2,ntyp,ntyp))
- allocate(alphasurcat(4,ntyp,ntyp),alphisocat(4,ntyp,ntyp))
- allocate(wqdipcat(2,ntyp,ntyp))
- allocate(wstatecat(4,ntyp,ntyp))
- allocate(dheadcat(2,2,ntyp,ntyp))
- allocate(nstatecat(ntyp,ntyp))
- allocate(debaykapcat(ntyp,ntyp))
-
- if (.not.allocated(epscat)) allocate (epscat(0:ntyp1,0:ntyp1))
- if (.not.allocated(sigmacat)) allocate(sigmacat(0:ntyp1,0:ntyp1))
+ allocate(alphapolcat(ntyp,-1:ntyp_molec(5)),epsheadcat(ntyp,-1:ntyp_molec(5)),sig0headcat(ntyp,-1:ntyp_molec(5)))
+ allocate(alphapolcat2(ntyp,-1:ntyp_molec(5)))
+ allocate(sigiso1cat(ntyp,-1:ntyp_molec(5)),rborn1cat(ntyp,-1:ntyp_molec(5)),rborn2cat(ntyp,-1:ntyp_molec(5)),sigmap1cat(ntyp,-1:ntyp_molec(5)))
+ allocate(sigmap2cat(ntyp,-1:ntyp_molec(5)),sigiso2cat(ntyp,-1:ntyp_molec(5)))
+ allocate(chis1cat(ntyp,-1:ntyp_molec(5)),chis2cat(ntyp,-1:ntyp_molec(5)),wquadcat(ntyp,-1:ntyp_molec(5)),chipp1cat(ntyp,-1:ntyp_molec(5)),chipp2cat(ntyp,-1:ntyp_molec(5)))
+ allocate(epsintabcat(ntyp,-1:ntyp_molec(5)))
+ allocate(dtailcat(2,ntyp,-1:ntyp_molec(5)))
+ allocate(alphasurcat(4,ntyp,-1:ntyp_molec(5)),alphisocat(4,ntyp,-1:ntyp_molec(5)))
+ allocate(wqdipcat(2,ntyp,-1:ntyp_molec(5)))
+ allocate(wstatecat(4,ntyp,-1:ntyp_molec(5)))
+ allocate(dheadcat(2,2,ntyp,-1:ntyp_molec(5)))
+ allocate(nstatecat(ntyp,-1:ntyp_molec(5)))
+ allocate(debaykapcat(ntyp,-1:ntyp_molec(5)))
+
+ if (.not.allocated(epscat)) allocate (epscat(0:ntyp1,-1:ntyp1))
+ if (.not.allocated(sigmacat)) allocate(sigmacat(0:ntyp1,-1:ntyp1))
! if (.not.allocated(chicat)) allocate(chicat(ntyp1,ntyp1)) !(ntyp,ntyp)
- if (.not.allocated(chi1cat)) allocate(chi1cat(ntyp1,ntyp1)) !(ntyp,ntyp)
- if (.not.allocated(chi2cat)) allocate(chi2cat(ntyp1,ntyp1)) !(ntyp,ntyp)
+ if (.not.allocated(chi1cat)) allocate(chi1cat(ntyp1,-1:ntyp1)) !(ntyp,ntyp)
+ if (.not.allocated(chi2cat)) allocate(chi2cat(ntyp1,-1:ntyp1)) !(ntyp,ntyp)
if (.not.allocated(ichargecat))&
! mp(5)=0.2
pstok(5)=3.0
!DIR$ NOUNROLL
- do j=1,ntyp_molec(5)-1 ! this is without Zn will be modified for ALL tranistion metals
+ do j=-1,ntyp_molec(5)-1 ! this is without Zn will be modified for ALL tranistion metals
+ if (j.eq.0) cycle
do i=1,ntyp
! do j=1,ntyp_molec(5)
! write (*,*) "Im in ALAB", i, " ", j
read(iion,*) &
epscat(i,j),sigmacat(i,j), &
! chicat(i,j),chicat(j,i),chippcat(i,j),chippcat(j,i), &
- chi1cat(i,j),chi2cat(i,j),chipp1cat(i,j),chipp2cat(i,j), &
+ chi1cat(i,j),chi2cat(i,j),chipp1cat(i,j),chipp2cat(i,j), & !6
- (alphasurcat(k,i,j),k=1,4),sigmap1cat(i,j),sigmap2cat(i,j),&
+ (alphasurcat(k,i,j),k=1,4),sigmap1cat(i,j),sigmap2cat(i,j),&!12
! chiscat(i,j),chiscat(j,i), &
chis1cat(i,j),chis2cat(i,j), &
- nstatecat(i,j),(wstatecat(k,i,j),k=1,4), & !5 w tej lini - 1 integer pierwszy
- dheadcat(1,1,i,j),dheadcat(1,2,i,j),dheadcat(2,1,i,j),dheadcat(2,2,i,j),&
+ nstatecat(i,j),(wstatecat(k,i,j),k=1,4), & !19 !5 w tej lini - 1 integer pierwszy
+ dheadcat(1,1,i,j),dheadcat(1,2,i,j),dheadcat(2,1,i,j),dheadcat(2,2,i,j),&!23
dtailcat(1,i,j),dtailcat(2,i,j), &
- epsheadcat(i,j),sig0headcat(i,j), &
+ epsheadcat(i,j),sig0headcat(i,j), &!27
!wdipcat = w1 , w2
! rborncat(i,j),rborncat(j,i),&
rborn1cat(i,j),rborn2cat(i,j),&
- (wqdipcat(k,i,j),k=1,2), &
- alphapolcat(i,j),alphapolcat2(j,i), &
+ (wqdipcat(k,i,j),k=1,2), &!31
+ alphapolcat(i,j),alphapolcat2(j,i), &!33
(alphisocat(k,i,j),k=1,4),sigiso1cat(i,j),sigiso2cat(i,j),epsintabcat(i,j),debaykapcat(i,j)
if (chi1cat(i,j).gt.0.9) write (*,*) "WTF ANISO", i,j, chi1cat(i,j)
END DO
END DO
- allocate(aa_aq_cat(-ntyp:ntyp,ntyp),bb_aq_cat(-ntyp:ntyp,ntyp))
+ allocate(aa_aq_cat(-ntyp:ntyp,-1:ntyp_molec(5)),&
+ bb_aq_cat(-ntyp:ntyp,-1:ntyp_molec(5)))
do i=1,ntyp
- do j=1,ntyp_molec(5)
+ do j=-1,ntyp_molec(5)
+ if (j.eq.0) cycle
epsij=epscat(i,j)
rrij=sigmacat(i,j)
rrij=rrij**expon
call exec_softreg
else if (modecalc.eq.12) then
call exec_MD
-! call exec_checkgrad
+ call exec_checkgrad
else if (modecalc.eq.14) then
call exec_MREMD
else