Add some Amber to NWChem comparisions

This commit is contained in:
Kurt Glaesmann 2009-11-03 18:10:52 +00:00
parent c28152a3de
commit 70d69d5e27
37 changed files with 17047 additions and 0 deletions

50
QA/amber-tests/README Normal file
View file

@ -0,0 +1,50 @@
The dotests.msub MOAB script runs a bunch of NWChem and Amber runs, you might have to
edit it to set the right path and account number. The resulting test.out.* file
should be empty when the job completes.
Each one of these subdirectories is a seperate test with its own README file.
big: amino acid sequence
big.shake: amino acid sequence, which shake enabled
ethane: single ethane molecule
ethane+: single ethane molecule with a +4 charge
small: ALA dimer
methanol: solvent only run (NWChem only)
Within each directory is a README file that describes what results you should get.
-------
Hints for generating your own amber tests:
First you need a PDB file of the molecule of interest (assume this is tiny.pdb)
Do this in an xterm
module purge
module load amber/10
./xleap &
Fix the PDB to match amber based upon these files:
${AMBERHOME}/dat/leap/lib/all_amino94.lib ("normal" amino acids)
${AMBERHOME}/dat/leap/lib/all_aminont94.lib (N terminating)
${AMBERHOME}/dat/leap/lib/all_aminoct94.lib (C terminating)
Until all atoms are correctly named when you type:
tiny = loadPdb tiny.pdb on the xleap command line
For example, I had to change CD to CD1 in ILE
You may have to "edit tiny" to set atom types, add bonds, set charges, etc. I
had to do that with ethane (there was nothing in the default libraries, because
it is not an amino acid).
To add bonds, you must be in "draw" mode
Once you think you are done, type this in the xleap command line:
check tiny
If all is well then save it as tiny2 (note this this is new name, since we 'fixed' things).
savePdb tiny tiny2.pdb
saveOff tiny tiny2.lib
saveAmberParm tiny tiny2.top tiny2.cor
-------
Hints for NWChem runs:
The atoms might need changed to match *.frg files for NWChem
This is why the pdb files might not be the same for amber and NWChem

View file

@ -0,0 +1,52 @@
This is a series of amino acids, with SHAKE ON
NWChem input files:
big.nw - nwchem input deck
big.pdb - pdb file with geometry
THR_N.sgm - amber95 segment file: one in SVN does not match
amber input files:
big_am.in - amber input deck
big_am.top - amber topology file
big_am.cor - coordinates
NWChem output files:
big.seq
big.nam
VAL.mat
THR.mat
SER.mat
PHE.mat
MET.mat
LEU.mat
ILE.mat
GLU.mat
big.top
big_md.rst
big_md.cmd
big.bsg
ASP.mat
ASN.mat
ARG.mat
ALA.mat
big_energy.out
big.out
big_md.qrs
big_md.out
big.db
amber output files:
mdinfo
big_am.rst
big_am.out
results from *.out files:
4.266015438E+03 kJ/mol for NWChem single point (10.1960E+02 kcal/mol)
-1.29173E+03 kJ/mol for 300 steps NWChem (-2.0873E+02 kcal/mol)
-2.70772E+03 kJ/mol for 1500 steps NWChem (-6.4716E+02 kcal/mol)
-2.70992E+03 kJ/mol for 3000 steps NWChem (-6.4769E+02 kcal/mol)
9.4289E+02 kcal/mol for amber single point
-4.5606E+02 kcal/mol for amber 300 steps
-6.4034E+02 kcal/mol for amber 1500 steps
-6.4192E+02 kcal/mol for amber 3000 steps
Note that they match, but I had to limit NWChem's maximum step size to get this

View file

@ -0,0 +1,185 @@
# This file matches the amber force field
# 07/31/09 13:21:28
4.600000
16 15 26 33 0 0 1 1
0.000000
1 N 0 0 0 1 1
N3 0.181200 0.000000
22H 0 0 0 1 1
H 0.193400 0.000000
33H 0 0 0 1 1
H 0.193400 0.000000
44H 0 0 0 1 1
H 0.193400 0.000000
5 CA 0 0 0 1 1
CT 0.003400 0.000000
6 HA 0 0 0 1 1
HP 0.108700 0.000000
7 CB 0 0 0 1 1
CT 0.451400 0.000000
8 HB 0 0 0 1 1
H1 -0.032300 0.000000
9 CG2 0 0 0 1 1
CT -0.255400 0.000000
102HG2 0 0 0 1 1
HC 0.062700 0.000000
113HG2 0 0 0 1 1
HC 0.062700 0.000000
124HG2 0 0 0 1 1
HC 0.062700 0.000000
13 OG1 0 0 0 1 1
OH -0.676400 0.000000
14 HG1 0 0 0 1 1
HO 0.407000 0.000000
15 C 2 1 0 1 1
C 0.616300 0.000000
16 O 0 0 0 1 1
O -0.572200 0.000000
1 1 2 0 0
0.000000 0.00000E+00
2 1 3 0 0
0.000000 0.00000E+00
3 1 4 0 0
0.000000 0.00000E+00
4 1 5 0 0
0.000000 0.00000E+00
5 5 6 0 0
0.000000 0.00000E+00
6 5 7 0 0
0.000000 0.00000E+00
7 5 15 0 0
0.000000 0.00000E+00
8 7 8 0 0
0.000000 0.00000E+00
9 7 9 0 0
0.000000 0.00000E+00
10 7 13 0 0
0.000000 0.00000E+00
11 9 10 0 0
0.000000 0.00000E+00
12 9 11 0 0
0.000000 0.00000E+00
13 9 12 0 0
0.000000 0.00000E+00
14 13 14 0 0
0.000000 0.00000E+00
15 15 16 0 0
0.000000 0.00000E+00
1 2 1 3 0 0
0.000000 0.00000E+00
2 2 1 4 0 0
0.000000 0.00000E+00
3 2 1 5 0 0
0.000000 0.00000E+00
4 3 1 4 0 0
0.000000 0.00000E+00
5 3 1 5 0 0
0.000000 0.00000E+00
6 4 1 5 0 0
0.000000 0.00000E+00
7 1 5 6 0 0
0.000000 0.00000E+00
8 1 5 7 0 0
0.000000 0.00000E+00
9 1 5 15 0 0
0.000000 0.00000E+00
10 6 5 7 0 0
0.000000 0.00000E+00
11 6 5 15 0 0
0.000000 0.00000E+00
12 7 5 15 0 0
0.000000 0.00000E+00
13 5 7 8 0 0
0.000000 0.00000E+00
14 5 7 9 0 0
0.000000 0.00000E+00
15 5 7 13 0 0
0.000000 0.00000E+00
16 8 7 9 0 0
0.000000 0.00000E+00
17 8 7 13 0 0
0.000000 0.00000E+00
18 9 7 13 0 0
0.000000 0.00000E+00
19 7 9 10 0 0
0.000000 0.00000E+00
20 7 9 11 0 0
0.000000 0.00000E+00
21 7 9 12 0 0
0.000000 0.00000E+00
22 10 9 11 0 0
0.000000 0.00000E+00
23 10 9 12 0 0
0.000000 0.00000E+00
24 11 9 12 0 0
0.000000 0.00000E+00
25 7 13 14 0 0
0.000000 0.00000E+00
26 5 15 16 0 0
0.000000 0.00000E+00
1 2 1 5 6 0 0
0 0.000000 0.00000E+00
2 2 1 5 7 0 0
0 0.000000 0.00000E+00
3 2 1 5 15 0 0
0 0.000000 0.00000E+00
4 3 1 5 6 0 0
0 0.000000 0.00000E+00
5 3 1 5 7 0 0
0 0.000000 0.00000E+00
6 3 1 5 15 0 0
0 0.000000 0.00000E+00
7 4 1 5 6 0 0
0 0.000000 0.00000E+00
8 4 1 5 7 0 0
0 0.000000 0.00000E+00
9 4 1 5 15 0 0
0 0.000000 0.00000E+00
10 1 5 7 8 0 0
0 0.000000 0.00000E+00
11 1 5 7 9 0 0
0 0.000000 0.00000E+00
12 1 5 7 13 0 0
0 0.000000 0.00000E+00
13 6 5 7 8 0 0
0 0.000000 0.00000E+00
14 6 5 7 9 0 0
0 0.000000 0.00000E+00
15 6 5 7 13 0 0
0 0.000000 0.00000E+00
16 15 5 7 8 0 0
0 0.000000 0.00000E+00
17 15 5 7 9 0 0
0 0.000000 0.00000E+00
18 15 5 7 13 0 0
0 0.000000 0.00000E+00
19 1 5 15 16 0 0
0 0.000000 0.00000E+00
20 6 5 15 16 0 0
0 0.000000 0.00000E+00
21 7 5 15 16 0 0
0 0.000000 0.00000E+00
22 5 7 9 10 0 0
0 0.000000 0.00000E+00
23 5 7 9 11 0 0
0 0.000000 0.00000E+00
24 5 7 9 12 0 0
0 0.000000 0.00000E+00
25 8 7 9 10 0 0
0 0.000000 0.00000E+00
26 8 7 9 11 0 0
0 0.000000 0.00000E+00
27 8 7 9 12 0 0
0 0.000000 0.00000E+00
28 13 7 9 10 0 0
0 0.000000 0.00000E+00
29 13 7 9 11 0 0
0 0.000000 0.00000E+00
30 13 7 9 12 0 0
0 0.000000 0.00000E+00
31 5 7 13 14 0 0
0 0.000000 0.00000E+00
32 8 7 13 14 0 0
0 0.000000 0.00000E+00
33 9 7 13 14 0 0
0 0.000000 0.00000E+00

View file

@ -0,0 +1,32 @@
#
# This is a AMBER test that uses a PDB geometry, and the *.frg files from NWChem to asign atom types
#
Title "big"
memory 3500 Mb
permanent_dir .
echo
prepare
system big
source big.pdb
new_top new_seq
new_rst
amber
vacuo
end
print debug
task prepare
md
system big_md
sd 3000 init 0.001 min 0.000000004 max 0.01
cutoff short 1500 long 1500
shake 1000 0.000001 # Match 0.00001 A in Amber
end
task md energy
task shell "cp big_md.out big_energy.out"
task md optimize

File diff suppressed because it is too large Load diff

View file

