Sat Oct 15 13:55:06 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 15 2011 at 13:55:06 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: 267 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 2668 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.6070942 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): 10.11200 60.97100 ymin,ymax (A): 9.696000 59.24400 zmin,zma (A): 8.711000 67.35100 x,y,z range (A): 50.85900 49.54800 58.64000 scale (grids/A): 0.6070942 object centre (A): 35.54150 34.47000 38.03100 number of atom coordinates read : 2668 total number of charged atoms : 2499 net assigned charge : 2.500012 assigned positive charge : 391.7929 centred at (gu) : 90.61850 88.94807 88.93442 assigned negative charge : -389.2935 centred at (gu) : 90.52785 89.01888 89.01552 start vw surface at 7.9999998E-02 fill in re-entrant regions at 9.0000004E-02 boundary points on vw surface= 6588 nv = 96 ne = 84 npr = 53847 no. pairs analyzed = 53847 no. exposed pairs = 3563 no. arc points = 24776 no. surface atoms = 1204 nbur = 1464 mkacc time = 9.0000004E-02 writing accessible surface arcs data to file ARCDAT initial cube size too small, in assigning accessible points to a grid therefore rescaling... grid for indexing accessible points = 1.456000 m= 300 mr = 3385 m= 1 mr = 26 m= 0 mr = 0 time to grow re-entrant surface = 0.1000000 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.2000000 number of dielectric boundary points 3478 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 2516 ## 280 charges are in solution ## iepsmp to db, and charging done at 0.2700000 number of grid points assigned charge 8090 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid 6.2207057E-09 -1.9323347E-09 midg,1,midg; midg,igrid,midg 5.6164309E-09 -1.5642676E-09 1,midg,midg; igrid,midg,midg -2.5285964E-09 6.1080900E-09 gauss-seidel spectral radius is 0.9859663 estimated iterations to convergence 69 setup time was (sec) 0.4000000 now iterating at: 13:55:06 rms-change max change #iterations 0.1166419 16.69517 at 10iterations 9.9601727E-03 1.584629 at 20iterations 8.7689125E-04 0.1378174 at 30iterations 8.6130953E-05 1.3562202E-02 at 40iterations 1.1502282E-05 2.2315979E-03 at 50iterations finished qdiffx linear iterations at : 13:55:09 total time elapsed so far: 3.410000 # loops : 50 mean,max change (kT/e) : 1.1502282E-05 2.2315979E-03 energy calculations done at 3.410000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 3.530000 DelPhi exited at 13:55:10 3.535u 0.217s 0:04.11 90.9% 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 15 2011 at 13:55:10 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: 267 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 2668 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: 1.821282 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): 10.11200 60.97100 ymin,ymax (A): 9.696000 59.24400 zmin,zma (A): 8.711000 67.35100 x,y,z range (A): 50.85900 49.54800 58.64000 scale (grids/A): 1.821282 object centre (A): 35.54150 34.47000 38.03100 number of atom coordinates read : 2668 total number of charged atoms : 2499 net assigned charge : 2.500012 assigned positive charge : 391.7929 centred at (gu) : 91.85507 86.84385 86.80269 assigned negative charge : -389.2935 centred at (gu) : 91.58322 87.05637 87.04620 start vw surface at 9.0000004E-02 fill in re-entrant regions at 0.1700000 boundary points on vw surface= 70981 reading accessible surface arcs data from file ARCDAT no. of arc points read = 24776 grid for indexing accessible points = 1.400000 m= 24307 mr = 47481 m= 3491 mr = 17250 m= 175 mr = 2007 m= 4 mr = 63 m= 0 mr = 0 time to grow re-entrant surface = 0.1700000 no. cavity mid-points inaccessible to solvent = 209 time to turn everything in is 0.3800000 number of dielectric boundary points 32157 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 172 ## 1 charges are in solution ## iepsmp to db, and charging done at 0.4700000 number of grid points assigned charge 19567 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.6070942 object centre at (A) : 35.54150 34.47000 38.03100 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 2.9854709E-02 -2.0705013E-02 midg,1,midg; midg,igrid,midg 1.7458742E-02 3.5929068E-03 1,midg,midg; igrid,midg,midg -9.6816849E-03 2.8966290E-03 gauss-seidel spectral radius is 0.9985238 estimated iterations to convergence 206 setup time was (sec) 0.6400000 now iterating at: 13:55:11 rms-change max change #iterations 3.895944 172.7562 at 10iterations 1.770876 66.08340 at 20iterations 0.8041475 30.27505 at 30iterations 0.3646057 11.91339 at 40iterations 0.1656106 5.560261 at 50iterations 7.5242050E-02 2.504967 at 60iterations 3.4183223E-02 1.146482 at 70iterations 1.5495269E-02 0.5224438 at 80iterations 7.0497193E-03 0.2220688 at 90iterations 3.1939591E-03 0.1034182 at 100iterations 1.4565043E-03 4.6741486E-02 at 110iterations 6.6207745E-04 2.3109436E-02 at 120iterations finished qdiffx linear iterations at : 13:55:19 total time elapsed so far: 8.880000 # loops : 120 mean,max change (kT/e) : 6.6207745E-04 2.3109436E-02 energy calculations done at 8.880000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 9.020000 DelPhi exited at 13:55:19 9.019u 0.279s 0:09.58 96.8% 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 15 2011 at 13:55:19 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: 267 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 2668 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.731924 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): 10.11200 60.97100 ymin,ymax (A): 9.696000 59.24400 zmin,zma (A): 8.711000 67.35100 x,y,z range (A): 50.85900 49.54800 58.64000 scale (grids/A): 2.731924 object centre (A): 35.54150 34.47000 38.03100 number of atom coordinates read : 2668 total number of charged atoms : 2499 net assigned charge : 2.500012 assigned positive charge : 391.7929 centred at (gu) : 92.78243 85.26565 85.20399 assigned negative charge : -389.2935 centred at (gu) : 92.37486 85.58453 85.56934 start vw surface at 9.0000004E-02 fill in re-entrant regions at 0.3400000 boundary points on vw surface= 164340 reading accessible surface arcs data from file ARCDAT no. of arc points read = 24776 grid for indexing accessible points = 1.400000 m= 79346 mr = 115891 m= 27320 mr = 63158 m= 6010 mr = 21370 m= 737 mr = 4308 m= 44 mr = 430 m= 0 mr = 12 time to grow re-entrant surface = 0.5200000 no. cavity mid-points inaccessible to solvent = 675 time to turn everything in is 0.9500000 number of dielectric boundary points 72628 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 2 iepsmp to db, and charging done at 1.060000 number of grid points assigned charge 19989 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.821282 object centre at (A) : 35.54150 34.47000 38.03100 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 0.6725084 -0.2575409 midg,1,midg; midg,igrid,midg 0.1448670 6.0570359E-02 1,midg,midg; igrid,midg,midg -0.1219840 4.4367447E-02 gauss-seidel spectral radius is 0.9992363 estimated iterations to convergence 286 setup time was (sec) 1.240000 now iterating at: 13:55:21 rms-change max change #iterations 9.279852 292.4873 at 10iterations 5.356729 115.8314 at 20iterations 3.070955 66.95775 at 30iterations 1.755273 43.77981 at 40iterations 1.003899 20.00493 at 50iterations 0.5722887 11.07297 at 60iterations 0.3263088 6.384615 at 70iterations 0.1857241 3.675601 at 80iterations 0.1058481 1.968338 at 90iterations 6.0199238E-02 1.245815 at 100iterations 3.4260914E-02 0.7034721 at 110iterations 1.9504229E-02 0.3467579 at 120iterations 1.1084517E-02 0.2024031 at 130iterations 6.3015381E-03 0.1207838 at 140iterations 3.5748533E-03 6.9576740E-02 at 150iterations 2.0329489E-03 4.0085793E-02 at 160iterations 1.1582882E-03 2.2303581E-02 at 170iterations 6.5789738E-04 1.2555122E-02 at 180iterations finished qdiffx linear iterations at : 13:55:34 total time elapsed so far: 14.58000 # loops : 180 mean,max change (kT/e) : 6.5789738E-04 1.2555122E-02 energy calculations done at 14.58000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 14.72000 DelPhi exited at 13:55:35 14.719u 0.288s 0:15.45 97.0% 0+0k 0+0io 0pf+0w