#  This is the control script for LAMMPS.

echo			both
log	    		2.1_Initialize.out

#-------------------------------------------------------------------------------
# Stage 2.1: Initialize LAMMPS run for 3-d periodic
#-------------------------------------------------------------------------------

units			real
boundary		p p p
atom_style		full

pair_style		lj/class2/coul/long 12
pair_modify             mix sixthpower
pair_modify		tail yes
bond_style		class2
angle_style		class2
dihedral_style		class2
improper_style		class2
special_bonds           lj 0.0 0.0 1.0 coul 0.0 0.0 1.0

box                     tilt large
read_data		structure.dat

include			parameters.dat

neighbor		2.0 bin
neigh_modify		delay 0 every 1 check yes
kspace_style		pppm 0.00001

variable		R		equal 0.00198722
variable		sysvol		equal vol
variable		sysmass		equal mass(all)/6.0221367e+23
variable		sysdensity	equal v_sysmass/v_sysvol/1.0e-24
variable		coulomb		equal ecoul+elong
variable		etotal		equal etotal
variable		pe		equal pe
variable		ke		equal ke
variable		evdwl		equal evdwl
variable		epair		equal epair
variable		ebond		equal ebond
variable		eangle		equal eangle
variable		edihed		equal edihed
variable		eimp		equal eimp
variable		lx		equal lx
variable		ly		equal ly
variable		lz		equal lz
variable		Nthermo		equal 0
variable		cella		equal lx
variable		cellb		equal sqrt(ly*ly+xy*xy)
variable		cellc		equal sqrt(lz*lz+xz*xz+yz*yz)
variable		cellalpha	equal acos((xy*xz+ly*yz)/(v_cellb*v_cellc))
variable		cellbeta	equal acos(xz/v_cellc)
variable		cellgamma	equal acos(xy/v_cellb)
variable		p		equal press
variable		pxx		equal pxx
variable		pyy		equal pyy
variable		pzz		equal pzz
variable		pyz		equal pyz
variable		pxz		equal pxz
variable		pxy		equal pxy
variable		sxx		equal -pxx
variable		syy		equal -pyy
variable		szz		equal -pzz
variable		syz		equal -pyz
variable		sxz		equal -pxz
variable		sxy		equal -pxy
variable		fmax		equal fmax
variable	        fnorm		equal fnorm
variable		time equal step*dt+0.000001
variable		surfacetension equal 0.5*v_lz*(0.5*(v_sxx+v_syy)-v_szz)

thermo_style		custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke
thermo_modify		flush yes


#
# Set up the fixed and movable groups
#

group		    	movable union all
group		    	fixed subtract all movable

compute 		sum_f1 movable reduce sum fx fy fz
variable 		sum_fx equal c_sum_f1[1]
variable 		sum_fy equal c_sum_f1[2]
variable  		sum_fz equal c_sum_f1[3]

#
# Subsets
#

group		    	subset_Li id 1 2 3 4 5 6 7

group		    	subset_Na id 8 9 10 11 12 13 14 15 16 17
group		    	subset_Na id 18 19 20 21 22 23 24 25 26 27
group		    	subset_Na id 28 29 30 31 32 33 34 35 36 37
group		    	subset_Na id 38 39 40 41 42 43 44 45 46 47
group		    	subset_Na id 48 49 50 51 52 53 54 55 56 57
group		    	subset_Na id 58 59 60 61 62 63

group		    	subset_K id 64 65 66 67 68 69 70 71 72 73
group		    	subset_K id 74 75 76 77 78 79 80 81 82 83
group		    	subset_K id 84 85 86 87 88 89 90 91 92 93
group		    	subset_K id 94 95 96 97 98 99 100 101 102 103

group		    	subset_Mg id 104 105 106 107 108 109 110 111 112 113
group		    	subset_Mg id 114 115 116 117 118 119 120 121 122 123
group		    	subset_Mg id 124 125 126 127 128 129 130 131 132 133
group		    	subset_Mg id 134 135 136 137 138 139 140 141 142 143
group		    	subset_Mg id 144 145 146 147 148 149 150 151 152 153
group		    	subset_Mg id 154 155 156 157 158 159 160 161 162 163
group		    	subset_Mg id 164 165 166 167 168 169 170 171 172 173
group		    	subset_Mg id 174 175 176 177 178 179 180 181 182 183
group		    	subset_Mg id 184 185 186 187 188 189 190 191 192 193
group		    	subset_Mg id 194 195 196 197 198 199 200 201 202 203
group		    	subset_Mg id 204 205 206 207 208 209 210 211 212 213
group		    	subset_Mg id 214 215 216 217 218 219 220 221 222 223
group		    	subset_Mg id 224 225 226 227 228 229 230 231 232 233
group		    	subset_Mg id 234 235 236 237 238 239 240 241 242 243
group		    	subset_Mg id 244 245 246 247 248 249 250 251 252 253
group		    	subset_Mg id 254

group		    	subset_Cl id 255 256 257 258 259 260 261 262 263 264
group		    	subset_Cl id 265 266 267 268 269 270 271 272 273 274
group		    	subset_Cl id 275 276 277 278 279 280 281 282 283 284
group		    	subset_Cl id 285 286 287 288 289 290 291 292 293 294
group		    	subset_Cl id 295 296 297 298 299 300 301 302 303 304
group		    	subset_Cl id 305 306 307 308 309 310 311 312 313 314
group		    	subset_Cl id 315 316 317 318 319 320 321 322 323 324
group		    	subset_Cl id 325 326 327 328 329 330 331 332 333 334
group		    	subset_Cl id 335 336 337 338 339 340 341 342 343 344
group		    	subset_Cl id 345 346 347 348 349 350 351 352 353 354
group		    	subset_Cl id 355 356 357 358 359 360 361 362 363 364
group		    	subset_Cl id 365 366 367 368 369 370 371 372 373 374
group		    	subset_Cl id 375 376 377 378 379 380 381 382 383 384
group		    	subset_Cl id 385 386 387 388 389 390 391 392 393 394
group		    	subset_Cl id 395 396 397 398 399 400 401 402 403 404
group		    	subset_Cl id 405 406 407 408 409 410 411 412 413 414
group		    	subset_Cl id 415 416 417 418 419 420 421 422 423 424
group		    	subset_Cl id 425 426 427 428 429 430 431 432 433 434
group		    	subset_Cl id 435 436 437 438 439 440 441 442 443 444
group		    	subset_Cl id 445 446 447 448 449 450 451 452 453 454
group		    	subset_Cl id 455 456 457 458 459 460 461 462 463 464
group		    	subset_Cl id 465 466 467 468 469 470 471 472 473 474
group		    	subset_Cl id 475 476 477 478 479 480 481 482 483 484
group		    	subset_Cl id 485 486 487 488 489 490 491 492 493 494
group		    	subset_Cl id 495 496 497 498 499 500 501 502 503 504
group		    	subset_Cl id 505 506 507 508 509 510 511 512 513 514
group		    	subset_Cl id 515 516 517 518 519 520 521 522 523 524
group		    	subset_Cl id 525 526 527 528 529 530 531 532 533 534
group		    	subset_Cl id 535 536 537 538 539 540 541 542 543 544
group		    	subset_Cl id 545 546 547 548 549 550 551 552 553 554
group		    	subset_Cl id 555 556 557 558 559 560 561 562 563 564
group		    	subset_Cl id 565 566 567 568 569 570 571 572 573 574
group		    	subset_Cl id 575 576 577 578 579 580 581 582 583 584
group		    	subset_Cl id 585 586 587 588 589 590 591 592 593 594
group		    	subset_Cl id 595 596 597

group		    	subset_O(SO4) id 599 600 601 602 604 605 606 607 609 610
group		    	subset_O(SO4) id 611 612 614 615 616 617 619 620 621 622
group		    	subset_O(SO4) id 624 625 626 627 629 630 631 632 634 635
group		    	subset_O(SO4) id 636 637 639 640 641 642 644 645 646 647
group		    	subset_O(SO4) id 649 650 651 652 654 655 656 657 659 660
group		    	subset_O(SO4) id 661 662 664 665 666 667 669 670 671 672
group		    	subset_O(SO4) id 674 675 676 677 679 680 681 682 684 685
group		    	subset_O(SO4) id 686 687 689 690 691 692 694 695 696 697
group		    	subset_O(SO4) id 699 700 701 702 704 705 706 707 709 710
group		    	subset_O(SO4) id 711 712 714 715 716 717 719 720 721 722
group		    	subset_O(SO4) id 724 725 726 727 729 730 731 732 734 735
group		    	subset_O(SO4) id 736 737 739 740 741 742 744 745 746 747
group		    	subset_O(SO4) id 749 750 751 752

group		    	subset_S(SO4) id 598 603 608 613 618 623 628 633 638 643
group		    	subset_S(SO4) id 648 653 658 663 668 673 678 683 688 693
group		    	subset_S(SO4) id 698 703 708 713 718 723 728 733 738 743
group		    	subset_S(SO4) id 748

