Sat Oct 8 19:18:45 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 Sat Oct 8 2011 at 19:18:45 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 = 2552 !!! WARNING: HIZ 94 has a net charge of 0.2000 !!! WARNING: HIZ 96 has a net charge of 0.2000 !!! WARNING: HIY 119 has a net charge of 0.2000 !!! WARNING: ZN5 502 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= 3.9999999E-02 grid size : 179 scale,in grids/A, set to be: 0.7263379 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): -13.63400 35.37900 ymin,ymax (A): -26.92400 21.02600 zmin,zma (A): 32.88400 81.58100 x,y,z range (A): 49.01300 47.95000 48.69700 scale (grids/A): 0.7263379 object centre (A): 10.87250 -2.948999 57.23250 number of atom coordinates read : 2552 total number of charged atoms : 2444 net assigned charge : 4.000010 assigned positive charge : 389.4371 centred at (gu) : 90.04939 89.30747 90.21904 assigned negative charge : -385.4376 centred at (gu) : 90.02779 89.48403 90.18217 start vw surface at 9.0000004E-02 fill in re-entrant regions at 9.0000004E-02 boundary points on vw surface= 9339 nv = 96 ne = 84 npr = 52290 no. pairs analyzed = 52290 no. exposed pairs = 3020 no. arc points = 21700 no. surface atoms = 1022 nbur = 1530 mkacc time = 8.9999996E-02 writing accessible surface arcs data to file ARCDAT grid for indexing accessible points = 1.400000 m= 653 mr = 5519 m= 5 mr = 117 m= 0 mr = 0 time to grow re-entrant surface = 0.1000000 no. cavity mid-points inaccessible to solvent = 208 time to turn everything in is 0.2000000 number of dielectric boundary points 4361 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 1880 ## 156 charges are in solution ## iepsmp to db, and charging done at 0.2700000 number of grid points assigned charge 10197 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid 1.6288839E-08 8.8372744E-08 midg,1,midg; midg,igrid,midg 2.4099856E-07 -9.9725909E-08 1,midg,midg; igrid,midg,midg 3.0748875E-08 7.1469678E-08 gauss-seidel spectral radius is 0.9901657 estimated iterations to convergence 82 setup time was (sec) 0.4100000 now iterating at: 19:18:46 rms-change max change #iterations 0.2508639 32.97848 at 10iterations 3.2562833E-02 4.912511 at 20iterations 4.2580124E-03 0.5200882 at 30iterations 5.6354865E-04 7.7663898E-02 at 40iterations 7.2926443E-05 8.6288452E-03 at 50iterations finished qdiffx linear iterations at : 19:18:49 total time elapsed so far: 3.530000 # loops : 50 mean,max change (kT/e) : 7.2926443E-05 8.6288452E-03 energy calculations done at 3.530000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 3.660000 DelPhi exited at 19:18:49 3.661u 0.224s 0:04.16 93.2% 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 Sat Oct 8 2011 at 19:18:49 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 = 2552 !!! WARNING: HIZ 94 has a net charge of 0.2000 !!! WARNING: HIZ 96 has a net charge of 0.2000 !!! WARNING: HIY 119 has a net charge of 0.2000 !!! WARNING: ZN5 502 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= 3.9999999E-02 grid size : 179 scale,in grids/A, set to be: 2.179014 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): -13.63400 35.37900 ymin,ymax (A): -26.92400 21.02600 zmin,zma (A): 32.88400 81.58100 x,y,z range (A): 49.01300 47.95000 48.69700 scale (grids/A): 2.179014 object centre (A): 10.87250 -2.948999 57.23250 number of atom coordinates read : 2552 total number of charged atoms : 2444 net assigned charge : 4.000010 assigned positive charge : 389.4371 centred at (gu) : 90.14785 87.92191 90.65677 assigned negative charge : -385.4376 centred at (gu) : 90.08321 88.45178 90.54627 start vw surface at 7.0000000E-02 fill in re-entrant regions at 0.1900000 boundary points on vw surface= 98223 reading accessible surface arcs data from file ARCDAT no. of arc points read = 21700 grid for indexing accessible points = 1.400000 m= 41353 mr = 70175 m= 8939 mr = 32366 m= 852 mr = 6212 m= 17 mr = 413 m= 0 mr = 3 time to grow re-entrant surface = 0.2600000 no. cavity mid-points inaccessible to solvent = 5699 time to turn everything in is 0.4900000 number of dielectric boundary points 40215 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 74 iepsmp to db, and charging done at 0.5700000 number of grid points assigned charge 19463 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.7263379 object centre at (A) : 10.87250 -2.948999 57.23250 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -3.0636932E-03 2.4293080E-02 midg,1,midg; midg,igrid,midg 6.7501739E-02 -2.7203448E-02 1,midg,midg; igrid,midg,midg -2.8505264E-02 9.6549029E-03 gauss-seidel spectral radius is 0.9989296 estimated iterations to convergence 242 setup time was (sec) 0.7200000 now iterating at: 19:18:50 rms-change max change #iterations 5.488468 201.1170 at 10iterations 2.795431 84.00446 at 20iterations 1.421817 42.91013 at 30iterations 0.7253547 20.04947 at 40iterations 0.3708071 11.91920 at 50iterations 0.1897209 4.997478 at 60iterations 9.6711114E-02 2.397546 at 70iterations 4.9460325E-02 1.395601 at 80iterations 2.5250509E-02 0.6566906 at 90iterations 1.2924959E-02 0.3789272 at 100iterations 6.5996102E-03 0.1831093 at 110iterations 3.3782856E-03 8.5714340E-02 at 120iterations finished qdiffx linear iterations at : 19:18:58 total time elapsed so far: 8.340000 # loops : 120 mean,max change (kT/e) : 3.3782856E-03 8.5714340E-02 energy calculations done at 8.340000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 8.480000 DelPhi exited at 19:18:58 8.482u 0.245s 0:09.08 96.0% 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 Sat Oct 8 2011 at 19:18:58 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 = 2552 !!! WARNING: HIZ 94 has a net charge of 0.2000 !!! WARNING: HIZ 96 has a net charge of 0.2000 !!! WARNING: HIY 119 has a net charge of 0.2000 !!! WARNING: ZN5 502 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= 3.9999999E-02 grid size : 179 scale,in grids/A, set to be: 3.268521 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): -13.63400 35.37900 ymin,ymax (A): -26.92400 21.02600 zmin,zma (A): 32.88400 81.58100 x,y,z range (A): 49.01300 47.95000 48.69700 scale (grids/A): 3.268521 object centre (A): 10.87250 -2.948999 57.23250 number of atom coordinates read : 2552 total number of charged atoms : 2444 net assigned charge : 4.000010 assigned positive charge : 389.4371 centred at (gu) : 90.22170 86.88275 90.98510 assigned negative charge : -385.4376 centred at (gu) : 90.12463 87.67757 90.81922 start vw surface at 7.9999998E-02 fill in re-entrant regions at 0.6400000 boundary points on vw surface= 226066 reading accessible surface arcs data from file ARCDAT no. of arc points read = 21700 grid for indexing accessible points = 1.400000 m= 126377 mr = 168120 m= 54967 mr = 107082 m= 17177 mr = 46338 m= 3604 mr = 13656 m= 396 mr = 2643 m= 15 mr = 232 m= 0 mr = 2 time to grow re-entrant surface = 0.8099999 no. cavity mid-points inaccessible to solvent = 18484 time to turn everything in is 1.560000 number of dielectric boundary points 90529 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 0 iepsmp to db, and charging done at 1.680000 number of grid points assigned charge 19552 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.179014 object centre at (A) : 10.87250 -2.948999 57.23250 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -4.9215447E-02 0.5163171 midg,1,midg; midg,igrid,midg 0.4107419 -0.2254488 1,midg,midg; igrid,midg,midg -0.6840577 0.2082736 gauss-seidel spectral radius is 0.9993255 estimated iterations to convergence 304 setup time was (sec) 1.850000 now iterating at: 19:19:01 rms-change max change #iterations 11.84352 342.2614 at 10iterations 7.035437 136.1355 at 20iterations 4.163350 70.25263 at 30iterations 2.463559 44.62391 at 40iterations 1.458767 24.90740 at 50iterations 0.8610149 12.94977 at 60iterations 0.5087388 10.41117 at 70iterations 0.3003420 4.817945 at 80iterations 0.1771510 2.688486 at 90iterations 0.1044530 1.784181 at 100iterations 6.1730254E-02 0.9747620 at 110iterations 3.6424994E-02 0.5558014 at 120iterations 2.1516811E-02 0.3346863 at 130iterations 1.2693195E-02 0.1890440 at 140iterations 7.4853371E-03 0.1102352 at 150iterations 4.4319849E-03 7.8515053E-02 at 160iterations 2.6228442E-03 4.0809631E-02 at 170iterations 1.5601603E-03 2.4174690E-02 at 180iterations finished qdiffx linear iterations at : 19:19:14 total time elapsed so far: 15.32000 # loops : 180 mean,max change (kT/e) : 1.5601603E-03 2.4174690E-02 energy calculations done at 15.32000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 15.45000 DelPhi exited at 19:19:14 15.456u 0.288s 0:16.01 98.2% 0+0k 0+0io 0pf+0w