Sat Oct 15 13:20:18 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+sph2.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 = amb2.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:20:18 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 = 3164 !!! WARNING: SEP 138 has a net charge of -0.4750 !!! WARNING: GLU 255 has a net charge of -0.2720 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.5611690 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): 8.805000 59.40600 ymin,ymax (A): -31.45600 31.98300 zmin,zma (A): 44.66800 95.74500 x,y,z range (A): 50.60100 63.43900 51.07700 scale (grids/A): 0.5611690 object centre (A): 34.10550 0.2635002 70.20650 number of atom coordinates read : 3164 total number of charged atoms : 3017 net assigned charge : -0.2469896 assigned positive charge : 478.1490 centred at (gu) : 89.79521 90.57281 91.22119 assigned negative charge : -478.3972 centred at (gu) : 89.85290 90.56130 91.10155 start vw surface at 0.1600000 fill in re-entrant regions at 0.1700000 boundary points on vw surface= 6467 nv = 96 ne = 84 npr = 62542 no. pairs analyzed = 62542 no. exposed pairs = 4132 no. arc points = 29281 no. surface atoms = 1377 nbur = 1787 mkacc time = 0.1400000 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.540000 m= 257 mr = 3260 m= 0 mr = 10 time to grow re-entrant surface = 0.1600000 no. cavity mid-points inaccessible to solvent = 35 time to turn everything in is 0.3300000 number of dielectric boundary points 3454 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 2958 ## 382 charges are in solution ## iepsmp to db, and charging done at 0.4300000 number of grid points assigned charge 8467 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid -1.1558139E-09 1.9139961E-09 midg,1,midg; midg,igrid,midg -1.6733764E-10 1.0788549E-10 1,midg,midg; igrid,midg,midg 6.9613121E-10 -7.3280010E-10 gauss-seidel spectral radius is 0.9836388 estimated iterations to convergence 64 setup time was (sec) 0.5600000 now iterating at: 13:20:18 rms-change max change #iterations 8.2942508E-02 15.45376 at 10iterations 5.9694005E-03 0.9255066 at 20iterations 4.4445836E-04 7.1061134E-02 at 30iterations 4.0159572E-05 6.6127777E-03 at 40iterations 6.2677359E-06 9.9277496E-04 at 50iterations finished qdiffx linear iterations at : 13:20:21 total time elapsed so far: 3.260000 # loops : 50 mean,max change (kT/e) : 6.2677359E-06 9.9277496E-04 energy calculations done at 3.260000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 3.440000 DelPhi exited at 13:20:22 3.439u 0.205s 0:04.01 90.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 Sat Oct 15 2011 at 13:20:22 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 = 3164 !!! WARNING: SEP 138 has a net charge of -0.4750 !!! WARNING: GLU 255 has a net charge of -0.2720 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 9.0000004E-02 grid size : 179 scale,in grids/A, set to be: 1.683507 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): 8.805000 59.40600 ymin,ymax (A): -31.45600 31.98300 zmin,zma (A): 44.66800 95.74500 x,y,z range (A): 50.60100 63.43900 51.07700 scale (grids/A): 1.683507 object centre (A): 34.10550 0.2635002 70.20650 number of atom coordinates read : 3164 total number of charged atoms : 3017 net assigned charge : -0.2469896 assigned positive charge : 478.1490 centred at (gu) : 89.38486 91.71758 93.66264 assigned negative charge : -478.3972 centred at (gu) : 89.55821 91.68374 93.30433 start vw surface at 0.1600000 fill in re-entrant regions at 0.2500000 boundary points on vw surface= 71579 reading accessible surface arcs data from file ARCDAT no. of arc points read = 29281 grid for indexing accessible points = 1.400000 m= 21941 mr = 47652 m= 2264 mr = 14898 m= 75 mr = 1092 m= 0 mr = 7 time to grow re-entrant surface = 0.2200000 no. cavity mid-points inaccessible to solvent = 1461 time to turn everything in is 0.5100000 number of dielectric boundary points 32210 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 290 ## 2 charges are in solution ## iepsmp to db, and charging done at 0.6500000 number of grid points assigned charge 23199 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.5611690 object centre at (A) : 34.10550 0.2635002 70.20650 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -3.0616417E-03 3.7847026E-03 midg,1,midg; midg,igrid,midg 2.1576687E-05 3.2121958E-03 1,midg,midg; igrid,midg,midg 8.0666477E-03 4.8168171E-03 gauss-seidel spectral radius is 0.9983293 estimated iterations to convergence 194 setup time was (sec) 0.8300000 now iterating at: 13:20:23 rms-change max change #iterations 3.668676 157.5415 at 10iterations 1.587159 62.23697 at 20iterations 0.6791557 22.64323 at 30iterations 0.2943295 10.14026 at 40iterations 0.1268899 4.566849 at 50iterations 5.4676760E-02 1.987932 at 60iterations 2.3623582E-02 0.8010674 at 70iterations 1.0188143E-02 0.3577671 at 80iterations 4.4052284E-03 0.1387959 at 90iterations 1.9073226E-03 6.8868637E-02 at 100iterations 8.2683121E-04 3.1081200E-02 at 110iterations 3.5710831E-04 1.1873245E-02 at 120iterations finished qdiffx linear iterations at : 13:20:30 total time elapsed so far: 7.610000 # loops : 120 mean,max change (kT/e) : 3.5710831E-04 1.1873245E-02 energy calculations done at 7.610000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 7.790000 DelPhi exited at 13:20:30 7.788u 0.220s 0:08.33 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 15 2011 at 13:20:30 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 = 3164 !!! WARNING: SEP 138 has a net charge of -0.4750 !!! WARNING: GLU 255 has a net charge of -0.2720 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: 2.525260 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): 8.805000 59.40600 ymin,ymax (A): -31.45600 31.98300 zmin,zma (A): 44.66800 95.74500 x,y,z range (A): 50.60100 63.43900 51.07700 scale (grids/A): 2.525260 object centre (A): 34.10550 0.2635002 70.20650 number of atom coordinates read : 3164 total number of charged atoms : 3017 net assigned charge : -0.2469896 assigned positive charge : 478.1490 centred at (gu) : 89.07733 92.57631 95.49383 assigned negative charge : -478.3972 centred at (gu) : 89.33733 92.52547 94.95647 start vw surface at 0.1500000 fill in re-entrant regions at 0.4200000 boundary points on vw surface= 166359 reading accessible surface arcs data from file ARCDAT no. of arc points read = 29281 grid for indexing accessible points = 1.400000 m= 75326 mr = 116759 m= 21617 mr = 59342 m= 3528 mr = 15972 m= 282 mr = 2210 m= 1 mr = 158 m= 0 mr = 0 time to grow re-entrant surface = 0.6800001 no. cavity mid-points inaccessible to solvent = 4703 time to turn everything in is 1.220000 number of dielectric boundary points 72672 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 13 iepsmp to db, and charging done at 1.370000 number of grid points assigned charge 24121 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.683507 object centre at (A) : 34.10550 0.2635002 70.20650 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid -5.7390952E-03 7.4873582E-02 midg,1,midg; midg,igrid,midg -0.1262636 0.1425445 1,midg,midg; igrid,midg,midg 0.1185832 9.3635127E-02 gauss-seidel spectral radius is 0.9991949 estimated iterations to convergence 278 setup time was (sec) 1.550000 now iterating at: 13:20:32 rms-change max change #iterations 9.477376 311.5768 at 10iterations 5.478401 135.2076 at 20iterations 3.145682 70.24051 at 30iterations 1.797041 48.55580 at 40iterations 1.023822 21.53813 at 50iterations 0.5807471 12.20033 at 60iterations 0.3289874 6.908325 at 70iterations 0.1860630 4.527874 at 80iterations 0.1048465 2.280426 at 90iterations 5.9223134E-02 1.329385 at 100iterations 3.3340767E-02 0.6709633 at 110iterations 1.8787244E-02 0.3887663 at 120iterations 1.0542946E-02 0.2176056 at 130iterations 5.9326557E-03 0.1213474 at 140iterations 3.3323013E-03 6.9713295E-02 at 150iterations 1.8737363E-03 4.1180044E-02 at 160iterations 1.0511278E-03 2.1319866E-02 at 170iterations 5.9099263E-04 1.2587130E-02 at 180iterations finished qdiffx linear iterations at : 13:20:43 total time elapsed so far: 12.48000 # loops : 180 mean,max change (kT/e) : 5.9099263E-04 1.2587130E-02 energy calculations done at 12.48000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 12.66000 DelPhi exited at 13:20:43 12.661u 0.241s 0:13.29 97.0% 0+0k 0+0io 0pf+0w