<div dir="ltr">Dear cp2krs<br><br>I am optimising some molecules using GEO_OPT.<br><br>is it possible read a geometry from a restart file? and the cell parameters?<br>I am trying but I cannot :-(<br><br>concerning the geometry this is what I put in the input:<br><br>    &TOPOLOGY<br>!!    COORD_FILE_NAME umsm-nh2.xyz<br>!!    COORD_FILE_FORMAT XYZ<br>      COORD_FILE_NAME umsm-nh2-1.restart<br>      COORD_FILE_FORMAT CP2K<br>    &END<br><br>and this is the error I had:<br><br> *******************************************************************************<br> *   ___                                                                       *<br> *  /   \                                                                      *<br> * [ABORT]                                                                     *<br> *  \___/    A floating point type object was expected, found end of the line, *<br> *    |         File: 'umsm-nh2-1.restart', Line: 1, Column: 1, Chunk: < >     *<br> *  O/|                                                                        *<br> * /| |                                                                        *<br> * / \                                          input/cp_parser_methods.F:1194 *<br> *******************************************************************************<br><br><br>Concerning the cell parameters it is even worse:<br><br>Input:<br>    &CELL<br>!!    ABC [angstrom] 41.868 41.880 17.564<br>!!    ALPHA_BETA_GAMMA 90. 90. 120.<br>    CELL_FILE_FORMAT CP2K<br>    CELL_FILE_NAME umsm-nh2-1.restart<br>      SYMMETRY NONE<br>    &END CELL<br><br>error:<br>forrtl: severe (59): list-directed I/O syntax error, unit -5, file Internal List-Directed Read<br><br>However in the restart file I can see the information are present:<br><br>[tiana@deneb1 tmp]$ grep -A3  COORD umsm-nh2-1.restart<br><br>     &COORD<br>Zn    7.0975619283291929E+00    1.7345103848396320E+01    1.3189045443483327E+01<br>Zn    1.3818777106160974E+01    1.8926557862281943E+01    4.4217589353253608E+00<br>Zn    2.3603220884596876E+00    3.3736335046100855E+01    1.3189876918006378E+01<br>--<br>     &END COORD<br><br>and <br><br>[tiana@deneb1 tmp]$ grep CELL umsm-nh2-1.restart <br>   RUN_TYPE  CELL_OPT<br>   &CELL_OPT<br>     TYPE  DIRECT_CELL_OPT<br>   &END CELL_OPT<br>     &CELL<br>       MULTIPLE_UNIT_CELL  1 1 1<br>     &END CELL<br><br>Any idea of/how can I parse that file?<br><br>Cheers<br>Davide<br></div>