abaqus子程序怎么提取主应力呀？ 100

在编写子程序的时候遇到了点问题…想用vsprinc来提取最大主应力的…但是得到的结果很怪 

dimension s（nblock, ndir+nshr), eigVal(nblock, 3) 

call vsprinc(nblock, s, eigVal, ndir, nshr) 

do k =1, nblock 

stateNew(k,2)=eigVal(k,1)

 刚开始接触子程序，还请各位大佬教教我叭，谢谢谢谢！