c c ----------------------------- c subroutine readsigmas real qq(3,3,50),oc(50) integer numsigs,i common /m2/ qq,oc,numsigs c open(4,file='misdat',status='old') do 50, i=1,8 read(4,*) read(4,*) read(4,*) c skip over the symmetry matrices! c 50 continue c do 100, i=1,50 write(*,*) ' sigma no. ', i read(4,*,end=200) ((qq(ii,jj,i),ii=1,3),jj=1,3),oc(i) c c NOTE the transposition of the input matrices c this is done so that when we form a Rodrigues vector from them c , if needed, the result is still in the fundamental zone! c but it means that one must (post-)multiply with the transpose c of the CSL matrix! c ADR vi 99 c 100 continue write(*,*) 'There may have been more than 50 Sigma matrices....' 200 numsigs=i-1 write(*,*) 'number of sigma values entered = ',numsigs return end c c ------------------------------------------------ c