Wed Dec 28 11:02:39 PST 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 Wed Dec 28 2011 at 11:02:39 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: 270 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 1966 !!! WARNING: ASH 25 has a net charge of -0.3190 !!! WARNING: ASH 125 has a net charge of -0.3190 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 5.9999999E-02 grid size : 179 scale,in grids/A, set to be: 0.6491494 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): -6.225000 34.29100 ymin,ymax (A): -3.427000 49.34200 zmin,zma (A): -23.96500 30.87600 x,y,z range (A): 40.51600 52.76900 54.84100 scale (grids/A): 0.6491494 object centre (A): 14.03300 22.95750 3.455500 number of atom coordinates read : 1966 total number of charged atoms : 1832 net assigned charge : 3.362005 assigned positive charge : 291.8019 centred at (gu) : 90.91340 92.20753 90.17400 assigned negative charge : -288.4402 centred at (gu) : 91.05152 92.36160 90.10315 start vw surface at 0.1400000 fill in re-entrant regions at 0.1500000 boundary points on vw surface= 5688 nv = 96 ne = 84 npr = 40226 no. pairs analyzed = 40226 no. exposed pairs = 2462 no. arc points = 17472 no. surface atoms = 835 nbur = 1131 mkacc time = 0.1300000 writing accessible surface arcs data to file ARCDAT grid for indexing accessible points = 1.400000 m= 276 mr = 3139 m= 1 mr = 37 m= 0 mr = 0 time to grow re-entrant surface = 0.1600000 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.3200000 number of dielectric boundary points 2789 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 1757 ## 172 charges are in solution ## iepsmp to db, and charging done at 0.4700000 number of grid points assigned charge 6598 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid -9.1348112E-10 1.5177616E-08 midg,1,midg; midg,igrid,midg 1.4846122E-08 -1.7038721E-08 1,midg,midg; igrid,midg,midg 1.8038621E-08 -1.0486286E-08 gauss-seidel spectral radius is 0.9876476 estimated iterations to convergence 74 setup time was (sec) 0.6400000 now iterating at: 11:02:40 rms-change max change #iterations 0.1354744 25.93972 at 10iterations 1.3561664E-02 2.334452 at 20iterations 1.4177690E-03 0.2629519 at 30iterations 1.4292920E-04 1.9981384E-02 at 40iterations 1.4344668E-05 2.5787354E-03 at 50iterations finished qdiffx linear iterations at : 11:02:44 total time elapsed so far: 4.210000 # loops : 50 mean,max change (kT/e) : 1.4344668E-05 2.5787354E-03 energy calculations done at 4.210000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 4.430000 DelPhi exited at 11:02:45 4.434u 0.475s 0:05.52 88.7% 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 Wed Dec 28 2011 at 11:02: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: 270 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 1966 !!! WARNING: ASH 25 has a net charge of -0.3190 !!! WARNING: ASH 125 has a net charge of -0.3190 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 5.9999999E-02 grid size : 179 scale,in grids/A, set to be: 1.947448 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): -6.225000 34.29100 ymin,ymax (A): -3.427000 49.34200 zmin,zma (A): -23.96500 30.87600 x,y,z range (A): 40.51600 52.76900 54.84100 scale (grids/A): 1.947448 object centre (A): 14.03300 22.95750 3.455500 number of atom coordinates read : 1966 total number of charged atoms : 1832 net assigned charge : 3.362005 assigned positive charge : 291.8019 centred at (gu) : 92.74001 96.62265 90.52214 assigned negative charge : -288.4402 centred at (gu) : 93.15468 97.08490 90.30968 start vw surface at 0.1400000 fill in re-entrant regions at 0.2700000 boundary points on vw surface= 60993 reading accessible surface arcs data from file ARCDAT no. of arc points read = 17472 grid for indexing accessible points = 1.400000 m= 22272 mr = 41925 m= 3530 mr = 17013 m= 201 mr = 2376 m= 1 mr = 90 m= 0 mr = 0 time to grow re-entrant surface = 0.3100000 no. cavity mid-points inaccessible to solvent = 2900 time to turn everything in is 0.6400000 number of dielectric boundary points 25593 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 140 iepsmp to db, and charging done at 0.7600000 number of grid points assigned charge 14454 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.6491494 object centre at (A) : 14.03300 22.95750 3.455500 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 4.3013226E-02 1.7676752E-02 midg,1,midg; midg,igrid,midg 3.5325345E-02 -1.3049771E-02 1,midg,midg; igrid,midg,midg 9.0309111E-03 -4.3388968E-03 gauss-seidel spectral radius is 0.9986170 estimated iterations to convergence 213 setup time was (sec) 0.9500000 now iterating at: 11:02:46 rms-change max change #iterations 3.782812 203.4365 at 10iterations 1.762953 64.17108 at 20iterations 0.8207853 30.27959 at 30iterations 0.3834387 13.26419 at 40iterations 0.1770710 6.412693 at 50iterations 8.2566880E-02 2.854403 at 60iterations 3.8634006E-02 1.360260 at 70iterations 1.7952323E-02 0.6468263 at 80iterations 8.3484799E-03 0.2774391 at 90iterations 3.8937468E-03 0.1414127 at 100iterations 1.8085164E-03 6.1851501E-02 at 110iterations 8.4545324E-04 2.9037476E-02 at 120iterations finished qdiffx linear iterations at : 11:02:55 total time elapsed so far: 9.930000 # loops : 120 mean,max change (kT/e) : 8.4545324E-04 2.9037476E-02 energy calculations done at 9.930000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 10.15000 DelPhi exited at 11:02:56 10.154u 0.473s 0:11.13 95.4% 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 Wed Dec 28 2011 at 11:02:56 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: 270 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 1966 !!! WARNING: ASH 25 has a net charge of -0.3190 !!! WARNING: ASH 125 has a net charge of -0.3190 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 5.9999999E-02 grid size : 179 scale,in grids/A, set to be: 2.921172 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): -6.225000 34.29100 ymin,ymax (A): -3.427000 49.34200 zmin,zma (A): -23.96500 30.87600 x,y,z range (A): 40.51600 52.76900 54.84100 scale (grids/A): 2.921172 object centre (A): 14.03300 22.95750 3.455500 number of atom coordinates read : 1966 total number of charged atoms : 1832 net assigned charge : 3.362005 assigned positive charge : 291.8019 centred at (gu) : 94.11012 99.93410 90.78301 assigned negative charge : -288.4402 centred at (gu) : 94.73217 100.6274 90.46455 start vw surface at 0.1400000 fill in re-entrant regions at 0.5500000 boundary points on vw surface= 140737 reading accessible surface arcs data from file ARCDAT no. of arc points read = 17472 grid for indexing accessible points = 1.400000 m= 71608 mr = 101255 m= 26772 mr = 59516 m= 6328 mr = 22183 m= 884 mr = 5023 m= 38 mr = 676 m= 0 mr = 20 time to grow re-entrant surface = 0.9600000 no. cavity mid-points inaccessible to solvent = 9424 time to turn everything in is 1.620000 number of dielectric boundary points 57694 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 0 iepsmp to db, and charging done at 1.760000 number of grid points assigned charge 14656 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.947448 object centre at (A) : 14.03300 22.95750 3.455500 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 1.011034 3.6932021E-02 midg,1,midg; midg,igrid,midg 0.3316632 -0.1567922 1,midg,midg; igrid,midg,midg 5.6536853E-02 -1.6238051E-02 gauss-seidel spectral radius is 0.9992406 estimated iterations to convergence 286 setup time was (sec) 1.960000 now iterating at: 11:02:58 rms-change max change #iterations 8.774607 458.4994 at 10iterations 5.012833 127.8131 at 20iterations 2.856112 63.56480 at 30iterations 1.628091 34.51715 at 40iterations 0.9289348 19.10796 at 50iterations 0.5295575 10.06346 at 60iterations 0.3014202 6.373438 at 70iterations 0.1708456 3.449316 at 80iterations 9.7425140E-02 1.948274 at 90iterations 5.5588592E-02 1.058990 at 100iterations 3.1559926E-02 0.5998374 at 110iterations 1.7961854E-02 0.3401794 at 120iterations 1.0211118E-02 0.2362785 at 130iterations 5.8067278E-03 0.1348286 at 140iterations 3.3090028E-03 6.2961578E-02 at 150iterations 1.8831195E-03 3.8929939E-02 at 160iterations 1.0741785E-03 2.1487921E-02 at 170iterations 6.1061449E-04 1.2019634E-02 at 180iterations finished qdiffx linear iterations at : 11:03:13 total time elapsed so far: 16.77000 # loops : 180 mean,max change (kT/e) : 6.1061449E-04 1.2019634E-02 energy calculations done at 16.77000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 16.98000 DelPhi exited at 11:03:14 16.986u 0.523s 0:18.02 97.1% 0+0k 0+0io 0pf+0w