group		    	subset_O(H2O) id 753 756 759 762 765 768 771 774 777 780
group		    	subset_O(H2O) id 783 786 789 792 795 798 801 804 807 810
group		    	subset_O(H2O) id 813 816 819 822 825 828 831 834 837 840
group		    	subset_O(H2O) id 843 846 849 852 855 858 861 864 867 870
group		    	subset_O(H2O) id 873 876 879 882 885 888 891 894 897 900
group		    	subset_O(H2O) id 903 906 909 912 915 918 921 924 927 930
group		    	subset_O(H2O) id 933 936 939 942 945 948 951 954 957 960
group		    	subset_O(H2O) id 963 966 969 972 975 978 981 984 987 990
group		    	subset_O(H2O) id 993 996 999 1002 1005 1008 1011 1014 1017 1020
group		    	subset_O(H2O) id 1023 1026 1029 1032 1035 1038 1041 1044 1047 1050
group		    	subset_O(H2O) id 1053 1056 1059 1062 1065 1068 1071 1074 1077 1080
group		    	subset_O(H2O) id 1083 1086 1089 1092 1095 1098 1101 1104 1107 1110
group		    	subset_O(H2O) id 1113 1116 1119 1122 1125 1128 1131 1134 1137 1140
group		    	subset_O(H2O) id 1143 1146 1149 1152 1155 1158 1161 1164 1167 1170
group		    	subset_O(H2O) id 1173 1176 1179 1182 1185 1188 1191 1194 1197 1200
group		    	subset_O(H2O) id 1203 1206 1209 1212 1215 1218 1221 1224 1227 1230
group		    	subset_O(H2O) id 1233 1236 1239 1242 1245 1248 1251 1254 1257 1260
group		    	subset_O(H2O) id 1263 1266 1269 1272 1275 1278 1281 1284 1287 1290
group		    	subset_O(H2O) id 1293 1296 1299 1302 1305 1308 1311 1314 1317 1320
group		    	subset_O(H2O) id 1323 1326 1329 1332 1335 1338 1341 1344 1347 1350
group		    	subset_O(H2O) id 1353 1356 1359 1362 1365 1368 1371 1374 1377 1380
group		    	subset_O(H2O) id 1383 1386 1389 1392 1395 1398 1401 1404 1407 1410
group		    	subset_O(H2O) id 1413 1416 1419 1422 1425 1428 1431 1434 1437 1440
group		    	subset_O(H2O) id 1443 1446 1449 1452 1455 1458 1461 1464 1467 1470
group		    	subset_O(H2O) id 1473 1476 1479 1482 1485 1488 1491 1494 1497 1500
group		    	subset_O(H2O) id 1503 1506 1509 1512 1515 1518 1521 1524 1527 1530
group		    	subset_O(H2O) id 1533 1536 1539 1542 1545 1548 1551 1554 1557 1560
group		    	subset_O(H2O) id 1563 1566 1569 1572 1575 1578 1581 1584 1587 1590
group		    	subset_O(H2O) id 1593 1596 1599 1602 1605 1608 1611 1614 1617 1620
group		    	subset_O(H2O) id 1623 1626 1629 1632 1635 1638 1641 1644 1647 1650
group		    	subset_O(H2O) id 1653 1656 1659 1662 1665 1668 1671 1674 1677 1680
group		    	subset_O(H2O) id 1683 1686 1689 1692 1695 1698 1701 1704 1707 1710
group		    	subset_O(H2O) id 1713 1716 1719 1722 1725 1728 1731 1734 1737 1740
group		    	subset_O(H2O) id 1743 1746 1749 1752 1755 1758 1761 1764 1767 1770
group		    	subset_O(H2O) id 1773 1776 1779 1782 1785 1788 1791 1794 1797 1800
group		    	subset_O(H2O) id 1803 1806 1809 1812 1815 1818 1821 1824 1827 1830
group		    	subset_O(H2O) id 1833 1836 1839 1842 1845 1848 1851 1854 1857 1860
group		    	subset_O(H2O) id 1863 1866 1869 1872 1875 1878 1881 1884 1887 1890
group		    	subset_O(H2O) id 1893 1896 1899 1902 1905 1908 1911 1914 1917 1920
group		    	subset_O(H2O) id 1923 1926 1929 1932 1935 1938 1941 1944 1947 1950
group		    	subset_O(H2O) id 1953 1956 1959 1962 1965 1968 1971 1974 1977 1980
group		    	subset_O(H2O) id 1983 1986 1989 1992 1995 1998 2001 2004 2007 2010
group		    	subset_O(H2O) id 2013 2016 2019 2022 2025 2028 2031 2034 2037 2040
group		    	subset_O(H2O) id 2043 2046 2049 2052 2055 2058 2061 2064 2067 2070
group		    	subset_O(H2O) id 2073 2076 2079 2082 2085 2088 2091 2094 2097 2100
group		    	subset_O(H2O) id 2103 2106 2109 2112 2115 2118 2121 2124 2127 2130
group		    	subset_O(H2O) id 2133 2136 2139 2142 2145 2148 2151 2154 2157 2160
group		    	subset_O(H2O) id 2163 2166 2169 2172 2175 2178 2181 2184 2187 2190
group		    	subset_O(H2O) id 2193 2196 2199 2202 2205 2208 2211 2214 2217 2220
group		    	subset_O(H2O) id 2223 2226 2229 2232 2235 2238 2241 2244 2247 2250
group		    	subset_O(H2O) id 2253 2256 2259 2262 2265 2268 2271 2274 2277 2280
group		    	subset_O(H2O) id 2283 2286 2289 2292 2295 2298 2301 2304 2307 2310
group		    	subset_O(H2O) id 2313 2316 2319 2322 2325 2328 2331 2334 2337 2340
group		    	subset_O(H2O) id 2343 2346 2349 2352 2355 2358 2361 2364 2367 2370
group		    	subset_O(H2O) id 2373 2376 2379 2382 2385 2388 2391 2394 2397 2400
group		    	subset_O(H2O) id 2403 2406 2409 2412 2415 2418 2421 2424 2427 2430
group		    	subset_O(H2O) id 2433 2436 2439 2442 2445 2448 2451 2454 2457 2460
group		    	subset_O(H2O) id 2463 2466 2469 2472 2475 2478 2481 2484 2487 2490
group		    	subset_O(H2O) id 2493 2496 2499 2502 2505 2508 2511 2514 2517 2520
group		    	subset_O(H2O) id 2523 2526 2529 2532 2535 2538 2541 2544 2547 2550
group		    	subset_O(H2O) id 2553 2556 2559 2562 2565 2568 2571 2574 2577 2580
group		    	subset_O(H2O) id 2583 2586 2589 2592 2595 2598 2601 2604 2607 2610
group		    	subset_O(H2O) id 2613 2616 2619 2622 2625 2628 2631 2634 2637 2640
group		    	subset_O(H2O) id 2643 2646 2649 2652 2655 2658 2661 2664 2667 2670
group		    	subset_O(H2O) id 2673 2676 2679 2682 2685 2688 2691 2694 2697 2700
group		    	subset_O(H2O) id 2703 2706 2709 2712 2715 2718 2721 2724 2727 2730
group		    	subset_O(H2O) id 2733 2736 2739 2742 2745 2748 2751 2754 2757 2760
group		    	subset_O(H2O) id 2763 2766 2769 2772 2775 2778 2781 2784 2787 2790
group		    	subset_O(H2O) id 2793 2796 2799 2802 2805 2808 2811 2814 2817 2820
group		    	subset_O(H2O) id 2823 2826 2829 2832 2835 2838 2841 2844 2847 2850
group		    	subset_O(H2O) id 2853 2856 2859 2862 2865 2868 2871 2874 2877 2880
group		    	subset_O(H2O) id 2883 2886 2889 2892 2895 2898 2901 2904 2907 2910
group		    	subset_O(H2O) id 2913 2916 2919 2922 2925 2928 2931 2934 2937 2940
group		    	subset_O(H2O) id 2943 2946 2949 2952 2955 2958 2961 2964 2967 2970
group		    	subset_O(H2O) id 2973 2976 2979 2982 2985 2988 2991 2994 2997 3000
group		    	subset_O(H2O) id 3003 3006 3009 3012 3015 3018 3021 3024 3027 3030
group		    	subset_O(H2O) id 3033 3036 3039 3042 3045 3048 3051 3054 3057 3060
group		    	subset_O(H2O) id 3063 3066 3069 3072 3075 3078 3081 3084 3087 3090
group		    	subset_O(H2O) id 3093 3096 3099 3102 3105 3108 3111 3114 3117 3120
group		    	subset_O(H2O) id 3123 3126 3129 3132 3135 3138 3141 3144 3147 3150
group		    	subset_O(H2O) id 3153 3156 3159 3162 3165 3168 3171 3174 3177 3180
group		    	subset_O(H2O) id 3183 3186 3189 3192 3195 3198 3201 3204 3207 3210
group		    	subset_O(H2O) id 3213 3216 3219 3222 3225 3228 3231 3234 3237 3240
group		    	subset_O(H2O) id 3243 3246 3249 3252 3255 3258 3261 3264 3267 3270
group		    	subset_O(H2O) id 3273 3276 3279 3282 3285 3288 3291 3294 3297 3300
group		    	subset_O(H2O) id 3303 3306 3309 3312 3315 3318 3321 3324 3327 3330
group		    	subset_O(H2O) id 3333 3336 3339 3342 3345 3348 3351 3354 3357 3360
group		    	subset_O(H2O) id 3363 3366 3369 3372 3375 3378 3381 3384 3387 3390
group		    	subset_O(H2O) id 3393 3396 3399 3402 3405 3408 3411 3414 3417 3420
group		    	subset_O(H2O) id 3423 3426 3429 3432 3435 3438 3441 3444 3447 3450
group		    	subset_O(H2O) id 3453 3456 3459 3462 3465 3468 3471 3474 3477 3480
group		    	subset_O(H2O) id 3483 3486 3489 3492 3495 3498 3501 3504 3507 3510
group		    	subset_O(H2O) id 3513 3516 3519 3522 3525 3528 3531 3534 3537 3540
group		    	subset_O(H2O) id 3543 3546 3549 3552 3555 3558 3561 3564 3567 3570
group		    	subset_O(H2O) id 3573 3576 3579 3582 3585 3588 3591 3594 3597 3600
group		    	subset_O(H2O) id 3603 3606 3609 3612 3615 3618 3621 3624 3627 3630
group		    	subset_O(H2O) id 3633 3636 3639 3642 3645 3648 3651 3654 3657 3660
group		    	subset_O(H2O) id 3663 3666 3669 3672 3675 3678 3681 3684 3687 3690
group		    	subset_O(H2O) id 3693 3696 3699 3702 3705 3708 3711 3714 3717 3720
group		    	subset_O(H2O) id 3723 3726 3729 3732 3735 3738 3741 3744 3747 3750
group		    	subset_O(H2O) id 3753 3756 3759 3762 3765 3768 3771 3774 3777 3780
group		    	subset_O(H2O) id 3783 3786 3789 3792 3795 3798 3801 3804 3807 3810
group		    	subset_O(H2O) id 3813 3816 3819 3822 3825 3828 3831 3834 3837 3840
group		    	subset_O(H2O) id 3843 3846 3849 3852 3855 3858 3861 3864 3867 3870
group		    	subset_O(H2O) id 3873 3876 3879 3882 3885 3888 3891 3894 3897 3900
group		    	subset_O(H2O) id 3903 3906 3909 3912 3915 3918 3921 3924 3927 3930
group		    	subset_O(H2O) id 3933 3936 3939 3942 3945 3948 3951 3954 3957 3960
group		    	subset_O(H2O) id 3963 3966 3969 3972 3975 3978 3981 3984 3987 3990
group		    	subset_O(H2O) id 3993 3996 3999 4002 4005 4008 4011 4014 4017 4020
group		    	subset_O(H2O) id 4023 4026 4029 4032 4035 4038 4041 4044 4047 4050
group		    	subset_O(H2O) id 4053 4056 4059 4062 4065 4068 4071 4074 4077 4080
group		    	subset_O(H2O) id 4083 4086 4089 4092 4095 4098 4101 4104 4107 4110
group		    	subset_O(H2O) id 4113 4116 4119 4122 4125 4128 4131 4134 4137 4140
group		    	subset_O(H2O) id 4143 4146 4149 4152 4155 4158 4161 4164 4167 4170
group		    	subset_O(H2O) id 4173 4176 4179 4182 4185 4188 4191 4194 4197 4200
group		    	subset_O(H2O) id 4203 4206 4209 4212 4215 4218 4221 4224 4227 4230
group		    	subset_O(H2O) id 4233 4236 4239 4242 4245 4248 4251 4254 4257 4260
group		    	subset_O(H2O) id 4263 4266 4269 4272 4275 4278 4281 4284 4287 4290
group		    	subset_O(H2O) id 4293 4296 4299 4302 4305 4308 4311 4314 4317 4320
group		    	subset_O(H2O) id 4323 4326 4329 4332 4335 4338 4341 4344 4347 4350
group		    	subset_O(H2O) id 4353 4356 4359 4362 4365 4368 4371 4374 4377 4380
group		    	subset_O(H2O) id 4383 4386 4389 4392 4395 4398 4401 4404 4407 4410
group		    	subset_O(H2O) id 4413 4416 4419 4422 4425 4428 4431 4434 4437 4440
group		    	subset_O(H2O) id 4443 4446 4449 4452 4455 4458 4461 4464 4467 4470
group		    	subset_O(H2O) id 4473 4476 4479 4482 4485 4488 4491 4494 4497 4500
group		    	subset_O(H2O) id 4503 4506 4509 4512 4515 4518 4521 4524 4527 4530
group		    	subset_O(H2O) id 4533 4536 4539 4542 4545 4548 4551 4554 4557 4560
group		    	subset_O(H2O) id 4563 4566 4569 4572 4575 4578 4581 4584 4587 4590
group		    	subset_O(H2O) id 4593 4596 4599 4602 4605 4608 4611 4614 4617 4620
group		    	subset_O(H2O) id 4623 4626 4629 4632 4635 4638 4641 4644 4647 4650
group		    	subset_O(H2O) id 4653 4656 4659 4662 4665 4668 4671 4674 4677 4680
group		    	subset_O(H2O) id 4683 4686 4689 4692 4695 4698 4701 4704 4707 4710
group		    	subset_O(H2O) id 4713 4716 4719 4722 4725 4728 4731 4734 4737 4740
group		    	subset_O(H2O) id 4743 4746 4749 4752 4755 4758 4761 4764 4767 4770
group		    	subset_O(H2O) id 4773 4776 4779 4782 4785 4788 4791 4794 4797 4800
group		    	subset_O(H2O) id 4803 4806 4809 4812 4815 4818 4821 4824 4827 4830
group		    	subset_O(H2O) id 4833 4836 4839 4842 4845 4848 4851 4854 4857 4860
group		    	subset_O(H2O) id 4863 4866 4869 4872 4875 4878 4881 4884 4887 4890
group		    	subset_O(H2O) id 4893 4896 4899 4902 4905 4908 4911 4914 4917 4920
group		    	subset_O(H2O) id 4923 4926 4929 4932 4935 4938 4941 4944 4947 4950
group		    	subset_O(H2O) id 4953 4956 4959 4962 4965 4968 4971 4974 4977 4980
group		    	subset_O(H2O) id 4983 4986 4989 4992 4995 4998 5001 5004 5007 5010
group		    	subset_O(H2O) id 5013 5016 5019 5022 5025 5028 5031 5034 5037 5040
group		    	subset_O(H2O) id 5043 5046 5049 5052 5055 5058 5061 5064 5067 5070
group		    	subset_O(H2O) id 5073 5076 5079 5082 5085 5088 5091 5094 5097 5100
group		    	subset_O(H2O) id 5103 5106 5109 5112 5115 5118 5121 5124 5127 5130
group		    	subset_O(H2O) id 5133 5136 5139 5142 5145 5148 5151 5154 5157 5160
group		    	subset_O(H2O) id 5163 5166 5169 5172 5175 5178 5181 5184 5187 5190
group		    	subset_O(H2O) id 5193 5196 5199 5202 5205 5208 5211 5214 5217 5220
group		    	subset_O(H2O) id 5223 5226 5229 5232 5235 5238 5241 5244 5247 5250
group		    	subset_O(H2O) id 5253 5256 5259 5262 5265 5268 5271 5274 5277 5280
group		    	subset_O(H2O) id 5283 5286 5289 5292 5295 5298 5301 5304 5307 5310
group		    	subset_O(H2O) id 5313 5316 5319 5322 5325 5328 5331 5334 5337 5340
group		    	subset_O(H2O) id 5343 5346 5349 5352 5355 5358 5361 5364 5367 5370
group		    	subset_O(H2O) id 5373 5376 5379 5382 5385 5388 5391 5394 5397 5400
group		    	subset_O(H2O) id 5403 5406 5409 5412 5415 5418 5421 5424 5427 5430
group		    	subset_O(H2O) id 5433 5436 5439 5442 5445 5448 5451 5454 5457 5460
group		    	subset_O(H2O) id 5463 5466 5469 5472 5475 5478 5481 5484 5487 5490
group		    	subset_O(H2O) id 5493 5496 5499 5502 5505 5508 5511 5514 5517 5520
group		    	subset_O(H2O) id 5523 5526 5529 5532 5535 5538 5541 5544 5547 5550
group		    	subset_O(H2O) id 5553 5556 5559 5562 5565 5568 5571 5574 5577 5580
group		    	subset_O(H2O) id 5583 5586 5589 5592 5595 5598 5601 5604 5607 5610
group		    	subset_O(H2O) id 5613 5616 5619 5622 5625 5628 5631 5634 5637 5640
group		    	subset_O(H2O) id 5643 5646 5649 5652 5655 5658 5661 5664 5667 5670
group		    	subset_O(H2O) id 5673 5676 5679 5682 5685 5688 5691 5694 5697 5700
group		    	subset_O(H2O) id 5703 5706 5709 5712 5715 5718 5721 5724 5727 5730
group		    	subset_O(H2O) id 5733 5736 5739 5742 5745 5748 5751 5754 5757 5760
group		    	subset_O(H2O) id 5763 5766 5769 5772 5775 5778 5781 5784 5787 5790
group		    	subset_O(H2O) id 5793 5796 5799 5802 5805 5808 5811 5814 5817 5820
group		    	subset_O(H2O) id 5823 5826 5829 5832 5835 5838 5841 5844 5847 5850
group		    	subset_O(H2O) id 5853 5856 5859 5862 5865 5868 5871 5874 5877 5880
group		    	subset_O(H2O) id 5883 5886 5889 5892 5895 5898 5901 5904 5907 5910
group		    	subset_O(H2O) id 5913 5916 5919 5922 5925 5928 5931 5934 5937 5940
group		    	subset_O(H2O) id 5943 5946 5949 5952 5955 5958 5961 5964 5967 5970
group		    	subset_O(H2O) id 5973 5976 5979 5982 5985 5988 5991 5994 5997 6000
group		    	subset_O(H2O) id 6003 6006 6009 6012 6015 6018 6021 6024 6027 6030
group		    	subset_O(H2O) id 6033 6036 6039 6042 6045 6048 6051 6054 6057 6060
group		    	subset_O(H2O) id 6063 6066 6069 6072 6075 6078 6081 6084 6087 6090
group		    	subset_O(H2O) id 6093 6096 6099 6102 6105 6108 6111 6114 6117 6120
group		    	subset_O(H2O) id 6123 6126 6129 6132 6135 6138 6141 6144 6147 6150
group		    	subset_O(H2O) id 6153 6156 6159 6162 6165 6168 6171 6174 6177 6180
group		    	subset_O(H2O) id 6183 6186 6189 6192 6195 6198 6201 6204 6207 6210
group		    	subset_O(H2O) id 6213 6216 6219 6222 6225 6228 6231 6234 6237 6240
group		    	subset_O(H2O) id 6243 6246 6249 6252 6255 6258 6261 6264 6267 6270
group		    	subset_O(H2O) id 6273 6276 6279 6282 6285 6288 6291 6294 6297 6300
group		    	subset_O(H2O) id 6303 6306 6309 6312 6315 6318 6321 6324 6327 6330
group		    	subset_O(H2O) id 6333 6336 6339 6342 6345 6348 6351 6354 6357 6360
group		    	subset_O(H2O) id 6363 6366 6369 6372 6375 6378 6381 6384 6387 6390
group		    	subset_O(H2O) id 6393 6396 6399 6402 6405 6408 6411 6414 6417 6420
group		    	subset_O(H2O) id 6423 6426 6429 6432 6435 6438 6441 6444 6447 6450
group		    	subset_O(H2O) id 6453 6456 6459 6462 6465 6468 6471 6474 6477 6480
group		    	subset_O(H2O) id 6483 6486 6489 6492 6495 6498 6501 6504 6507 6510
group		    	subset_O(H2O) id 6513 6516 6519 6522 6525 6528 6531 6534 6537 6540
group		    	subset_O(H2O) id 6543 6546 6549 6552 6555 6558 6561 6564 6567 6570
group		    	subset_O(H2O) id 6573 6576 6579 6582 6585 6588 6591 6594 6597 6600
group		    	subset_O(H2O) id 6603 6606 6609 6612 6615 6618 6621 6624 6627 6630
group		    	subset_O(H2O) id 6633 6636 6639 6642 6645 6648 6651 6654 6657 6660
group		    	subset_O(H2O) id 6663 6666 6669 6672 6675 6678 6681 6684 6687 6690
group		    	subset_O(H2O) id 6693 6696 6699 6702 6705 6708 6711 6714 6717 6720
group		    	subset_O(H2O) id 6723 6726 6729 6732 6735 6738 6741 6744 6747 6750
group		    	subset_O(H2O) id 6753 6756 6759 6762 6765 6768 6771 6774 6777 6780
group		    	subset_O(H2O) id 6783 6786 6789 6792 6795 6798 6801 6804 6807 6810
group		    	subset_O(H2O) id 6813 6816 6819 6822 6825 6828 6831 6834 6837 6840
group		    	subset_O(H2O) id 6843 6846 6849 6852 6855 6858 6861 6864 6867 6870
group		    	subset_O(H2O) id 6873 6876 6879 6882 6885 6888 6891 6894 6897 6900
group		    	subset_O(H2O) id 6903 6906 6909 6912 6915 6918 6921 6924 6927 6930
group		    	subset_O(H2O) id 6933 6936 6939 6942 6945 6948 6951 6954 6957 6960
group		    	subset_O(H2O) id 6963 6966 6969 6972 6975 6978 6981 6984 6987 6990
group		    	subset_O(H2O) id 6993 6996 6999 7002 7005 7008 7011 7014 7017 7020
group		    	subset_O(H2O) id 7023 7026 7029 7032 7035 7038 7041 7044 7047 7050
group		    	subset_O(H2O) id 7053 7056 7059 7062 7065 7068 7071 7074 7077 7080
group		    	subset_O(H2O) id 7083 7086 7089 7092 7095 7098 7101 7104 7107 7110
group		    	subset_O(H2O) id 7113 7116 7119 7122 7125 7128 7131 7134 7137 7140
group		    	subset_O(H2O) id 7143 7146 7149 7152 7155 7158 7161 7164 7167 7170
group		    	subset_O(H2O) id 7173 7176 7179 7182 7185 7188 7191 7194 7197 7200
group		    	subset_O(H2O) id 7203 7206 7209 7212 7215 7218 7221 7224 7227 7230
group		    	subset_O(H2O) id 7233 7236 7239 7242 7245 7248 7251 7254 7257 7260
group		    	subset_O(H2O) id 7263 7266 7269 7272 7275 7278 7281 7284 7287 7290
group		    	subset_O(H2O) id 7293 7296 7299 7302 7305 7308 7311 7314 7317 7320
group		    	subset_O(H2O) id 7323 7326 7329 7332 7335 7338 7341 7344 7347 7350
group		    	subset_O(H2O) id 7353 7356 7359 7362 7365 7368 7371 7374 7377 7380
group		    	subset_O(H2O) id 7383 7386 7389 7392 7395 7398 7401 7404 7407 7410
group		    	subset_O(H2O) id 7413 7416 7419 7422 7425 7428 7431 7434 7437 7440
group		    	subset_O(H2O) id 7443 7446 7449 7452 7455 7458 7461 7464 7467 7470
group		    	subset_O(H2O) id 7473 7476 7479 7482 7485 7488 7491 7494 7497 7500
group		    	subset_O(H2O) id 7503 7506 7509 7512 7515 7518 7521 7524 7527 7530
group		    	subset_O(H2O) id 7533 7536 7539 7542 7545 7548 7551 7554 7557 7560
group		    	subset_O(H2O) id 7563 7566 7569 7572 7575 7578 7581 7584 7587 7590
group		    	subset_O(H2O) id 7593 7596 7599 7602 7605 7608 7611 7614 7617 7620
group		    	subset_O(H2O) id 7623 7626 7629 7632 7635 7638 7641 7644 7647 7650
group		    	subset_O(H2O) id 7653 7656 7659 7662 7665 7668 7671 7674 7677 7680
group		    	subset_O(H2O) id 7683 7686 7689 7692 7695 7698 7701 7704 7707 7710
group		    	subset_O(H2O) id 7713 7716 7719 7722 7725 7728 7731 7734 7737 7740
group		    	subset_O(H2O) id 7743 7746 7749 7752 7755 7758 7761 7764 7767 7770
group		    	subset_O(H2O) id 7773 7776 7779 7782 7785 7788 7791 7794 7797 7800
group		    	subset_O(H2O) id 7803 7806 7809 7812 7815 7818 7821 7824 7827 7830
group		    	subset_O(H2O) id 7833 7836 7839 7842 7845 7848 7851 7854 7857 7860
group		    	subset_O(H2O) id 7863 7866 7869 7872 7875 7878 7881 7884 7887 7890
group		    	subset_O(H2O) id 7893 7896 7899 7902 7905 7908 7911 7914 7917 7920
group		    	subset_O(H2O) id 7923 7926 7929 7932 7935 7938 7941 7944 7947 7950
group		    	subset_O(H2O) id 7953 7956 7959 7962 7965 7968 7971 7974 7977 7980
group		    	subset_O(H2O) id 7983 7986 7989 7992 7995 7998 8001 8004 8007 8010
group		    	subset_O(H2O) id 8013 8016 8019 8022 8025 8028 8031 8034 8037 8040
group		    	subset_O(H2O) id 8043 8046 8049 8052 8055 8058 8061 8064 8067 8070
group		    	subset_O(H2O) id 8073 8076 8079 8082 8085 8088 8091 8094 8097 8100
group		    	subset_O(H2O) id 8103 8106 8109 8112 8115 8118 8121 8124 8127 8130
group		    	subset_O(H2O) id 8133 8136 8139 8142 8145 8148 8151 8154 8157 8160
group		    	subset_O(H2O) id 8163 8166 8169 8172 8175 8178 8181 8184 8187 8190
group		    	subset_O(H2O) id 8193 8196 8199 8202 8205 8208 8211 8214 8217 8220
group		    	subset_O(H2O) id 8223 8226 8229 8232 8235 8238 8241 8244 8247 8250
group		    	subset_O(H2O) id 8253 8256 8259 8262 8265 8268 8271 8274 8277 8280
group		    	subset_O(H2O) id 8283 8286 8289 8292 8295 8298 8301 8304 8307 8310
group		    	subset_O(H2O) id 8313 8316 8319 8322 8325 8328 8331 8334 8337 8340
group		    	subset_O(H2O) id 8343 8346 8349 8352 8355 8358 8361 8364 8367 8370
group		    	subset_O(H2O) id 8373 8376 8379 8382 8385 8388 8391 8394 8397 8400
group		    	subset_O(H2O) id 8403 8406 8409 8412 8415 8418 8421 8424 8427 8430
group		    	subset_O(H2O) id 8433 8436 8439 8442 8445 8448 8451 8454 8457 8460
group		    	subset_O(H2O) id 8463 8466 8469 8472 8475 8478 8481 8484 8487 8490
group		    	subset_O(H2O) id 8493 8496 8499 8502 8505 8508 8511 8514 8517 8520
group		    	subset_O(H2O) id 8523 8526 8529 8532 8535 8538 8541 8544 8547 8550
group		    	subset_O(H2O) id 8553 8556 8559 8562 8565 8568 8571 8574 8577 8580
group		    	subset_O(H2O) id 8583 8586 8589 8592 8595 8598 8601 8604 8607 8610
group		    	subset_O(H2O) id 8613 8616 8619 8622 8625 8628 8631 8634 8637 8640
group		    	subset_O(H2O) id 8643 8646 8649 8652 8655 8658 8661 8664 8667 8670
group		    	subset_O(H2O) id 8673 8676 8679 8682 8685 8688 8691 8694 8697 8700
group		    	subset_O(H2O) id 8703 8706 8709 8712 8715 8718 8721 8724 8727 8730
group		    	subset_O(H2O) id 8733 8736 8739 8742 8745 8748 8751 8754 8757 8760
group		    	subset_O(H2O) id 8763 8766 8769 8772 8775 8778 8781 8784 8787 8790
group		    	subset_O(H2O) id 8793 8796 8799 8802 8805 8808 8811 8814 8817 8820
group		    	subset_O(H2O) id 8823 8826 8829 8832 8835 8838 8841 8844 8847 8850
group		    	subset_O(H2O) id 8853 8856 8859 8862 8865 8868 8871 8874 8877 8880
group		    	subset_O(H2O) id 8883 8886 8889 8892 8895 8898 8901 8904 8907 8910
group		    	subset_O(H2O) id 8913 8916 8919 8922 8925 8928 8931 8934 8937 8940
group		    	subset_O(H2O) id 8943 8946 8949 8952 8955 8958 8961 8964 8967 8970
group		    	subset_O(H2O) id 8973 8976 8979 8982 8985 8988 8991 8994 8997 9000
group		    	subset_O(H2O) id 9003 9006 9009 9012 9015 9018 9021 9024 9027 9030
group		    	subset_O(H2O) id 9033 9036 9039 9042 9045 9048 9051 9054 9057 9060
group		    	subset_O(H2O) id 9063 9066 9069 9072 9075 9078 9081 9084 9087 9090
group		    	subset_O(H2O) id 9093 9096 9099 9102 9105 9108 9111 9114 9117 9120
group		    	subset_O(H2O) id 9123 9126 9129 9132 9135 9138 9141 9144 9147 9150
group		    	subset_O(H2O) id 9153 9156 9159 9162 9165 9168 9171 9174 9177 9180
group		    	subset_O(H2O) id 9183 9186

