Sun Oct 9 03:56:43 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 03:56:43 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: 291 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 4627 !!! WARNING: HIZ 357 has a net charge of 0.2000 !!! WARNING: HEM 2209 has a net charge of 1.4000 !!! WARNING: FE 5 has a net charge of 1.4000 !!! WARNING: FE 6 has a net charge of 1.4000 !!! WARNING: FE 7 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= 0.1000000 grid size : 179 scale,in grids/A, set to be: 0.3075727 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): 3.763000 84.72600 ymin,ymax (A): -56.14600 11.93500 zmin,zma (A): -42.01500 73.73000 x,y,z range (A): 80.96300 68.08100 115.7450 scale (grids/A): 0.3075727 object centre (A): 44.24450 -22.10550 15.85750 number of atom coordinates read : 4627 total number of charged atoms : 4310 net assigned charge : 21.80001 assigned positive charge : 673.5243 centred at (gu) : 87.05880 88.86136 82.01704 assigned negative charge : -651.7288 centred at (gu) : 86.94517 88.84684 81.95345 start vw surface at 0.1900000 fill in re-entrant regions at 0.2000000 boundary points on vw surface= 2266 nv = 96 ne = 84 npr = 93066 no. pairs analyzed = 93066 no. exposed pairs = 6011 no. arc points = 42335 no. surface atoms = 2010 nbur = 2617 mkacc time = 0.2100000 writing accessible surface arcs data to file ARCDAT initial cube size too small, in assigning accessible points to a grid therefore rescaling... initial cube size too small, in assigning accessible points to a grid therefore rescaling... grid for indexing accessible points = 2.766400 m= 33 mr = 895 m= 0 mr = 0 time to grow re-entrant surface = 0.2200000 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.4300000 number of dielectric boundary points 1404 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 1977 ## 674 charges are in solution ## iepsmp to db, and charging done at 0.5600000 number of grid points assigned charge 3140 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid 7.2996564E-16 6.7955538E-18 midg,1,midg; midg,igrid,midg 1.3521391E-16 7.2374585E-17 1,midg,midg; igrid,midg,midg -1.1401601E-16 6.5696590E-17 gauss-seidel spectral radius is 0.9470770 estimated iterations to convergence 37 setup time was (sec) 0.7100000 now iterating at: 03:56:44 rms-change max change #iterations 5.8493102E-03 2.187862 at 10iterations 6.2465074E-04 0.2142601 at 20iterations 8.8970446E-05 3.5074234E-02 at 30iterations 1.3304557E-05 5.5007935E-03 at 40iterations 2.0168095E-06 8.5067749E-04 at 50iterations finished qdiffx linear iterations at : 03:56:47 total time elapsed so far: 4.420000 # loops : 50 mean,max change (kT/e) : 2.0168095E-06 8.5067749E-04 energy calculations done at 4.420000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 4.600000 DelPhi exited at 03:56:48 4.600u 0.247s 0:05.19 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 Sun Oct 9 2011 at 03:56:48 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: 291 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 4627 !!! WARNING: HIZ 357 has a net charge of 0.2000 !!! WARNING: HEM 2209 has a net charge of 1.4000 !!! WARNING: FE 5 has a net charge of 1.4000 !!! WARNING: FE 6 has a net charge of 1.4000 !!! WARNING: FE 7 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= 0.1300000 grid size : 179 scale,in grids/A, set to be: 0.9227180 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): 3.763000 84.72600 ymin,ymax (A): -56.14600 11.93500 zmin,zma (A): -42.01500 73.73000 x,y,z range (A): 80.96300 68.08100 115.7450 scale (grids/A): 0.9227180 object centre (A): 44.24450 -22.10550 15.85750 number of atom coordinates read : 4627 total number of charged atoms : 4310 net assigned charge : 21.80001 assigned positive charge : 673.5243 centred at (gu) : 81.17487 86.58236 66.04937 assigned negative charge : -651.7288 centred at (gu) : 80.83502 86.54022 65.86034 start vw surface at 0.2300000 fill in re-entrant regions at 0.2700000 boundary points on vw surface= 28968 reading accessible surface arcs data from file ARCDAT no. of arc points read = 42335 initial cube size too small, in assigning accessible points to a grid therefore rescaling... initial cube size too small, in assigning accessible points to a grid therefore rescaling... grid for indexing accessible points = 2.591680 m= 3476 mr = 17337 m= 49 mr = 1312 m= 0 mr = 0 time to grow re-entrant surface = 0.1300000 no. cavity mid-points inaccessible to solvent = 317 time to turn everything in is 0.4300000 number of dielectric boundary points 13844 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 1672 ## 142 charges are in solution ## iepsmp to db, and charging done at 0.5600000 number of grid points assigned charge 23472 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.3075727 object centre at (A) : 44.24450 -22.10550 15.85750 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 2.8326022E-03 3.6630209E-04 midg,1,midg; midg,igrid,midg 8.1440048E-05 5.6548492E-05 1,midg,midg; igrid,midg,midg -1.1357245E-04 8.1456819E-05 gauss-seidel spectral radius is 0.9940647 estimated iterations to convergence 105 setup time was (sec) 0.7500000 now iterating at: 03:56:49 rms-change max change #iterations 0.8398941 86.29123 at 10iterations 0.1747698 11.03018 at 20iterations 3.6655352E-02 2.875378 at 30iterations 8.7348968E-03 0.6213226 at 40iterations 2.7836496E-03 0.1704807 at 50iterations 1.2375028E-03 6.9444299E-02 at 60iterations 6.0975994E-04 3.0348778E-02 at 70iterations 3.1007419E-04 1.5465736E-02 at 80iterations 1.5870886E-04 7.9498291E-03 at 90iterations 8.1321872E-05 4.0912628E-03 at 100iterations 4.1732801E-05 2.1066666E-03 at 110iterations 2.1427328E-05 1.0852814E-03 at 120iterations finished qdiffx linear iterations at : 03:56:57 total time elapsed so far: 9.150000 # loops : 120 mean,max change (kT/e) : 2.1427328E-05 1.0852814E-03 energy calculations done at 9.150000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 9.330000 DelPhi exited at 03:56:58 9.328u 0.332s 0:09.99 96.5% 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 03:56: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: 291 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 4627 !!! WARNING: HIZ 357 has a net charge of 0.2000 !!! WARNING: HEM 2209 has a net charge of 1.4000 !!! WARNING: FE 5 has a net charge of 1.4000 !!! WARNING: FE 6 has a net charge of 1.4000 !!! WARNING: FE 7 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= 0.1300000 grid size : 179 scale,in grids/A, set to be: 1.384077 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): 3.763000 84.72600 ymin,ymax (A): -56.14600 11.93500 zmin,zma (A): -42.01500 73.73000 x,y,z range (A): 80.96300 68.08100 115.7450 scale (grids/A): 1.384077 object centre (A): 44.24450 -22.10550 15.85750 number of atom coordinates read : 4627 total number of charged atoms : 4310 net assigned charge : 21.80001 assigned positive charge : 673.5243 centred at (gu) : 76.76196 84.87318 54.07364 assigned negative charge : -651.7288 centred at (gu) : 76.25234 84.81027 53.79034 start vw surface at 0.2500000 fill in re-entrant regions at 0.3500000 boundary points on vw surface= 69364 reading accessible surface arcs data from file ARCDAT no. of arc points read = 42335 initial cube size too small, in assigning accessible points to a grid therefore rescaling... initial cube size too small, in assigning accessible points to a grid therefore rescaling... grid for indexing accessible points = 2.562560 m= 16763 mr = 45105 m= 1010 mr = 10327 m= 7 mr = 281 m= 0 mr = 0 time to grow re-entrant surface = 0.3400000 no. cavity mid-points inaccessible to solvent = 1106 time to turn everything in is 0.7800000 number of dielectric boundary points 31431 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 688 ## 18 charges are in solution ## iepsmp to db, and charging done at 0.9300000 number of grid points assigned charge 31087 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.9227180 object centre at (A) : 44.24450 -22.10550 15.85750 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 0.2608012 1.0945798E-02 midg,1,midg; midg,igrid,midg 4.7333352E-03 3.1523057E-03 1,midg,midg; igrid,midg,midg -7.9138670E-03 9.0662204E-03 gauss-seidel spectral radius is 0.9974021 estimated iterations to convergence 156 setup time was (sec) 1.120000 now iterating at: 03:56:59 rms-change max change #iterations 3.281251 148.5878 at 10iterations 1.210652 56.68987 at 20iterations 0.4405806 19.12888 at 30iterations 0.1584791 6.877455 at 40iterations 5.7315566E-02 2.490238 at 50iterations 2.0433735E-02 0.8214855 at 60iterations 7.3245917E-03 0.3021545 at 70iterations 2.6232121E-03 0.1081485 at 80iterations 9.4430853E-04 3.9688110E-02 at 90iterations 3.4327497E-04 1.3652802E-02 at 100iterations 1.2663699E-04 5.2185059E-03 at 110iterations 4.8358110E-05 2.3612976E-03 at 120iterations 1.9339637E-05 1.0833740E-03 at 130iterations 8.2632232E-06 4.9591064E-04 at 140iterations 3.7736127E-06 2.6702881E-04 at 150iterations 1.9085767E-06 1.4877319E-04 at 160iterations 1.1477654E-06 9.9182129E-05 at 170iterations 8.7265647E-07 6.1035156E-05 at 180iterations finished qdiffx linear iterations at : 03:57:12 total time elapsed so far: 13.87000 # loops : 180 mean,max change (kT/e) : 8.7265647E-07 6.1035156E-05 energy calculations done at 13.87000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 14.05000 DelPhi exited at 03:57:13 14.053u 0.308s 0:14.63 98.0% 0+0k 0+0io 0pf+0w