在編寫子程序的時候遇到了點問題…想用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)

 剛開始接觸子程序，還請各位大佬教教我叭，謝謝謝謝！