@ -0,0 +1,540 @@
1076
-7.1470000 23.8710000 25.0130000 -7.2952944 23.9901834 26.0049199
-6.2946862 23.3475465 24.8727927 -8.0087968 23.5550476 24.5915977
-6.8910000 25.1930000 24.4280000 -5.8850703 25.2745585 24.0162320
-7.9230000 25.6260000 23.3230000 -7.8345081 26.7056384 23.2019616
-7.7040000 24.9430000 21.9740000 -7.7928053 23.8632894 22.0941593
-8.4534963 25.2918616 21.2636128 -6.7092326 25.1858804 21.6004395
-9.2380000 25.3860000 23.9080000 -9.9185752 25.6325634 23.2774255
-6.8010000 26.2280000 25.5430000 -5.8290000 26.9990000 25.5610000
-7.7810000 26.2170000 26.4190000 -8.5466857 25.5741995 26.2753607
-7.8580000 27.1540000 27.5570000 -8.1329647 28.1400214 27.1824931
-8.8810000 26.7260000 28.5920000 -9.7307785 26.3972520 27.9937521
-8.4181838 25.8626578 29.0700601 -9.3580000 27.6990000 29.6570000
-8.4615931 28.0538152 30.1655880 -9.8244452 28.5273104 29.1236802
-10.3220000 27.1490000 30.6680000 -10.7750000 26.0340000 30.3340000
-10.6480000 27.7010000 31.7060000 -6.5020000 27.3430000 28.2240000
-6.0400000 28.4830000 28.4230000 -5.8780000 26.2240000 28.5550000
-6.3422524 25.3435754 28.3834728 -4.5540000 26.2010000 29.1940000
-4.5434862 26.8959116 30.0336943 -4.2200000 24.8180000 29.7070000
-4.9839982 24.6501070 30.4660902 -4.4211444 24.1841216 28.8433802
-2.8590000 24.4770000 30.2770000 -2.1496331 24.5751268 29.4552515
-2.6308417 25.2174871 31.0436302 -2.7550000 23.0960000 30.8670000
-3.0960000 22.0690000 30.3040000 -2.2840000 23.1150000 32.0290000
-3.4880000 26.6540000 28.1940000 -2.4920000 27.2840000 28.5750000
-3.7420000 26.3350000 26.9340000 -4.5632789 25.7725947 26.7628262
-2.8740000 26.6630000 25.7970000 -1.8532639 26.3764321 26.0501323
-3.2730000 25.9550000 24.5100000 -4.3625372 25.9292851 24.4913794
-2.9105382 26.5718993 23.6877125 -2.7310000 24.5440000 24.3610000
-1.6569867 24.5496612 24.5469148 -3.2233721 23.9090851 25.0975832
-2.9800000 23.9700000 22.9850000 -4.0820000 23.8580000 22.4520000
-1.8930000 23.5810000 22.3350000 -1.9761311 23.1864487 21.4089763
-0.9862449 23.6808062 22.7685151 -2.8380000 28.1870000 25.5880000
-1.7320000 28.7490000 25.5180000 -4.0170000 28.7720000 25.4880000
-4.8578319 28.2179494 25.4097118 -4.1180000 30.2260000 25.2840000
-3.4600921 30.5740649 24.4876915 -5.5690000 30.7200000 25.0110000
-5.4951060 31.6448304 24.4388734 -6.4410000 29.6970000 24.2310000
-6.5157420 28.7720090 24.8027568 -7.4378377 30.1105520 24.0780621
-5.9818748 29.4899270 23.2643439 -6.2490000 31.1530000 26.3400000
-6.2198611 30.2967482 27.0138579 -5.6550142 31.9650681 26.7593157
-7.7190000 31.6240000 26.1650000 -8.3137428 30.8122554 25.7461306
-8.1272405 31.9101133 27.1343184 -7.7488960 32.4805745 25.4915861
-3.4050000 30.9160000 26.4440000 -2.6860000 31.9090000 26.2400000
-3.5680000 30.3860000 27.6410000 -4.2157133 29.6183928 27.7475173
-2.9330000 30.9330000 28.8460000 -3.3394863 31.9323099 29.0017169
-3.2810000 30.0700000 30.0580000 -2.8727698 29.0686695 29.9209405
-2.8545399 30.5158083 30.9566028 -4.3642872 30.0090282 30.1622613
-1.4270000 31.0940000 28.6770000 -0.8490000 32.1200000 29.0600000
-0.7870000 30.0840000 28.1070000 -1.2993264 29.2576866 27.8334544
0.6700000 30.1020000 27.9000000 1.1638343 30.3981376 28.8254887
1.1910000 28.7180000 27.4870000 0.7371114 28.4606583 26.5299944
2.2722547 28.7883350 27.3685037 0.8590000 27.6300000 28.5200000
1.3569297 27.9177537 29.4459399 -0.2209283 27.6562455 28.6654830
1.2770000 26.2450000 28.1390000 2.4510000 26.0280000 28.5170000
0.5840000 25.4400000 27.5470000 1.0980000 31.1160000 26.8520000
2.1830000 31.7040000 26.9570000 0.2230000 31.2830000 25.8790000
-0.5972086 30.6953200 25.8343969 0.4620000 32.2360000 24.7770000
1.5117698 32.1719031 24.4906886 -0.3930000 31.9690000 23.5500000
-1.3988517 31.7051427 23.8767123 -0.4338604 32.8751424 22.9455693
0.0660000 30.8670000 22.6410000 1.3280000 30.9200000 22.0460000
1.9859561 31.7643528 22.2515773 1.7520000 29.8910000 21.1850000
2.7360664 29.9255917 20.7175415 0.8930000 28.8130000 20.9310000
1.2186363 28.0152960 20.2633705 -0.3690000 28.7550000 21.5240000
-1.0396694 27.9196850 21.3226260 -0.7670000 29.7740000 22.3770000
-1.7451538 29.7327705 22.8561818 0.2710000 33.6720000 25.2990000
1.0020000 34.5320000 24.8260000 -0.6590000 33.8570000 26.2030000
-1.1889504 33.0459457 26.4883819 -0.9870000 35.1590000 26.8170000
-1.2154939 35.9157294 26.0664999 -2.2170000 34.9300000 27.6610000
-2.8372107 34.2830112 27.0406426 -1.8416226 34.3639684 28.5135257
-3.0410000 36.0940000 28.1540000 -2.5671267 36.5361902 29.0303634
-3.1245080 36.8453666 27.3687760 -4.4420000 35.5890000 28.5300000
-4.3468270 34.6503177 29.0758174 -4.9284546 36.3298969 29.1644551
-5.2760000 35.3660000 27.2830000 -5.6469602 36.3283130 26.9302356
-4.6477231 34.9155513 26.5145839 -6.4240000 34.4690000 27.5720000
-7.0065357 34.8864127 28.2837011 -6.9656059 34.3335555 26.7303244
-6.0804657 33.5770978 27.8984870 0.1920000 35.6210000 27.6560000
0.6130000 36.7880000 27.7760000 0.7740000 34.6060000 28.3030000
0.3905739 33.6728934 28.2540534 1.9480000 34.8320000 29.1540000
1.7002357 35.6194119 29.8658248 2.4020000 33.6050000 29.9260000
2.0895570 32.7513289 29.3245649 3.4905710 33.6583681 29.9422819
1.8810000 33.4420000 31.3470000 2.4357002 34.1636184 31.9467295
0.8264939 33.7174548 31.3314828 2.0320000 32.0700000 31.9470000
3.1850000 31.6130000 31.7950000 1.1110000 31.4900000 32.4970000
3.0920000 35.3460000 28.2880000 3.7460000 36.3130000 28.6840000
3.3120000 34.7090000 27.1510000 2.7882604 33.8653449 26.9664945
4.3750000 35.0890000 26.2100000 5.3128376 35.1331221 26.7637265
4.4990000 34.0200000 25.1230000 3.5706476 33.9702961 24.5539734
5.3212291 34.2747119 24.4542987 4.6936105 33.0521390 25.5850294
4.1460000 36.5060000 25.6650000 5.1080000 37.2850000 25.6350000
2.9340000 36.8130000 25.2870000 2.2645321 36.0592810 25.3488043
2.4820000 38.1060000 24.7440000 3.0646038 38.3254547 23.8492868
1.0020000 38.1230000 24.3780000 0.8294572 37.3041294 23.6795881
0.4132467 37.9600628 25.2807298 0.4630000 39.3530000 23.7340000
0.5340000 39.5400000 22.3650000 0.9821541 38.7505500 21.7616538
0.0500000 40.7030000 21.7500000 0.0914942 40.8622831 20.6724993
-0.5060000 41.6860000 22.5840000 -0.8822032 42.5909965 22.1069756
-0.5890000 41.5410000 23.9750000 -1.0098747 42.3213205 24.6090857
-0.1150000 40.3590000 24.5450000 -0.1884435 40.2054087 25.6216224
2.7710000 39.2120000 25.7550000 3.3410000 40.2480000 25.3760000
2.3980000 38.9680000 27.0170000 1.8849722 38.1377501 27.2769747
2.6200000 40.0030000 28.0440000 2.1258918 40.9265118 27.7422018
2.0540000 39.5650000 29.4030000 2.6079684 38.7161387 29.8038152
2.0773809 40.3850272 30.1207108 0.6940000 39.1780000 29.1640000
0.2925768 38.8938291 29.9884429 4.0650000 40.4240000 28.1640000
4.2810000 41.5960000 28.5420000 5.0170000 39.5370000 27.8630000
4.7950387 38.5902288 27.5901321 6.4290000 39.9160000 27.9720000
6.6568737 40.2769909 28.9749249 7.2920000 38.6690000 27.6890000
7.2337959 38.0745335 28.6007681 6.7955835 38.1375481 26.8770705
8.7570000 38.9000000 27.3250000 8.8029770 39.7503783 26.6446738
9.6280000 39.1780000 28.5450000 9.5828944 38.3278105 29.2256189
10.6590326 39.3346095 28.2279102 9.2656828 40.0706203 29.0549565
9.2880000 37.6530000 26.6110000 8.7055910 37.4756998 25.7068629
10.3341965 37.8038745 26.3448966 9.2025703 36.7908889 27.2724871
6.7600000 41.0540000 26.9940000 7.5800000 41.9470000 27.2900000
6.0910000 41.0070000 25.8520000 5.4307525 40.2667025 25.6619154
6.3260000 42.0200000 24.8060000 7.3684392 42.3377761 24.8269289
6.0180000 41.4410000 23.4200000 5.0388943 40.9626546 23.4452378
6.0136512 42.2373766 22.6757787 6.9850000 40.4100000 22.9370000
8.1600000 40.7720000 22.2780000 8.3954353 41.8236562 22.1146311
9.0360000 39.7820000 21.8270000 9.9567870 40.0267817 21.2975344
8.7050000 38.4490000 22.0710000 9.4067756 37.6949061 21.7146948
7.5370000 38.0590000 22.7380000 7.3047022 37.0115918 22.9305459
6.6500000 39.0670000 23.1640000 5.7185119 38.7983747 23.6622672
5.5010000 43.3130000 24.9660000 6.0550000 44.3900000 24.6970000
4.2660000 43.1210000 25.3760000 3.9263426 42.1758088 25.4825166
3.3040000 44.2450000 25.5730000 3.4293272 45.0030645 24.7998696
1.8920000 43.6880000 25.4830000 1.8101074 43.3303314 24.4566151
1.7314459 42.8575425 26.1705047 0.8380000 44.7660000 25.7320000
1.3620000 45.9000000 25.6800000 -0.3270000 44.4680000 26.0380000
3.6700000 44.9290000 26.9000000 2.8890000 44.7430000 27.8230000
4.7280000 45.6930000 26.9420000 5.2193341 45.7357518 26.0606020
5.2930000 46.4000000 28.0820000 5.5879603 45.6460336 28.8118166
6.5030000 47.2920000 27.7830000 6.1870956 47.9567022 26.9789635
6.6711418 47.8734652 28.6894913 7.7990000 46.5880000 27.3820000
7.7433488 46.2989768 26.3324909 8.6390003 47.2672771 27.5272040
7.9970000 45.3450000 28.2410000 7.6903950 45.5856862 29.2589207
7.3566998 44.5562019 27.8461422 9.4320000 44.8730000 28.2460000
9.8710370 45.0204001 27.2592780 10.0019493 45.4366840 28.9845915
9.4660000 43.4200000 28.5940000 8.9382879 42.8974492 27.9094854
10.4230759 43.0973880 28.5992206 9.0596122 43.2832510 29.5084658
4.2940000 47.2800000 28.8280000 4.3150000 47.3110000 30.0600000
3.4440000 47.9600000 28.0720000 3.5341250 48.0617198 27.0711854
2.4670000 48.8110000 28.7760000 2.7003921 48.8893755 29.8378308
2.5490000 50.2110000 28.1510000 2.1072477 50.9568183 28.8118397
3.5933687 50.4627538 27.9665711 1.7910000 50.1910000 26.8300000
1.4800000 49.1260000 26.2870000 1.4870000 51.2940000 26.3370000
1.1140000 48.1780000 28.8810000 0.1890000 48.8260000 29.3960000
0.9010000 46.9280000 28.4700000 1.6406840 46.3673524 28.0717020
-0.4170000 46.3410000 28.5890000 -0.2759187 45.3145075 28.2506091
-0.6079840 46.3462481 29.6621248 -1.6030000 46.9320000 27.8780000
-2.7250000 46.7820000 28.4030000 -1.4540000 47.5460000 26.7100000
-0.5276198 47.7736193 26.3781741 -2.6160000 48.0790000 25.9660000
-3.4504435 48.2676213 26.6414445 -2.2130000 49.4070000 25.2990000
-3.0941511 50.0152164 25.0946844 -1.5244527 49.9669124 25.9318513
-1.5220000 49.0540000 23.9830000 -0.6470000 48.1980000 24.1440000
-1.7960000 49.5540000 22.9200000 -3.2190000 47.0540000 25.0400000
-4.2130000 47.3330000 24.3510000 -2.6310000 45.8600000 24.9310000
-1.7937363 45.6465941 25.4540169 -3.1850000 44.8220000 24.0470000
-3.0486580 43.8557301 24.5326259 -4.2498607 45.0177026 23.9210175
-2.5100000 44.7850000 22.6510000 -3.0400000 44.0510000 21.7800000
-1.4680000 45.5220000 22.5170000 -1.1664003 46.1374373 23.2588719
-0.7420000 45.5700000 21.2250000 -0.8399244 44.6081876 20.7215848
-1.2500000 46.7020000 20.2520000 -0.7118239 46.5492103 19.3165204
-2.7520000 46.7520000 19.9730000 -3.2908520 46.9052421 20.9080160
-2.9681714 47.5741142 19.2907264 -3.0695956 45.8123482 19.5210176
-0.8240000 47.9600000 20.8890000 -1.1013569 48.7027124 20.3476575
0.7320000 45.7030000 21.4940000 1.2300000 46.2590000 22.4720000
1.4870000 45.1760000 20.4980000 1.0487052 44.7120341 19.7152276
2.9480000 45.2260000 20.5450000 3.2927821 45.6756274 21.4761604
3.6310000 43.8280000 20.3560000 3.4213691 43.4301437 19.3630923
5.1700000 44.0380000 20.5180000 5.3804681 44.4354312 21.5109009
5.6825375 43.0841965 20.3928404 5.5226892 44.7411541 19.7634910
3.0610000 42.7720000 21.3220000 3.4990538 42.9604110 22.3021587
1.9830975 42.9264891 21.3705680 3.3570000 41.2980000 20.8940000
4.4347835 41.1426106 20.8456600 2.9240844 40.6134118 21.6233986
2.9188288 41.1086894 19.9140677 3.3730000 46.1960000 19.4300000
2.9680000 46.0580000 18.2830000 4.1300000 47.1870000 19.7840000
4.4578927 47.2108641 20.7389954 4.6140000 48.2310000 18.8480000
4.4416071 48.6362990 17.8509480 4.3970000 49.4070000 19.9260000
3.6976466 49.0691218 20.6907505 5.6970000 49.7590000 20.5630000
6.3967706 50.0974334 19.7988772 5.5416789 50.5552591 21.2909739
6.1050029 48.8822004 21.0658487 3.3280000 50.2610000 19.4650000
3.1872212 50.9679928 20.0989896 6.0660000 48.0390000 18.4870000
6.7090000 47.1940000 19.1490000 6.6250000 48.7930000 17.5510000
6.0615290 49.4637404 17.0483002 8.0420000 48.6610000 17.1860000
8.2457151 47.6409129 16.8603904 8.4750000 49.6540000 16.0360000
9.5427527 49.5091862 15.8715818 7.7300000 49.4510000 14.7290000
6.6621536 49.5961173 14.8925376 8.0853423 50.1711737 13.9919920
7.9068130 48.4397407 14.3626772 8.2360000 51.0000000 16.5430000
8.4857601 51.6426002 15.8749542 8.9720000 48.8980000 18.3700000
10.0100000 48.2160000 18.4770000 8.6200000 49.8650000 19.2040000
7.7717407 50.3747804 19.0023085 9.3950000 50.2470000 20.3860000
10.4308228 50.3724346 20.0706643 8.9370000 51.5280000 21.0760000
8.2371571 51.2283151 21.8560693 9.8227235 51.9691485 21.5331455
8.2710000 52.5550000 20.1620000 8.6228345 53.5554102 20.4139779
8.5146751 52.3358367 19.1224367 6.7550000 52.4840000 20.3560000
6.4870885 51.4546526 20.5942568 6.4848006 53.1315829 21.1901032
6.0090000 52.9250000 19.1120000 4.9614904 52.6374860 19.2023324
6.0809586 54.0084672 19.0170267 6.5970000 52.2820000 17.9070000
6.5300323 51.2780381 17.9946099 6.0927912 52.5826221 17.0851120
7.5675309 52.5486303 17.8228739 9.3590000 49.1630000 21.4490000
10.3710000 48.9180000 22.1160000 8.2200000 48.5160000 21.6010000
7.3706318 48.7590600 21.1115157 8.1600000 47.4340000 22.6290000
8.6362109 47.7829374 23.5452778 6.7230000 47.0390000 22.9470000
6.2081555 46.8944370 21.9971916 6.7582628 46.0934837 23.4881612
5.9710000 48.0750000 23.7810000 6.4722400 48.0394485 24.7482618
6.1727312 49.0335229 23.3028395 4.5170000 47.9190000 23.9750000
3.8300000 47.4500000 23.0670000 3.9000000 48.2470000 24.9960000
8.9910000 46.2690000 22.1360000 9.8140000 45.6820000 22.8800000
8.7980000 45.8960000 20.8860000 8.0650008 46.3224684 20.3374235
9.5710000 44.7760000 20.2820000 9.3974131 43.8553891 20.8391735
9.1200000 44.5610000 18.8210000 8.0825974 44.2294929 18.8657011
9.1707349 45.5275830 18.3197583 9.9530000 43.5330000 18.0370000
10.9942312 43.8549877 18.0208013 9.8560000 42.1680000 18.6940000
8.8150463 41.8451171 18.7101507 10.4500498 41.4493893 18.1293725
10.2331529 42.2279525 19.7149122 9.3550000 43.5290000 16.6210000
9.4375806 44.5260385 16.1883406 9.8977874 42.8174042 15.9988147
8.3049703 43.2407174 16.6702918 11.0360000 45.1150000 20.4430000
11.8210000 44.2260000 20.7670000 11.3990000 46.3700000 20.2520000
10.6638686 47.0183286 20.0083790 12.7700000 46.8830000 20.3660000
13.3936257 46.3743074 19.6308660 12.7595595 47.9523569 20.1551299
13.3560000 46.6570000 21.7550000 14.4850000 46.1500000 21.9040000
12.6240000 47.0050000 22.8050000 11.8363664 47.6260262 22.6864285
13.1530000 46.7970000 24.1690000 14.2381405 46.8567258 24.2526785
12.5100000 47.7260000 25.2720000 13.3168588 48.2756585 25.7567213
11.5050000 48.7040000 24.7020000 10.6975457 48.1550331 24.2174874
11.0962617 49.3158320 25.5061727 11.9964922 49.3463652 23.9713134
12.0100000 46.9180000 26.3790000 11.6213860 47.4913887 27.0436841
13.2680000 45.3120000 24.4460000 14.2690000 44.9430000 25.1200000
12.3610000 44.4970000 23.9230000 11.4984580 44.7926706 23.4886044
12.5570000 43.0500000 24.1840000 12.8219255 42.8653490 25.2250656
11.3560000 42.1830000 23.8260000 11.0758154 42.3954076 22.7942642
11.6660000 40.6940000 23.9730000 11.9458376 40.4807247 25.0046518
10.7839525 40.1108527 23.7083646 12.4896642 40.4267414 23.3109941
10.1300000 42.5500000 24.6370000 9.8668269 43.5907198 24.4479385
9.2978421 41.9073533 24.3495663 10.3420545 42.4163368 25.6977864
13.8590000 42.5480000 23.5470000 14.6440000 41.8510000 24.2410000
14.0890000 42.8560000 22.2950000 13.3741552 43.3462853 21.7766213
15.2910000 42.4030000 21.5710000 15.3568925 41.3178980 21.6504419
15.2710000 42.7830000 20.0950000 15.3761435 43.8657943 20.0271723
16.1249066 42.3047946 19.6151425 14.0100000 42.3600000 19.4060000
13.1516776 42.7160588 19.9757417 13.9849165 42.7906195 18.4049807
13.9330000 40.5450000 19.2790000 14.8980000 40.2740000 17.7850000
15.9071919 40.6599851 17.9286878 14.9449912 39.2062711 17.5708795
14.4273058 40.7913226 16.9489824 16.5740000 42.8760000 22.2350000
17.5090000 42.0580000 22.3820000 16.6210000 44.1340000 22.5990000
15.8433992 44.7279029 22.3485689 17.7840000 44.7480000 23.2480000
18.6734757 44.5904851 22.6379827 17.5700000 46.2440000 23.4180000
16.5887329 46.3411187 23.8825235 18.3380191 46.5641260 24.1221066
17.6290000 47.1080000 22.1830000 18.4612880 46.7466522 21.5789922
16.6951838 46.9566544 21.6415355 17.8100000 48.5480000 22.4590000
18.6153312 48.6879582 23.1800772 18.0541516 49.0736381 21.5358564
16.5770000 49.1200000 23.0140000 16.4304543 49.1638046 24.0123511
15.6340000 49.6200000 22.2080000 15.8000000 49.6480000 20.8890000
16.6510771 49.2867985 20.4824274 15.0742499 50.0303423 20.2997625
14.5040000 50.0690000 22.7510000 14.3764623 50.0261604 23.7519990
13.7796800 50.4503191 22.1593455 18.1270000 44.0710000 24.5770000
19.3100000 44.0620000 24.9760000 17.1300000 43.5320000 25.2370000
16.1819017 43.6285591 24.9025075 17.3510000 42.8510000 26.5310000
18.1006952 43.4146220 27.0863257 16.1090000 42.7340000 27.3800000
16.4075238 42.5390385 28.4100353 15.5637809 43.6766452 27.3324994
15.2770000 41.6850000 26.9240000 14.4958081 41.6317037 27.4794264
17.9600000 41.4760000 26.2330000 18.4830000 40.8140000 27.1430000
17.8850000 41.1050000 24.9610000 17.3702110 41.6741689 24.3043886
18.4260000 39.8200000 24.4950000 18.5448978 39.1311775 25.3313528
17.4560000 39.2020000 23.4870000 17.3976802 39.9632300 22.7090366
17.9710446 38.3238721 23.0974869 16.0590000 38.8160000 23.9200000
15.5462519 39.7144276 24.2635375 15.2690000 38.1850000 22.7830000
15.7809236 37.2861968 22.4392182 14.2714616 37.9221799 23.1350561
15.1878641 38.8939399 21.9590330 16.1370000 37.8380000 25.0900000
16.6571006 38.3091710 25.9240220 15.1295347 37.5610792 25.4005314
16.6801254 36.9446040 24.7818485 19.8410000 39.9720000 23.9450000
20.4100000 38.9480000 23.5300000 20.4050000 41.1680000 23.9430000
19.9064009 41.9715179 24.2977637 21.7740000 41.3660000 23.4250000
22.2237548 42.2067299 23.9531979 22.3509995 40.4626543 23.6228318
21.7970000 41.6520000 21.9240000 22.8160000 41.5940000 21.2130000
20.6380000 42.0150000 21.4270000 19.8481890 42.1000016 22.0507569
20.4420000 42.3450000 19.9980000 21.3541481 42.2113277 19.4164360
19.4030000 41.3740000 19.4950000 19.3941143 40.6024764 20.2649166
18.4800343 41.9513167 19.5492195 19.5110000 40.7330000 18.1490000
18.8811776 41.2645339 17.4356303 20.5382028 40.7172244 17.7847096
19.0140000 39.3010000 18.2630000 18.3870000 38.7610000 17.3670000
19.3000000 38.7140000 19.4280000 19.0091555 37.7612667 19.5947639
19.8058293 39.2254902 20.1369541 19.9200000 43.7750000 19.9210000
18.8960000 44.0340000 20.5780000 20.6000000 44.6250000 19.1790000
21.5120573 44.3847775 18.8176912 20.1420000 46.0260000 19.0520000
19.1721218 46.1491509 19.5339433 21.1380000 46.9550000 19.7220000
22.1011591 46.8244575 19.2286619 20.7690715 47.9637383 19.5364228
21.3310000 46.7760000 21.2200000 21.6970000 45.7290000 21.7750000
21.1160000 47.8420000 22.0060000 21.2344381 47.7605592 23.0057198
20.8367051 48.7209974 21.5943462 19.7970000 46.2680000 17.5760000
20.5720000 46.8110000 16.7650000 18.5870000 45.8500000 17.2230000
17.5830000 45.2170000 18.0910000 17.4766527 45.8107647 18.9988734
17.9227872 44.2144524 18.3508883 16.3110000 45.2000000 17.2500000
15.7123333 46.0838319 17.4703143 15.7324264 44.3043085 17.4760323
16.7780000 45.1930000 15.8070000 16.0396029 45.7066193 15.1913184
16.9778681 44.2071427 15.3871933 18.0690000 46.0090000 15.8590000
18.7986935 45.6190448 15.1493659 17.8160000 47.4800000 15.5680000
17.6570000 48.3130000 16.4860000 17.7860000 47.7800000 14.2760000
17.9730368 47.0890682 13.5634483 17.5410000 49.1730000 13.8530000
17.8684358 49.8699509 14.6244561 18.3210000 49.4900000 12.5260000
18.1763448 50.5517187 12.3261790 19.8200000 49.1770000 12.5610000
19.9654504 48.1152974 12.7603272 20.2671238 49.4302702 11.5997333
20.2956525 49.7623796 13.3478837 17.6680000 48.6750000 11.4960000
18.0992062 48.8216280 10.6509191 16.0270000 49.3370000 13.7010000
15.3340000 48.3170000 13.6440000 15.5580000 50.5700000 13.6140000
16.1570402 51.3709428 13.7545017 14.1130000 50.8020000 13.4220000
13.5625563 50.3241590 14.2324190 13.7770000 52.2870000 13.4150000
14.2846932 52.6807685 14.2955074 14.2621650 52.6672836 12.5160565
12.3300000 52.7250000 13.4300000 12.3871283 53.7720854 13.1325957
11.9101289 52.1381779 12.6130255 11.4570000 52.5970000 14.6340000
12.0590000 52.6470000 15.7350000 10.2340000 52.4680000 14.5570000
13.6740000 50.0870000 12.1480000 12.6040000 49.4510000 12.1180000
14.4720000 50.0990000 11.0930000 15.3251579 50.6375123 11.1401751
14.1410000 49.4270000 9.8280000 13.2237600 49.8850254 9.4578882
15.2480000 49.6410000 8.7910000 16.1706867 49.1799022 9.1433322
14.9533634 49.1863599 7.8451713 15.4087574 50.7092258 8.6455696
13.8810000 47.9430000 9.9980000 12.9390000 47.3380000 9.4580000
14.7710000 47.3070000 10.7520000 15.5649499 47.8044131 11.1292574
14.6580000 45.8690000 11.0230000 14.4990698 45.3468547 10.0794934
15.9040000 45.2970000 11.6930000 16.1750154 45.9802643 12.4978601
15.6240451 44.3318838 12.1152273 17.1020000 45.1150000 10.7520000
16.9489246 44.1276919 10.3162346 17.0300843 45.8715698 9.9706362
18.4610000 45.1770000 11.3840000 18.7270000 45.7390000 12.4370000
19.3280000 44.6160000 10.6720000 13.4320000 45.6430000 11.9120000
12.7530000 44.6550000 11.6780000 13.1810000 46.5440000 12.8450000
13.8397939 47.2946893 12.9951848 12.0240000 46.4290000 13.7660000
12.0167921 45.4484474 14.2419872 12.1250000 47.4920000 14.8460000
12.9304203 48.0954125 14.4273225 11.1870079 48.0433789 14.7807916
12.4280000 47.1830000 16.2950000 11.4905644 47.2865091 16.8414441
13.0070000 45.8060000 16.5430000 13.9447512 45.7018990 15.9972121
13.1914111 45.6754787 17.6093290 12.3022181 45.0479910 16.2012323
13.3730000 48.2470000 16.8600000 12.9010405 49.2270166 16.7898490
13.5904257 48.0229286 17.9043272 14.3012308 48.2495884 16.2886158
10.7550000 46.4920000 12.9150000 9.8490000 45.6880000 13.0890000
10.7170000 47.3630000 11.9090000 11.4366556 48.0652258 11.8137426
9.5290000 47.4600000 11.0460000 8.6186195 47.4735814 11.6452679
9.6070000 48.7640000 10.2440000 9.8682419 49.5820301 10.9153260
10.3738501 48.6662895 9.4755642 8.2650000 49.0560000 9.5870000
8.3820560 49.9046832 8.9131042 7.9359251 48.1851840 9.0199980
7.2040000 49.3950000 10.6120000 6.1410000 48.7870000 10.6450000
7.5300000 50.3940000 11.4340000 6.8805975 50.6872866 12.1497938
8.4247036 50.8524155 11.3367265 9.2830000 46.2080000 10.2240000
8.1210000 45.7510000 10.0850000 10.3200000 45.5840000 9.6620000
11.2351532 46.0059162 9.7296800 10.1790000 44.3410000 8.8800000
9.4481679 44.4873952 8.0846686 11.5180000 43.9460000 8.2480000
12.3148984 44.1143648 8.9723655 11.4794088 42.8875505 7.9905212
11.7920000 44.7630000 6.9940000 10.9000000 45.1790000 6.2500000
13.0030000 44.9620000 6.8170000 9.6290000 43.1990000 9.7460000
8.7950000 42.4180000 9.2930000 10.1050000 43.1390000 10.9780000
10.8714537 43.7252866 11.2761881 9.6130000 42.1060000 11.9070000
9.7114786 41.1246179 11.4430062 10.4040000 42.1110000 13.2060000
10.2110669 43.0548754 13.7158777 10.0470936 41.2860137 13.8225351
11.8820000 41.9600000 12.9740000 12.0771689 41.1626214 12.2569245
12.3084779 42.8925898 12.6045478 12.6620000 41.5340000 14.5750000
14.2680000 42.3040000 14.3430000 14.1385066 43.3741826 14.1816322
14.8805628 42.1423283 15.2299779 14.7609542 41.8634437 13.4763931
8.1170000 42.2530000 12.1280000 7.3760000 41.2260000 12.0980000
7.6360000 43.4550000 12.3720000 8.2755430 44.2310828 12.4656995
6.2060000 43.6970000 12.6140000 5.7085959 42.9745032 13.2610601
6.0110000 45.1450000 13.1980000 6.7264250 45.8207170 12.7293061
4.5730000 45.6400000 12.9150000 3.8569736 44.9649959 13.3838030
4.4466363 46.6429979 13.3225882 4.4006938 45.6612915 11.8389162
6.4510000 45.0960000 14.6820000 6.5661231 44.0402357 14.9273757
7.4235170 45.5864366 14.7242157 5.4990000 45.7570000 15.6950000
4.5263461 45.2667574 15.6536901 5.9118232 45.6600198 16.6991270
5.3837382 46.8129584 15.4505284 5.4020000 43.4650000 11.3440000
4.3480000 42.8280000 11.3270000 5.9490000 43.9350000 10.2320000
6.8160469 44.4514524 10.2720732 5.2950000 43.7750000 8.9300000
4.3500404 44.3177030 8.9549789 6.1310000 44.2980000 7.7500000
7.1774104 44.0403974 7.9136013 5.7924755 43.8613362 6.8104118
6.0120000 45.8160000 7.6530000 6.7040000 46.4240000 6.8340000
5.1550000 46.3250000 8.5170000 4.9978723 47.3224405 8.5398697
4.6599149 45.7149285 9.1516687 4.9440000 42.3130000 8.6540000
3.8500000 41.9930000 8.2010000 5.8990000 41.4510000 8.8980000
6.8066044 41.7526926 9.2225537 5.6800000 40.0170000 8.6790000
5.4796299 39.7952959 7.6307642 6.9800000 39.3060000 9.0920000
7.7488843 39.7837080 8.4847841 7.1162693 39.5864363 10.1364550
7.0910000 37.8030000 8.9550000 6.2081742 37.3902016 9.4431771
7.0426139 37.5994319 7.8852720 8.3240000 37.1670000 9.5300000
9.2130000 38.0030000 9.7940000 8.4390000 35.9600000 9.7210000
4.5200000 39.3970000 9.4430000 3.8960000 38.4410000 8.9500000
4.2200000 39.9020000 10.6270000 4.8023075 40.6555000 10.9635340
3.1560000 39.3980000 11.5060000 2.7589606 38.4240464 11.2198649
3.7070000 39.3030000 12.9520000 3.0643343 38.6382584 13.5292336
5.1220000 38.7660000 12.9570000 5.7653952 39.4305070 12.3803078
5.4858682 38.7088708 13.9828834 5.1352635 37.7713254 12.5114126
3.6490000 40.6010000 13.7140000 2.6139648 40.9353169 13.7849684
4.0519974 40.4530630 14.7159019 4.2390049 41.3552542 13.1933316
1.8800000 40.2070000 11.5010000 0.8100000 39.7220000 11.8920000
1.9650000 41.4610000 11.0750000 2.8523991 41.8844408 10.8440887
0.7380000 42.3290000 11.0680000 0.0683888 42.1517921 11.9096159
1.3380000 43.7300000 11.2070000 1.8966859 43.7434195 12.1428361
2.0156868 43.9297695 10.3769802 0.2560000 44.8000000 11.2470000
-0.8340000 44.3780000 11.6580000 0.6210000 45.9080000 10.8290000
-0.1710000 42.1080000 9.8880000 -0.0650000 42.6850000 8.7910000
-1.1350000 41.2070000 10.0620000 -1.0932079 40.6783382 10.9215752
-2.1140000 40.8190000 9.0620000 -1.4995110 40.5714357 8.1964283
-2.8290000 39.5520000 9.5320000 -3.4554460 39.7857931 10.3928172
-3.4510659 39.1646698 8.7250864 -2.0909727 38.8009458 9.8136601
-3.0850000 41.8920000 8.6090000 -3.3640000 41.8760000 7.3920000
-3.5420000 42.7830000 9.4580000 -3.2521300 42.6884042 10.4208740
-4.4980000 43.8380000 9.1220000 -5.0481160 43.5664304 8.2210441
-5.5010000 44.0350000 10.2690000 -6.2196930 44.7624702 9.8916850
-6.0097314 43.0916045 10.4672357 -4.8870000 44.5510000 11.5430000
-3.6690000 44.7510000 11.6950000 -5.6420000 44.7420000 12.5210000
-3.8650000 45.1840000 8.7630000 -4.6090000 46.1550000 8.5770000
-2.5630000 45.2240000 8.7310000 -2.1168595 44.3495499 8.9684764
-1.7330000 46.3420000 8.4050000 -0.7063512 46.0069440 8.5527463
-1.9006548 46.5419150 7.3466876 -1.9500000 47.6070000 9.1810000
-1.6660000 48.6770000 8.6140000 -2.3820000 47.4820000 10.4280000
-2.7084892 46.5800094 10.7440966 -2.6080000 48.6450000 11.2920000
-2.8114059 49.5076689 10.6575493 -3.7720000 48.4670000 12.2600000
-4.0728627 49.4481712 12.6272654 -4.6136767 47.9958820 11.7523306
-3.3770000 47.6060000 13.4450000 -2.4290000 46.7930000 13.3250000
-4.0740000 47.7460000 14.5520000 -3.8432575 47.1914023 15.3639600
-4.8365902 48.4074265 14.5847403 -1.3250000 49.1020000 11.9560000
-1.3830000 50.1180000 12.6880000 -0.2240000 48.4060000 11.7800000
-0.2665681 47.5582621 11.2326173 1.0530000 48.8000000 12.3370000
1.7893039 48.5123646 11.5865191 1.0086723 49.8870455 12.4038314
1.4740000 48.2400000 13.6810000 2.6040000 48.5630000 14.0890000
0.5990000 47.4510000 14.3060000 -0.2939121 47.1959440 13.9088178
0.9640000 46.8930000 15.6180000 2.0433582 46.9844996 15.7392963
0.2560000 47.6550000 16.8070000 0.4982517 47.1048902 17.7162819
0.6630000 49.1290000 16.9760000 0.4206141 49.6798859 16.0672242
0.1221278 49.5616090 17.8176685 1.7350280 49.1913148 17.1630118
-1.1570000 47.6240000 16.5260000 -1.6324216 48.0779383 17.2256533
0.4520000 45.4390000 15.5580000 -0.3910000 45.0990000 14.7240000
0.9940000 44.7340000 16.5580000 1.8161315 44.9826462 17.0893890
0.5130000 43.3310000 16.6210000 -0.3816466 43.0236638 16.0794733
1.6520000 42.3180000 16.3520000 1.3189825 41.7845216 15.4617189
3.0940000 42.8030000 16.1480000 3.4279009 43.3362088 17.0381119
3.7446931 41.9462510 15.9728418 3.1348291 43.4714304 15.2879808
1.5630000 41.0660000 17.2800000 2.5489704 40.6065641 17.3498694
1.2392446 41.3827483 18.2714396 0.5560000 40.0400000 16.7130000
0.8793576 39.7224152 15.7216980 0.5077359 39.1742922 17.3735451
-0.4303684 40.4986006 16.6432660 -0.3670000 43.2000000 17.8720000
0.0150000 43.6420000 18.9550000 -1.4990000 42.5920000 17.6430000
-1.7110279 42.3732543 16.6800390 -2.5170000 42.3090000 18.6610000
-2.2613917 42.8830088 19.5516612 -3.9290000 42.6510000 18.2560000
-4.5588318 42.4626868 19.1254534 -3.8766979 43.7239870 18.0714375
-4.5160000 41.9640000 17.0540000 -4.0240000 40.9170000 16.5730000
-5.5540000 42.4400000 16.5610000 -2.3650000 40.8170000 19.0550000
-1.5510000 40.1260000 18.4370000 -3.1150000 40.4010000 20.0580000
-3.6505254 41.0491054 20.6177064 -3.0210000 38.9840000 20.5130000
-1.9854985 38.7282991 20.7376172 -3.8620000 38.7930000 21.7870000
-3.5535846 39.4917336 22.5646571 -4.9306561 38.9035154 21.6029876
-3.6610000 37.4020000 22.3500000 -2.3950000 36.9960000 22.7580000
-1.5389113 37.6684638 22.7031870 -2.2320000 35.6910000 23.2480000
-1.2558779 35.3292977 23.5711993 -3.3480000 34.8500000 23.3180000
-3.1908190 33.8473042 23.7154860 -4.6230000 35.2460000 22.9050000
-5.4903283 34.5870160 22.9447715 -4.7580000 36.5450000 22.4270000
-5.7337679 36.9084922 22.1047416 -3.3070000 37.9470000 19.4450000
-2.5520000 36.9610000 19.2560000 -4.4260000 38.0420000 18.7630000
-5.4080000 39.1300000 18.9450000 -5.0310577 40.0512827 18.5008683
-5.5959046 39.2894819 20.0067704 -6.6180000 38.5900000 18.2210000
-7.1664490 39.4224391 17.7801509 -7.2568894 38.0732794 18.9371853
-6.1210000 37.6480000 17.1730000 -5.9461591 38.1792382 16.2374169
-6.8203147 36.8298953 17.0004797 -4.8150000 37.1060000 17.7080000
-5.0190985 36.1363826 18.1621851 -3.7130000 36.9560000 16.6650000
-3.4270000 35.8250000 16.2350000 -3.0710000 38.0410000 16.2750000
-3.3898602 38.9199052 16.6570380 -1.9810000 38.0650000 15.2770000
-2.2984313 37.5007550 14.4000940 -1.6120000 39.4810000 14.8390000
-1.5756776 40.0513214 15.7671775 -0.6052708 39.3855059 14.4322112
-2.5040000 40.2080000 13.8280000 -2.3438978 39.7325559 12.8603115
-3.5312388 40.0377872 14.1503466 -2.2530000 41.7000000 13.7060000
-1.4700000 42.2170000 14.4380000 -2.9200000 42.2490000 12.8000000
-0.7470000 37.3510000 15.8140000 -0.1260000 36.5320000 15.1100000
-0.3870000 37.6120000 17.0610000 -0.9224354 38.3152746 17.5496851
0.7440000 36.9900000 17.7580000 1.6556967 37.1629756 17.1861716
0.9300000 37.5630000 19.1540000 1.2431540 38.5990928 19.0253702
-0.0315200 37.5355703 19.6666664 1.9230000 36.9270000 20.0730000
1.5360000 35.8880000 20.9290000 0.5065296 35.5298248 20.9300188
2.4780000 35.3100000 21.7850000 2.1850407 34.4811522 22.4294276
3.7900000 35.7900000 21.8170000 4.5125121 35.3533056 22.5064738
4.1760000 36.8290000 20.9660000 5.1977016 37.2087209 20.9721751
3.2320000 37.3810000 20.1000000 3.5177123 38.1876079 19.4248318
0.5460000 35.4670000 17.8550000 1.5200000 34.7200000 17.6360000
-0.6500000 35.1010000 18.2310000 -1.2808998 35.8462356 18.4892415
-1.0990000 33.7060000 18.4200000 -0.5070283 33.1804220 19.1692906
-2.5330000 33.8040000 18.9150000 -2.3633403 34.5386643 19.7021361
-3.0231382 34.2855066 18.0688227 -3.4060000 32.7050000 19.4380000
-4.3416278 33.1465858 19.7810793 -3.7000000 31.6570000 18.3750000
-2.7648195 31.2146090 18.0317413 -4.3359716 30.8789406 18.7972138
-4.2101408 32.1255564 17.5333874 -2.7450000 32.0620000 20.6570000
-2.6047522 32.8137549 21.4337204 -3.3812843 31.2625026 21.0365358
-1.7769239 31.6505581 20.3712671 -0.8830000 32.9180000 17.1330000
-0.1700000 31.8970000 17.1730000 -1.4580000 33.3930000 16.0430000
-2.0675224 34.1953011 16.1129624 -1.3350000 32.7230000 14.7410000
-1.5927128 31.6664337 14.8141499 -2.1650000 33.3640000 13.5620000
-2.3259343 32.5694003 12.8334318 -3.5040000 33.9720000 13.9530000
-3.3438244 34.7670454 14.6812487 -3.9892716 34.3831644 13.0678139
-4.1394211 33.2019319 14.3904194 -1.2710000 34.3300000 12.9300000
-1.7250251 34.7555836 12.1990143 0.1170000 32.5850000 14.2910000
0.4640000 31.5270000 13.7210000 0.9290000 33.5940000 14.5170000
0.5897961 34.4573011 14.9166890 2.3390000 33.5240000 14.1020000
2.3949025 33.0981081 13.1002065 2.9950000 34.8970000 14.0840000
2.6590331 35.4383892 13.1996227 2.6971849 35.4434870 14.9788509
4.5050000 34.7650000 14.0520000 4.8019601 33.7381626 13.8386445
4.9386143 35.4306087 13.3056450 5.1190000 35.2360000 15.7050000
6.5170000 34.1250000 15.9040000 7.2410915 34.3053971 15.1094893
6.9877587 34.3021433 16.8710094 6.1714507 33.0925108 15.8523961
3.1160000 32.5570000 14.9980000 3.9650000 31.7880000 14.5170000
2.8130000 32.6440000 16.2750000 2.1268704 33.3357519 16.5410918
3.4450000 31.8500000 17.3320000 4.5185715 32.0380105 17.3459848
2.8660000 32.2260000 18.6940000 2.7363349 33.3071781 18.7423751
1.8991247 31.7386539 18.8194831 3.7970000 31.7830000 19.7910000
3.4266478 32.0702633 20.7750833 3.9656480 30.7063951 19.7668774
5.3740000 32.6450000 19.4680000 6.3060000 32.1630000 20.9360000
6.3951179 31.0771864 20.9701447 7.3004570 32.6076368 20.8978690
5.7862161 32.5126027 21.8280224 3.2600000 30.3610000 17.0490000
4.1820000 29.5530000 17.2000000 2.0410000 30.0690000 16.6600000
1.3495806 30.8049244 16.6386987 1.6130000 28.7010000 16.3340000
1.8130267 28.0417278 17.1786589 0.1220000 28.6610000 16.0280000
-0.0926107 29.3093637 15.1784903 -0.1729689 27.6394045 15.7883414
-0.4373155 29.0055565 16.8977967 2.4600000 28.1840000 15.1700000
3.2140000 27.2150000 15.3210000 2.3150000 28.8640000 14.0470000
1.6195853 29.5952346 14.0046481 3.0370000 28.5170000 12.8150000
2.7308555 27.5211251 12.4946752 2.7550000 29.5040000 11.6800000
2.7912544 30.5115590 12.0942580 3.5406484 29.3903393 10.9330477
1.3930000 29.2710000 11.0260000 1.1896752 28.2001737 11.0355335
0.6382752 29.7888838 11.6178521 1.3530000 29.7710000 9.6240000
2.2853781 29.4823606 9.1387605 0.5165736 29.2871750 9.1196156
1.1950000 31.2140000 9.5100000 2.0402889 31.7559472 9.4010079
0.0250000 31.8560000 9.5300000 -1.1240000 31.1950000 9.6440000
-1.1184480 30.1876771 9.7172792 -1.9965222 31.7035690 9.6567368
-0.0110000 33.1900000 9.4450000 0.8503955 33.7113265 9.3655232
-0.8994213 33.6701677 9.4606653 4.5230000 28.3770000 13.0940000
5.1770000 27.4790000 12.5260000 5.0380000 29.2340000 13.9550000
4.4689045 29.9577865 14.3701659 6.4590000 29.1810000 14.3060000
7.0273061 29.1604095 13.3761057 6.9530000 30.3450000 15.1520000
6.6315203 31.2648667 14.6635378 6.4750346 30.2665599 16.1284713
8.4760000 30.3660000 15.3290000 8.7022319 30.6903673 16.3447292
8.8562798 29.3557734 15.1775682 9.1340000 31.3060000 14.3410000
9.1644809 32.2898693 14.8091572 10.1511710 30.9456415 14.1873632
8.4300000 31.4040000 13.0120000 7.3749346 31.6073838 13.1952283
8.8666652 32.2302269 12.4509547 8.5630000 30.1510000 12.2300000
8.1580512 29.3853497 12.7495150 8.0809714 30.2517421 11.3481840
9.5405535 29.9625850 12.0597412 6.7500000 27.8840000 15.0730000
7.7960000 27.2570000 14.8140000 6.0009031 27.5098717 15.9726230