group		    	subset_H(H2O) id 754 755 757 758 760 761 763 764 766 767
group		    	subset_H(H2O) id 769 770 772 773 775 776 778 779 781 782
group		    	subset_H(H2O) id 784 785 787 788 790 791 793 794 796 797
group		    	subset_H(H2O) id 799 800 802 803 805 806 808 809 811 812
group		    	subset_H(H2O) id 814 815 817 818 820 821 823 824 826 827
group		    	subset_H(H2O) id 829 830 832 833 835 836 838 839 841 842
group		    	subset_H(H2O) id 844 845 847 848 850 851 853 854 856 857
group		    	subset_H(H2O) id 859 860 862 863 865 866 868 869 871 872
group		    	subset_H(H2O) id 874 875 877 878 880 881 883 884 886 887
group		    	subset_H(H2O) id 889 890 892 893 895 896 898 899 901 902
group		    	subset_H(H2O) id 904 905 907 908 910 911 913 914 916 917
group		    	subset_H(H2O) id 919 920 922 923 925 926 928 929 931 932
group		    	subset_H(H2O) id 934 935 937 938 940 941 943 944 946 947
group		    	subset_H(H2O) id 949 950 952 953 955 956 958 959 961 962
group		    	subset_H(H2O) id 964 965 967 968 970 971 973 974 976 977
group		    	subset_H(H2O) id 979 980 982 983 985 986 988 989 991 992
group		    	subset_H(H2O) id 994 995 997 998 1000 1001 1003 1004 1006 1007
group		    	subset_H(H2O) id 1009 1010 1012 1013 1015 1016 1018 1019 1021 1022
group		    	subset_H(H2O) id 1024 1025 1027 1028 1030 1031 1033 1034 1036 1037
group		    	subset_H(H2O) id 1039 1040 1042 1043 1045 1046 1048 1049 1051 1052
group		    	subset_H(H2O) id 1054 1055 1057 1058 1060 1061 1063 1064 1066 1067
group		    	subset_H(H2O) id 1069 1070 1072 1073 1075 1076 1078 1079 1081 1082
group		    	subset_H(H2O) id 1084 1085 1087 1088 1090 1091 1093 1094 1096 1097
group		    	subset_H(H2O) id 1099 1100 1102 1103 1105 1106 1108 1109 1111 1112
group		    	subset_H(H2O) id 1114 1115 1117 1118 1120 1121 1123 1124 1126 1127
group		    	subset_H(H2O) id 1129 1130 1132 1133 1135 1136 1138 1139 1141 1142
group		    	subset_H(H2O) id 1144 1145 1147 1148 1150 1151 1153 1154 1156 1157
group		    	subset_H(H2O) id 1159 1160 1162 1163 1165 1166 1168 1169 1171 1172
group		    	subset_H(H2O) id 1174 1175 1177 1178 1180 1181 1183 1184 1186 1187
group		    	subset_H(H2O) id 1189 1190 1192 1193 1195 1196 1198 1199 1201 1202
group		    	subset_H(H2O) id 1204 1205 1207 1208 1210 1211 1213 1214 1216 1217
group		    	subset_H(H2O) id 1219 1220 1222 1223 1225 1226 1228 1229 1231 1232
group		    	subset_H(H2O) id 1234 1235 1237 1238 1240 1241 1243 1244 1246 1247
group		    	subset_H(H2O) id 1249 1250 1252 1253 1255 1256 1258 1259 1261 1262
group		    	subset_H(H2O) id 1264 1265 1267 1268 1270 1271 1273 1274 1276 1277
group		    	subset_H(H2O) id 1279 1280 1282 1283 1285 1286 1288 1289 1291 1292
group		    	subset_H(H2O) id 1294 1295 1297 1298 1300 1301 1303 1304 1306 1307
group		    	subset_H(H2O) id 1309 1310 1312 1313 1315 1316 1318 1319 1321 1322
group		    	subset_H(H2O) id 1324 1325 1327 1328 1330 1331 1333 1334 1336 1337
group		    	subset_H(H2O) id 1339 1340 1342 1343 1345 1346 1348 1349 1351 1352
group		    	subset_H(H2O) id 1354 1355 1357 1358 1360 1361 1363 1364 1366 1367
group		    	subset_H(H2O) id 1369 1370 1372 1373 1375 1376 1378 1379 1381 1382
group		    	subset_H(H2O) id 1384 1385 1387 1388 1390 1391 1393 1394 1396 1397
group		    	subset_H(H2O) id 1399 1400 1402 1403 1405 1406 1408 1409 1411 1412
group		    	subset_H(H2O) id 1414 1415 1417 1418 1420 1421 1423 1424 1426 1427
group		    	subset_H(H2O) id 1429 1430 1432 1433 1435 1436 1438 1439 1441 1442
group		    	subset_H(H2O) id 1444 1445 1447 1448 1450 1451 1453 1454 1456 1457
group		    	subset_H(H2O) id 1459 1460 1462 1463 1465 1466 1468 1469 1471 1472
group		    	subset_H(H2O) id 1474 1475 1477 1478 1480 1481 1483 1484 1486 1487
group		    	subset_H(H2O) id 1489 1490 1492 1493 1495 1496 1498 1499 1501 1502
group		    	subset_H(H2O) id 1504 1505 1507 1508 1510 1511 1513 1514 1516 1517
group		    	subset_H(H2O) id 1519 1520 1522 1523 1525 1526 1528 1529 1531 1532
group		    	subset_H(H2O) id 1534 1535 1537 1538 1540 1541 1543 1544 1546 1547
group		    	subset_H(H2O) id 1549 1550 1552 1553 1555 1556 1558 1559 1561 1562
group		    	subset_H(H2O) id 1564 1565 1567 1568 1570 1571 1573 1574 1576 1577
group		    	subset_H(H2O) id 1579 1580 1582 1583 1585 1586 1588 1589 1591 1592
group		    	subset_H(H2O) id 1594 1595 1597 1598 1600 1601 1603 1604 1606 1607
group		    	subset_H(H2O) id 1609 1610 1612 1613 1615 1616 1618 1619 1621 1622
group		    	subset_H(H2O) id 1624 1625 1627 1628 1630 1631 1633 1634 1636 1637
group		    	subset_H(H2O) id 1639 1640 1642 1643 1645 1646 1648 1649 1651 1652
group		    	subset_H(H2O) id 1654 1655 1657 1658 1660 1661 1663 1664 1666 1667
group		    	subset_H(H2O) id 1669 1670 1672 1673 1675 1676 1678 1679 1681 1682
group		    	subset_H(H2O) id 1684 1685 1687 1688 1690 1691 1693 1694 1696 1697
group		    	subset_H(H2O) id 1699 1700 1702 1703 1705 1706 1708 1709 1711 1712
group		    	subset_H(H2O) id 1714 1715 1717 1718 1720 1721 1723 1724 1726 1727
group		    	subset_H(H2O) id 1729 1730 1732 1733 1735 1736 1738 1739 1741 1742
group		    	subset_H(H2O) id 1744 1745 1747 1748 1750 1751 1753 1754 1756 1757
group		    	subset_H(H2O) id 1759 1760 1762 1763 1765 1766 1768 1769 1771 1772
group		    	subset_H(H2O) id 1774 1775 1777 1778 1780 1781 1783 1784 1786 1787
group		    	subset_H(H2O) id 1789 1790 1792 1793 1795 1796 1798 1799 1801 1802
group		    	subset_H(H2O) id 1804 1805 1807 1808 1810 1811 1813 1814 1816 1817
group		    	subset_H(H2O) id 1819 1820 1822 1823 1825 1826 1828 1829 1831 1832
group		    	subset_H(H2O) id 1834 1835 1837 1838 1840 1841 1843 1844 1846 1847
group		    	subset_H(H2O) id 1849 1850 1852 1853 1855 1856 1858 1859 1861 1862
group		    	subset_H(H2O) id 1864 1865 1867 1868 1870 1871 1873 1874 1876 1877
group		    	subset_H(H2O) id 1879 1880 1882 1883 1885 1886 1888 1889 1891 1892
group		    	subset_H(H2O) id 1894 1895 1897 1898 1900 1901 1903 1904 1906 1907
group		    	subset_H(H2O) id 1909 1910 1912 1913 1915 1916 1918 1919 1921 1922
group		    	subset_H(H2O) id 1924 1925 1927 1928 1930 1931 1933 1934 1936 1937
group		    	subset_H(H2O) id 1939 1940 1942 1943 1945 1946 1948 1949 1951 1952
group		    	subset_H(H2O) id 1954 1955 1957 1958 1960 1961 1963 1964 1966 1967
group		    	subset_H(H2O) id 1969 1970 1972 1973 1975 1976 1978 1979 1981 1982
group		    	subset_H(H2O) id 1984 1985 1987 1988 1990 1991 1993 1994 1996 1997
group		    	subset_H(H2O) id 1999 2000 2002 2003 2005 2006 2008 2009 2011 2012
group		    	subset_H(H2O) id 2014 2015 2017 2018 2020 2021 2023 2024 2026 2027
group		    	subset_H(H2O) id 2029 2030 2032 2033 2035 2036 2038 2039 2041 2042
group		    	subset_H(H2O) id 2044 2045 2047 2048 2050 2051 2053 2054 2056 2057
group		    	subset_H(H2O) id 2059 2060 2062 2063 2065 2066 2068 2069 2071 2072
group		    	subset_H(H2O) id 2074 2075 2077 2078 2080 2081 2083 2084 2086 2087
group		    	subset_H(H2O) id 2089 2090 2092 2093 2095 2096 2098 2099 2101 2102
group		    	subset_H(H2O) id 2104 2105 2107 2108 2110 2111 2113 2114 2116 2117
group		    	subset_H(H2O) id 2119 2120 2122 2123 2125 2126 2128 2129 2131 2132
group		    	subset_H(H2O) id 2134 2135 2137 2138 2140 2141 2143 2144 2146 2147
group		    	subset_H(H2O) id 2149 2150 2152 2153 2155 2156 2158 2159 2161 2162
group		    	subset_H(H2O) id 2164 2165 2167 2168 2170 2171 2173 2174 2176 2177
group		    	subset_H(H2O) id 2179 2180 2182 2183 2185 2186 2188 2189 2191 2192
group		    	subset_H(H2O) id 2194 2195 2197 2198 2200 2201 2203 2204 2206 2207
group		    	subset_H(H2O) id 2209 2210 2212 2213 2215 2216 2218 2219 2221 2222
group		    	subset_H(H2O) id 2224 2225 2227 2228 2230 2231 2233 2234 2236 2237
group		    	subset_H(H2O) id 2239 2240 2242 2243 2245 2246 2248 2249 2251 2252
group		    	subset_H(H2O) id 2254 2255 2257 2258 2260 2261 2263 2264 2266 2267
group		    	subset_H(H2O) id 2269 2270 2272 2273 2275 2276 2278 2279 2281 2282
group		    	subset_H(H2O) id 2284 2285 2287 2288 2290 2291 2293 2294 2296 2297
group		    	subset_H(H2O) id 2299 2300 2302 2303 2305 2306 2308 2309 2311 2312
group		    	subset_H(H2O) id 2314 2315 2317 2318 2320 2321 2323 2324 2326 2327
group		    	subset_H(H2O) id 2329 2330 2332 2333 2335 2336 2338 2339 2341 2342
group		    	subset_H(H2O) id 2344 2345 2347 2348 2350 2351 2353 2354 2356 2357
group		    	subset_H(H2O) id 2359 2360 2362 2363 2365 2366 2368 2369 2371 2372
group		    	subset_H(H2O) id 2374 2375 2377 2378 2380 2381 2383 2384 2386 2387
group		    	subset_H(H2O) id 2389 2390 2392 2393 2395 2396 2398 2399 2401 2402
group		    	subset_H(H2O) id 2404 2405 2407 2408 2410 2411 2413 2414 2416 2417
group		    	subset_H(H2O) id 2419 2420 2422 2423 2425 2426 2428 2429 2431 2432
group		    	subset_H(H2O) id 2434 2435 2437 2438 2440 2441 2443 2444 2446 2447
group		    	subset_H(H2O) id 2449 2450 2452 2453 2455 2456 2458 2459 2461 2462
group		    	subset_H(H2O) id 2464 2465 2467 2468 2470 2471 2473 2474 2476 2477
group		    	subset_H(H2O) id 2479 2480 2482 2483 2485 2486 2488 2489 2491 2492
group		    	subset_H(H2O) id 2494 2495 2497 2498 2500 2501 2503 2504 2506 2507
group		    	subset_H(H2O) id 2509 2510 2512 2513 2515 2516 2518 2519 2521 2522
group		    	subset_H(H2O) id 2524 2525 2527 2528 2530 2531 2533 2534 2536 2537
group		    	subset_H(H2O) id 2539 2540 2542 2543 2545 2546 2548 2549 2551 2552
group		    	subset_H(H2O) id 2554 2555 2557 2558 2560 2561 2563 2564 2566 2567
group		    	subset_H(H2O) id 2569 2570 2572 2573 2575 2576 2578 2579 2581 2582
group		    	subset_H(H2O) id 2584 2585 2587 2588 2590 2591 2593 2594 2596 2597
group		    	subset_H(H2O) id 2599 2600 2602 2603 2605 2606 2608 2609 2611 2612
group		    	subset_H(H2O) id 2614 2615 2617 2618 2620 2621 2623 2624 2626 2627
group		    	subset_H(H2O) id 2629 2630 2632 2633 2635 2636 2638 2639 2641 2642
group		    	subset_H(H2O) id 2644 2645 2647 2648 2650 2651 2653 2654 2656 2657
group		    	subset_H(H2O) id 2659 2660 2662 2663 2665 2666 2668 2669 2671 2672
group		    	subset_H(H2O) id 2674 2675 2677 2678 2680 2681 2683 2684 2686 2687
group		    	subset_H(H2O) id 2689 2690 2692 2693 2695 2696 2698 2699 2701 2702
group		    	subset_H(H2O) id 2704 2705 2707 2708 2710 2711 2713 2714 2716 2717
group		    	subset_H(H2O) id 2719 2720 2722 2723 2725 2726 2728 2729 2731 2732
group		    	subset_H(H2O) id 2734 2735 2737 2738 2740 2741 2743 2744 2746 2747
group		    	subset_H(H2O) id 2749 2750 2752 2753 2755 2756 2758 2759 2761 2762
group		    	subset_H(H2O) id 2764 2765 2767 2768 2770 2771 2773 2774 2776 2777
group		    	subset_H(H2O) id 2779 2780 2782 2783 2785 2786 2788 2789 2791 2792
group		    	subset_H(H2O) id 2794 2795 2797 2798 2800 2801 2803 2804 2806 2807
group		    	subset_H(H2O) id 2809 2810 2812 2813 2815 2816 2818 2819 2821 2822
group		    	subset_H(H2O) id 2824 2825 2827 2828 2830 2831 2833 2834 2836 2837
group		    	subset_H(H2O) id 2839 2840 2842 2843 2845 2846 2848 2849 2851 2852
group		    	subset_H(H2O) id 2854 2855 2857 2858 2860 2861 2863 2864 2866 2867
group		    	subset_H(H2O) id 2869 2870 2872 2873 2875 2876 2878 2879 2881 2882
group		    	subset_H(H2O) id 2884 2885 2887 2888 2890 2891 2893 2894 2896 2897
group		    	subset_H(H2O) id 2899 2900 2902 2903 2905 2906 2908 2909 2911 2912
group		    	subset_H(H2O) id 2914 2915 2917 2918 2920 2921 2923 2924 2926 2927
group		    	subset_H(H2O) id 2929 2930 2932 2933 2935 2936 2938 2939 2941 2942
group		    	subset_H(H2O) id 2944 2945 2947 2948 2950 2951 2953 2954 2956 2957
group		    	subset_H(H2O) id 2959 2960 2962 2963 2965 2966 2968 2969 2971 2972
group		    	subset_H(H2O) id 2974 2975 2977 2978 2980 2981 2983 2984 2986 2987
group		    	subset_H(H2O) id 2989 2990 2992 2993 2995 2996 2998 2999 3001 3002
group		    	subset_H(H2O) id 3004 3005 3007 3008 3010 3011 3013 3014 3016 3017
group		    	subset_H(H2O) id 3019 3020 3022 3023 3025 3026 3028 3029 3031 3032
group		    	subset_H(H2O) id 3034 3035 3037 3038 3040 3041 3043 3044 3046 3047
group		    	subset_H(H2O) id 3049 3050 3052 3053 3055 3056 3058 3059 3061 3062
group		    	subset_H(H2O) id 3064 3065 3067 3068 3070 3071 3073 3074 3076 3077
group		    	subset_H(H2O) id 3079 3080 3082 3083 3085 3086 3088 3089 3091 3092
group		    	subset_H(H2O) id 3094 3095 3097 3098 3100 3101 3103 3104 3106 3107
group		    	subset_H(H2O) id 3109 3110 3112 3113 3115 3116 3118 3119 3121 3122
group		    	subset_H(H2O) id 3124 3125 3127 3128 3130 3131 3133 3134 3136 3137
group		    	subset_H(H2O) id 3139 3140 3142 3143 3145 3146 3148 3149 3151 3152
group		    	subset_H(H2O) id 3154 3155 3157 3158 3160 3161 3163 3164 3166 3167
group		    	subset_H(H2O) id 3169 3170 3172 3173 3175 3176 3178 3179 3181 3182
group		    	subset_H(H2O) id 3184 3185 3187 3188 3190 3191 3193 3194 3196 3197
group		    	subset_H(H2O) id 3199 3200 3202 3203 3205 3206 3208 3209 3211 3212
group		    	subset_H(H2O) id 3214 3215 3217 3218 3220 3221 3223 3224 3226 3227
group		    	subset_H(H2O) id 3229 3230 3232 3233 3235 3236 3238 3239 3241 3242
group		    	subset_H(H2O) id 3244 3245 3247 3248 3250 3251 3253 3254 3256 3257
group		    	subset_H(H2O) id 3259 3260 3262 3263 3265 3266 3268 3269 3271 3272
group		    	subset_H(H2O) id 3274 3275 3277 3278 3280 3281 3283 3284 3286 3287
group		    	subset_H(H2O) id 3289 3290 3292 3293 3295 3296 3298 3299 3301 3302
group		    	subset_H(H2O) id 3304 3305 3307 3308 3310 3311 3313 3314 3316 3317
group		    	subset_H(H2O) id 3319 3320 3322 3323 3325 3326 3328 3329 3331 3332
group		    	subset_H(H2O) id 3334 3335 3337 3338 3340 3341 3343 3344 3346 3347
group		    	subset_H(H2O) id 3349 3350 3352 3353 3355 3356 3358 3359 3361 3362
group		    	subset_H(H2O) id 3364 3365 3367 3368 3370 3371 3373 3374 3376 3377
group		    	subset_H(H2O) id 3379 3380 3382 3383 3385 3386 3388 3389 3391 3392
group		    	subset_H(H2O) id 3394 3395 3397 3398 3400 3401 3403 3404 3406 3407
group		    	subset_H(H2O) id 3409 3410 3412 3413 3415 3416 3418 3419 3421 3422
group		    	subset_H(H2O) id 3424 3425 3427 3428 3430 3431 3433 3434 3436 3437
group		    	subset_H(H2O) id 3439 3440 3442 3443 3445 3446 3448 3449 3451 3452
group		    	subset_H(H2O) id 3454 3455 3457 3458 3460 3461 3463 3464 3466 3467
group		    	subset_H(H2O) id 3469 3470 3472 3473 3475 3476 3478 3479 3481 3482
group		    	subset_H(H2O) id 3484 3485 3487 3488 3490 3491 3493 3494 3496 3497
group		    	subset_H(H2O) id 3499 3500 3502 3503 3505 3506 3508 3509 3511 3512
group		    	subset_H(H2O) id 3514 3515 3517 3518 3520 3521 3523 3524 3526 3527
group		    	subset_H(H2O) id 3529 3530 3532 3533 3535 3536 3538 3539 3541 3542
group		    	subset_H(H2O) id 3544 3545 3547 3548 3550 3551 3553 3554 3556 3557
group		    	subset_H(H2O) id 3559 3560 3562 3563 3565 3566 3568 3569 3571 3572
group		    	subset_H(H2O) id 3574 3575 3577 3578 3580 3581 3583 3584 3586 3587
group		    	subset_H(H2O) id 3589 3590 3592 3593 3595 3596 3598 3599 3601 3602
group		    	subset_H(H2O) id 3604 3605 3607 3608 3610 3611 3613 3614 3616 3617
group		    	subset_H(H2O) id 3619 3620 3622 3623 3625 3626 3628 3629 3631 3632
group		    	subset_H(H2O) id 3634 3635 3637 3638 3640 3641 3643 3644 3646 3647
group		    	subset_H(H2O) id 3649 3650 3652 3653 3655 3656 3658 3659 3661 3662
group		    	subset_H(H2O) id 3664 3665 3667 3668 3670 3671 3673 3674 3676 3677
group		    	subset_H(H2O) id 3679 3680 3682 3683 3685 3686 3688 3689 3691 3692
group		    	subset_H(H2O) id 3694 3695 3697 3698 3700 3701 3703 3704 3706 3707
group		    	subset_H(H2O) id 3709 3710 3712 3713 3715 3716 3718 3719 3721 3722
group		    	subset_H(H2O) id 3724 3725 3727 3728 3730 3731 3733 3734 3736 3737
group		    	subset_H(H2O) id 3739 3740 3742 3743 3745 3746 3748 3749 3751 3752
group		    	subset_H(H2O) id 3754 3755 3757 3758 3760 3761 3763 3764 3766 3767
group		    	subset_H(H2O) id 3769 3770 3772 3773 3775 3776 3778 3779 3781 3782
group		    	subset_H(H2O) id 3784 3785 3787 3788 3790 3791 3793 3794 3796 3797
group		    	subset_H(H2O) id 3799 3800 3802 3803 3805 3806 3808 3809 3811 3812
group		    	subset_H(H2O) id 3814 3815 3817 3818 3820 3821 3823 3824 3826 3827
group		    	subset_H(H2O) id 3829 3830 3832 3833 3835 3836 3838 3839 3841 3842
group		    	subset_H(H2O) id 3844 3845 3847 3848 3850 3851 3853 3854 3856 3857
group		    	subset_H(H2O) id 3859 3860 3862 3863 3865 3866 3868 3869 3871 3872
group		    	subset_H(H2O) id 3874 3875 3877 3878 3880 3881 3883 3884 3886 3887
group		    	subset_H(H2O) id 3889 3890 3892 3893 3895 3896 3898 3899 3901 3902
group		    	subset_H(H2O) id 3904 3905 3907 3908 3910 3911 3913 3914 3916 3917
group		    	subset_H(H2O) id 3919 3920 3922 3923 3925 3926 3928 3929 3931 3932
group		    	subset_H(H2O) id 3934 3935 3937 3938 3940 3941 3943 3944 3946 3947
group		    	subset_H(H2O) id 3949 3950 3952 3953 3955 3956 3958 3959 3961 3962
group		    	subset_H(H2O) id 3964 3965 3967 3968 3970 3971 3973 3974 3976 3977
group		    	subset_H(H2O) id 3979 3980 3982 3983 3985 3986 3988 3989 3991 3992
group		    	subset_H(H2O) id 3994 3995 3997 3998 4000 4001 4003 4004 4006 4007
group		    	subset_H(H2O) id 4009 4010 4012 4013 4015 4016 4018 4019 4021 4022
group		    	subset_H(H2O) id 4024 4025 4027 4028 4030 4031 4033 4034 4036 4037
group		    	subset_H(H2O) id 4039 4040 4042 4043 4045 4046 4048 4049 4051 4052
group		    	subset_H(H2O) id 4054 4055 4057 4058 4060 4061 4063 4064 4066 4067
group		    	subset_H(H2O) id 4069 4070 4072 4073 4075 4076 4078 4079 4081 4082
group		    	subset_H(H2O) id 4084 4085 4087 4088 4090 4091 4093 4094 4096 4097
group		    	subset_H(H2O) id 4099 4100 4102 4103 4105 4106 4108 4109 4111 4112
group		    	subset_H(H2O) id 4114 4115 4117 4118 4120 4121 4123 4124 4126 4127
group		    	subset_H(H2O) id 4129 4130 4132 4133 4135 4136 4138 4139 4141 4142
group		    	subset_H(H2O) id 4144 4145 4147 4148 4150 4151 4153 4154 4156 4157
group		    	subset_H(H2O) id 4159 4160 4162 4163 4165 4166 4168 4169 4171 4172
group		    	subset_H(H2O) id 4174 4175 4177 4178 4180 4181 4183 4184 4186 4187
group		    	subset_H(H2O) id 4189 4190 4192 4193 4195 4196 4198 4199 4201 4202
group		    	subset_H(H2O) id 4204 4205 4207 4208 4210 4211 4213 4214 4216 4217
group		    	subset_H(H2O) id 4219 4220 4222 4223 4225 4226 4228 4229 4231 4232
group		    	subset_H(H2O) id 4234 4235 4237 4238 4240 4241 4243 4244 4246 4247
group		    	subset_H(H2O) id 4249 4250 4252 4253 4255 4256 4258 4259 4261 4262
group		    	subset_H(H2O) id 4264 4265 4267 4268 4270 4271 4273 4274 4276 4277
group		    	subset_H(H2O) id 4279 4280 4282 4283 4285 4286 4288 4289 4291 4292
group		    	subset_H(H2O) id 4294 4295 4297 4298 4300 4301 4303 4304 4306 4307
group		    	subset_H(H2O) id 4309 4310 4312 4313 4315 4316 4318 4319 4321 4322
group		    	subset_H(H2O) id 4324 4325 4327 4328 4330 4331 4333 4334 4336 4337
group		    	subset_H(H2O) id 4339 4340 4342 4343 4345 4346 4348 4349 4351 4352
group		    	subset_H(H2O) id 4354 4355 4357 4358 4360 4361 4363 4364 4366 4367
group		    	subset_H(H2O) id 4369 4370 4372 4373 4375 4376 4378 4379 4381 4382
group		    	subset_H(H2O) id 4384 4385 4387 4388 4390 4391 4393 4394 4396 4397
group		    	subset_H(H2O) id 4399 4400 4402 4403 4405 4406 4408 4409 4411 4412
group		    	subset_H(H2O) id 4414 4415 4417 4418 4420 4421 4423 4424 4426 4427
group		    	subset_H(H2O) id 4429 4430 4432 4433 4435 4436 4438 4439 4441 4442
group		    	subset_H(H2O) id 4444 4445 4447 4448 4450 4451 4453 4454 4456 4457
group		    	subset_H(H2O) id 4459 4460 4462 4463 4465 4466 4468 4469 4471 4472
group		    	subset_H(H2O) id 4474 4475 4477 4478 4480 4481 4483 4484 4486 4487
group		    	subset_H(H2O) id 4489 4490 4492 4493 4495 4496 4498 4499 4501 4502
group		    	subset_H(H2O) id 4504 4505 4507 4508 4510 4511 4513 4514 4516 4517
group		    	subset_H(H2O) id 4519 4520 4522 4523 4525 4526 4528 4529 4531 4532
group		    	subset_H(H2O) id 4534 4535 4537 4538 4540 4541 4543 4544 4546 4547
group		    	subset_H(H2O) id 4549 4550 4552 4553 4555 4556 4558 4559 4561 4562
group		    	subset_H(H2O) id 4564 4565 4567 4568 4570 4571 4573 4574 4576 4577
group		    	subset_H(H2O) id 4579 4580 4582 4583 4585 4586 4588 4589 4591 4592
group		    	subset_H(H2O) id 4594 4595 4597 4598 4600 4601 4603 4604 4606 4607
group		    	subset_H(H2O) id 4609 4610 4612 4613 4615 4616 4618 4619 4621 4622
group		    	subset_H(H2O) id 4624 4625 4627 4628 4630 4631 4633 4634 4636 4637
group		    	subset_H(H2O) id 4639 4640 4642 4643 4645 4646 4648 4649 4651 4652
group		    	subset_H(H2O) id 4654 4655 4657 4658 4660 4661 4663 4664 4666 4667
group		    	subset_H(H2O) id 4669 4670 4672 4673 4675 4676 4678 4679 4681 4682
group		    	subset_H(H2O) id 4684 4685 4687 4688 4690 4691 4693 4694 4696 4697
group		    	subset_H(H2O) id 4699 4700 4702 4703 4705 4706 4708 4709 4711 4712
group		    	subset_H(H2O) id 4714 4715 4717 4718 4720 4721 4723 4724 4726 4727
group		    	subset_H(H2O) id 4729 4730 4732 4733 4735 4736 4738 4739 4741 4742
group		    	subset_H(H2O) id 4744 4745 4747 4748 4750 4751 4753 4754 4756 4757
group		    	subset_H(H2O) id 4759 4760 4762 4763 4765 4766 4768 4769 4771 4772
group		    	subset_H(H2O) id 4774 4775 4777 4778 4780 4781 4783 4784 4786 4787
group		    	subset_H(H2O) id 4789 4790 4792 4793 4795 4796 4798 4799 4801 4802
group		    	subset_H(H2O) id 4804 4805 4807 4808 4810 4811 4813 4814 4816 4817
group		    	subset_H(H2O) id 4819 4820 4822 4823 4825 4826 4828 4829 4831 4832
group		    	subset_H(H2O) id 4834 4835 4837 4838 4840 4841 4843 4844 4846 4847
group		    	subset_H(H2O) id 4849 4850 4852 4853 4855 4856 4858 4859 4861 4862
group		    	subset_H(H2O) id 4864 4865 4867 4868 4870 4871 4873 4874 4876 4877
group		    	subset_H(H2O) id 4879 4880 4882 4883 4885 4886 4888 4889 4891 4892
group		    	subset_H(H2O) id 4894 4895 4897 4898 4900 4901 4903 4904 4906 4907
group		    	subset_H(H2O) id 4909 4910 4912 4913 4915 4916 4918 4919 4921 4922
group		    	subset_H(H2O) id 4924 4925 4927 4928 4930 4931 4933 4934 4936 4937
group		    	subset_H(H2O) id 4939 4940 4942 4943 4945 4946 4948 4949 4951 4952
group		    	subset_H(H2O) id 4954 4955 4957 4958 4960 4961 4963 4964 4966 4967
group		    	subset_H(H2O) id 4969 4970 4972 4973 4975 4976 4978 4979 4981 4982
group		    	subset_H(H2O) id 4984 4985 4987 4988 4990 4991 4993 4994 4996 4997
group		    	subset_H(H2O) id 4999 5000 5002 5003 5005 5006 5008 5009 5011 5012
group		    	subset_H(H2O) id 5014 5015 5017 5018 5020 5021 5023 5024 5026 5027
group		    	subset_H(H2O) id 5029 5030 5032 5033 5035 5036 5038 5039 5041 5042
group		    	subset_H(H2O) id 5044 5045 5047 5048 5050 5051 5053 5054 5056 5057
group		    	subset_H(H2O) id 5059 5060 5062 5063 5065 5066 5068 5069 5071 5072
group		    	subset_H(H2O) id 5074 5075 5077 5078 5080 5081 5083 5084 5086 5087
group		    	subset_H(H2O) id 5089 5090 5092 5093 5095 5096 5098 5099 5101 5102
group		    	subset_H(H2O) id 5104 5105 5107 5108 5110 5111 5113 5114 5116 5117
group		    	subset_H(H2O) id 5119 5120 5122 5123 5125 5126 5128 5129 5131 5132
group		    	subset_H(H2O) id 5134 5135 5137 5138 5140 5141 5143 5144 5146 5147
group		    	subset_H(H2O) id 5149 5150 5152 5153 5155 5156 5158 5159 5161 5162
group		    	subset_H(H2O) id 5164 5165 5167 5168 5170 5171 5173 5174 5176 5177
group		    	subset_H(H2O) id 5179 5180 5182 5183 5185 5186 5188 5189 5191 5192
group		    	subset_H(H2O) id 5194 5195 5197 5198 5200 5201 5203 5204 5206 5207
group		    	subset_H(H2O) id 5209 5210 5212 5213 5215 5216 5218 5219 5221 5222
group		    	subset_H(H2O) id 5224 5225 5227 5228 5230 5231 5233 5234 5236 5237
group		    	subset_H(H2O) id 5239 5240 5242 5243 5245 5246 5248 5249 5251 5252
group		    	subset_H(H2O) id 5254 5255 5257 5258 5260 5261 5263 5264 5266 5267
group		    	subset_H(H2O) id 5269 5270 5272 5273 5275 5276 5278 5279 5281 5282
group		    	subset_H(H2O) id 5284 5285 5287 5288 5290 5291 5293 5294 5296 5297
group		    	subset_H(H2O) id 5299 5300 5302 5303 5305 5306 5308 5309 5311 5312
group		    	subset_H(H2O) id 5314 5315 5317 5318 5320 5321 5323 5324 5326 5327
group		    	subset_H(H2O) id 5329 5330 5332 5333 5335 5336 5338 5339 5341 5342
group		    	subset_H(H2O) id 5344 5345 5347 5348 5350 5351 5353 5354 5356 5357
group		    	subset_H(H2O) id 5359 5360 5362 5363 5365 5366 5368 5369 5371 5372
group		    	subset_H(H2O) id 5374 5375 5377 5378 5380 5381 5383 5384 5386 5387
group		    	subset_H(H2O) id 5389 5390 5392 5393 5395 5396 5398 5399 5401 5402
group		    	subset_H(H2O) id 5404 5405 5407 5408 5410 5411 5413 5414 5416 5417
group		    	subset_H(H2O) id 5419 5420 5422 5423 5425 5426 5428 5429 5431 5432
group		    	subset_H(H2O) id 5434 5435 5437 5438 5440 5441 5443 5444 5446 5447
group		    	subset_H(H2O) id 5449 5450 5452 5453 5455 5456 5458 5459 5461 5462
group		    	subset_H(H2O) id 5464 5465 5467 5468 5470 5471 5473 5474 5476 5477
group		    	subset_H(H2O) id 5479 5480 5482 5483 5485 5486 5488 5489 5491 5492
group		    	subset_H(H2O) id 5494 5495 5497 5498 5500 5501 5503 5504 5506 5507
group		    	subset_H(H2O) id 5509 5510 5512 5513 5515 5516 5518 5519 5521 5522
group		    	subset_H(H2O) id 5524 5525 5527 5528 5530 5531 5533 5534 5536 5537
group		    	subset_H(H2O) id 5539 5540 5542 5543 5545 5546 5548 5549 5551 5552
group		    	subset_H(H2O) id 5554 5555 5557 5558 5560 5561 5563 5564 5566 5567
group		    	subset_H(H2O) id 5569 5570 5572 5573 5575 5576 5578 5579 5581 5582
group		    	subset_H(H2O) id 5584 5585 5587 5588 5590 5591 5593 5594 5596 5597
group		    	subset_H(H2O) id 5599 5600 5602 5603 5605 5606 5608 5609 5611 5612
group		    	subset_H(H2O) id 5614 5615 5617 5618 5620 5621 5623 5624 5626 5627
group		    	subset_H(H2O) id 5629 5630 5632 5633 5635 5636 5638 5639 5641 5642
group		    	subset_H(H2O) id 5644 5645 5647 5648 5650 5651 5653 5654 5656 5657
group		    	subset_H(H2O) id 5659 5660 5662 5663 5665 5666 5668 5669 5671 5672
group		    	subset_H(H2O) id 5674 5675 5677 5678 5680 5681 5683 5684 5686 5687
group		    	subset_H(H2O) id 5689 5690 5692 5693 5695 5696 5698 5699 5701 5702
group		    	subset_H(H2O) id 5704 5705 5707 5708 5710 5711 5713 5714 5716 5717
group		    	subset_H(H2O) id 5719 5720 5722 5723 5725 5726 5728 5729 5731 5732
group		    	subset_H(H2O) id 5734 5735 5737 5738 5740 5741 5743 5744 5746 5747
group		    	subset_H(H2O) id 5749 5750 5752 5753 5755 5756 5758 5759 5761 5762
group		    	subset_H(H2O) id 5764 5765 5767 5768 5770 5771 5773 5774 5776 5777
group		    	subset_H(H2O) id 5779 5780 5782 5783 5785 5786 5788 5789 5791 5792
group		    	subset_H(H2O) id 5794 5795 5797 5798 5800 5801 5803 5804 5806 5807
group		    	subset_H(H2O) id 5809 5810 5812 5813 5815 5816 5818 5819 5821 5822
group		    	subset_H(H2O) id 5824 5825 5827 5828 5830 5831 5833 5834 5836 5837
group		    	subset_H(H2O) id 5839 5840 5842 5843 5845 5846 5848 5849 5851 5852
group		    	subset_H(H2O) id 5854 5855 5857 5858 5860 5861 5863 5864 5866 5867
group		    	subset_H(H2O) id 5869 5870 5872 5873 5875 5876 5878 5879 5881 5882
group		    	subset_H(H2O) id 5884 5885 5887 5888 5890 5891 5893 5894 5896 5897
group		    	subset_H(H2O) id 5899 5900 5902 5903 5905 5906 5908 5909 5911 5912
group		    	subset_H(H2O) id 5914 5915 5917 5918 5920 5921 5923 5924 5926 5927
group		    	subset_H(H2O) id 5929 5930 5932 5933 5935 5936 5938 5939 5941 5942
group		    	subset_H(H2O) id 5944 5945 5947 5948 5950 5951 5953 5954 5956 5957
group		    	subset_H(H2O) id 5959 5960 5962 5963 5965 5966 5968 5969 5971 5972
group		    	subset_H(H2O) id 5974 5975 5977 5978 5980 5981 5983 5984 5986 5987
group		    	subset_H(H2O) id 5989 5990 5992 5993 5995 5996 5998 5999 6001 6002
group		    	subset_H(H2O) id 6004 6005 6007 6008 6010 6011 6013 6014 6016 6017
group		    	subset_H(H2O) id 6019 6020 6022 6023 6025 6026 6028 6029 6031 6032
group		    	subset_H(H2O) id 6034 6035 6037 6038 6040 6041 6043 6044 6046 6047
group		    	subset_H(H2O) id 6049 6050 6052 6053 6055 6056 6058 6059 6061 6062
group		    	subset_H(H2O) id 6064 6065 6067 6068 6070 6071 6073 6074 6076 6077
group		    	subset_H(H2O) id 6079 6080 6082 6083 6085 6086 6088 6089 6091 6092
group		    	subset_H(H2O) id 6094 6095 6097 6098 6100 6101 6103 6104 6106 6107
group		    	subset_H(H2O) id 6109 6110 6112 6113 6115 6116 6118 6119 6121 6122
group		    	subset_H(H2O) id 6124 6125 6127 6128 6130 6131 6133 6134 6136 6137
group		    	subset_H(H2O) id 6139 6140 6142 6143 6145 6146 6148 6149 6151 6152
group		    	subset_H(H2O) id 6154 6155 6157 6158 6160 6161 6163 6164 6166 6167
group		    	subset_H(H2O) id 6169 6170 6172 6173 6175 6176 6178 6179 6181 6182
group		    	subset_H(H2O) id 6184 6185 6187 6188 6190 6191 6193 6194 6196 6197
group		    	subset_H(H2O) id 6199 6200 6202 6203 6205 6206 6208 6209 6211 6212
group		    	subset_H(H2O) id 6214 6215 6217 6218 6220 6221 6223 6224 6226 6227
group		    	subset_H(H2O) id 6229 6230 6232 6233 6235 6236 6238 6239 6241 6242
group		    	subset_H(H2O) id 6244 6245 6247 6248 6250 6251 6253 6254 6256 6257
group		    	subset_H(H2O) id 6259 6260 6262 6263 6265 6266 6268 6269 6271 6272
group		    	subset_H(H2O) id 6274 6275 6277 6278 6280 6281 6283 6284 6286 6287
group		    	subset_H(H2O) id 6289 6290 6292 6293 6295 6296 6298 6299 6301 6302
group		    	subset_H(H2O) id 6304 6305 6307 6308 6310 6311 6313 6314 6316 6317
group		    	subset_H(H2O) id 6319 6320 6322 6323 6325 6326 6328 6329 6331 6332
group		    	subset_H(H2O) id 6334 6335 6337 6338 6340 6341 6343 6344 6346 6347
group		    	subset_H(H2O) id 6349 6350 6352 6353 6355 6356 6358 6359 6361 6362
group		    	subset_H(H2O) id 6364 6365 6367 6368 6370 6371 6373 6374 6376 6377
group		    	subset_H(H2O) id 6379 6380 6382 6383 6385 6386 6388 6389 6391 6392
group		    	subset_H(H2O) id 6394 6395 6397 6398 6400 6401 6403 6404 6406 6407
group		    	subset_H(H2O) id 6409 6410 6412 6413 6415 6416 6418 6419 6421 6422
group		    	subset_H(H2O) id 6424 6425 6427 6428 6430 6431 6433 6434 6436 6437
group		    	subset_H(H2O) id 6439 6440 6442 6443 6445 6446 6448 6449 6451 6452
group		    	subset_H(H2O) id 6454 6455 6457 6458 6460 6461 6463 6464 6466 6467
group		    	subset_H(H2O) id 6469 6470 6472 6473 6475 6476 6478 6479 6481 6482
group		    	subset_H(H2O) id 6484 6485 6487 6488 6490 6491 6493 6494 6496 6497
group		    	subset_H(H2O) id 6499 6500 6502 6503 6505 6506 6508 6509 6511 6512
group		    	subset_H(H2O) id 6514 6515 6517 6518 6520 6521 6523 6524 6526 6527
group		    	subset_H(H2O) id 6529 6530 6532 6533 6535 6536 6538 6539 6541 6542
group		    	subset_H(H2O) id 6544 6545 6547 6548 6550 6551 6553 6554 6556 6557
group		    	subset_H(H2O) id 6559 6560 6562 6563 6565 6566 6568 6569 6571 6572
group		    	subset_H(H2O) id 6574 6575 6577 6578 6580 6581 6583 6584 6586 6587
group		    	subset_H(H2O) id 6589 6590 6592 6593 6595 6596 6598 6599 6601 6602
group		    	subset_H(H2O) id 6604 6605 6607 6608 6610 6611 6613 6614 6616 6617
group		    	subset_H(H2O) id 6619 6620 6622 6623 6625 6626 6628 6629 6631 6632
group		    	subset_H(H2O) id 6634 6635 6637 6638 6640 6641 6643 6644 6646 6647
group		    	subset_H(H2O) id 6649 6650 6652 6653 6655 6656 6658 6659 6661 6662
group		    	subset_H(H2O) id 6664 6665 6667 6668 6670 6671 6673 6674 6676 6677
group		    	subset_H(H2O) id 6679 6680 6682 6683 6685 6686 6688 6689 6691 6692
group		    	subset_H(H2O) id 6694 6695 6697 6698 6700 6701 6703 6704 6706 6707
group		    	subset_H(H2O) id 6709 6710 6712 6713 6715 6716 6718 6719 6721 6722
group		    	subset_H(H2O) id 6724 6725 6727 6728 6730 6731 6733 6734 6736 6737
group		    	subset_H(H2O) id 6739 6740 6742 6743 6745 6746 6748 6749 6751 6752
group		    	subset_H(H2O) id 6754 6755 6757 6758 6760 6761 6763 6764 6766 6767
group		    	subset_H(H2O) id 6769 6770 6772 6773 6775 6776 6778 6779 6781 6782
group		    	subset_H(H2O) id 6784 6785 6787 6788 6790 6791 6793 6794 6796 6797
group		    	subset_H(H2O) id 6799 6800 6802 6803 6805 6806 6808 6809 6811 6812
group		    	subset_H(H2O) id 6814 6815 6817 6818 6820 6821 6823 6824 6826 6827
group		    	subset_H(H2O) id 6829 6830 6832 6833 6835 6836 6838 6839 6841 6842
group		    	subset_H(H2O) id 6844 6845 6847 6848 6850 6851 6853 6854 6856 6857
group		    	subset_H(H2O) id 6859 6860 6862 6863 6865 6866 6868 6869 6871 6872
group		    	subset_H(H2O) id 6874 6875 6877 6878 6880 6881 6883 6884 6886 6887
group		    	subset_H(H2O) id 6889 6890 6892 6893 6895 6896 6898 6899 6901 6902
group		    	subset_H(H2O) id 6904 6905 6907 6908 6910 6911 6913 6914 6916 6917
group		    	subset_H(H2O) id 6919 6920 6922 6923 6925 6926 6928 6929 6931 6932
group		    	subset_H(H2O) id 6934 6935 6937 6938 6940 6941 6943 6944 6946 6947
group		    	subset_H(H2O) id 6949 6950 6952 6953 6955 6956 6958 6959 6961 6962
group		    	subset_H(H2O) id 6964 6965 6967 6968 6970 6971 6973 6974 6976 6977
group		    	subset_H(H2O) id 6979 6980 6982 6983 6985 6986 6988 6989 6991 6992
group		    	subset_H(H2O) id 6994 6995 6997 6998 7000 7001 7003 7004 7006 7007
group		    	subset_H(H2O) id 7009 7010 7012 7013 7015 7016 7018 7019 7021 7022
group		    	subset_H(H2O) id 7024 7025 7027 7028 7030 7031 7033 7034 7036 7037
group		    	subset_H(H2O) id 7039 7040 7042 7043 7045 7046 7048 7049 7051 7052
group		    	subset_H(H2O) id 7054 7055 7057 7058 7060 7061 7063 7064 7066 7067
group		    	subset_H(H2O) id 7069 7070 7072 7073 7075 7076 7078 7079 7081 7082
group		    	subset_H(H2O) id 7084 7085 7087 7088 7090 7091 7093 7094 7096 7097
group		    	subset_H(H2O) id 7099 7100 7102 7103 7105 7106 7108 7109 7111 7112
group		    	subset_H(H2O) id 7114 7115 7117 7118 7120 7121 7123 7124 7126 7127
group		    	subset_H(H2O) id 7129 7130 7132 7133 7135 7136 7138 7139 7141 7142
group		    	subset_H(H2O) id 7144 7145 7147 7148 7150 7151 7153 7154 7156 7157
group		    	subset_H(H2O) id 7159 7160 7162 7163 7165 7166 7168 7169 7171 7172
group		    	subset_H(H2O) id 7174 7175 7177 7178 7180 7181 7183 7184 7186 7187
group		    	subset_H(H2O) id 7189 7190 7192 7193 7195 7196 7198 7199 7201 7202
group		    	subset_H(H2O) id 7204 7205 7207 7208 7210 7211 7213 7214 7216 7217
group		    	subset_H(H2O) id 7219 7220 7222 7223 7225 7226 7228 7229 7231 7232
group		    	subset_H(H2O) id 7234 7235 7237 7238 7240 7241 7243 7244 7246 7247
group		    	subset_H(H2O) id 7249 7250 7252 7253 7255 7256 7258 7259 7261 7262
group		    	subset_H(H2O) id 7264 7265 7267 7268 7270 7271 7273 7274 7276 7277
group		    	subset_H(H2O) id 7279 7280 7282 7283 7285 7286 7288 7289 7291 7292
group		    	subset_H(H2O) id 7294 7295 7297 7298 7300 7301 7303 7304 7306 7307
group		    	subset_H(H2O) id 7309 7310 7312 7313 7315 7316 7318 7319 7321 7322
group		    	subset_H(H2O) id 7324 7325 7327 7328 7330 7331 7333 7334 7336 7337
group		    	subset_H(H2O) id 7339 7340 7342 7343 7345 7346 7348 7349 7351 7352
group		    	subset_H(H2O) id 7354 7355 7357 7358 7360 7361 7363 7364 7366 7367
group		    	subset_H(H2O) id 7369 7370 7372 7373 7375 7376 7378 7379 7381 7382
group		    	subset_H(H2O) id 7384 7385 7387 7388 7390 7391 7393 7394 7396 7397
group		    	subset_H(H2O) id 7399 7400 7402 7403 7405 7406 7408 7409 7411 7412
group		    	subset_H(H2O) id 7414 7415 7417 7418 7420 7421 7423 7424 7426 7427
group		    	subset_H(H2O) id 7429 7430 7432 7433 7435 7436 7438 7439 7441 7442
group		    	subset_H(H2O) id 7444 7445 7447 7448 7450 7451 7453 7454 7456 7457
group		    	subset_H(H2O) id 7459 7460 7462 7463 7465 7466 7468 7469 7471 7472
group		    	subset_H(H2O) id 7474 7475 7477 7478 7480 7481 7483 7484 7486 7487
group		    	subset_H(H2O) id 7489 7490 7492 7493 7495 7496 7498 7499 7501 7502
group		    	subset_H(H2O) id 7504 7505 7507 7508 7510 7511 7513 7514 7516 7517
group		    	subset_H(H2O) id 7519 7520 7522 7523 7525 7526 7528 7529 7531 7532
group		    	subset_H(H2O) id 7534 7535 7537 7538 7540 7541 7543 7544 7546 7547
group		    	subset_H(H2O) id 7549 7550 7552 7553 7555 7556 7558 7559 7561 7562
group		    	subset_H(H2O) id 7564 7565 7567 7568 7570 7571 7573 7574 7576 7577
group		    	subset_H(H2O) id 7579 7580 7582 7583 7585 7586 7588 7589 7591 7592
group		    	subset_H(H2O) id 7594 7595 7597 7598 7600 7601 7603 7604 7606 7607
group		    	subset_H(H2O) id 7609 7610 7612 7613 7615 7616 7618 7619 7621 7622
group		    	subset_H(H2O) id 7624 7625 7627 7628 7630 7631 7633 7634 7636 7637
group		    	subset_H(H2O) id 7639 7640 7642 7643 7645 7646 7648 7649 7651 7652
group		    	subset_H(H2O) id 7654 7655 7657 7658 7660 7661 7663 7664 7666 7667
group		    	subset_H(H2O) id 7669 7670 7672 7673 7675 7676 7678 7679 7681 7682
group		    	subset_H(H2O) id 7684 7685 7687 7688 7690 7691 7693 7694 7696 7697
group		    	subset_H(H2O) id 7699 7700 7702 7703 7705 7706 7708 7709 7711 7712
group		    	subset_H(H2O) id 7714 7715 7717 7718 7720 7721 7723 7724 7726 7727
group		    	subset_H(H2O) id 7729 7730 7732 7733 7735 7736 7738 7739 7741 7742
group		    	subset_H(H2O) id 7744 7745 7747 7748 7750 7751 7753 7754 7756 7757
group		    	subset_H(H2O) id 7759 7760 7762 7763 7765 7766 7768 7769 7771 7772
group		    	subset_H(H2O) id 7774 7775 7777 7778 7780 7781 7783 7784 7786 7787
group		    	subset_H(H2O) id 7789 7790 7792 7793 7795 7796 7798 7799 7801 7802
group		    	subset_H(H2O) id 7804 7805 7807 7808 7810 7811 7813 7814 7816 7817
group		    	subset_H(H2O) id 7819 7820 7822 7823 7825 7826 7828 7829 7831 7832
group		    	subset_H(H2O) id 7834 7835 7837 7838 7840 7841 7843 7844 7846 7847
group		    	subset_H(H2O) id 7849 7850 7852 7853 7855 7856 7858 7859 7861 7862
group		    	subset_H(H2O) id 7864 7865 7867 7868 7870 7871 7873 7874 7876 7877
group		    	subset_H(H2O) id 7879 7880 7882 7883 7885 7886 7888 7889 7891 7892
group		    	subset_H(H2O) id 7894 7895 7897 7898 7900 7901 7903 7904 7906 7907
group		    	subset_H(H2O) id 7909 7910 7912 7913 7915 7916 7918 7919 7921 7922
group		    	subset_H(H2O) id 7924 7925 7927 7928 7930 7931 7933 7934 7936 7937
group		    	subset_H(H2O) id 7939 7940 7942 7943 7945 7946 7948 7949 7951 7952
group		    	subset_H(H2O) id 7954 7955 7957 7958 7960 7961 7963 7964 7966 7967
group		    	subset_H(H2O) id 7969 7970 7972 7973 7975 7976 7978 7979 7981 7982
group		    	subset_H(H2O) id 7984 7985 7987 7988 7990 7991 7993 7994 7996 7997
group		    	subset_H(H2O) id 7999 8000 8002 8003 8005 8006 8008 8009 8011 8012
group		    	subset_H(H2O) id 8014 8015 8017 8018 8020 8021 8023 8024 8026 8027
group		    	subset_H(H2O) id 8029 8030 8032 8033 8035 8036 8038 8039 8041 8042
group		    	subset_H(H2O) id 8044 8045 8047 8048 8050 8051 8053 8054 8056 8057
group		    	subset_H(H2O) id 8059 8060 8062 8063 8065 8066 8068 8069 8071 8072
group		    	subset_H(H2O) id 8074 8075 8077 8078 8080 8081 8083 8084 8086 8087
group		    	subset_H(H2O) id 8089 8090 8092 8093 8095 8096 8098 8099 8101 8102
group		    	subset_H(H2O) id 8104 8105 8107 8108 8110 8111 8113 8114 8116 8117
group		    	subset_H(H2O) id 8119 8120 8122 8123 8125 8126 8128 8129 8131 8132
group		    	subset_H(H2O) id 8134 8135 8137 8138 8140 8141 8143 8144 8146 8147
group		    	subset_H(H2O) id 8149 8150 8152 8153 8155 8156 8158 8159 8161 8162
group		    	subset_H(H2O) id 8164 8165 8167 8168 8170 8171 8173 8174 8176 8177
group		    	subset_H(H2O) id 8179 8180 8182 8183 8185 8186 8188 8189 8191 8192
group		    	subset_H(H2O) id 8194 8195 8197 8198 8200 8201 8203 8204 8206 8207
group		    	subset_H(H2O) id 8209 8210 8212 8213 8215 8216 8218 8219 8221 8222
group		    	subset_H(H2O) id 8224 8225 8227 8228 8230 8231 8233 8234 8236 8237
group		    	subset_H(H2O) id 8239 8240 8242 8243 8245 8246 8248 8249 8251 8252
group		    	subset_H(H2O) id 8254 8255 8257 8258 8260 8261 8263 8264 8266 8267
group		    	subset_H(H2O) id 8269 8270 8272 8273 8275 8276 8278 8279 8281 8282
group		    	subset_H(H2O) id 8284 8285 8287 8288 8290 8291 8293 8294 8296 8297
group		    	subset_H(H2O) id 8299 8300 8302 8303 8305 8306 8308 8309 8311 8312
group		    	subset_H(H2O) id 8314 8315 8317 8318 8320 8321 8323 8324 8326 8327
group		    	subset_H(H2O) id 8329 8330 8332 8333 8335 8336 8338 8339 8341 8342
group		    	subset_H(H2O) id 8344 8345 8347 8348 8350 8351 8353 8354 8356 8357
group		    	subset_H(H2O) id 8359 8360 8362 8363 8365 8366 8368 8369 8371 8372
group		    	subset_H(H2O) id 8374 8375 8377 8378 8380 8381 8383 8384 8386 8387
group		    	subset_H(H2O) id 8389 8390 8392 8393 8395 8396 8398 8399 8401 8402
group		    	subset_H(H2O) id 8404 8405 8407 8408 8410 8411 8413 8414 8416 8417
group		    	subset_H(H2O) id 8419 8420 8422 8423 8425 8426 8428 8429 8431 8432
group		    	subset_H(H2O) id 8434 8435 8437 8438 8440 8441 8443 8444 8446 8447
group		    	subset_H(H2O) id 8449 8450 8452 8453 8455 8456 8458 8459 8461 8462
group		    	subset_H(H2O) id 8464 8465 8467 8468 8470 8471 8473 8474 8476 8477
group		    	subset_H(H2O) id 8479 8480 8482 8483 8485 8486 8488 8489 8491 8492
group		    	subset_H(H2O) id 8494 8495 8497 8498 8500 8501 8503 8504 8506 8507
group		    	subset_H(H2O) id 8509 8510 8512 8513 8515 8516 8518 8519 8521 8522
group		    	subset_H(H2O) id 8524 8525 8527 8528 8530 8531 8533 8534 8536 8537
group		    	subset_H(H2O) id 8539 8540 8542 8543 8545 8546 8548 8549 8551 8552
group		    	subset_H(H2O) id 8554 8555 8557 8558 8560 8561 8563 8564 8566 8567
group		    	subset_H(H2O) id 8569 8570 8572 8573 8575 8576 8578 8579 8581 8582
group		    	subset_H(H2O) id 8584 8585 8587 8588 8590 8591 8593 8594 8596 8597
group		    	subset_H(H2O) id 8599 8600 8602 8603 8605 8606 8608 8609 8611 8612
group		    	subset_H(H2O) id 8614 8615 8617 8618 8620 8621 8623 8624 8626 8627
group		    	subset_H(H2O) id 8629 8630 8632 8633 8635 8636 8638 8639 8641 8642
group		    	subset_H(H2O) id 8644 8645 8647 8648 8650 8651 8653 8654 8656 8657
group		    	subset_H(H2O) id 8659 8660 8662 8663 8665 8666 8668 8669 8671 8672
group		    	subset_H(H2O) id 8674 8675 8677 8678 8680 8681 8683 8684 8686 8687
group		    	subset_H(H2O) id 8689 8690 8692 8693 8695 8696 8698 8699 8701 8702
group		    	subset_H(H2O) id 8704 8705 8707 8708 8710 8711 8713 8714 8716 8717
group		    	subset_H(H2O) id 8719 8720 8722 8723 8725 8726 8728 8729 8731 8732
group		    	subset_H(H2O) id 8734 8735 8737 8738 8740 8741 8743 8744 8746 8747
group		    	subset_H(H2O) id 8749 8750 8752 8753 8755 8756 8758 8759 8761 8762
group		    	subset_H(H2O) id 8764 8765 8767 8768 8770 8771 8773 8774 8776 8777
group		    	subset_H(H2O) id 8779 8780 8782 8783 8785 8786 8788 8789 8791 8792
group		    	subset_H(H2O) id 8794 8795 8797 8798 8800 8801 8803 8804 8806 8807
group		    	subset_H(H2O) id 8809 8810 8812 8813 8815 8816 8818 8819 8821 8822
group		    	subset_H(H2O) id 8824 8825 8827 8828 8830 8831 8833 8834 8836 8837
group		    	subset_H(H2O) id 8839 8840 8842 8843 8845 8846 8848 8849 8851 8852
group		    	subset_H(H2O) id 8854 8855 8857 8858 8860 8861 8863 8864 8866 8867
group		    	subset_H(H2O) id 8869 8870 8872 8873 8875 8876 8878 8879 8881 8882
group		    	subset_H(H2O) id 8884 8885 8887 8888 8890 8891 8893 8894 8896 8897
group		    	subset_H(H2O) id 8899 8900 8902 8903 8905 8906 8908 8909 8911 8912
group		    	subset_H(H2O) id 8914 8915 8917 8918 8920 8921 8923 8924 8926 8927
group		    	subset_H(H2O) id 8929 8930 8932 8933 8935 8936 8938 8939 8941 8942
group		    	subset_H(H2O) id 8944 8945 8947 8948 8950 8951 8953 8954 8956 8957
group		    	subset_H(H2O) id 8959 8960 8962 8963 8965 8966 8968 8969 8971 8972
group		    	subset_H(H2O) id 8974 8975 8977 8978 8980 8981 8983 8984 8986 8987
group		    	subset_H(H2O) id 8989 8990 8992 8993 8995 8996 8998 8999 9001 9002
group		    	subset_H(H2O) id 9004 9005 9007 9008 9010 9011 9013 9014 9016 9017
group		    	subset_H(H2O) id 9019 9020 9022 9023 9025 9026 9028 9029 9031 9032
group		    	subset_H(H2O) id 9034 9035 9037 9038 9040 9041 9043 9044 9046 9047
group		    	subset_H(H2O) id 9049 9050 9052 9053 9055 9056 9058 9059 9061 9062
group		    	subset_H(H2O) id 9064 9065 9067 9068 9070 9071 9073 9074 9076 9077
group		    	subset_H(H2O) id 9079 9080 9082 9083 9085 9086 9088 9089 9091 9092
group		    	subset_H(H2O) id 9094 9095 9097 9098 9100 9101 9103 9104 9106 9107
group		    	subset_H(H2O) id 9109 9110 9112 9113 9115 9116 9118 9119 9121 9122
group		    	subset_H(H2O) id 9124 9125 9127 9128 9130 9131 9133 9134 9136 9137
group		    	subset_H(H2O) id 9139 9140 9142 9143 9145 9146 9148 9149 9151 9152
group		    	subset_H(H2O) id 9154 9155 9157 9158 9160 9161 9163 9164 9166 9167
group		    	subset_H(H2O) id 9169 9170 9172 9173 9175 9176 9178 9179 9181 9182
group		    	subset_H(H2O) id 9184 9185 9187 9188

