通常の結晶構造であれば、データベースか何かから引っ張って来れるが、それをベースにした超周期構造を作ろうと思うと、対称性（空間群）が変わるため、手でやるには歯が立たない。 そこで、入力した原子座標に対して対称性を見つけてくれる spglib を使って…

何も考えずにプログラムを作っていたら、配列の次元が7個を超えるとコンパイルできないらしい。 program test_array7 real( kind( 0d0 ) ) :: r( 1, 2, 3, 4, 5, 6, 7, 8 ) end program test_array7 これを gfortran でコンパイルすると test_array7.f90:3:4…