View file

@ -0,0 +1,17 @@
initial minimization prior to MD
&cntrl
imin = 1,
maxcyc = 3000,
ncyc = 2000,
ntb = 0,
igb = 6,
cut = 500,
ntmin = 1,
drms = 0.000001,
scee = 1.2,
scnb = 2.0,
nsnb = 1,
ntc=3,
ntf=1
/

File diff suppressed because it is too large Load diff

52
QA/amber-tests/big/README Normal file
View file

@ -0,0 +1,52 @@
This is a series of amino acids
NWChem input files:
big.nw - nwchem input deck
big.pdb - pdb file with geometry
THR_N.sgm - amber95 segment file: one in SVN does not match
amber input files:
big_am.in - amber input deck
big_am.top - amber topology file
big_am.cor - coordinates
NWChem output files:
big.seq
big.nam
VAL.mat
THR.mat
SER.mat
PHE.mat
MET.mat
LEU.mat
ILE.mat
GLU.mat
big.top
big_md.rst
big_md.cmd
big.bsg
ASP.mat
ASN.mat
ARG.mat
ALA.mat
big_energy.out
big.out
big_md.qrs
big_md.out
big.db
amber output files:
mdinfo
big_am.rst
big_am.out
results from *.out files:
4.269883063E+03 kJ/mol for nwchem single point (1.0205d+03 kcal/mol)
-1.36273E+03 for step 500 of nwchem (-3.2570E+02 kcal/mol)
-2.19845E+03 for step 1500 of nwchem (-5.2544E+02 kcal/mol)
-2.63110E+03 for step 3000 of nwchem (-6.2885E+02 kcal/mol)
1.0211E+03 kcal/mol for amber single point
-4.6434E+02 kcal/mol for amber 500 steps
-6.4649E+02 kcal/mol for amber 1500 steps
-7.5345E+02 kcal/mol for amber 3000 steps
Note that this is an unzipping of the amino acid chain, so the energy just keeps dropping