log			2.2_Minimization.out
#-------------------------------------------------------------------------------
# Stage 2.2: Minimization
#-------------------------------------------------------------------------------

min_style		cg
min_modify		dmax 0.05 line fast
reset_timestep		0
thermo_style		custom step fmax fnorm press vol v_sysdensity v_sxx v_syy v_szz v_syz v_sxz v_sxy pe v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_sum_fx v_sum_fy v_sum_fz
dump			sci all custom 3000 2.2.xyz id mol type q xs ys zs
thermo			100
minimize		0.0 1.0 3000 10000

undump			sci

log		    	2.3_Velocities.out
#-------------------------------------------------------------------------------
# Stage 2.3: Set the initial velocities for $T
#-------------------------------------------------------------------------------

velocity     	    	all create 298.2 72489 dist gaussian mom yes rot no


log			2.4_NPT.out
#-------------------------------------------------------------------------------
# Stage 2.4: NPT integration for 2 ns with a timestep of 1 fs
#             Temperature 298.2 K
#             Pressure 1 bar
#-------------------------------------------------------------------------------

# Fix shake is needed for water in pcff+
fix			shaken all shake .0001 20 50000 a 2 b 2

change_box 		all triclinic
kspace_style        pppm 0.00001
reset_timestep		0
thermo_style		custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke
thermo			${Nthermo}
change_box 		all triclinic
fix			1 movable npt temp 298.2 298.2 100 z 0.9869233 0.9869233 100  drag 0 mtk yes nreset 40000
fix			2 movable ave/time 1 199999 200000 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.4_averages.txt off 1
fix			3 movable ave/time 200  1     200 v_time  c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.4_instantaneous.txt
restart 		2000000 2.4.restart
dump 			sci all custom 200000 2.4.xyz id mol type q xs ys zs

