<div dir="ltr">Dear Molpro Users,<br><br>I would like to do Diatomic potential curve analysis. I have gathered the bond distance and the corresponding energies in a data file from a previously calculated PEC. How could I direct Molpro to read these data and perform Diatomic potential curve analysis. I am doing something like this:<br><br>Sample of my input file: <a href="http://input.in" target="_blank">input.in</a><br><br>  Â  Â  ANG<br>  Â  Â  readvar,data.dat<br>  Â  Â  table,distance, energy<br>  Â  Â  diatomic, Mass=14<br><br>My data file: data.dat<br><br>  Â  Â  Â distance Â  Â  energy<br>  Â  Â  Â  1.0 Â  Â  Â  Â  Â  Â -10.00<br>  Â  Â  Â  2.0 Â  Â  Â  Â  Â  Â -20.00<br><br>If run the input file, I get result such as:<br><br><font size="1">DISTANCE Â  ENERGY </font><div><font size="1">DISTANCE Â  ENERGY<br><br> Polynomial fit for diatomic using reduced mass Â 14.00000 u<br><br> CALCULATED SPECTROSCOPIC CONSTANTS FOR Â NPOINTS= 1 DEGREE= 0 Â RMIN= 0.00 Â RMAX= 0.00<br><br> =============================<wbr>==============================<wbr>==============================<wbr>===============<br><br> TITLE Â  Â  Â  Â  ENRE Â  Â  RE(AU) Â  RE(A) Â  BE Â  Â  Â  AE Â  Â  Â WE Â  Â  WEXE Â  WEYE Â  DE Â  Â D0 Â  Â  Â SIG<br><br> =============================<wbr>==============================<wbr>==============================<wbr>===============<br><br> ENERGY Â  Â  Â  Â 0.0000000 Â  Â  NaN Â  Â  NaN Â  Â  NaN Â  Â  NaN Â  Â  Â NaN Â  Â NaN Â  NaN Â 0.00 Â 0.00 Â 0.0000000E+00<br><br> =============================<wbr>==============================<wbr>==============================<wbr>===============</font><div><font size="1"><br></font></div><div><br></div><div>The program does not the read the data.</div><div><br></div><div>Thanks in advance.</div><div><br></div><div>Somnath Bhowmick</div></div></div>