Sun Oct 9 07:03:07 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 Sun Oct 9 2011 at 07:03:07 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: 284 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 2001 !!! WARNING: ASP 186 has a net charge of -2.0000 !!! WARNING: NDP 187 has a net charge of -3.0000 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 5.0000001E-02 grid size : 179 scale,in grids/A, set to be: 0.6621285 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): -24.26800 26.21300 ymin,ymax (A): -45.88700 -1.027000 zmin,zma (A): -56.02800 -2.262000 x,y,z range (A): 50.48100 44.86000 53.76600 scale (grids/A): 0.6621285 object centre (A): 0.9725008 -23.45700 -29.14500 number of atom coordinates read : 2001 total number of charged atoms : 1862 net assigned charge : -2.499985 assigned positive charge : 299.6560 centred at (gu) : 89.11776 90.39911 91.50404 assigned negative charge : -302.1561 centred at (gu) : 89.05792 90.20154 91.55605 start vw surface at 0.1500000 fill in re-entrant regions at 0.1600000 boundary points on vw surface= 5789 nv = 96 ne = 84 npr = 42749 no. pairs analyzed = 42749 no. exposed pairs = 2410 no. arc points = 17303 no. surface atoms = 825 nbur = 1176 mkacc time = 0.1200000 writing accessible surface arcs data to file ARCDAT grid for indexing accessible points = 1.400000 m= 288 mr = 3171 m= 1 mr = 36 m= 0 mr = 0 time to grow re-entrant surface = 0.1400000 no. cavity mid-points inaccessible to solvent = 39 time to turn everything in is 0.3000000 number of dielectric boundary points 2871 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 1634 ## 143 charges are in solution ## iepsmp to db, and charging done at 0.4300000 number of grid points assigned charge 6774 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid 1.2350391E-09 -2.0596234E-08 midg,1,midg; midg,igrid,midg -3.6132228E-08 2.7019244E-08 1,midg,midg; igrid,midg,midg -1.9628491E-08 2.4882318E-09 gauss-seidel spectral radius is 0.9881204 estimated iterations to convergence 75 setup time was (sec) 0.5800000 now iterating at: 07:03:07 rms-change max change #iterations 0.1530971 24.93024 at 10iterations 1.6310772E-02 2.507682 at 20iterations 1.7158425E-03 0.2642140 at 30iterations 1.8167357E-04 2.8351307E-02 at 40iterations 1.9399975E-05 3.0117035E-03 at 50iterations finished qdiffx linear iterations at : 07:03:11 total time elapsed so far: 3.840000 # loops : 50 mean,max change (kT/e) : 1.9399975E-05 3.0117035E-03 energy calculations done at 3.840000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 4.010000 DelPhi exited at 07:03:11 4.018u 0.256s 0:04.53 94.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 Sun Oct 9 2011 at 07:03:11 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: 284 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 2001 !!! WARNING: ASP 186 has a net charge of -2.0000 !!! WARNING: NDP 187 has a net charge of -3.0000 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 5.0000001E-02 grid size : 179 scale,in grids/A, set to be: 1.986385 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): -24.26800 26.21300 ymin,ymax (A): -45.88700 -1.027000 zmin,zma (A): -56.02800 -2.262000 x,y,z range (A): 50.48100 44.86000 53.76600 scale (grids/A): 1.986385 object centre (A): 0.9725008 -23.45700 -29.14500 number of atom coordinates read : 2001 total number of charged atoms : 1862 net assigned charge : -2.499985 assigned positive charge : 299.6560 centred at (gu) : 87.35352 91.19740 94.51225 assigned negative charge : -302.1561 centred at (gu) : 87.17401 90.60471 94.66832 start vw surface at 0.1500000 fill in re-entrant regions at 0.2700000 boundary points on vw surface= 62518 reading accessible surface arcs data from file ARCDAT no. of arc points read = 17303 grid for indexing accessible points = 1.400000 m= 21809 mr = 43157 m= 3321 mr = 16159 m= 207 mr = 2032 m= 2 mr = 73 m= 1 mr = 1 m= 0 mr = 0 time to grow re-entrant surface = 0.2600000 no. cavity mid-points inaccessible to solvent = 1314 time to turn everything in is 0.5700000 number of dielectric boundary points 26436 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 106 iepsmp to db, and charging done at 0.7000000 number of grid points assigned charge 14718 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.6621285 object centre at (A) : 0.9725008 -23.45700 -29.14500 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 2.4078688E-02 -5.8194634E-02 midg,1,midg; midg,igrid,midg -3.1191964E-02 4.1486476E-02 1,midg,midg; igrid,midg,midg -4.3312475E-02 -1.0399003E-02 gauss-seidel spectral radius is 0.9986753 estimated iterations to convergence 218 setup time was (sec) 0.8700000 now iterating at: 07:03:12 rms-change max change #iterations 4.057086 192.5409 at 10iterations 1.912395 71.78093 at 20iterations 0.9031904 31.22960 at 30iterations 0.4278265 17.10281 at 40iterations 0.2027030 6.812346 at 50iterations 9.5534936E-02 3.540417 at 60iterations 4.5369130E-02 1.602968 at 70iterations 2.1526812E-02 0.7544470 at 80iterations 1.0126814E-02 0.3624735 at 90iterations 4.8172176E-03 0.1675541 at 100iterations 2.2839941E-03 7.3838234E-02 at 110iterations 1.0867256E-03 3.6954880E-02 at 120iterations finished qdiffx linear iterations at : 07:03:21 total time elapsed so far: 9.050000 # loops : 120 mean,max change (kT/e) : 1.0867256E-03 3.6954880E-02 energy calculations done at 9.050000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 9.240000 DelPhi exited at 07:03:21 9.239u 0.288s 0:09.79 97.1% 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 Sun Oct 9 2011 at 07:03:21 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: 284 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 2001 !!! WARNING: ASP 186 has a net charge of -2.0000 !!! WARNING: NDP 187 has a net charge of -3.0000 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 5.0000001E-02 grid size : 179 scale,in grids/A, set to be: 2.979578 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): -24.26800 26.21300 ymin,ymax (A): -45.88700 -1.027000 zmin,zma (A): -56.02800 -2.262000 x,y,z range (A): 50.48100 44.86000 53.76600 scale (grids/A): 2.979578 object centre (A): 0.9725008 -23.45700 -29.14500 number of atom coordinates read : 2001 total number of charged atoms : 1862 net assigned charge : -2.499985 assigned positive charge : 299.6560 centred at (gu) : 86.03024 91.79594 96.76836 assigned negative charge : -302.1561 centred at (gu) : 85.76102 90.90704 97.00253 start vw surface at 0.1400000 fill in re-entrant regions at 0.5100000 boundary points on vw surface= 144545 reading accessible surface arcs data from file ARCDAT no. of arc points read = 17303 grid for indexing accessible points = 1.400000 m= 71055 mr = 104410 m= 25000 mr = 58258 m= 6102 mr = 20175 m= 919 mr = 4585 m= 58 mr = 599 m= 1 mr = 25 m= 0 mr = 0 time to grow re-entrant surface = 0.8300000 no. cavity mid-points inaccessible to solvent = 4527 time to turn everything in is 1.440000 number of dielectric boundary points 59628 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 0 iepsmp to db, and charging done at 1.590000 number of grid points assigned charge 14892 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) : 1.986385 object centre at (A) : 0.9725008 -23.45700 -29.14500 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 0.8059314 -0.8427657 midg,1,midg; midg,igrid,midg -0.3425573 0.4603940 1,midg,midg; igrid,midg,midg -0.4081455 -0.1166536 gauss-seidel spectral radius is 0.9992569 estimated iterations to convergence 290 setup time was (sec) 1.750000 now iterating at: 07:03:23 rms-change max change #iterations 9.305408 382.6743 at 10iterations 5.331792 120.5898 at 20iterations 3.071851 63.22128 at 30iterations 1.760153 35.44805 at 40iterations 1.008456 21.08786 at 50iterations 0.5782892 11.28810 at 60iterations 0.3305440 6.434589 at 70iterations 0.1899942 3.880623 at 80iterations 0.1087737 2.070055 at 90iterations 6.2281296E-02 1.174694 at 100iterations 3.5633940E-02 0.7297058 at 110iterations 2.0432314E-02 0.3629036 at 120iterations 1.1706850E-02 0.2271852 at 130iterations 6.6969758E-03 0.1240005 at 140iterations 3.8396912E-03 7.4668884E-02 at 150iterations 2.1944123E-03 4.2169571E-02 at 160iterations 1.2584218E-03 2.5285721E-02 at 170iterations 7.2146655E-04 1.4377594E-02 at 180iterations finished qdiffx linear iterations at : 07:03:36 total time elapsed so far: 15.00000 # loops : 180 mean,max change (kT/e) : 7.2146655E-04 1.4377594E-02 energy calculations done at 15.00000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 15.18000 DelPhi exited at 07:03:37 15.178u 0.310s 0:15.75 98.2% 0+0k 0+0io 0pf+0w