View file

@ -0,0 +1,185 @@
# This file matches the amber force field
# 07/31/09 13:21:28
4.600000
16 15 26 33 0 0 1 1
0.000000
1 N 0 0 0 1 1
N3 0.181200 0.000000
22H 0 0 0 1 1
H 0.193400 0.000000
33H 0 0 0 1 1
H 0.193400 0.000000
44H 0 0 0 1 1
H 0.193400 0.000000
5 CA 0 0 0 1 1
CT 0.003400 0.000000
6 HA 0 0 0 1 1
HP 0.108700 0.000000
7 CB 0 0 0 1 1
CT 0.451400 0.000000
8 HB 0 0 0 1 1
H1 -0.032300 0.000000
9 CG2 0 0 0 1 1
CT -0.255400 0.000000
102HG2 0 0 0 1 1
HC 0.062700 0.000000
113HG2 0 0 0 1 1
HC 0.062700 0.000000
124HG2 0 0 0 1 1
HC 0.062700 0.000000
13 OG1 0 0 0 1 1
OH -0.676400 0.000000
14 HG1 0 0 0 1 1
HO 0.407000 0.000000
15 C 2 1 0 1 1
C 0.616300 0.000000
16 O 0 0 0 1 1
O -0.572200 0.000000
1 1 2 0 0
0.000000 0.00000E+00
2 1 3 0 0
0.000000 0.00000E+00
3 1 4 0 0
0.000000 0.00000E+00
4 1 5 0 0
0.000000 0.00000E+00
5 5 6 0 0
0.000000 0.00000E+00
6 5 7 0 0
0.000000 0.00000E+00
7 5 15 0 0
0.000000 0.00000E+00
8 7 8 0 0
0.000000 0.00000E+00
9 7 9 0 0
0.000000 0.00000E+00
10 7 13 0 0
0.000000 0.00000E+00
11 9 10 0 0
0.000000 0.00000E+00
12 9 11 0 0
0.000000 0.00000E+00
13 9 12 0 0
0.000000 0.00000E+00
14 13 14 0 0
0.000000 0.00000E+00
15 15 16 0 0
0.000000 0.00000E+00
1 2 1 3 0 0
0.000000 0.00000E+00
2 2 1 4 0 0
0.000000 0.00000E+00
3 2 1 5 0 0
0.000000 0.00000E+00
4 3 1 4 0 0
0.000000 0.00000E+00
5 3 1 5 0 0
0.000000 0.00000E+00
6 4 1 5 0 0
0.000000 0.00000E+00
7 1 5 6 0 0
0.000000 0.00000E+00
8 1 5 7 0 0
0.000000 0.00000E+00
9 1 5 15 0 0
0.000000 0.00000E+00
10 6 5 7 0 0
0.000000 0.00000E+00
11 6 5 15 0 0
0.000000 0.00000E+00
12 7 5 15 0 0
0.000000 0.00000E+00
13 5 7 8 0 0
0.000000 0.00000E+00
14 5 7 9 0 0
0.000000 0.00000E+00
15 5 7 13 0 0
0.000000 0.00000E+00
16 8 7 9 0 0
0.000000 0.00000E+00
17 8 7 13 0 0
0.000000 0.00000E+00
18 9 7 13 0 0
0.000000 0.00000E+00
19 7 9 10 0 0
0.000000 0.00000E+00
20 7 9 11 0 0
0.000000 0.00000E+00
21 7 9 12 0 0
0.000000 0.00000E+00
22 10 9 11 0 0
0.000000 0.00000E+00
23 10 9 12 0 0
0.000000 0.00000E+00
24 11 9 12 0 0
0.000000 0.00000E+00
25 7 13 14 0 0
0.000000 0.00000E+00
26 5 15 16 0 0
0.000000 0.00000E+00
1 2 1 5 6 0 0
0 0.000000 0.00000E+00
2 2 1 5 7 0 0
0 0.000000 0.00000E+00
3 2 1 5 15 0 0
0 0.000000 0.00000E+00
4 3 1 5 6 0 0
0 0.000000 0.00000E+00
5 3 1 5 7 0 0
0 0.000000 0.00000E+00
6 3 1 5 15 0 0
0 0.000000 0.00000E+00
7 4 1 5 6 0 0
0 0.000000 0.00000E+00
8 4 1 5 7 0 0
0 0.000000 0.00000E+00
9 4 1 5 15 0 0
0 0.000000 0.00000E+00
10 1 5 7 8 0 0
0 0.000000 0.00000E+00
11 1 5 7 9 0 0
0 0.000000 0.00000E+00
12 1 5 7 13 0 0
0 0.000000 0.00000E+00
13 6 5 7 8 0 0
0 0.000000 0.00000E+00
14 6 5 7 9 0 0
0 0.000000 0.00000E+00
15 6 5 7 13 0 0
0 0.000000 0.00000E+00
16 15 5 7 8 0 0
0 0.000000 0.00000E+00
17 15 5 7 9 0 0
0 0.000000 0.00000E+00
18 15 5 7 13 0 0
0 0.000000 0.00000E+00
19 1 5 15 16 0 0
0 0.000000 0.00000E+00
20 6 5 15 16 0 0
0 0.000000 0.00000E+00
21 7 5 15 16 0 0
0 0.000000 0.00000E+00
22 5 7 9 10 0 0
0 0.000000 0.00000E+00
23 5 7 9 11 0 0
0 0.000000 0.00000E+00
24 5 7 9 12 0 0
0 0.000000 0.00000E+00
25 8 7 9 10 0 0
0 0.000000 0.00000E+00
26 8 7 9 11 0 0
0 0.000000 0.00000E+00
27 8 7 9 12 0 0
0 0.000000 0.00000E+00
28 13 7 9 10 0 0
0 0.000000 0.00000E+00
29 13 7 9 11 0 0
0 0.000000 0.00000E+00
30 13 7 9 12 0 0
0 0.000000 0.00000E+00
31 5 7 13 14 0 0
0 0.000000 0.00000E+00
32 8 7 13 14 0 0
0 0.000000 0.00000E+00
33 9 7 13 14 0 0
0 0.000000 0.00000E+00

32
QA/amber-tests/big/big.nw Normal file
View file

@ -0,0 +1,32 @@
#
# This is a AMBER test that uses a PDB geometry, and the *.frg files from NWChem to asign atom types
#
Title "big"
memory 3500 Mb
permanent_dir .
echo
prepare
system big
source big.pdb
new_top new_seq
new_rst
amber
vacuo
end
print debug
task prepare
md
system big_md
sd 3000 init 0.001 min 0.000000004 max 0.05
noshake solvent solute
cutoff short 500 long 500
end
task md energy
task shell "cp big_md.out big_energy.out"
task md optimize

1078
QA/amber-tests/big/big.pdb Normal file

File diff suppressed because it is too large Load diff

View file

