<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello All!<div><br></div><div>Regarding my problem below I have discovered a new threshold in Molpro which puzzles me.</div><div><br></div><div>I can only calculate the electrostatic potential at 30 points per calculation before it hits the maximum value 30 of the MXOPER variable in propa.f. I then determined and saved the integrals in one job and ran individual jobs re-using these data with an evaluation of 30 points per run.</div><div><br></div><div>This works for 148 points (as I found before) and then it hits the maximum value of the variable MAXOPER (note the additional A, it is actually set to 150 what is reached during a loop). This happens in readop.f (routine write_op) and Molpro dies with the same error as before, "Too many operators".</div><div><br></div><div>Thanks a lot for any ideas or hints!</div><div><br></div><div>Cheers,</div><div><br></div><div>Ben</div><div><br></div><div><span class="Apple-style-span" style="font-family: Helvetica; "><div><br class="khtml-block-placeholder"></div><div>======================================</div><div><br class="khtml-block-placeholder"></div><div>Benjamin Bulheller</div><div><br class="khtml-block-placeholder"></div><div>Research Group Prof. J. D. Hirst</div><div><a href="http://comp.chem.nottingham.ac.uk">http://comp.chem.nottingham.ac.uk</a></div><div>School of Physical Chemistry</div><div>University of Nottingham</div><div><br class="khtml-block-placeholder"></div><div>======================================</div><div><br class="khtml-block-placeholder"></div></span></div><div><br></div><div><div><br><div>Begin forwarded message:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" color="#000000" style="font: 12.0px Helvetica; color: #000000"><b>From: </b></font><font face="Helvetica" size="3" style="font: 12.0px Helvetica">Benjamin Bulheller <<a href="mailto:pcxbb2@nottingham.ac.uk">pcxbb2@nottingham.ac.uk</a>></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" color="#000000" style="font: 12.0px Helvetica; color: #000000"><b>Date: </b></font><font face="Helvetica" size="3" style="font: 12.0px Helvetica">3 September 2008 17:44:06 BST</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" color="#000000" style="font: 12.0px Helvetica; color: #000000"><b>To: </b></font><font face="Helvetica" size="3" style="font: 12.0px Helvetica"><a href="mailto:molpro-user@molpro.net">molpro-user@molpro.net</a></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" color="#000000" style="font: 12.0px Helvetica; color: #000000"><b>Subject: </b></font><font face="Helvetica" size="3" style="font: 12.0px Helvetica"><b>Calculation of Potential</b></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; min-height: 14px; "><br></div> </div><div>Hi!<br><br>I need to calculate the electrostatic potential of a compound at several thousand points around the molecule. This has been done before (with Molpro 2002.6) using the following code:<br><br> property<br> density,2100.2<br> pot,,-12.28322, -2.83459, -0.94486<br> pot,,-12.24522, -0.94486, 1.88973<br> pot,,-12.67322, 1.88973, 0.94486<br> pot,,-11.33836, -2.83459, -3.77945<br> pot,,-10.12736, 2.83459, 2.83459<br> etc...<br><br>In 2002.6 this worked with up to 2500 "pot" statements while in 2006.1 only 30 points are possible and I get a "Too many operators" error with more than that. This happens in propa.f caused by the MXOPER variable.<br><br>Using a slightly crude workaround I can wrap the pot statement in nested do loops and for some reason Molpro will then compute 148 points and fail thereafter.<br><br>Am I doing something wrong there in general? What puzzles me is that it obviously must have worked in 2002.6 but stopped working in 2006.1 (without changes in the concerned code section), is there any explanation for that?<br><br>Thank you very much in advance,<br><br>Ben<br><br><br><br>=================================================<br>Example Input for H2<br>=================================================<br><br>print,basis<br>gprint,orbitals<br>geomtyp=xyz<br><br>geometry={<br>nosym<br>2<br>Title<br>H 0.0 0.0 0.0<br>H 0.0 0.0 1.0<br>}<br><br>basis={<br> ! HYDROGEN [7s,3p] -> (2s,1p) ANO-S<br> s,H ,82.636374,12.409558,2.8238540,0.79767000,0.25805300,0.08989100,0.03146200<br> c,1.7, 0.0025641, 0.0192126, 0.0950512, 0.3018678, 0.4874459, 0.23547943, 0.01577776<br> c,1.7,-0.0046008,-0.0324447,-0.1710737,-0.5048672,-0.0948920, 0.37085551, 0.65565134<br> p,H ,1.6625000,0.41560000,0.10390000<br> c,1.3, 0.3012771, 0.6800574, 0.1985627<br>}<br><br>int;<br>hf;<br>property<br>density,2100.2<br>dm,,0.0,0.0,0.0;<br><br>property<br>density,2100.2<br>pot,,-12.28322, -2.83459, -0.94486<br>pot,,-12.24522, -0.94486, 1.88973<br>pot,,-12.67322, 1.88973, 0.94486<br>pot,,-11.33836, -2.83459, -3.77945<br>pot,,-10.12736, 2.83459, 2.83459<br>pot,, -9.33836, 2.83459, 3.77945<br>pot,, -8.25736, 3.77945, 1.88973<br>! plus 1000 and more additional lines for other points<br><br><br>=================================================<br></div></blockquote></div><br></div><br><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><div><br></div></span> </div><br></body></html>