DMAP修改sol103输出特征向量的问题？？？求助，急急急

big-dish

各位大侠，我需要修改sol103的求解序列，使得特征向量输出到一个文件中，
比如：
ASSIGN OUTPUT4='ssplate2_eg.dat', UNIT=11 , FORMATTED delete
SOL 103
DIAG 8,15
compile xread $
alter 'IF ( READAPP'(,-1) $
OUTPUT4 phix//-1/11/0/ $
我用的是这个alter语句，但是得到的结果中的特征向量，不包含零的项（DMAP中要求的是输出第一个非零元素到最后一个非零元素之间所有的元素，其中包含第一个非零元素和最后一个非零元素之间的零元素，即上面的输出结果中不包含第一个非零元素与最后一个非零元素间的零元素）。我需要的是输出整个自由度中的特征向量（包含零元素的特征向量），课题很紧，真心的希望见到帖子的大侠伸出您的援助之手，小弟在这里先谢您了！！！！！！