EMAP reading

Posted by: pmj

EMAP reading - 05/06/19 07:50 AM

Dear all,

can anyone help me with emap reading? I have segmentation fault error as before the emap I am reading a lot of files so when I want to use the following script the emap is opened as the unit number used previously.

define fit sele .byres. ALL .and. .not. HYDROGEN end

set map "/home//test/test_1/EMmap/emd_8149.map"

emap read emd name @map
emap refe emd sele fit end
emap assi emd as RIG0 sele fit end
emap cons 0.5 RIG0

I am aware I can open the emap as a unit, but I am not sure how to write the following lines in that case:

emap refe emd sele fit end
emap assi emd as RIG0 sele fit end
emap cons 0.5 RIG0

Thank you!