do sp =41,61,2 震源位置のループ作成

出力ファイル名
open(33,file=trim(adjustl(filename))//“filename2.d”)