@ -0,0 +1,540 @@
1076
-7.1470000 23.8710000 25.0130000 -7.2952944 23.9901834 26.0049199
-6.2946862 23.3475465 24.8727927 -8.0087968 23.5550476 24.5915977
-6.8910000 25.1930000 24.4280000 -5.8850703 25.2745585 24.0162320
-7.9230000 25.6260000 23.3230000 -7.8345081 26.7056384 23.2019616
-7.7040000 24.9430000 21.9740000 -7.7928053 23.8632894 22.0941593
-8.4534963 25.2918616 21.2636128 -6.7092326 25.1858804 21.6004395
-9.2380000 25.3860000 23.9080000 -9.9185752 25.6325634 23.2774255
-6.8010000 26.2280000 25.5430000 -5.8290000 26.9990000 25.5610000
-7.7810000 26.2170000 26.4190000 -8.5466857 25.5741995 26.2753607
-7.8580000 27.1540000 27.5570000 -8.1329647 28.1400214 27.1824931
-8.8810000 26.7260000 28.5920000 -9.7307785 26.3972520 27.9937521
-8.4181838 25.8626578 29.0700601 -9.3580000 27.6990000 29.6570000
-8.4615931 28.0538152 30.1655880 -9.8244452 28.5273104 29.1236802
-10.3220000 27.1490000 30.6680000 -10.7750000 26.0340000 30.3340000
-10.6480000 27.7010000 31.7060000 -6.5020000 27.3430000 28.2240000
-6.0400000 28.4830000 28.4230000 -5.8780000 26.2240000 28.5550000
-6.3422524 25.3435754 28.3834728 -4.5540000 26.2010000 29.1940000
-4.5434862 26.8959116 30.0336943 -4.2200000 24.8180000 29.7070000
-4.9839982 24.6501070 30.4660902 -4.4211444 24.1841216 28.8433802
-2.8590000 24.4770000 30.2770000 -2.1496331 24.5751268 29.4552515
-2.6308417 25.2174871 31.0436302 -2.7550000 23.0960000 30.8670000
-3.0960000 22.0690000 30.3040000 -2.2840000 23.1150000 32.0290000
-3.4880000 26.6540000 28.1940000 -2.4920000 27.2840000 28.5750000
-3.7420000 26.3350000 26.9340000 -4.5632789 25.7725947 26.7628262
-2.8740000 26.6630000 25.7970000 -1.8532639 26.3764321 26.0501323
-3.2730000 25.9550000 24.5100000 -4.3625372 25.9292851 24.4913794
-2.9105382 26.5718993 23.6877125 -2.7310000 24.5440000 24.3610000
-1.6569867 24.5496612 24.5469148 -3.2233721 23.9090851 25.0975832
-2.9800000 23.9700000 22.9850000 -4.0820000 23.8580000 22.4520000
-1.8930000 23.5810000 22.3350000 -1.9761311 23.1864487 21.4089763
-0.9862449 23.6808062 22.7685151 -2.8380000 28.1870000 25.5880000
-1.7320000 28.7490000 25.5180000 -4.0170000 28.7720000 25.4880000
-4.8578319 28.2179494 25.4097118 -4.1180000 30.2260000 25.2840000
-3.4600921 30.5740649 24.4876915 -5.5690000 30.7200000 25.0110000
-5.4951060 31.6448304 24.4388734 -6.4410000 29.6970000 24.2310000
-6.5157420 28.7720090 24.8027568 -7.4378377 30.1105520 24.0780621
-5.9818748 29.4899270 23.2643439 -6.2490000 31.1530000 26.3400000
-6.2198611 30.2967482 27.0138579 -5.6550142 31.9650681 26.7593157
-7.7190000 31.6240000 26.1650000 -8.3137428 30.8122554 25.7461306
-8.1272405 31.9101133 27.1343184 -7.7488960 32.4805745 25.4915861
-3.4050000 30.9160000 26.4440000 -2.6860000 31.9090000 26.2400000
-3.5680000 30.3860000 27.6410000 -4.2157133 29.6183928 27.7475173
-2.9330000 30.9330000 28.8460000 -3.3394863 31.9323099 29.0017169
-3.2810000 30.0700000 30.0580000 -2.8727698 29.0686695 29.9209405
-2.8545399 30.5158083 30.9566028 -4.3642872 30.0090282 30.1622613
-1.4270000 31.0940000 28.6770000 -0.8490000 32.1200000 29.0600000
-0.7870000 30.0840000 28.1070000 -1.2993264 29.2576866 27.8334544
0.6700000 30.1020000 27.9000000 1.1638343 30.3981376 28.8254887
1.1910000 28.7180000 27.4870000 0.7371114 28.4606583 26.5299944
2.2722547 28.7883350 27.3685037 0.8590000 27.6300000 28.5200000
1.3569297 27.9177537 29.4459399 -0.2209283 27.6562455 28.6654830
1.2770000 26.2450000 28.1390000 2.4510000 26.0280000 28.5170000
0.5840000 25.4400000 27.5470000 1.0980000 31.1160000 26.8520000
2.1830000 31.7040000 26.9570000 0.2230000 31.2830000 25.8790000
-0.5972086 30.6953200 25.8343969 0.4620000 32.2360000 24.7770000
1.5117698 32.1719031 24.4906886 -0.3930000 31.9690000 23.5500000
-1.3988517 31.7051427 23.8767123 -0.4338604 32.8751424 22.9455693
0.0660000 30.8670000 22.6410000 1.3280000 30.9200000 22.0460000
1.9859561 31.7643528 22.2515773 1.7520000 29.8910000 21.1850000
2.7360664 29.9255917 20.7175415 0.8930000 28.8130000 20.9310000
1.2186363 28.0152960 20.2633705 -0.3690000 28.7550000 21.5240000
-1.0396694 27.9196850 21.3226260 -0.7670000 29.7740000 22.3770000
-1.7451538 29.7327705 22.8561818 0.2710000 33.6720000 25.2990000
1.0020000 34.5320000 24.8260000 -0.6590000 33.8570000 26.2030000
-1.1889504 33.0459457 26.4883819 -0.9870000 35.1590000 26.8170000
-1.2154939 35.9157294 26.0664999 -2.2170000 34.9300000 27.6610000
-2.8372107 34.2830112 27.0406426 -1.8416226 34.3639684 28.5135257
-3.0410000 36.0940000 28.1540000 -2.5671267 36.5361902 29.0303634
-3.1245080 36.8453666 27.3687760 -4.4420000 35.5890000 28.5300000
-4.3468270 34.6503177 29.0758174 -4.9284546 36.3298969 29.1644551
-5.2760000 35.3660000 27.2830000 -5.6469602 36.3283130 26.9302356
-4.6477231 34.9155513 26.5145839 -6.4240000 34.4690000 27.5720000
-7.0065357 34.8864127 28.2837011 -6.9656059 34.3335555 26.7303244
-6.0804657 33.5770978 27.8984870 0.1920000 35.6210000 27.6560000
0.6130000 36.7880000 27.7760000 0.7740000 34.6060000 28.3030000
0.3905739 33.6728934 28.2540534 1.9480000 34.8320000 29.1540000
1.7002357 35.6194119 29.8658248 2.4020000 33.6050000 29.9260000
2.0895570 32.7513289 29.3245649 3.4905710 33.6583681 29.9422819
1.8810000 33.4420000 31.3470000 2.4357002 34.1636184 31.9467295
0.8264939 33.7174548 31.3314828 2.0320000 32.0700000 31.9470000
3.1850000 31.6130000 31.7950000 1.1110000 31.4900000 32.4970000
3.0920000 35.3460000 28.2880000 3.7460000 36.3130000 28.6840000
3.3120000 34.7090000 27.1510000 2.7882604 33.8653449 26.9664945
4.3750000 35.0890000 26.2100000 5.3128376 35.1331221 26.7637265
4.4990000 34.0200000 25.1230000 3.5706476 33.9702961 24.5539734
5.3212291 34.2747119 24.4542987 4.6936105 33.0521390 25.5850294
4.1460000 36.5060000 25.6650000 5.1080000 37.2850000 25.6350000
2.9340000 36.8130000 25.2870000 2.2645321 36.0592810 25.3488043
2.4820000 38.1060000 24.7440000 3.0646038 38.3254547 23.8492868
1.0020000 38.1230000 24.3780000 0.8294572 37.3041294 23.6795881
0.4132467 37.9600628 25.2807298 0.4630000 39.3530000 23.7340000
0.5340000 39.5400000 22.3650000 0.9821541 38.7505500 21.7616538
0.0500000 40.7030000 21.7500000 0.0914942 40.8622831 20.6724993
-0.5060000 41.6860000 22.5840000 -0.8822032 42.5909965 22.1069756
-0.5890000 41.5410000 23.9750000 -1.0098747 42.3213205 24.6090857
-0.1150000 40.3590000 24.5450000 -0.1884435 40.2054087 25.6216224
2.7710000 39.2120000 25.7550000 3.3410000 40.2480000 25.3760000
2.3980000 38.9680000 27.0170000 1.8849722 38.1377501 27.2769747
2.6200000 40.0030000 28.0440000 2.1258918 40.9265118 27.7422018
2.0540000 39.5650000 29.4030000 2.6079684 38.7161387 29.8038152
2.0773809 40.3850272 30.1207108 0.6940000 39.1780000 29.1640000
0.2925768 38.8938291 29.9884429 4.0650000 40.4240000 28.1640000
4.2810000 41.5960000 28.5420000 5.0170000 39.5370000 27.8630000
4.7950387 38.5902288 27.5901321 6.4290000 39.9160000 27.9720000
6.6568737 40.2769909 28.9749249 7.2920000 38.6690000 27.6890000
7.2337959 38.0745335 28.6007681 6.7955835 38.1375481 26.8770705
8.7570000 38.9000000 27.3250000 8.8029770 39.7503783 26.6446738
9.6280000 39.1780000 28.5450000 9.5828944 38.3278105 29.2256189
10.6590326 39.3346095 28.2279102 9.2656828 40.0706203 29.0549565
9.2880000 37.6530000 26.6110000 8.7055910 37.4756998 25.7068629
10.3341965 37.8038745 26.3448966 9.2025703 36.7908889 27.2724871
6.7600000 41.0540000 26.9940000 7.5800000 41.9470000 27.2900000
6.0910000 41.0070000 25.8520000 5.4307525 40.2667025 25.6619154
6.3260000 42.0200000 24.8060000 7.3684392 42.3377761 24.8269289
6.0180000 41.4410000 23.4200000 5.0388943 40.9626546 23.4452378
6.0136512 42.2373766 22.6757787 6.9850000 40.4100000 22.9370000
8.1600000 40.7720000 22.2780000 8.3954353 41.8236562 22.1146311
9.0360000 39.7820000 21.8270000 9.9567870 40.0267817 21.2975344
8.7050000 38.4490000 22.0710000 9.4067756 37.6949061 21.7146948
7.5370000 38.0590000 22.7380000 7.3047022 37.0115918 22.9305459
6.6500000 39.0670000 23.1640000 5.7185119 38.7983747 23.6622672
5.5010000 43.3130000 24.9660000 6.0550000 44.3900000 24.6970000
4.2660000 43.1210000 25.3760000 3.9263426 42.1758088 25.4825166
3.3040000 44.2450000 25.5730000 3.4293272 45.0030645 24.7998696
1.8920000 43.6880000 25.4830000 1.8101074 43.3303314 24.4566151
1.7314459 42.8575425 26.1705047 0.8380000 44.7660000 25.7320000
1.3620000 45.9000000 25.6800000 -0.3270000 44.4680000 26.0380000
3.6700000 44.9290000 26.9000000 2.8890000 44.7430000 27.8230000
4.7280000 45.6930000 26.9420000 5.2193341 45.7357518 26.0606020
5.2930000 46.4000000 28.0820000 5.5879603 45.6460336 28.8118166
6.5030000 47.2920000 27.7830000 6.1870956 47.9567022 26.9789635
6.6711418 47.8734652 28.6894913 7.7990000 46.5880000 27.3820000
7.7433488 46.2989768 26.3324909 8.6390003 47.2672771 27.5272040
7.9970000 45.3450000 28.2410000 7.6903950 45.5856862 29.2589207
7.3566998 44.5562019 27.8461422 9.4320000 44.8730000 28.2460000
9.8710370 45.0204001 27.2592780 10.0019493 45.4366840 28.9845915
9.4660000 43.4200000 28.5940000 8.9382879 42.8974492 27.9094854
10.4230759 43.0973880 28.5992206 9.0596122 43.2832510 29.5084658
4.2940000 47.2800000 28.8280000 4.3150000 47.3110000 30.0600000
3.4440000 47.9600000 28.0720000 3.5341250 48.0617198 27.0711854
2.4670000 48.8110000 28.7760000 2.7003921 48.8893755 29.8378308
2.5490000 50.2110000 28.1510000 2.1072477 50.9568183 28.8118397
3.5933687 50.4627538 27.9665711 1.7910000 50.1910000 26.8300000
1.4800000 49.1260000 26.2870000 1.4870000 51.2940000 26.3370000
1.1140000 48.1780000 28.8810000 0.1890000 48.8260000 29.3960000
0.9010000 46.9280000 28.4700000 1.6406840 46.3673524 28.0717020
-0.4170000 46.3410000 28.5890000 -0.2759187 45.3145075 28.2506091
-0.6079840 46.3462481 29.6621248 -1.6030000 46.9320000 27.8780000
-2.7250000 46.7820000 28.4030000 -1.4540000 47.5460000 26.7100000
-0.5276198 47.7736193 26.3781741 -2.6160000 48.0790000 25.9660000
-3.4504435 48.2676213 26.6414445 -2.2130000 49.4070000 25.2990000
-3.0941511 50.0152164 25.0946844 -1.5244527 49.9669124 25.9318513
-1.5220000 49.0540000 23.9830000 -0.6470000 48.1980000 24.1440000
-1.7960000 49.5540000 22.9200000 -3.2190000 47.0540000 25.0400000
-4.2130000 47.3330000 24.3510000 -2.6310000 45.8600000 24.9310000
-1.7937363 45.6465941 25.4540169 -3.1850000 44.8220000 24.0470000
-3.0486580 43.8557301 24.5326259 -4.2498607 45.0177026 23.9210175
-2.5100000 44.7850000 22.6510000 -3.0400000 44.0510000 21.7800000
-1.4680000 45.5220000 22.5170000 -1.1664003 46.1374373 23.2588719
-0.7420000 45.5700000 21.2250000 -0.8399244 44.6081876 20.7215848
-1.2500000 46.7020000 20.2520000 -0.7118239 46.5492103 19.3165204
-2.7520000 46.7520000 19.9730000 -3.2908520 46.9052421 20.9080160
-2.9681714 47.5741142 19.2907264 -3.0695956 45.8123482 19.5210176
-0.8240000 47.9600000 20.8890000 -1.1013569 48.7027124 20.3476575
0.7320000 45.7030000 21.4940000 1.2300000 46.2590000 22.4720000
1.4870000 45.1760000 20.4980000 1.0487052 44.7120341 19.7152276
2.9480000 45.2260000 20.5450000 3.2927821 45.6756274 21.4761604
3.6310000 43.8280000 20.3560000 3.4213691 43.4301437 19.3630923
5.1700000 44.0380000 20.5180000 5.3804681 44.4354312 21.5109009
5.6825375 43.0841965 20.3928404 5.5226892 44.7411541 19.7634910
3.0610000 42.7720000 21.3220000 3.4990538 42.9604110 22.3021587
1.9830975 42.9264891 21.3705680 3.3570000 41.2980000 20.8940000
4.4347835 41.1426106 20.8456600 2.9240844 40.6134118 21.6233986
2.9188288 41.1086894 19.9140677 3.3730000 46.1960000 19.4300000
2.9680000 46.0580000 18.2830000 4.1300000 47.1870000 19.7840000
4.4578927 47.2108641 20.7389954 4.6140000 48.2310000 18.8480000
4.4416071 48.6362990 17.8509480 4.3970000 49.4070000 19.9260000
3.6976466 49.0691218 20.6907505 5.6970000 49.7590000 20.5630000
6.3967706 50.0974334 19.7988772 5.5416789 50.5552591 21.2909739
6.1050029 48.8822004 21.0658487 3.3280000 50.2610000 19.4650000
3.1872212 50.9679928 20.0989896 6.0660000 48.0390000 18.4870000
6.7090000 47.1940000 19.1490000 6.6250000 48.7930000 17.5510000
6.0615290 49.4637404 17.0483002 8.0420000 48.6610000 17.1860000
8.2457151 47.6409129 16.8603904 8.4750000 49.6540000 16.0360000
9.5427527 49.5091862 15.8715818 7.7300000 49.4510000 14.7290000
6.6621536 49.5961173 14.8925376 8.0853423 50.1711737 13.9919920
7.9068130 48.4397407 14.3626772 8.2360000 51.0000000 16.5430000
8.4857601 51.6426002 15.8749542 8.9720000 48.8980000 18.3700000
10.0100000 48.2160000 18.4770000 8.6200000 49.8650000 19.2040000
7.7717407 50.3747804 19.0023085 9.3950000 50.2470000 20.3860000
10.4308228 50.3724346 20.0706643 8.9370000 51.5280000 21.0760000
8.2371571 51.2283151 21.8560693 9.8227235 51.9691485 21.5331455
8.2710000 52.5550000 20.1620000 8.6228345 53.5554102 20.4139779
8.5146751 52.3358367 19.1224367 6.7550000 52.4840000 20.3560000
6.4870885 51.4546526 20.5942568 6.4848006 53.1315829 21.1901032
6.0090000 52.9250000 19.1120000 4.9614904 52.6374860 19.2023324
6.0809586 54.0084672 19.0170267 6.5970000 52.2820000 17.9070000
6.5300323 51.2780381 17.9946099 6.0927912 52.5826221 17.0851120
7.5675309 52.5486303 17.8228739 9.3590000 49.1630000 21.4490000
10.3710000 48.9180000 22.1160000 8.2200000 48.5160000 21.6010000
7.3706318 48.7590600 21.1115157 8.1600000 47.4340000 22.6290000
8.6362109 47.7829374 23.5452778 6.7230000 47.0390000 22.9470000
6.2081555 46.8944370 21.9971916 6.7582628 46.0934837 23.4881612
5.9710000 48.0750000 23.7810000 6.4722400 48.0394485 24.7482618
6.1727312 49.0335229 23.3028395 4.5170000 47.9190000 23.9750000
3.8300000 47.4500000 23.0670000 3.9000000 48.2470000 24.9960000
8.9910000 46.2690000 22.1360000 9.8140000 45.6820000 22.8800000
8.7980000 45.8960000 20.8860000 8.0650008 46.3224684 20.3374235
9.5710000 44.7760000 20.2820000 9.3974131 43.8553891 20.8391735
9.1200000 44.5610000 18.8210000 8.0825974 44.2294929 18.8657011
9.1707349 45.5275830 18.3197583 9.9530000 43.5330000 18.0370000
10.9942312 43.8549877 18.0208013 9.8560000 42.1680000 18.6940000
8.8150463 41.8451171 18.7101507 10.4500498 41.4493893 18.1293725
10.2331529 42.2279525 19.7149122 9.3550000 43.5290000 16.6210000
9.4375806 44.5260385 16.1883406 9.8977874 42.8174042 15.9988147
8.3049703 43.2407174 16.6702918 11.0360000 45.1150000 20.4430000
11.8210000 44.2260000 20.7670000 11.3990000 46.3700000 20.2520000
10.6638686 47.0183286 20.0083790 12.7700000 46.8830000 20.3660000
13.3936257 46.3743074 19.6308660 12.7595595 47.9523569 20.1551299
13.3560000 46.6570000 21.7550000 14.4850000 46.1500000 21.9040000
12.6240000 47.0050000 22.8050000 11.8363664 47.6260262 22.6864285
13.1530000 46.7970000 24.1690000 14.2381405 46.8567258 24.2526785
12.5100000 47.7260000 25.2720000 13.3168588 48.2756585 25.7567213
11.5050000 48.7040000 24.7020000 10.6975457 48.1550331 24.2174874
11.0962617 49.3158320 25.5061727 11.9964922 49.3463652 23.9713134
12.0100000 46.9180000 26.3790000 11.6213860 47.4913887 27.0436841
13.2680000 45.3120000 24.4460000 14.2690000 44.9430000 25.1200000
12.3610000 44.4970000 23.9230000 11.4984580 44.7926706 23.4886044
12.5570000 43.0500000 24.1840000 12.8219255 42.8653490 25.2250656
11.3560000 42.1830000 23.8260000 11.0758154 42.3954076 22.7942642
11.6660000 40.6940000 23.9730000 11.9458376 40.4807247 25.0046518
10.7839525 40.1108527 23.7083646 12.4896642 40.4267414 23.3109941
10.1300000 42.5500000 24.6370000 9.8668269 43.5907198 24.4479385
9.2978421 41.9073533 24.3495663 10.3420545 42.4163368 25.6977864
13.8590000 42.5480000 23.5470000 14.6440000 41.8510000 24.2410000
14.0890000 42.8560000 22.2950000 13.3741552 43.3462853 21.7766213
15.2910000 42.4030000 21.5710000 15.3568925 41.3178980 21.6504419
15.2710000 42.7830000 20.0950000 15.3761435 43.8657943 20.0271723
16.1249066 42.3047946 19.6151425 14.0100000 42.3600000 19.4060000
13.1516776 42.7160588 19.9757417 13.9849165 42.7906195 18.4049807
13.9330000 40.5450000 19.2790000 14.8980000 40.2740000 17.7850000
15.9071919 40.6599851 17.9286878 14.9449912 39.2062711 17.5708795
14.4273058 40.7913226 16.9489824 16.5740000 42.8760000 22.2350000
17.5090000 42.0580000 22.3820000 16.6210000 44.1340000 22.5990000
15.8433992 44.7279029 22.3485689 17.7840000 44.7480000 23.2480000
18.6734757 44.5904851 22.6379827 17.5700000 46.2440000 23.4180000
16.5887329 46.3411187 23.8825235 18.3380191 46.5641260 24.1221066
17.6290000 47.1080000 22.1830000 18.4612880 46.7466522 21.5789922
16.6951838 46.9566544 21.6415355 17.8100000 48.5480000 22.4590000
18.6153312 48.6879582 23.1800772 18.0541516 49.0736381 21.5358564
16.5770000 49.1200000 23.0140000 16.4304543 49.1638046 24.0123511
15.6340000 49.6200000 22.2080000 15.8000000 49.6480000 20.8890000
16.6510771 49.2867985 20.4824274 15.0742499 50.0303423 20.2997625
14.5040000 50.0690000 22.7510000 14.3764623 50.0261604 23.7519990
13.7796800 50.4503191 22.1593455 18.1270000 44.0710000 24.5770000
19.3100000 44.0620000 24.9760000 17.1300000 43.5320000 25.2370000
16.1819017 43.6285591 24.9025075 17.3510000 42.8510000 26.5310000
18.1006952 43.4146220 27.0863257 16.1090000 42.7340000 27.3800000
16.4075238 42.5390385 28.4100353 15.5637809 43.6766452 27.3324994
15.2770000 41.6850000 26.9240000 14.4958081 41.6317037 27.4794264
17.9600000 41.4760000 26.2330000 18.4830000 40.8140000 27.1430000
17.8850000 41.1050000 24.9610000 17.3702110 41.6741689 24.3043886
18.4260000 39.8200000 24.4950000 18.5448978 39.1311775 25.3313528
17.4560000 39.2020000 23.4870000 17.3976802 39.9632300 22.7090366
17.9710446 38.3238721 23.0974869 16.0590000 38.8160000 23.9200000
15.5462519 39.7144276 24.2635375 15.2690000 38.1850000 22.7830000
15.7809236 37.2861968 22.4392182 14.2714616 37.9221799 23.1350561
15.1878641 38.8939399 21.9590330 16.1370000 37.8380000 25.0900000
16.6571006 38.3091710 25.9240220 15.1295347 37.5610792 25.4005314
16.6801254 36.9446040 24.7818485 19.8410000 39.9720000 23.9450000
20.4100000 38.9480000 23.5300000 20.4050000 41.1680000 23.9430000
19.9064009 41.9715179 24.2977637 21.7740000 41.3660000 23.4250000
22.2237548 42.2067299 23.9531979 22.3509995 40.4626543 23.6228318
21.7970000 41.6520000 21.9240000 22.8160000 41.5940000 21.2130000
20.6380000 42.0150000 21.4270000 19.8481890 42.1000016 22.0507569
20.4420000 42.3450000 19.9980000 21.3541481 42.2113277 19.4164360
19.4030000 41.3740000 19.4950000 19.3941143 40.6024764 20.2649166
18.4800343 41.9513167 19.5492195 19.5110000 40.7330000 18.1490000
18.8811776 41.2645339 17.4356303 20.5382028 40.7172244 17.7847096
19.0140000 39.3010000 18.2630000 18.3870000 38.7610000 17.3670000
19.3000000 38.7140000 19.4280000 19.0091555 37.7612667 19.5947639
19.8058293 39.2254902 20.1369541 19.9200000 43.7750000 19.9210000
18.8960000 44.0340000 20.5780000 20.6000000 44.6250000 19.1790000
21.5120573 44.3847775 18.8176912 20.1420000 46.0260000 19.0520000
19.1721218 46.1491509 19.5339433 21.1380000 46.9550000 19.7220000
22.1011591 46.8244575 19.2286619 20.7690715 47.9637383 19.5364228
21.3310000 46.7760000 21.2200000 21.6970000 45.7290000 21.7750000
21.1160000 47.8420000 22.0060000 21.2344381 47.7605592 23.0057198
20.8367051 48.7209974 21.5943462 19.7970000 46.2680000 17.5760000
20.5720000 46.8110000 16.7650000 18.5870000 45.8500000 17.2230000
17.5830000 45.2170000 18.0910000 17.4766527 45.8107647 18.9988734
17.9227872 44.2144524 18.3508883 16.3110000 45.2000000 17.2500000
15.7123333 46.0838319 17.4703143 15.7324264 44.3043085 17.4760323
16.7780000 45.1930000 15.8070000 16.0396029 45.7066193 15.1913184
16.9778681 44.2071427 15.3871933 18.0690000 46.0090000 15.8590000
18.7986935 45.6190448 15.1493659 17.8160000 47.4800000 15.5680000
17.6570000 48.3130000 16.4860000 17.7860000 47.7800000 14.2760000
17.9730368 47.0890682 13.5634483 17.5410000 49.1730000 13.8530000
17.8684358 49.8699509 14.6244561 18.3210000 49.4900000 12.5260000
18.1763448 50.5517187 12.3261790 19.8200000 49.1770000 12.5610000
19.9654504 48.1152974 12.7603272 20.2671238 49.4302702 11.5997333
20.2956525 49.7623796 13.3478837 17.6680000 48.6750000 11.4960000
18.0992062 48.8216280 10.6509191 16.0270000 49.3370000 13.7010000
15.3340000 48.3170000 13.6440000 15.5580000 50.5700000 13.6140000
16.1570402 51.3709428 13.7545017 14.1130000 50.8020000 13.4220000
13.5625563 50.3241590 14.2324190 13.7770000 52.2870000 13.4150000
14.2846932 52.6807685 14.2955074 14.2621650 52.6672836 12.5160565
12.3300000 52.7250000 13.4300000 12.3871283 53.7720854 13.1325957
11.9101289 52.1381779 12.6130255 11.4570000 52.5970000 14.6340000
12.0590000 52.6470000 15.7350000 10.2340000 52.4680000 14.5570000
13.6740000 50.0870000 12.1480000 12.6040000 49.4510000 12.1180000
14.4720000 50.0990000 11.0930000 15.3251579 50.6375123 11.1401751
14.1410000 49.4270000 9.8280000 13.2237600 49.8850254 9.4578882
15.2480000 49.6410000 8.7910000 16.1706867 49.1799022 9.1433322
14.9533634 49.1863599 7.8451713 15.4087574 50.7092258 8.6455696
13.8810000 47.9430000 9.9980000 12.9390000 47.3380000 9.4580000
14.7710000 47.3070000 10.7520000 15.5649499 47.8044131 11.1292574
14.6580000 45.8690000 11.0230000 14.4990698 45.3468547 10.0794934
15.9040000 45.2970000 11.6930000 16.1750154 45.9802643 12.4978601
15.6240451 44.3318838 12.1152273 17.1020000 45.1150000 10.7520000
16.9489246 44.1276919 10.3162346 17.0300843 45.8715698 9.9706362
18.4610000 45.1770000 11.3840000 18.7270000 45.7390000 12.4370000
19.3280000 44.6160000 10.6720000 13.4320000 45.6430000 11.9120000
12.7530000 44.6550000 11.6780000 13.1810000 46.5440000 12.8450000
13.8397939 47.2946893 12.9951848 12.0240000 46.4290000 13.7660000
12.0167921 45.4484474 14.2419872 12.1250000 47.4920000 14.8460000
12.9304203 48.0954125 14.4273225 11.1870079 48.0433789 14.7807916
12.4280000 47.1830000 16.2950000 11.4905644 47.2865091 16.8414441
13.0070000 45.8060000 16.5430000 13.9447512 45.7018990 15.9972121
13.1914111 45.6754787 17.6093290 12.3022181 45.0479910 16.2012323
13.3730000 48.2470000 16.8600000 12.9010405 49.2270166 16.7898490
13.5904257 48.0229286 17.9043272 14.3012308 48.2495884 16.2886158
10.7550000 46.4920000 12.9150000 9.8490000 45.6880000 13.0890000
10.7170000 47.3630000 11.9090000 11.4366556 48.0652258 11.8137426
9.5290000 47.4600000 11.0460000 8.6186195 47.4735814 11.6452679
9.6070000 48.7640000 10.2440000 9.8682419 49.5820301 10.9153260
10.3738501 48.6662895 9.4755642 8.2650000 49.0560000 9.5870000
8.3820560 49.9046832 8.9131042 7.9359251 48.1851840 9.0199980
7.2040000 49.3950000 10.6120000 6.1410000 48.7870000 10.6450000
7.5300000 50.3940000 11.4340000 6.8805975 50.6872866 12.1497938
8.4247036 50.8524155 11.3367265 9.2830000 46.2080000 10.2240000
8.1210000 45.7510000 10.0850000 10.3200000 45.5840000 9.6620000
11.2351532 46.0059162 9.7296800 10.1790000 44.3410000 8.8800000
9.4481679 44.4873952 8.0846686 11.5180000 43.9460000 8.2480000
12.3148984 44.1143648 8.9723655 11.4794088 42.8875505 7.9905212
11.7920000 44.7630000 6.9940000 10.9000000 45.1790000 6.2500000
13.0030000 44.9620000 6.8170000 9.6290000 43.1990000 9.7460000
8.7950000 42.4180000 9.2930000 10.1050000 43.1390000 10.9780000
10.8714537 43.7252866 11.2761881 9.6130000 42.1060000 11.9070000
9.7114786 41.1246179 11.4430062 10.4040000 42.1110000 13.2060000
10.2110669 43.0548754 13.7158777 10.0470936 41.2860137 13.8225351
11.8820000 41.9600000 12.9740000 12.0771689 41.1626214 12.2569245
12.3084779 42.8925898 12.6045478 12.6620000 41.5340000 14.5750000
14.2680000 42.3040000 14.3430000 14.1385066 43.3741826 14.1816322
14.8805628 42.1423283 15.2299779 14.7609542 41.8634437 13.4763931
8.1170000 42.2530000 12.1280000 7.3760000 41.2260000 12.0980000
7.6360000 43.4550000 12.3720000 8.2755430 44.2310828 12.4656995
6.2060000 43.6970000 12.6140000 5.7085959 42.9745032 13.2610601
6.0110000 45.1450000 13.1980000 6.7264250 45.8207170 12.7293061
4.5730000 45.6400000 12.9150000 3.8569736 44.9649959 13.3838030
4.4466363 46.6429979 13.3225882 4.4006938 45.6612915 11.8389162
6.4510000 45.0960000 14.6820000 6.5661231 44.0402357 14.9273757
7.4235170 45.5864366 14.7242157 5.4990000 45.7570000 15.6950000
4.5263461 45.2667574 15.6536901 5.9118232 45.6600198 16.6991270
5.3837382 46.8129584 15.4505284 5.4020000 43.4650000 11.3440000
4.3480000 42.8280000 11.3270000 5.9490000 43.9350000 10.2320000
6.8160469 44.4514524 10.2720732 5.2950000 43.7750000 8.9300000
4.3500404 44.3177030 8.9549789 6.1310000 44.2980000 7.7500000
7.1774104 44.0403974 7.9136013 5.7924755 43.8613362 6.8104118
6.0120000 45.8160000 7.6530000 6.7040000 46.4240000 6.8340000
5.1550000 46.3250000 8.5170000 4.9978723 47.3224405 8.5398697
4.6599149 45.7149285 9.1516687 4.9440000 42.3130000 8.6540000
3.8500000 41.9930000 8.2010000 5.8990000 41.4510000 8.8980000
6.8066044 41.7526926 9.2225537 5.6800000 40.0170000 8.6790000
5.4796299 39.7952959 7.6307642 6.9800000 39.3060000 9.0920000
7.7488843 39.7837080 8.4847841 7.1162693 39.5864363 10.1364550
7.0910000 37.8030000 8.9550000 6.2081742 37.3902016 9.4431771
7.0426139 37.5994319 7.8852720 8.3240000 37.1670000 9.5300000
9.2130000 38.0030000 9.7940000 8.4390000 35.9600000 9.7210000
4.5200000 39.3970000 9.4430000 3.8960000 38.4410000 8.9500000
4.2200000 39.9020000 10.6270000 4.8023075 40.6555000 10.9635340
3.1560000 39.3980000 11.5060000 2.7589606 38.4240464 11.2198649
3.7070000 39.3030000 12.9520000 3.0643343 38.6382584 13.5292336
5.1220000 38.7660000 12.9570000 5.7653952 39.4305070 12.3803078
5.4858682 38.7088708 13.9828834 5.1352635 37.7713254 12.5114126
3.6490000 40.6010000 13.7140000 2.6139648 40.9353169 13.7849684
4.0519974 40.4530630 14.7159019 4.2390049 41.3552542 13.1933316
1.8800000 40.2070000 11.5010000 0.8100000 39.7220000 11.8920000
1.9650000 41.4610000 11.0750000 2.8523991 41.8844408 10.8440887
0.7380000 42.3290000 11.0680000 0.0683888 42.1517921 11.9096159
1.3380000 43.7300000 11.2070000 1.8966859 43.7434195 12.1428361
2.0156868 43.9297695 10.3769802 0.2560000 44.8000000 11.2470000
-0.8340000 44.3780000 11.6580000 0.6210000 45.9080000 10.8290000
-0.1710000 42.1080000 9.8880000 -0.0650000 42.6850000 8.7910000
-1.1350000 41.2070000 10.0620000 -1.0932079 40.6783382 10.9215752
-2.1140000 40.8190000 9.0620000 -1.4995110 40.5714357 8.1964283
-2.8290000 39.5520000 9.5320000 -3.4554460 39.7857931 10.3928172
-3.4510659 39.1646698 8.7250864 -2.0909727 38.8009458 9.8136601
-3.0850000 41.8920000 8.6090000 -3.3640000 41.8760000 7.3920000
-3.5420000 42.7830000 9.4580000 -3.2521300 42.6884042 10.4208740
-4.4980000 43.8380000 9.1220000 -5.0481160 43.5664304 8.2210441
-5.5010000 44.0350000 10.2690000 -6.2196930 44.7624702 9.8916850
-6.0097314 43.0916045 10.4672357 -4.8870000 44.5510000 11.5430000
-3.6690000 44.7510000 11.6950000 -5.6420000 44.7420000 12.5210000
-3.8650000 45.1840000 8.7630000 -4.6090000 46.1550000 8.5770000
-2.5630000 45.2240000 8.7310000 -2.1168595 44.3495499 8.9684764
-1.7330000 46.3420000 8.4050000 -0.7063512 46.0069440 8.5527463
-1.9006548 46.5419150 7.3466876 -1.9500000 47.6070000 9.1810000
-1.6660000 48.6770000 8.6140000 -2.3820000 47.4820000 10.4280000
-2.7084892 46.5800094 10.7440966 -2.6080000 48.6450000 11.2920000
-2.8114059 49.5076689 10.6575493 -3.7720000 48.4670000 12.2600000
-4.0728627 49.4481712 12.6272654 -4.6136767 47.9958820 11.7523306
-3.3770000 47.6060000 13.4450000 -2.4290000 46.7930000 13.3250000
-4.0740000 47.7460000 14.5520000 -3.8432575 47.1914023 15.3639600
-4.8365902 48.4074265 14.5847403 -1.3250000 49.1020000 11.9560000
-1.3830000 50.1180000 12.6880000 -0.2240000 48.4060000 11.7800000
-0.2665681 47.5582621 11.2326173 1.0530000 48.8000000 12.3370000
1.7893039 48.5123646 11.5865191 1.0086723 49.8870455 12.4038314
1.4740000 48.2400000 13.6810000 2.6040000 48.5630000 14.0890000
0.5990000 47.4510000 14.3060000 -0.2939121 47.1959440 13.9088178
0.9640000 46.8930000 15.6180000 2.0433582 46.9844996 15.7392963
0.2560000 47.6550000 16.8070000 0.4982517 47.1048902 17.7162819
0.6630000 49.1290000 16.9760000 0.4206141 49.6798859 16.0672242
0.1221278 49.5616090 17.8176685 1.7350280 49.1913148 17.1630118
-1.1570000 47.6240000 16.5260000 -1.6324216 48.0779383 17.2256533
0.4520000 45.4390000 15.5580000 -0.3910000 45.0990000 14.7240000
0.9940000 44.7340000 16.5580000 1.8161315 44.9826462 17.0893890
0.5130000 43.3310000 16.6210000 -0.3816466 43.0236638 16.0794733
1.6520000 42.3180000 16.3520000 1.3189825 41.7845216 15.4617189
3.0940000 42.8030000 16.1480000 3.4279009 43.3362088 17.0381119
3.7446931 41.9462510 15.9728418 3.1348291 43.4714304 15.2879808
1.5630000 41.0660000 17.2800000 2.5489704 40.6065641 17.3498694
1.2392446 41.3827483 18.2714396 0.5560000 40.0400000 16.7130000
0.8793576 39.7224152 15.7216980 0.5077359 39.1742922 17.3735451
-0.4303684 40.4986006 16.6432660 -0.3670000 43.2000000 17.8720000
0.0150000 43.6420000 18.9550000 -1.4990000 42.5920000 17.6430000
-1.7110279 42.3732543 16.6800390 -2.5170000 42.3090000 18.6610000
-2.2613917 42.8830088 19.5516612 -3.9290000 42.6510000 18.2560000
-4.5588318 42.4626868 19.1254534 -3.8766979 43.7239870 18.0714375
-4.5160000 41.9640000 17.0540000 -4.0240000 40.9170000 16.5730000
-5.5540000 42.4400000 16.5610000 -2.3650000 40.8170000 19.0550000
-1.5510000 40.1260000 18.4370000 -3.1150000 40.4010000 20.0580000
-3.6505254 41.0491054 20.6177064 -3.0210000 38.9840000 20.5130000
-1.9854985 38.7282991 20.7376172 -3.8620000 38.7930000 21.7870000
-3.5535846 39.4917336 22.5646571 -4.9306561 38.9035154 21.6029876
-3.6610000 37.4020000 22.3500000 -2.3950000 36.9960000 22.7580000
-1.5389113 37.6684638 22.7031870 -2.2320000 35.6910000 23.2480000
-1.2558779 35.3292977 23.5711993 -3.3480000 34.8500000 23.3180000
-3.1908190 33.8473042 23.7154860 -4.6230000 35.2460000 22.9050000
-5.4903283 34.5870160 22.9447715 -4.7580000 36.5450000 22.4270000
-5.7337679 36.9084922 22.1047416 -3.3070000 37.9470000 19.4450000
-2.5520000 36.9610000 19.2560000 -4.4260000 38.0420000 18.7630000
-5.4080000 39.1300000 18.9450000 -5.0310577 40.0512827 18.5008683
-5.5959046 39.2894819 20.0067704 -6.6180000 38.5900000 18.2210000
-7.1664490 39.4224391 17.7801509 -7.2568894 38.0732794 18.9371853
-6.1210000 37.6480000 17.1730000 -5.9461591 38.1792382 16.2374169
-6.8203147 36.8298953 17.0004797 -4.8150000 37.1060000 17.7080000
-5.0190985 36.1363826 18.1621851 -3.7130000 36.9560000 16.6650000
-3.4270000 35.8250000 16.2350000 -3.0710000 38.0410000 16.2750000
-3.3898602 38.9199052 16.6570380 -1.9810000 38.0650000 15.2770000
-2.2984313 37.5007550 14.4000940 -1.6120000 39.4810000 14.8390000
-1.5756776 40.0513214 15.7671775 -0.6052708 39.3855059 14.4322112
-2.5040000 40.2080000 13.8280000 -2.3438978 39.7325559 12.8603115
-3.5312388 40.0377872 14.1503466 -2.2530000 41.7000000 13.7060000
-1.4700000 42.2170000 14.4380000 -2.9200000 42.2490000 12.8000000
-0.7470000 37.3510000 15.8140000 -0.1260000 36.5320000 15.1100000
-0.3870000 37.6120000 17.0610000 -0.9224354 38.3152746 17.5496851
0.7440000 36.9900000 17.7580000 1.6556967 37.1629756 17.1861716
0.9300000 37.5630000 19.1540000 1.2431540 38.5990928 19.0253702
-0.0315200 37.5355703 19.6666664 1.9230000 36.9270000 20.0730000
1.5360000 35.8880000 20.9290000 0.5065296 35.5298248 20.9300188
2.4780000 35.3100000 21.7850000 2.1850407 34.4811522 22.4294276
3.7900000 35.7900000 21.8170000 4.5125121 35.3533056 22.5064738
4.1760000 36.8290000 20.9660000 5.1977016 37.2087209 20.9721751
3.2320000 37.3810000 20.1000000 3.5177123 38.1876079 19.4248318
0.5460000 35.4670000 17.8550000 1.5200000 34.7200000 17.6360000
-0.6500000 35.1010000 18.2310000 -1.2808998 35.8462356 18.4892415
-1.0990000 33.7060000 18.4200000 -0.5070283 33.1804220 19.1692906
-2.5330000 33.8040000 18.9150000 -2.3633403 34.5386643 19.7021361
-3.0231382 34.2855066 18.0688227 -3.4060000 32.7050000 19.4380000
-4.3416278 33.1465858 19.7810793 -3.7000000 31.6570000 18.3750000
-2.7648195 31.2146090 18.0317413 -4.3359716 30.8789406 18.7972138
-4.2101408 32.1255564 17.5333874 -2.7450000 32.0620000 20.6570000
-2.6047522 32.8137549 21.4337204 -3.3812843 31.2625026 21.0365358
-1.7769239 31.6505581 20.3712671 -0.8830000 32.9180000 17.1330000
-0.1700000 31.8970000 17.1730000 -1.4580000 33.3930000 16.0430000
-2.0675224 34.1953011 16.1129624 -1.3350000 32.7230000 14.7410000
-1.5927128 31.6664337 14.8141499 -2.1650000 33.3640000 13.5620000
-2.3259343 32.5694003 12.8334318 -3.5040000 33.9720000 13.9530000
-3.3438244 34.7670454 14.6812487 -3.9892716 34.3831644 13.0678139
-4.1394211 33.2019319 14.3904194 -1.2710000 34.3300000 12.9300000
-1.7250251 34.7555836 12.1990143 0.1170000 32.5850000 14.2910000
0.4640000 31.5270000 13.7210000 0.9290000 33.5940000 14.5170000
0.5897961 34.4573011 14.9166890 2.3390000 33.5240000 14.1020000
2.3949025 33.0981081 13.1002065 2.9950000 34.8970000 14.0840000
2.6590331 35.4383892 13.1996227 2.6971849 35.4434870 14.9788509
4.5050000 34.7650000 14.0520000 4.8019601 33.7381626 13.8386445
4.9386143 35.4306087 13.3056450 5.1190000 35.2360000 15.7050000
6.5170000 34.1250000 15.9040000 7.2410915 34.3053971 15.1094893
6.9877587 34.3021433 16.8710094 6.1714507 33.0925108 15.8523961
3.1160000 32.5570000 14.9980000 3.9650000 31.7880000 14.5170000
2.8130000 32.6440000 16.2750000 2.1268704 33.3357519 16.5410918
3.4450000 31.8500000 17.3320000 4.5185715 32.0380105 17.3459848
2.8660000 32.2260000 18.6940000 2.7363349 33.3071781 18.7423751
1.8991247 31.7386539 18.8194831 3.7970000 31.7830000 19.7910000
3.4266478 32.0702633 20.7750833 3.9656480 30.7063951 19.7668774
5.3740000 32.6450000 19.4680000 6.3060000 32.1630000 20.9360000
6.3951179 31.0771864 20.9701447 7.3004570 32.6076368 20.8978690
5.7862161 32.5126027 21.8280224 3.2600000 30.3610000 17.0490000
4.1820000 29.5530000 17.2000000 2.0410000 30.0690000 16.6600000
1.3495806 30.8049244 16.6386987 1.6130000 28.7010000 16.3340000
1.8130267 28.0417278 17.1786589 0.1220000 28.6610000 16.0280000
-0.0926107 29.3093637 15.1784903 -0.1729689 27.6394045 15.7883414
-0.4373155 29.0055565 16.8977967 2.4600000 28.1840000 15.1700000
3.2140000 27.2150000 15.3210000 2.3150000 28.8640000 14.0470000
1.6195853 29.5952346 14.0046481 3.0370000 28.5170000 12.8150000
2.7308555 27.5211251 12.4946752 2.7550000 29.5040000 11.6800000
2.7912544 30.5115590 12.0942580 3.5406484 29.3903393 10.9330477
1.3930000 29.2710000 11.0260000 1.1896752 28.2001737 11.0355335
0.6382752 29.7888838 11.6178521 1.3530000 29.7710000 9.6240000
2.2853781 29.4823606 9.1387605 0.5165736 29.2871750 9.1196156
1.1950000 31.2140000 9.5100000 2.0402889 31.7559472 9.4010079
0.0250000 31.8560000 9.5300000 -1.1240000 31.1950000 9.6440000
-1.1184480 30.1876771 9.7172792 -1.9965222 31.7035690 9.6567368
-0.0110000 33.1900000 9.4450000 0.8503955 33.7113265 9.3655232
-0.8994213 33.6701677 9.4606653 4.5230000 28.3770000 13.0940000
5.1770000 27.4790000 12.5260000 5.0380000 29.2340000 13.9550000
4.4689045 29.9577865 14.3701659 6.4590000 29.1810000 14.3060000
7.0273061 29.1604095 13.3761057 6.9530000 30.3450000 15.1520000
6.6315203 31.2648667 14.6635378 6.4750346 30.2665599 16.1284713
8.4760000 30.3660000 15.3290000 8.7022319 30.6903673 16.3447292
8.8562798 29.3557734 15.1775682 9.1340000 31.3060000 14.3410000
9.1644809 32.2898693 14.8091572 10.1511710 30.9456415 14.1873632
8.4300000 31.4040000 13.0120000 7.3749346 31.6073838 13.1952283
8.8666652 32.2302269 12.4509547 8.5630000 30.1510000 12.2300000
8.1580512 29.3853497 12.7495150 8.0809714 30.2517421 11.3481840
9.5405535 29.9625850 12.0597412 6.7500000 27.8840000 15.0730000
7.7960000 27.2570000 14.8140000 6.0009031 27.5098717 15.9726230

