Mon Oct 10 17:00:56 PDT 2011 #!/bin/csh -f # input: # $1 receptor coordinates file; # $2 grid size; # $3 # of iterations for the first run; # $4 # of iterations for the second run; # $5 # of iterations for the third run; # output: # delphi map: $1.phi # log file: delphi.out set RECCRG = rec+sph.crg set GRID = 179 set n1 = 50 set n2 = 120 set n3 = 180 set DELPHI = $DOCK_BASE/private/delphi set VDW = vdw.siz #set VDW = vdw.hydren set CRG = amb.crg.oxt set PRM = genric.prm # three step focussing calc. if (-e ARCDAT) /bin/rm ARCDAT cp $VDW fort.11 cp $CRG fort.12 cp $RECCRG fort.13 echo "gsize=$GRID, perfil=20, bndcon=2, linit=$n1\n" >! fort.10 cat $PRM >> fort.10 time $DELPHI if ($status != 0) exit(100) echo "gsize=$GRID, perfil=60, bndcon=3, linit=$n2\n" >! fort.10 cat $PRM >> fort.10 mv fort.14 fort.18 time $DELPHI if ($status != 0) exit(100) echo "gsize=$GRID, perfil=90, bndcon=3, linit=$n3\n" >! fort.10 cat $PRM >> fort.10 mv fort.14 fort.18 time $DELPHI if ($status != 0) exit(100) /bin/rm fort.1[0-3] fort.18 ARCDAT mv fort.14 $RECCRG:r.phi ___________________DelPhi II____________________ / | A program to solve the PB equation | | in 3D, using non-linear form, incorporating | | 2 dielectric regions, ionic strength, periodic | | and focussing boundary conditions, utilizing | | stripped optimum successive over-relaxation | | and an improved algorithm for mapping the | | Mol. Surface to the finite-Difference grid | __________________ ____________________/ DelPhi II program started on Mon Oct 10 2011 at 17:00:56 opening parameter file fort.10 atom radii read from file fort.11 ! default extended atom radii based loosely ! on mike connolly's MS program- note H's are 0 # of radius parameter records: 18 atomic charges read from file fort.12 # of charge parameter records: 254 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 1966 !!! WARNING: HIZ 69 has a net charge of 0.2000 !!! WARNING: HIZ 84 has a net charge of 0.2000 !!! WARNING: HIY 97 has a net charge of 0.2000 !!! WARNING: HIZ 119 has a net charge of 0.2000 !!! WARNING: HIZ 123 has a net charge of 0.2000 !!! WARNING: HIZ 129 has a net charge of 0.2000 !!! WARNING: HIZ 230 has a net charge of 0.2000 !!! WARNING: HIZ 245 has a net charge of 0.2000 !!! WARNING: ZN6 647 has a net charge of 1.4000 !!! WARNING: ZN6 651 has a net charge of 1.4000 !!! WARNING: ZN6 652 has a net charge of 1.4000 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 7.9999998E-02 grid size : 179 scale,in grids/A, set to be: 0.6719770 object centred at (gu) : 0.0000000E+00 0.0000000E+00 0.0000000E+00 inner,outer dielectrics : 2.060000 78.50000 ionic strength (M) : 0.1450000 debye length (A) : 8.001814 ion exclusion radius (A) : 2.000000 probe radius (A) : 1.400000 boundary conditions : dipolar x,y,z periodic bc. flags : F F F # of linear iterations : 50 # of non-linear iterations : 0 concentration map output : F spherical charge distbn. : F INSIGHT format output : F site potential output : F modified atom file output : T map file label : qdiffxas: qdiffxs4 with an improved surfacing routine convergence graph turned off potential listings turned off box fill (%): 20.00000 xmin,xmax (A): -18.14600 28.39100 ymin,ymax (A): 2.243000 55.22100 zmin,zma (A): 111.7510 161.3650 x,y,z range (A): 46.53700 52.97800 49.61401 scale (grids/A): 0.6719770 object centre (A): 5.122499 28.73200 136.5580 number of atom coordinates read : 1966 total number of charged atoms : 1835 net assigned charge : -8.699987 assigned positive charge : 279.4951 centred at (gu) : 88.27779 91.45943 90.56959 assigned negative charge : -288.1950 centred at (gu) : 88.41496 91.17426 90.45670 start vw surface at 0.1800000 fill in re-entrant regions at 0.1900000 boundary points on vw surface= 6006 nv = 96 ne = 84 npr = 41512 no. pairs analyzed = 41512 no. exposed pairs = 2670 no. arc points = 19017 no. surface atoms = 894 nbur = 1072 mkacc time = 0.1200000 writing accessible surface arcs data to file ARCDAT grid for indexing accessible points = 1.400000 m= 360 mr = 3043 m= 0 mr = 47 time to grow re-entrant surface = 0.1400000 no. cavity mid-points inaccessible to solvent = 15 time to turn everything in is 0.3300000 number of dielectric boundary points 3276 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 1760 ## 145 charges are in solution ## iepsmp to db, and charging done at 0.4700000 number of grid points assigned charge 6896 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid -4.8405148E-08 -1.0632497E-08 midg,1,midg; midg,igrid,midg -6.6750651E-08 3.2038660E-08 1,midg,midg; igrid,midg,midg -6.2526544E-09 -4.1903167E-08 gauss-seidel spectral radius is 0.9884690 estimated iterations to convergence 76 setup time was (sec) 0.6000000 now iterating at: 17:00:57 rms-change max change #iterations 0.1610271 26.62651 at 10iterations 1.7524192E-02 3.618614 at 20iterations 1.9414951E-03 0.2797365 at 30iterations 2.1269695E-04 3.5049438E-02 at 40iterations 2.3403079E-05 3.4828186E-03 at 50iterations finished qdiffx linear iterations at : 17:01:00 total time elapsed so far: 3.770000 # loops : 50 mean,max change (kT/e) : 2.3403079E-05 3.4828186E-03 energy calculations done at 3.770000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 4.000000 DelPhi exited at 17:01:02 4.006u 0.162s 0:05.45 76.3% 0+0k 0+0io 0pf+0w ___________________DelPhi II____________________ / | A program to solve the PB equation | | in 3D, using non-linear form, incorporating | | 2 dielectric regions, ionic strength, periodic | | and focussing boundary conditions, utilizing | | stripped optimum successive over-relaxation | | and an improved algorithm for mapping the | | Mol. Surface to the finite-Difference grid | __________________ ____________________/ DelPhi II program started on Mon Oct 10 2011 at 17:01:02 opening parameter file fort.10 atom radii read from file fort.11 ! default extended atom radii based loosely ! on mike connolly's MS program- note H's are 0 # of radius parameter records: 18 atomic charges read from file fort.12 # of charge parameter records: 254 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 1966 !!! WARNING: HIZ 69 has a net charge of 0.2000 !!! WARNING: HIZ 84 has a net charge of 0.2000 !!! WARNING: HIY 97 has a net charge of 0.2000 !!! WARNING: HIZ 119 has a net charge of 0.2000 !!! WARNING: HIZ 123 has a net charge of 0.2000 !!! WARNING: HIZ 129 has a net charge of 0.2000 !!! WARNING: HIZ 230 has a net charge of 0.2000 !!! WARNING: HIZ 245 has a net charge of 0.2000 !!! WARNING: ZN6 647 has a net charge of 1.4000 !!! WARNING: ZN6 651 has a net charge of 1.4000 !!! WARNING: ZN6 652 has a net charge of 1.4000 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 7.0000000E-02 grid size : 179 scale,in grids/A, set to be: 2.015931 object centred at (gu) : 0.0000000E+00 0.0000000E+00 0.0000000E+00 inner,outer dielectrics : 2.060000 78.50000 ionic strength (M) : 0.1450000 debye length (A) : 8.001814 ion exclusion radius (A) : 2.000000 probe radius (A) : 1.400000 boundary conditions : focussing x,y,z periodic bc. flags : F F F # of linear iterations : 120 # of non-linear iterations : 0 concentration map output : F spherical charge distbn. : F INSIGHT format output : F site potential output : F modified atom file output : T map file label : qdiffxas: qdiffxs4 with an improved surfacing routine convergence graph turned off potential listings turned off box fill (%): 60.00000 xmin,xmax (A): -18.14600 28.39100 ymin,ymax (A): 2.243000 55.22100 zmin,zma (A): 111.7510 161.3650 x,y,z range (A): 46.53700 52.97800 49.61401 scale (grids/A): 2.015931 object centre (A): 5.122499 28.73200 136.5580 number of atom coordinates read : 1966 total number of charged atoms : 1835 net assigned charge : -8.699987 assigned positive charge : 279.4951 centred at (gu) : 84.83344 94.37835 91.70905 assigned negative charge : -288.1950 centred at (gu) : 85.24502 93.52280 91.37012 start vw surface at 0.1600000 fill in re-entrant regions at 0.2800000 boundary points on vw surface= 63639 reading accessible surface arcs data from file ARCDAT no. of arc points read = 19017 grid for indexing accessible points = 1.400000 m= 22198 mr = 41377 m= 3658 mr = 15846 m= 205 mr = 2204 m= 4 mr = 87 m= 0 mr = 0 time to grow re-entrant surface = 0.2800000 no. cavity mid-points inaccessible to solvent = 27 time to turn everything in is 0.6000000 number of dielectric boundary points 30190 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 73 iepsmp to db, and charging done at 0.7600000 number of grid points assigned charge 14545 setting boundary conditions focussing boundary condition read from file fort.18 focussing potential map: qdiffxas: qdiffxs4 with an improved surfacing routine original scale (grids/A) : 0.6719770 object centre at (A) : 5.122499 28.73200 136.5580 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -4.5054521E-02 -2.7265705E-02 midg,1,midg; midg,igrid,midg -5.5189520E-02 -1.1325704E-02 1,midg,midg; igrid,midg,midg 2.6385186E-02 -2.0467658E-02 gauss-seidel spectral radius is 0.9987136 estimated iterations to convergence 221 setup time was (sec) 0.9500000 now iterating at: 17:01:03 rms-change max change #iterations 4.026398 202.1111 at 10iterations 1.905417 70.00694 at 20iterations 0.9135787 33.73885 at 30iterations 0.4370669 15.88625 at 40iterations 0.2079193 6.942419 at 50iterations 9.9465631E-02 3.343196 at 60iterations 4.7371004E-02 1.541466 at 70iterations 2.2690998E-02 0.7780828 at 80iterations 1.0835934E-02 0.3416519 at 90iterations 5.1724757E-03 0.1578889 at 100iterations 2.4798880E-03 8.5129738E-02 at 110iterations 1.1808577E-03 3.7803650E-02 at 120iterations finished qdiffx linear iterations at : 17:01:11 total time elapsed so far: 8.700000 # loops : 120 mean,max change (kT/e) : 1.1808577E-03 3.7803650E-02 energy calculations done at 8.700000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 8.930000 DelPhi exited at 17:01:12 8.936u 0.205s 0:10.32 88.4% 0+0k 0+0io 0pf+0w ___________________DelPhi II____________________ / | A program to solve the PB equation | | in 3D, using non-linear form, incorporating | | 2 dielectric regions, ionic strength, periodic | | and focussing boundary conditions, utilizing | | stripped optimum successive over-relaxation | | and an improved algorithm for mapping the | | Mol. Surface to the finite-Difference grid | __________________ ____________________/ DelPhi II program started on Mon Oct 10 2011 at 17:01:13 opening parameter file fort.10 atom radii read from file fort.11 ! default extended atom radii based loosely ! on mike connolly's MS program- note H's are 0 # of radius parameter records: 18 atomic charges read from file fort.12 # of charge parameter records: 254 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 1966 !!! WARNING: HIZ 69 has a net charge of 0.2000 !!! WARNING: HIZ 84 has a net charge of 0.2000 !!! WARNING: HIY 97 has a net charge of 0.2000 !!! WARNING: HIZ 119 has a net charge of 0.2000 !!! WARNING: HIZ 123 has a net charge of 0.2000 !!! WARNING: HIZ 129 has a net charge of 0.2000 !!! WARNING: HIZ 230 has a net charge of 0.2000 !!! WARNING: HIZ 245 has a net charge of 0.2000 !!! WARNING: ZN6 647 has a net charge of 1.4000 !!! WARNING: ZN6 651 has a net charge of 1.4000 !!! WARNING: ZN6 652 has a net charge of 1.4000 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 7.9999998E-02 grid size : 179 scale,in grids/A, set to be: 3.023897 object centred at (gu) : 0.0000000E+00 0.0000000E+00 0.0000000E+00 inner,outer dielectrics : 2.060000 78.50000 ionic strength (M) : 0.1450000 debye length (A) : 8.001814 ion exclusion radius (A) : 2.000000 probe radius (A) : 1.400000 boundary conditions : focussing x,y,z periodic bc. flags : F F F # of linear iterations : 180 # of non-linear iterations : 0 concentration map output : F spherical charge distbn. : F INSIGHT format output : F site potential output : F modified atom file output : T map file label : qdiffxas: qdiffxs4 with an improved surfacing routine convergence graph turned off potential listings turned off box fill (%): 90.00000 xmin,xmax (A): -18.14600 28.39100 ymin,ymax (A): 2.243000 55.22100 zmin,zma (A): 111.7510 161.3650 x,y,z range (A): 46.53700 52.97800 49.61401 scale (grids/A): 3.023897 object centre (A): 5.122499 28.73200 136.5580 number of atom coordinates read : 1966 total number of charged atoms : 1835 net assigned charge : -8.699987 assigned positive charge : 279.4951 centred at (gu) : 82.25008 96.56752 92.56365 assigned negative charge : -288.1950 centred at (gu) : 82.86744 95.28418 92.05520 start vw surface at 0.1800000 fill in re-entrant regions at 0.5600000 boundary points on vw surface= 146858 reading accessible surface arcs data from file ARCDAT no. of arc points read = 19017 grid for indexing accessible points = 1.400000 m= 70927 mr = 100813 m= 25976 mr = 56451 m= 6359 mr = 20448 m= 979 mr = 4758 m= 78 mr = 670 m= 0 mr = 37 time to grow re-entrant surface = 0.8900000 no. cavity mid-points inaccessible to solvent = 81 time to turn everything in is 1.600000 number of dielectric boundary points 68000 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 0 iepsmp to db, and charging done at 1.760000 number of grid points assigned charge 14680 setting boundary conditions focussing boundary condition read from file fort.18 focussing potential map: qdiffxas: qdiffxs4 with an improved surfacing routine original scale (grids/A) : 2.015931 object centre at (A) : 5.122499 28.73200 136.5580 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -0.4666076 -0.2561345 midg,1,midg; midg,igrid,midg -0.4337623 -0.1928344 1,midg,midg; igrid,midg,midg 0.4047585 -0.1805411 gauss-seidel spectral radius is 0.9992748 estimated iterations to convergence 293 setup time was (sec) 1.990000 now iterating at: 17:01:15 rms-change max change #iterations 9.208879 317.3104 at 10iterations 5.390988 147.9424 at 20iterations 3.138286 74.81004 at 30iterations 1.822250 39.55180 at 40iterations 1.056025 25.04329 at 50iterations 0.6131252 13.73545 at 60iterations 0.3531457 6.729979 at 70iterations 0.2040289 3.748611 at 80iterations 0.1178775 2.555168 at 90iterations 6.7782663E-02 1.418740 at 100iterations 3.9186139E-02 0.7708057 at 110iterations 2.2572249E-02 0.4564409 at 120iterations 1.3041264E-02 0.2572441 at 130iterations 7.5161764E-03 0.1558646 at 140iterations 4.3335473E-03 8.0142975E-02 at 150iterations 2.5058175E-03 5.1759243E-02 at 160iterations 1.4480912E-03 2.5707245E-02 at 170iterations 8.3920313E-04 1.5762329E-02 at 180iterations finished qdiffx linear iterations at : 17:01:29 total time elapsed so far: 15.80000 # loops : 180 mean,max change (kT/e) : 8.3920313E-04 1.5762329E-02 energy calculations done at 15.80000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 16.05000 DelPhi exited at 17:01:30 16.049u 0.298s 0:16.99 96.1% 0+0k 0+0io 0pf+0w