<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello All!<div><br></div><div>I have found a workaround for the problem and since similar thresholds may occur in other parts of Molpro as well I'd just like to post the solution. Maybe it's helpful for somebody stumbling over the question in the archive later on.</div><div><br></div><div>Two limits apply to PROPERTY:</div><div>- Per invocation in an execution block only 30 points can be calculated ("Too many operators" in propa.f thereafter)</div><div>- For each job at most 148 points can be calculated in all combined PROPERTY invocations ("Too many operators" in read_op.f afterwards) </div><div><br></div><div>The trick is to calculate the wavefunction and integrals, save the restart data (.wfu, .int & .aux) and calculate all needed points in portions of at most 148 using always a *fresh* copy of the scratch files. That is, copy the three files to new names, run a "restarted" job with e.g. 5 PROPERTY blocks with 29 points each, then delete the scratch files and create a new copy for the next 5 * 29 points. When restarting the job the files have to be in the same directory, copying the scratch files to a sub folder is not enough ("Recomputing integrals since basis changed" and "Illegal geometry input - no atoms").</div><div><br></div><div>Example input follows below. Hope this helps anyone!</div><div><br></div><div><br></div><div>Best regards,</div><div><br></div><div>Ben</div><div><br></div><div><br></div><div><br></div><div>=================================================<br>Example Input for the Property Calculations<br>=================================================<br><br></div><div><div>***, ELECPOT</div><div>file, 1, scratch.copy.int</div><div>file, 2, scratch.copy.wfu</div><div>file, 3, scratch.copy.aux</div><div><br></div><div>{matrop</div><div>read, savmat, DEN, charge, file=density.mat</div><div>save, savmat, 2140.2, DENSITY}</div><div><br></div><div>{property</div><div>density,2140.2</div><div>pot,, -5.66918, -7.55890, -4.72431</div><div>pot,, -5.66918, -7.55890, -3.77945</div><div>pot,, -5.66918, -7.55890, -2.83459</div><div> [ etc., etc., 23 other points here ]</div><div>pot,, -5.66918, -5.66918, 5.66918</div><div>pot,, -5.66918, -5.66918, 6.61404</div><div>pot,, -5.66918, -4.72431, -6.61404}</div><div><br></div><div>[etc, 5 blocks like the above]</div></div><div><br></div><div>=================================================</div><div><br></div><div><br></div><div><br><div><div>On 10 Sep 2008, at 10:52, Benjamin Bulheller wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div 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><br></div></span></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></div></blockquote></div><br></div></body></html>