View file

@ -0,0 +1,15 @@
initial minimization prior to MD
&cntrl
imin = 1,
maxcyc = 3000,
ncyc = 3500,
ntb = 0,
igb = 6,
cut = 500,
ntmin = 1,
drms = 0.000001,
scee = 1.2,
scnb = 2.0,
nsnb = 1
/

File diff suppressed because it is too large Load diff

View file

@ -0,0 +1,68 @@
#!/bin/csh -f
#MSUB -A mscfcons
#MSUB -l nodes=1:ppn=8,walltime=30:00
#MSUB -e test.out.%j
#MSUB -o test.out.%j
#MSUB -N NWChemAmber
### If you do this as a bash shell then NWChem gives bogus errors at completion related the the task system calls
source /etc/profile.d/modules.csh
source /home/scicons/bin/set_modulepath.csh
module purge
module load amber/10
setenv AMBEREXE ${AMBERHOME}/bin/sander.MPI
setenv DOAMBER "srun --mpi=none -n 1 -N 1 ${AMBEREXE}"
cd big
${DOAMBER} -O -i big_am.in -o big_am.out -p big_am.top -c big_am.cor -r big_am.rst
cd ..
cd big.shake
${DOAMBER} -O -i big_am.in -o big_am.out -p big_am.top -c big_am.cor -r big_am.rst
cd ..
cd ethane
${DOAMBER} -O -i eth_am.in -o eth_am.out -p eth_am.top -c eth_am.cor -r eth_am.rst
cd ..
cd ethane+
${DOAMBER} -O -i ethplus_am.in -o ethplus_am.out -p ethplus_am.top -c ethplus_am.cor -r ethplus_am.rst
cd ..
cd small
${DOAMBER} -O -i small_am.in -o small_am.out -p small_am.top -c small_am.cor -r small_am.rst
cd ..
module purge
module load nwchem
setenv NWCHEMEXE ${NWCHEM_TOP}/bin/${NWCHEM_TARGET}/nwchem
setenv DONWCHEM "mpirun -srun -n 1 -N 1 ${NWCHEMEXE}"
cd big
${DONWCHEM} big.nw >& big.out
cd ..
cd big.shake
${DONWCHEM} big.nw >& big.out
cd ..
cd ethane
${DONWCHEM} eth.nw >& eth.out
cd ..
cd ethane+
${DONWCHEM} ethplus.nw >& ethplus.out
cd ../
cd small
${DONWCHEM} small.nw >& small.out
cd ..
cd methanol
${DONWCHEM} meoh.nw >& meoh.out
cd ..

