Sun Oct 9 08:21: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 Sun Oct 9 2011 at 08:21: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: 236 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 5761 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 0.1500000 grid size : 179 scale,in grids/A, set to be: 0.5179086 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): 14.92700 81.19800 ymin,ymax (A): 17.88800 84.45100 zmin,zma (A): 0.2020000 68.93999 x,y,z range (A): 66.27100 66.56300 68.73799 scale (grids/A): 0.5179086 object centre (A): 48.06250 51.16950 34.57100 number of atom coordinates read : 5761 total number of charged atoms : 5595 net assigned charge : 2.000055 assigned positive charge : 868.2553 centred at (gu) : 90.54848 89.34621 90.05988 assigned negative charge : -866.2629 centred at (gu) : 90.51460 89.39357 90.24268 start vw surface at 0.2200000 fill in re-entrant regions at 0.2400000 boundary points on vw surface= 9869 nv = 96 ne = 84 npr = 115127 no. pairs analyzed = 115127 no. exposed pairs = 7428 no. arc points = 50681 no. surface atoms = 2480 nbur = 3281 mkacc time = 0.3300000 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.652000 m= 409 mr = 5100 m= 2 mr = 13 m= 0 mr = 0 time to grow re-entrant surface = 0.3700000 no. cavity mid-points inaccessible to solvent = 0 time to turn everything in is 0.6200000 number of dielectric boundary points 5167 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 5084 ## 736 charges are in solution ## iepsmp to db, and charging done at 0.7500000 number of grid points assigned charge 13255 setting boundary conditions some initial phi values: midg,midg,1; midg,midg,igrid 7.8967988E-10 -7.7620815E-10 midg,1,midg; midg,igrid,midg 2.8456773E-10 -1.3952599E-10 1,midg,midg; igrid,midg,midg -9.3005638E-11 2.1074523E-10 gauss-seidel spectral radius is 0.9810045 estimated iterations to convergence 60 setup time was (sec) 0.9000000 now iterating at: 08:21:07 rms-change max change #iterations 7.9945609E-02 11.37967 at 10iterations 4.9385219E-03 0.7020476 at 20iterations 5.7664176E-04 0.1183071 at 30iterations 1.6427881E-04 3.5169601E-02 at 40iterations 5.4464061E-05 1.1186600E-02 at 50iterations finished qdiffx linear iterations at : 08:21:10 total time elapsed so far: 4.200000 # loops : 50 mean,max change (kT/e) : 5.4464061E-05 1.1186600E-02 energy calculations done at 4.200000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 4.380000 DelPhi exited at 08:21:11 4.383u 0.280s 0:05.07 91.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 Sun Oct 9 2011 at 08:21: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: 236 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 5761 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 0.1500000 grid size : 179 scale,in grids/A, set to be: 1.553726 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): 14.92700 81.19800 ymin,ymax (A): 17.88800 84.45100 zmin,zma (A): 0.2020000 68.93999 x,y,z range (A): 66.27100 66.56300 68.73799 scale (grids/A): 1.553726 object centre (A): 48.06250 51.16950 34.57100 number of atom coordinates read : 5761 total number of charged atoms : 5595 net assigned charge : 2.000055 assigned positive charge : 868.2553 centred at (gu) : 91.64376 88.03696 90.17796 assigned negative charge : -866.2629 centred at (gu) : 91.54389 88.18044 90.72791 start vw surface at 0.2600000 fill in re-entrant regions at 0.4400000 boundary points on vw surface= 109806 reading accessible surface arcs data from file ARCDAT no. of arc points read = 50681 initial cube size too small, in assigning accessible points to a grid therefore rescaling... grid for indexing accessible points = 1.540000 m= 34312 mr = 74824 m= 3257 mr = 22952 m= 75 mr = 1406 m= 0 mr = 3 time to grow re-entrant surface = 0.5000000 no. cavity mid-points inaccessible to solvent = 787 time to turn everything in is 1.000000 number of dielectric boundary points 48265 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 755 ## 8 charges are in solution ## iepsmp to db, and charging done at 1.150000 number of grid points assigned charge 42261 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.5179086 object centre at (A) : 48.06250 51.16950 34.57100 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 3.6763784E-02 -4.2107172E-02 midg,1,midg; midg,igrid,midg 9.2785375E-04 -1.7192632E-02 1,midg,midg; igrid,midg,midg 5.7883891E-03 2.6804775E-03 gauss-seidel spectral radius is 0.9982657 estimated iterations to convergence 191 setup time was (sec) 1.310000 now iterating at: 08:21:12 rms-change max change #iterations 4.359704 182.6132 at 10iterations 1.851086 70.21503 at 20iterations 0.7867420 22.62394 at 30iterations 0.3342115 11.59094 at 40iterations 0.1422365 4.231629 at 50iterations 6.0656868E-02 1.691328 at 60iterations 2.5909940E-02 0.7839508 at 70iterations 1.1069249E-02 0.3197823 at 80iterations 4.7431118E-03 0.1361122 at 90iterations 2.0585831E-03 6.0970306E-02 at 100iterations 9.0896606E-04 3.0068159E-02 at 110iterations 4.1985780E-04 1.6763687E-02 at 120iterations finished qdiffx linear iterations at : 08:21:21 total time elapsed so far: 9.910000 # loops : 120 mean,max change (kT/e) : 4.1985780E-04 1.6763687E-02 energy calculations done at 9.910000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 10.09000 DelPhi exited at 08:21:22 10.089u 0.291s 0:10.64 97.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 Sun Oct 9 2011 at 08:21: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: 236 assigning charges and radii... opening formatted file:fort.13 number of atoms read in = 5761 atomic coordinates, charges and radii written to file fort.19 time to read in and/or assign rad/chrg= 0.1500000 grid size : 179 scale,in grids/A, set to be: 2.330589 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): 14.92700 81.19800 ymin,ymax (A): 17.88800 84.45100 zmin,zma (A): 0.2020000 68.93999 x,y,z range (A): 66.27100 66.56300 68.73799 scale (grids/A): 2.330589 object centre (A): 48.06250 51.16950 34.57100 number of atom coordinates read : 5761 total number of charged atoms : 5595 net assigned charge : 2.000055 assigned positive charge : 868.2553 centred at (gu) : 92.46527 87.05482 90.26653 assigned negative charge : -866.2629 centred at (gu) : 92.31570 87.27058 91.09174 start vw surface at 0.2400000 fill in re-entrant regions at 0.7700000 boundary points on vw surface= 254952 reading accessible surface arcs data from file ARCDAT no. of arc points read = 50681 initial cube size too small, in assigning accessible points to a grid therefore rescaling... grid for indexing accessible points = 1.540000 m= 119433 mr = 183580 m= 33210 mr = 92881 m= 4858 mr = 24264 m= 242 mr = 2973 m= 4 mr = 57 m= 0 mr = 0 time to grow re-entrant surface = 1.530000 no. cavity mid-points inaccessible to solvent = 2545 time to turn everything in is 2.460000 number of dielectric boundary points 108944 no. dielectric boundary points in salt = 0 no. grid points charged and at boundary= 71 iepsmp to db, and charging done at 2.640000 number of grid points assigned charge 44638 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.553726 object centre at (A) : 48.06250 51.16950 34.57100 pulling boundary values out of old potential map... some initial phi values: midg,midg,1; midg,midg,igrid 0.4265940 -0.9199491 midg,1,midg; midg,igrid,midg -2.8531158E-02 -0.3401276 1,midg,midg; igrid,midg,midg 0.1674294 0.1572611 gauss-seidel spectral radius is 0.9992376 estimated iterations to convergence 286 setup time was (sec) 2.820000 now iterating at: 08:21:25 rms-change max change #iterations 11.79216 435.2331 at 10iterations 6.996517 145.4960 at 20iterations 4.121662 92.46027 at 30iterations 2.408570 53.60278 at 40iterations 1.402158 28.53492 at 50iterations 0.8145837 15.98348 at 60iterations 0.4717790 8.238562 at 70iterations 0.2720054 5.003338 at 80iterations 0.1571301 2.625206 at 90iterations 9.0252310E-02 1.419411 at 100iterations 5.2000362E-02 0.8892365 at 110iterations 2.9849185E-02 0.5015879 at 120iterations 1.7187482E-02 0.3003845 at 130iterations 9.8496908E-03 0.1586285 at 140iterations 5.6654294E-03 9.4909668E-02 at 150iterations 3.2597717E-03 5.8046341E-02 at 160iterations 1.8770536E-03 3.3823013E-02 at 170iterations 1.0877399E-03 1.8741608E-02 at 180iterations finished qdiffx linear iterations at : 08:21:39 total time elapsed so far: 17.46000 # loops : 180 mean,max change (kT/e) : 1.0877399E-03 1.8741608E-02 energy calculations done at 17.46000 writing potential map in DELPHI format potential map written to file fort.14 total cpu time was (sec) 17.64000 DelPhi exited at 08:21:40 17.639u 0.312s 0:18.21 98.5% 0+0k 0+0io 0pf+0w