timestep		1
run			2000000
undump 			sci

restart 			0
dump			sci all custom 2000000 2.4.xyz id mol type q xs ys zs
run 			0
undump			sci


unfix			1
unfix			2
unfix			3
unfix			shaken

log		2.5_Custom.out
#-------------------------------------------------------------------------------
# Stage 2.5: Custom code
#-------------------------------------------------------------------------------

thermo_style      custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke
thermo            100
fix               1 all nvt temp 298.2 298.2 100 drag 0.0

compute myADF1 all adf 200 4 7 7 0 2.7 0 2.7
fix ad1 all ave/time 1 3000000 3000000 c_myADF1[*] file tmp1.adf mode vector

compute myADF2 all adf 200 6 7 7 0 3.3 0 3.3
fix ad2 all ave/time 1 3000000 3000000 c_myADF2[*] file tmp2.adf mode vector

compute myADF3 all adf 200 3 7 7 0 3.5 0 3.5
fix ad3 all ave/time 1 3000000 3000000 c_myADF3[*] file tmp3.adf mode vector

compute myADF4 all adf 200 5 7 7 0 3.0 0 3.0
fix ad4 all ave/time 1 3000000 3000000 c_myADF4[*] file tmp4.adf mode vector

compute myADF5 all adf 200 1 7 7 0 3.9 0 3.9
fix ad5 all ave/time 1 3000000 3000000 c_myADF5[*] file tmp5.adf mode vector

compute myADF6 all adf 200 9 7 7 0 4.8 0 4.8
fix ad6 all ave/time 1 3000000 3000000 c_myADF6[*] file tmp6.adf mode vector

timestep             0.2
run                  3000000

unfix 1
unfix ad1
unfix ad2
unfix ad3
unfix ad4
unfix ad5
unfix ad6