View file

@ -0,0 +1,47 @@
This is an Ethane 4 plus ion, that NWChem generates new data for
NWChem input file:
ethplus.nw - nwchem input deck
Amber input files:
ethplus_am.in - amber input deck
ethplus_am.cor - coordinates
ethplus_am.top - topology file
NWChem output files:
ethplus.zmat
ethplus.p
ethplus.nam
ethplus.movecs
ethplus.c
ethplus.b^-1
ethplus.b
ethplus.xyz
ethplus.top
ethplus.seq
ethplus.q
ethplus.plt
ethplus_md.rst
ethplus_md.cmd
ethplus.grid
ethplus_energy.out
ethplus.bsg
ethpl.sgm
ethpl.frg
ethplus.out
ethplus_md.qrs
ethplus_md.out
ethplus.db
Amber output files:
mdinfo
ethplus_am.rst
ethplus_am.out
Results from *.out files:
1.471673291E+03 kJ/mol NWChem single point (3.5174E+02 kcal/mol)
1.38277E+03 kJ/mol NWChem step 25 (3.3049E+02 kcal/mol)
3.5181E+02 kcal/mol Amber single point
3.3048E+02 kcal/mol Amber step 400
Note that NWChem and Amber match both optimized and initial

View file

@ -0,0 +1,57 @@
#
# This is an AMBER test that uses a NWChem geometry
# It is uses the automatic atom type determination at the end of the Amber force file
#
Title "ethane plus"
memory 3000 Mb
permanent_dir .
start ethplus
echo
charge +4
geometry autosym
H -0.97624726 -0.26158466 -1.15438725
H 0.71466259 -0.71466259 -1.15438725
H 0.26158466 0.97624726 -1.15438725
C 0.00000000 0.00000000 -0.76221612
C 0.00000000 0.00000000 0.76221612
H -0.71466259 0.71466259 1.15438725
H -0.26158466 -0.97624726 1.15438725
H 0.97624726 0.26158466 1.15438725
end
basis
* library sto-3g
end
task scf
esp
constrain equal 1 2 3 6 7 8
constrain equal 4 5
end
prepare
source rtdb
amber
new_top new_seq
new_rst
# write pdb ethplus.pdb
vacuo
end
task prepare
md
system ethplus_md
noshake solvent solute
end
task md energy
task shell "cp ethplus_md.out ethplus_energy.out"
task md optimize

View file

@ -0,0 +1,6 @@
8
-0.7150000 0.7150000 1.1540000 0.9760000 0.2620000 1.1540000
-0.2620000 -0.9760000 1.1540000 0.0000000 0.0000000 0.7620000
0.0000000 0.0000000 -0.7620000 -0.9760000 -0.2620000 -1.1540000
0.2620000 0.9760000 -1.1540000 0.7150000 -0.7150000 -1.1540000

View file

@ -0,0 +1,12 @@
initial minimization prior to MD
&cntrl
imin = 1,
maxcyc = 400,
ncyc = 400,
ntb = 0,
igb = 0,
cut = 50
ntmin = 1
drms = 0.000001
/

View file

@ -0,0 +1,130 @@
%VERSION VERSION_STAMP = V0001.000 DATE = 07/31/09 09:02:17
%FLAG TITLE
%FORMAT(20a4)
%FLAG POINTERS
%FORMAT(10I8)
8 2 6 1 12 0 9 0 0 0
29 1 1 0 0 2 2 1 2 0
0 0 0 0 0 0 0 0 8 0
0
%FLAG ATOM_NAME
%FORMAT(20a4)
H12 H13 H14 C1 C5 H22 H23 H24
%FLAG CHARGE
%FORMAT(5E16.8)
1.12008834E+01 1.12008834E+01 1.12008834E+01 2.84196813E+00 2.84196813E+00
1.12008834E+01 1.12008834E+01 1.12008834E+01
%FLAG MASS
%FORMAT(5E16.8)
1.00800000E+00 1.00800000E+00 1.00800000E+00 1.20100000E+01 1.20100000E+01
1.00800000E+00 1.00800000E+00 1.00800000E+00
%FLAG ATOM_TYPE_INDEX
%FORMAT(10I8)
1 1 1 2 2 1 1 1
%FLAG NUMBER_EXCLUDED_ATOMS
%FORMAT(10I8)
7 6 5 4 3 2 1 1
%FLAG NONBONDED_PARM_INDEX
%FORMAT(10I8)
1 2 2 3
%FLAG RESIDUE_LABEL
%FORMAT(20a4)
ETH
%FLAG RESIDUE_POINTER
%FORMAT(10I8)
1
%FLAG BOND_FORCE_CONSTANT
%FORMAT(5E16.8)
3.40000000E+02 3.10000000E+02
%FLAG BOND_EQUIL_VALUE
%FORMAT(5E16.8)
1.09000000E+00 1.52600000E+00
%FLAG ANGLE_FORCE_CONSTANT
%FORMAT(5E16.8)
5.00000000E+01 3.50000000E+01
%FLAG ANGLE_EQUIL_VALUE
%FORMAT(5E16.8)
1.91113635E+00 1.91113635E+00
%FLAG DIHEDRAL_FORCE_CONSTANT
%FORMAT(5E16.8)
1.50000000E-01
%FLAG DIHEDRAL_PERIODICITY
%FORMAT(5E16.8)
3.00000000E+00
%FLAG DIHEDRAL_PHASE
%FORMAT(5E16.8)
0.00000000E+00
%FLAG SOLTY
%FORMAT(5E16.8)
0.00000000E+00 0.00000000E+00
%FLAG LENNARD_JONES_ACOEF
%FORMAT(5E16.8)
7.51607703E+03 9.71708117E+04 1.04308023E+06
%FLAG LENNARD_JONES_BCOEF
%FORMAT(5E16.8)
2.17257828E+01 1.26919150E+02 6.75612247E+02
%FLAG BONDS_INC_HYDROGEN
%FORMAT(10I8)
0 9 1 3 9 1 6 9 1 12
21 1 12 18 1 12 15 1
%FLAG BONDS_WITHOUT_HYDROGEN
%FORMAT(10I8)
9 12 2
%FLAG ANGLES_INC_HYDROGEN
%FORMAT(10I8)
0 9 12 1 0 9 6 2 0 9
3 2 3 9 12 1 3 9 6 2
6 9 12 1 9 12 21 1 9 12
18 1 9 12 15 1 15 12 21 2
15 12 18 2 18 12 21 2
%FLAG ANGLES_WITHOUT_HYDROGEN
%FORMAT(10I8)
%FLAG DIHEDRALS_INC_HYDROGEN
%FORMAT(10I8)
0 9 12 21 1 0 9 12 18 1
0 9 12 15 1 3 9 12 21 1
3 9 12 18 1 3 9 12 15 1
6 9 12 21 1 6 9 12 18 1
6 9 12 15 1
%FLAG DIHEDRALS_WITHOUT_HYDROGEN
%FORMAT(10I8)
%FLAG EXCLUDED_ATOMS_LIST
%FORMAT(10I8)
2 3 4 5 6 7 8 3 4 5
6 7 8 4 5 6 7 8 5 6
7 8 6 7 8 7 8 8 0
%FLAG HBOND_ACOEF
%FORMAT(5E16.8)
%FLAG HBOND_BCOEF
%FORMAT(5E16.8)
%FLAG HBCUT
%FORMAT(5E16.8)
%FLAG AMBER_ATOM_TYPE
%FORMAT(20a4)
HC HC HC CT CT HC HC HC
%FLAG TREE_CHAIN_CLASSIFICATION
%FORMAT(20a4)
BLA BLA BLA BLA BLA BLA BLA BLA
%FLAG JOIN_ARRAY
%FORMAT(10I8)
0 0 0 0 0 0 0 0
%FLAG IROTAT
%FORMAT(10I8)
0 0 0 0 0 0 0 0
%FLAG RADIUS_SET
%FORMAT(1a80)
modified Bondi radii (mbondi)
%FLAG RADII
%FORMAT(5E16.8)
1.30000000E+00 1.30000000E+00 1.30000000E+00 1.70000000E+00 1.70000000E+00
1.30000000E+00 1.30000000E+00 1.30000000E+00
%FLAG SCREEN
%FORMAT(5E16.8)
8.50000000E-01 8.50000000E-01 8.50000000E-01 7.20000000E-01 7.20000000E-01
8.50000000E-01 8.50000000E-01 8.50000000E-01

View file

@ -0,0 +1,47 @@
This is an Ethane molecule, that NWChem generates new data for
NWChem input file:
eth.nw - nwchem input deck
Amber input files:
eth_am.in - amber input deck
eth_am.cor - coordinates
eth_am.top - topology file
NWChem output files:
-rw-rw-r-- 1 d3p852 d3p852 248 Oct 19 14:12 eth.zmat
-rw-rw-r-- 1 d3p852 d3p852 6.2K Oct 19 14:12 eth.p
-rw-rw-r-- 1 d3p852 d3p852 6.2K Oct 19 14:12 eth.c
-rw-rw-r-- 1 d3p852 d3p852 5.3K Oct 19 14:12 eth.b^-1
-rw-rw-r-- 1 d3p852 d3p852 5.3K Oct 19 14:12 eth.b
-rw-rw-r-- 1 d3p852 d3p852 803 Oct 19 14:12 eth.q
-rw-rw-r-- 1 d3p852 d3p852 104 Oct 19 14:12 eth.nam
-rw-rw-r-- 1 d3p852 d3p852 2.7K Oct 19 14:12 eth.movecs
-rw-rw-r-- 1 d3p852 d3p852 167K Oct 19 14:12 eth.grid
-rw-rw-r-- 1 d3p852 d3p852 442 Oct 19 14:12 eth.xyz
-rw-rw-r-- 1 d3p852 d3p852 4.7K Oct 19 14:12 eth.top
-rw-rw-r-- 1 d3p852 d3p852 2.3K Oct 19 14:12 eth.sgm
-rw-rw-r-- 1 d3p852 d3p852 40 Oct 19 14:12 eth.seq
-rw-rw-r-- 1 d3p852 d3p852 169K Oct 19 14:12 eth.plt
-rw-rw-r-- 1 d3p852 d3p852 41K Oct 19 14:12 eth.out
-rw-rw-r-- 1 d3p852 d3p852 1.3K Oct 19 14:12 eth_md.rst
-rw-rw-r-- 1 d3p852 d3p852 20K Oct 19 14:12 eth_md.qrs
-rw-rw-r-- 1 d3p852 d3p852 9.6K Oct 19 14:12 eth_md.out
-rw-rw-r-- 1 d3p852 d3p852 0 Oct 19 14:12 eth_md.cmd
-rw-rw-r-- 1 d3p852 d3p852 824 Oct 19 14:12 eth.frg
-rw-rw-r-- 1 d3p852 d3p852 12K Oct 19 14:12 eth_energy.out
-rw-rw-r-- 1 d3p852 d3p852 1.1M Oct 19 14:12 eth.db
-rw-rw-r-- 1 d3p852 d3p852 426 Oct 19 14:12 eth.bsg
Amber output files:
-rw-rw-r-- 1 d3p852 d3p852 384 Oct 19 14:11 mdinfo
-rw-rw-r-- 1 d3p852 d3p852 303 Oct 19 14:11 eth_am.rst
-rw-rw-r-- 1 d3p852 d3p852 13K Oct 19 14:11 eth_am.out
Results from *.out files:
3.714703721E+00 kJ/mol Single point energy NWChem (8.8784E-01 kcal/mol)
1.63128E+00 kJ/mol after 19 steps in NWChem (3.8989E-01 kcal/mol)
8.8371E-01 kcal/mol single point energy Amber
3.8987E-01 kcal/mol after 500 steps in Amber
Note that NWChem and Amber match both optimized and initial

View file

@ -0,0 +1,54 @@
#
# This is an AMBER test that uses a NWChem geometry
# It is uses the automatic atom type determination at the end of the Amber force file
#
Title "ethane"
memory 3000 Mb
permanent_dir .
start eth
echo
geometry autosym
H -0.97624726 -0.26158466 -1.15438725
H 0.71466259 -0.71466259 -1.15438725
H 0.26158466 0.97624726 -1.15438725
C 0.00000000 0.00000000 -0.76221612
C 0.00000000 0.00000000 0.76221612
H -0.71466259 0.71466259 1.15438725
H -0.26158466 -0.97624726 1.15438725
H 0.97624726 0.26158466 1.15438725
end
basis
* library sto-3g
end
task scf
esp
constrain equal 1 2 3 6 7 8
constrain equal 4 5
end
prepare
source rtdb
amber
new_top new_seq
new_rst
# write pdb eth.pdb
end
task prepare
md
system eth_md
noshake solvent solute
end
task md energy
task shell "cp eth_md.out eth_energy.out"
task md optimize

View file

@ -0,0 +1,14 @@
HEADER
TITLE THIS IS ETHANE
REMARK 4 XXXX COMPLIES WITH FORMAT V. 2.1, 25-OCT-1996
CRYST1 1.952 1.952 2.309 90.00 90.00 90.00
ATOM 1 2H1 ETH 1 -0.715 0.715 1.154 0.00 H
ATOM 2 3H1 ETH 1 0.976 0.262 1.154 0.00 H
ATOM 3 4H1 ETH 1 -0.262 -0.976 1.154 0.00 H
ATOM 4 C1 ETH 1 0.000 0.000 0.762 0.00 C
ATOM 5 C2 ETH 1 0.000 0.000 -0.762 0.00 C
ATOM 6 2H2 ETH 1 -0.976 -0.262 -1.154 0.00 H
ATOM 7 3H2 ETH 1 0.262 0.976 -1.154 0.00 H
ATOM 8 4H2 ETH 1 0.715 -0.715 -1.154 0.00 H
TER
END

View file

@ -0,0 +1,6 @@
8
-0.7150000 0.7150000 1.1540000 0.9760000 0.2620000 1.1540000
-0.2620000 -0.9760000 1.1540000 0.0000000 0.0000000 0.7620000
0.0000000 0.0000000 -0.7620000 -0.9760000 -0.2620000 -1.1540000
0.2620000 0.9760000 -1.1540000 0.7150000 -0.7150000 -1.1540000

View file

@ -0,0 +1,12 @@
initial minimization prior to MD
&cntrl
imin = 1,
maxcyc = 700,
ncyc = 500,
ntb = 0,
igb = 0,
cut = 50
ntmin = 1
drms = 0.000001
/

View file

@ -0,0 +1,130 @@
%VERSION VERSION_STAMP = V0001.000 DATE = 07/14/09 15:01:10
%FLAG TITLE
%FORMAT(20a4)
%FLAG POINTERS
%FORMAT(10I8)
8 2 6 1 12 0 9 0 0 0
29 1 1 0 0 2 2 1 2 0
0 0 0 0 0 0 0 0 8 0
0
%FLAG ATOM_NAME
%FORMAT(20a4)
H12 H13 H14 C1 C5 H22 H23 H24
%FLAG CHARGE
%FORMAT(5E16.8)
3.11674219E-01 3.11674219E-01 3.11674219E-01 -9.35040880E-01 -9.35040880E-01
3.11674219E-01 3.11674219E-01 3.11674219E-01
%FLAG MASS
%FORMAT(5E16.8)
1.00800000E+00 1.00800000E+00 1.00800000E+00 1.20100000E+01 1.20100000E+01
1.00800000E+00 1.00800000E+00 1.00800000E+00
%FLAG ATOM_TYPE_INDEX
%FORMAT(10I8)
1 1 1 2 2 1 1 1
%FLAG NUMBER_EXCLUDED_ATOMS
%FORMAT(10I8)
7 6 5 4 3 2 1 1
%FLAG NONBONDED_PARM_INDEX
%FORMAT(10I8)
1 2 2 3
%FLAG RESIDUE_LABEL
%FORMAT(20a4)
ETH
%FLAG RESIDUE_POINTER
%FORMAT(10I8)
1
%FLAG BOND_FORCE_CONSTANT
%FORMAT(5E16.8)
3.40000000E+02 3.10000000E+02
%FLAG BOND_EQUIL_VALUE
%FORMAT(5E16.8)
1.09000000E+00 1.52600000E+00
%FLAG ANGLE_FORCE_CONSTANT
%FORMAT(5E16.8)
5.00000000E+01 3.50000000E+01
%FLAG ANGLE_EQUIL_VALUE
%FORMAT(5E16.8)
1.91113635E+00 1.91113635E+00
%FLAG DIHEDRAL_FORCE_CONSTANT
%FORMAT(5E16.8)
1.50000000E-01
%FLAG DIHEDRAL_PERIODICITY
%FORMAT(5E16.8)
3.00000000E+00
%FLAG DIHEDRAL_PHASE
%FORMAT(5E16.8)
0.00000000E+00
%FLAG SOLTY
%FORMAT(5E16.8)
0.00000000E+00 0.00000000E+00
%FLAG LENNARD_JONES_ACOEF
%FORMAT(5E16.8)
7.51607703E+03 9.71708117E+04 1.04308023E+06
%FLAG LENNARD_JONES_BCOEF
%FORMAT(5E16.8)
2.17257828E+01 1.26919150E+02 6.75612247E+02
%FLAG BONDS_INC_HYDROGEN
%FORMAT(10I8)
0 9 1 3 9 1 6 9 1 12
21 1 12 15 1 12 18 1
%FLAG BONDS_WITHOUT_HYDROGEN
%FORMAT(10I8)
9 12 2
%FLAG ANGLES_INC_HYDROGEN
%FORMAT(10I8)
0 9 12 1 0 9 6 2 0 9
3 2 3 9 12 1 3 9 6 2
6 9 12 1 9 12 21 1 9 12
15 1 9 12 18 1 15 12 21 2
15 12 18 2 18 12 21 2
%FLAG ANGLES_WITHOUT_HYDROGEN
%FORMAT(10I8)
%FLAG DIHEDRALS_INC_HYDROGEN
%FORMAT(10I8)
0 9 12 21 1 0 9 12 15 1
0 9 12 18 1 3 9 12 21 1
3 9 12 15 1 3 9 12 18 1
6 9 12 21 1 6 9 12 15 1
6 9 12 18 1
%FLAG DIHEDRALS_WITHOUT_HYDROGEN
%FORMAT(10I8)
%FLAG EXCLUDED_ATOMS_LIST
%FORMAT(10I8)
2 3 4 5 6 7 8 3 4 5
6 7 8 4 5 6 7 8 5 6
7 8 6 7 8 7 8 8 0
%FLAG HBOND_ACOEF
%FORMAT(5E16.8)
%FLAG HBOND_BCOEF
%FORMAT(5E16.8)
%FLAG HBCUT
%FORMAT(5E16.8)
%FLAG AMBER_ATOM_TYPE
%FORMAT(20a4)
HC HC HC CT CT HC HC HC
%FLAG TREE_CHAIN_CLASSIFICATION
%FORMAT(20a4)
BLA BLA BLA BLA BLA BLA BLA BLA
%FLAG JOIN_ARRAY
%FORMAT(10I8)
0 0 0 0 0 0 0 0
%FLAG IROTAT
%FORMAT(10I8)
0 0 0 0 0 0 0 0
%FLAG RADIUS_SET
%FORMAT(1a80)
modified Bondi radii (mbondi)
%FLAG RADII
%FORMAT(5E16.8)
1.30000000E+00 1.30000000E+00 1.30000000E+00 1.70000000E+00 1.70000000E+00
1.30000000E+00 1.30000000E+00 1.30000000E+00
%FLAG SCREEN
%FORMAT(5E16.8)
8.50000000E-01 8.50000000E-01 8.50000000E-01 7.20000000E-01 7.20000000E-01
8.50000000E-01 8.50000000E-01 8.50000000E-01

View file

@ -0,0 +1,21 @@
This is a solvent only NWChem run
NWChem input files:
meoh.nw - nwchem input deck
meoh.pdb - EMPTY pdb file, to fool nwchem
meoh.rst - restart file - same as scr/data/solvents/meoh.rst
NWChem output files:
meoh.top
meoh.seq
meoh_out.pdb - This is a PDB file of the resulting structure
meoh.out
meoh.nam
meoh_em.rst
meoh.db
meoh.bsg
Result:
Expect a 5780 line PDB file that starts like
CRYST1 42.476 42.476 42.476 90.00 90.00 90.00
Note that this is not a test of a result, it is a functionality test

View file

@ -0,0 +1,33 @@
#
# This is a test of a Solvent Only run
# The meoh.pdb file is empty, because NWChem needs a solute
# This uses the segment file provided with NWChem
#
Title "solvent methanol"
memory 1800 Mb
permanent_dir .
start meoh
echo
basis
* library sto-3g
end
prepare
amber
source pdb
read rst meoh.rst
system meoh_em
solvent name MOH model meoh
touch 0.25
expand 0.1
center
solvate
new_top new_seq
new_rst
write pdb meoh_out.pdb
end
task prepare

View file

@ -0,0 +1,2 @@
REMARK EMPTY PDB FILE
END

File diff suppressed because it is too large Load diff

View file

@ -0,0 +1,37 @@
This is an ALA dimer
NWChem input files:
small.nw - nwchem job file
small.pdb - geometry as a PDB file, used by small.nw
NWChem output file:
small_energy.out
small_md.out
small.top
small.seq
small.out
small.nam
small_md.rst
small_md.qrs
small_md.cmd
small.db
small.bsg
Amber input files:
small_am.in - amber job file
small_am.cor - amber coordinate file
small_am.top - amber topology file
Amber output files:
small_am.out
small_am.rst
mdinfo
Results from *.out files:
-1.374134919E+02 kJ/mol for single point NWChem (-3.2843E+01 kcal/mol)
-1.83007E+02 kJ/mol for 500 steps NWChem (-4.3740E+01 kcal/mol)
-1.90732E+02 kj/mol for 7000 steps NWChem (-4.5586E+01 kcal/mol)
-3.2847E+01 kcal/mol for single point Amber
-4.5025E+01 kcal/mol for 500 steps Amber
-4.5867E+01 kcal/mol for 7000 steps Amber
Note that the initial and optimized energies match well

View file

@ -0,0 +1,33 @@
#
# This is a amber test that uses a PDB geometry
#
Title "small"
memory 3500 Mb
permanent_dir .
# use the next line to verify that all atoms are really found
# print debug
echo
prepare
system small
source small.pdb
new_top new_seq
new_rst
amber
vacuo
end
print debug
task prepare
md
system small_md
sd 7000 init 0.001 min 0.000000001 max 0.05
noshake solvent solute
end
task md energy
task shell "cp small_md.out small_energy.out"
task md optimize

View file

@ -0,0 +1,25 @@
ATOM 1 N ALA 1 -0.491 0.404 1.207 1.00 0.00
ATOM 2 2H ALA 1 -1.404 0.049 1.455 1.00 0.00
ATOM 3 3H ALA 1 0.141 0.360 1.994 1.00 0.00
ATOM 4 4H ALA 1 -0.544 1.373 0.924 1.00 0.00
ATOM 5 CA ALA 1 0.112 -0.360 0.098 1.00 0.00
ATOM 6 HA ALA 1 -0.506 -0.224 -0.791 1.00 0.00
ATOM 7 CB ALA 1 0.185 -1.856 0.423 1.00 0.00
ATOM 8 2HB ALA 1 0.769 -2.016 1.330 1.00 0.00
ATOM 9 3HB ALA 1 0.650 -2.395 -0.404 1.00 0.00
ATOM 10 4HB ALA 1 -0.819 -2.256 0.576 1.00 0.00
ATOM 11 C ALA 1 1.501 0.207 -0.184 1.00 0.00
ATOM 12 O ALA 1 1.862 1.202 0.441 1.00 0.00
ATOM 13 N ALA 2 2.246 -0.417 -1.094 1.00 0.00
ATOM 14 H ALA 2 1.869 -1.205 -1.607 1.00 0.00
ATOM 15 CA ALA 2 3.541 0.027 -1.590 1.00 0.00
ATOM 16 HA ALA 2 4.313 -0.337 -0.910 1.00 0.00
ATOM 17 CB ALA 2 3.608 1.556 -1.663 1.00 0.00
ATOM 18 2HB ALA 2 2.830 1.929 -2.330 1.00 0.00
ATOM 19 3HB ALA 2 4.585 1.870 -2.032 1.00 0.00
ATOM 20 4HB ALA 2 3.454 1.974 -0.668 1.00 0.00
ATOM 21 C ALA 2 3.746 -0.574 -2.985 1.00 0.00
ATOM 22 O ALA 2 4.786 -0.260 -3.601 1.00 0.00
ATOM 23 OXT ALA 2 2.844 -1.340 -3.399 1.00 0.00
TER
END

View file

@ -0,0 +1,14 @@
23
-0.4910000 0.4040000 1.2070000 -1.4040000 0.0490000 1.4550000
0.1410000 0.3600000 1.9940000 -0.5440000 1.3730000 0.9240000
0.1120000 -0.3600000 0.0980000 -0.5060000 -0.2240000 -0.7910000
0.1850000 -1.8560000 0.4230000 0.7690000 -2.0160000 1.3300000
0.6500000 -2.3950000 -0.4040000 -0.8190000 -2.2560000 0.5760000
1.5010000 0.2070000 -0.1840000 1.8620000 1.2020000 0.4410000
2.2460000 -0.4170000 -1.0940000 1.8690000 -1.2050000 -1.6070000
3.5410000 0.0270000 -1.5900000 4.3130000 -0.3370000 -0.9100000
3.6080000 1.5560000 -1.6630000 2.8300000 1.9290000 -2.3300000
4.5850000 1.8700000 -2.0320000 3.4540000 1.9740000 -0.6680000
3.7460000 -0.5740000 -2.9850000 4.7860000 -0.2600000 -3.6010000
2.8440000 -1.3400000 -3.3990000

View file

@ -0,0 +1,12 @@
initial minimization prior to MD
&cntrl
imin = 1,
maxcyc = 7000,
ncyc = 3500,
ntb = 0,
igb = 0,
cut = 50
ntmin = 1
drms = 0.000001
/

View file

@ -0,0 +1,240 @@
%VERSION VERSION_STAMP = V0001.000 DATE = 07/13/09 12:33:19
%FLAG TITLE
%FORMAT(20a4)
%FLAG POINTERS
%FORMAT(10I8)
23 8 12 10 26 13 39 19 0 0
111 2 10 13 19 12 23 13 10 0
0 0 0 0 0 0 0 0 12 0
0
%FLAG ATOM_NAME
%FORMAT(20a4)
N H1 H2 H3 CA HA CB HB1 HB2 HB3 C O N H CA HA CB HB1 HB2 HB3
C O OXT
%FLAG CHARGE
%FORMAT(5E16.8)
2.57663322E+00 3.63899331E+00 3.63899331E+00 3.63899331E+00 1.75298526E+00
1.61996247E+00 -1.08787131E+00 5.46669000E-01 5.46669000E-01 5.46669000E-01
1.12304035E+01 -1.04268001E+01 -6.96274083E+00 4.88539863E+00 -3.18343581E+00
1.94431941E+00 -3.81392739E+00 1.39218372E+00 1.39218372E+00 1.39218372E+00
1.40876601E+01 -1.46780627E+01 -1.46780627E+01
%FLAG MASS
%FORMAT(5E16.8)
1.40100000E+01 1.00800000E+00 1.00800000E+00 1.00800000E+00 1.20100000E+01
1.00800000E+00 1.20100000E+01 1.00800000E+00 1.00800000E+00 1.00800000E+00
1.20100000E+01 1.60000000E+01 1.40100000E+01 1.00800000E+00 1.20100000E+01
1.00800000E+00 1.20100000E+01 1.00800000E+00 1.00800000E+00 1.00800000E+00
1.20100000E+01 1.60000000E+01 1.60000000E+01
%FLAG ATOM_TYPE_INDEX
%FORMAT(10I8)
1 2 2 2 3 4 3 5 5 5
6 7 1 2 3 8 3 5 5 5
6 7 7
%FLAG NUMBER_EXCLUDED_ATOMS
%FORMAT(10I8)
12 6 5 4 10 7 6 3 2 1
7 3 10 4 8 7 6 3 2 1
2 1 1
%FLAG NONBONDED_PARM_INDEX
%FORMAT(10I8)
1 2 4 7 11 16 22 29 2 3
5 8 12 17 23 30 4 5 6 9
13 18 24 31 7 8 9 10 14 19
25 32 11 12 13 14 15 20 26 33
16 17 18 19 20 21 27 34 22 23
24 25 26 27 28 35 29 30 31 32
33 34 35 36
%FLAG RESIDUE_LABEL
%FORMAT(20a4)
ALA ALA
%FLAG RESIDUE_POINTER
%FORMAT(10I8)
1 13
%FLAG BOND_FORCE_CONSTANT
%FORMAT(5E16.8)
5.70000000E+02 4.90000000E+02 3.40000000E+02 3.40000000E+02 3.10000000E+02
3.17000000E+02 4.34000000E+02 3.67000000E+02 6.56000000E+02 3.40000000E+02
4.34000000E+02 3.37000000E+02
%FLAG BOND_EQUIL_VALUE
%FORMAT(5E16.8)
1.22900000E+00 1.33500000E+00 1.09000000E+00 1.09000000E+00 1.52600000E+00
1.52200000E+00 1.01000000E+00 1.47100000E+00 1.25000000E+00 1.09000000E+00
1.01000000E+00 1.44900000E+00
%FLAG ANGLE_FORCE_CONSTANT
%FORMAT(5E16.8)
8.00000000E+01 5.00000000E+01 5.00000000E+01 3.50000000E+01 6.30000000E+01
5.00000000E+01 5.00000000E+01 5.00000000E+01 8.00000000E+01 7.00000000E+01
5.00000000E+01 3.50000000E+01 5.00000000E+01 8.00000000E+01 8.00000000E+01
8.00000000E+01 5.00000000E+01 5.00000000E+01 7.00000000E+01 5.00000000E+01
5.00000000E+01 8.00000000E+01 6.30000000E+01
%FLAG ANGLE_EQUIL_VALUE
%FORMAT(5E16.8)
2.14501057E+00 2.09439600E+00 2.12755727E+00 1.91113635E+00 1.93906163E+00
1.91113635E+00 1.91113635E+00 1.91113635E+00 2.10137732E+00 2.03505478E+00
1.91113635E+00 1.91113635E+00 1.91113635E+00 1.94080696E+00 1.94080696E+00
2.19911580E+00 1.91113635E+00 1.91113635E+00 2.04203610E+00 2.06018753E+00
1.91113635E+00 1.91462701E+00 1.92160833E+00
%FLAG DIHEDRAL_FORCE_CONSTANT
%FORMAT(5E16.8)
2.00000000E+00 2.50000000E+00 0.00000000E+00 5.30000000E-01 1.50000000E-01
5.00000000E-01 8.00000000E-01 8.50000000E-01 1.55555556E-01 7.00000000E-02
1.00000000E-01 1.05000000E+01 1.10000000E+00
%FLAG DIHEDRAL_PERIODICITY
%FORMAT(5E16.8)
1.00000000E+00 2.00000000E+00 2.00000000E+00 1.00000000E+00 3.00000000E+00
4.00000000E+00 1.00000000E+00 2.00000000E+00 3.00000000E+00 2.00000000E+00
4.00000000E+00 2.00000000E+00 2.00000000E+00
%FLAG DIHEDRAL_PHASE
%FORMAT(5E16.8)
0.00000000E+00 3.14159400E+00 0.00000000E+00 0.00000000E+00 3.14159400E+00
3.14159400E+00 0.00000000E+00 3.14159400E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 3.14159400E+00 3.14159400E+00
%FLAG SOLTY
%FORMAT(5E16.8)
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
%FLAG LENNARD_JONES_ACOEF
%FORMAT(5E16.8)
9.44293233E+05 2.12601181E+03 1.39982777E-01 9.95480466E+05 2.56678134E+03
1.04308023E+06 2.01791425E+04 9.14716912E+00 2.27401052E+04 2.01823541E+02
8.96776989E+04 1.07193646E+02 9.71708117E+04 1.41077189E+03 7.51607703E+03
8.82619071E+05 2.27577561E+03 9.24822270E+05 2.01619733E+04 8.61541883E+04
8.19971662E+05 6.06829342E+05 1.02595236E+03 6.47841731E+05 1.12780457E+04
5.44261042E+04 5.74393458E+05 3.79876399E+05 6.20665997E+04 5.94667300E+01
6.78771368E+04 8.79040886E+02 4.98586848E+03 6.01816484E+04 3.69471530E+04
3.25969625E+03
%FLAG LENNARD_JONES_BCOEF
%FORMAT(5E16.8)
8.01323529E+02 2.09604198E+01 9.37598976E-02 7.36907417E+02 2.06278363E+01
6.75612247E+02 6.45756063E+01 7.57919667E-01 6.13981767E+01 3.56012899E+00
1.36131731E+02 2.59456373E+00 1.26919150E+02 9.41257003E+00 2.17257828E+01
6.53361429E+02 1.82891803E+01 5.99015525E+02 5.44372326E+01 1.12529845E+02
5.31102864E+02 6.77220874E+02 1.53505284E+01 6.26720080E+02 5.08951803E+01
1.11805549E+02 5.55666448E+02 5.64885984E+02 1.13252061E+02 1.93248820E+00
1.06076943E+02 7.42992380E+00 1.76949863E+01 9.40505980E+01 9.21192136E+01
1.43076527E+01
%FLAG BONDS_INC_HYDROGEN
%FORMAT(10I8)
18 21 3 18 24 3 18 27 3 12
15 4 0 3 7 0 6 7 0 9
7 48 51 3 48 54 3 48 57 3
42 45 10 36 39 11
%FLAG BONDS_WITHOUT_HYDROGEN
%FORMAT(10I8)
30 33 1 30 36 2 12 18 5 12
30 6 0 12 8 60 63 9 60 66
9 42 48 5 42 60 6 36 42 12
%FLAG ANGLES_INC_HYDROGEN
%FORMAT(10I8)
30 36 39 2 24 18 27 4 21 18
24 4 21 18 27 4 15 12 18 6
15 12 30 7 12 18 21 8 12 18
24 8 12 18 27 8 9 0 12 11
6 0 9 12 6 0 12 11 3 0
6 12 3 0 9 12 3 0 12 11
0 12 15 13 54 48 57 4 51 48
54 4 51 48 57 4 45 42 48 17
45 42 60 18 42 48 51 8 42 48
54 8 42 48 57 8 39 36 42 20
36 42 45 21
%FLAG ANGLES_WITHOUT_HYDROGEN
%FORMAT(10I8)
33 30 36 1 30 36 42 3 18 12
30 5 12 30 33 9 12 30 36 10
0 12 18 14 0 12 30 15 63 60
66 16 48 42 60 5 42 60 63 19
42 60 66 19 36 42 48 22 36 42
60 23
%FLAG DIHEDRALS_INC_HYDROGEN
%FORMAT(10I8)
33 30 36 39 1 33 30 -36 39 2
30 36 42 45 3 27 18 12 30 9
24 18 12 30 9 21 18 12 30 9
15 12 18 21 9 15 12 18 24 9
15 12 18 27 9 15 12 30 33 3
15 12 30 36 3 12 30 36 39 2
9 0 12 15 9 9 0 12 18 9
9 0 12 30 9 6 0 12 15 9
6 0 12 18 9 6 0 12 30 9
3 0 12 15 9 3 0 12 18 9
3 0 12 30 9 0 12 18 21 9
0 12 18 24 9 0 12 18 27 9
57 48 42 60 9 54 48 42 60 9
51 48 42 60 9 45 42 48 51 9
45 42 48 54 9 45 42 48 57 9
45 42 60 63 3 45 42 60 66 3
39 36 42 45 3 39 36 42 48 3
39 36 42 60 3 36 42 48 51 9
36 42 48 54 9 36 42 48 57 9
30 42 -36 -39 13
%FLAG DIHEDRALS_WITHOUT_HYDROGEN
%FORMAT(10I8)
33 30 36 42 2 30 36 42 48 4
30 36 -42 48 3 30 36 -42 48 5
30 36 -42 48 6 30 36 42 60 7
30 36 -42 60 8 18 12 30 33 3
18 12 30 36 10 18 12 -30 36 11
12 30 36 42 2 0 12 30 33 3
0 12 30 36 3 48 42 60 63 3
48 42 60 66 3 36 42 60 63 3
36 42 60 66 3 12 36 -30 -33 12
42 63 -60 -66 12
%FLAG EXCLUDED_ATOMS_LIST
%FORMAT(10I8)
2 3 4 5 6 7 8 9 10 11
12 13 3 4 5 6 7 11 4 5
6 7 11 5 6 7 11 6 7 8
9 10 11 12 13 14 15 7 8 9
10 11 12 13 8 9 10 11 12 13
9 10 11 10 11 11 12 13 14 15
16 17 21 13 14 15 14 15 16 17
18 19 20 21 22 23 15 16 17 21
16 17 18 19 20 21 22 23 17 18
19 20 21 22 23 18 19 20 21 22
23 19 20 21 20 21 21 22 23 23
0
%FLAG HBOND_ACOEF
%FORMAT(5E16.8)
%FLAG HBOND_BCOEF
%FORMAT(5E16.8)
%FLAG HBCUT
%FORMAT(5E16.8)
%FLAG AMBER_ATOM_TYPE
%FORMAT(20a4)
N3 H H H CT HP CT HC HC HC C O N H CT H1 CT HC HC HC
C O2 O2
%FLAG TREE_CHAIN_CLASSIFICATION
%FORMAT(20a4)
BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA BLA
BLA BLA BLA
%FLAG JOIN_ARRAY
%FORMAT(10I8)
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0
%FLAG IROTAT
%FORMAT(10I8)
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0
%FLAG RADIUS_SET
%FORMAT(1a80)
modified Bondi radii (mbondi)
%FLAG RADII
%FORMAT(5E16.8)
1.55000000E+00 1.30000000E+00 1.30000000E+00 1.30000000E+00 1.70000000E+00
1.30000000E+00 1.70000000E+00 1.30000000E+00 1.30000000E+00 1.30000000E+00
1.70000000E+00 1.50000000E+00 1.55000000E+00 1.30000000E+00 1.70000000E+00
1.30000000E+00 1.70000000E+00 1.30000000E+00 1.30000000E+00 1.30000000E+00
1.70000000E+00 1.50000000E+00 1.50000000E+00
%FLAG SCREEN
%FORMAT(5E16.8)
7.90000000E-01 8.50000000E-01 8.50000000E-01 8.50000000E-01 7.20000000E-01
8.50000000E-01 7.20000000E-01 8.50000000E-01 8.50000000E-01 8.50000000E-01
7.20000000E-01 8.50000000E-01 7.90000000E-01 8.50000000E-01 7.20000000E-01
8.50000000E-01 7.20000000E-01 8.50000000E-01 8.50000000E-01 8.50000000E-01
7.20000000E-01 8.50000000E-01 8.50000000E-01