!让MSC.MARC程序批量执行
  !by 陆新征,清华大学
subroutine Run_Program(EQNo) ! 进行地震运算
  use DFLIB
  integer EQNo
  character*60 :: RunCommand ! 程序运行命令
 RunCommand(1: 29)="run_marc -b no -j in.dat -u THUFIBER"
  RunCommand(30:60)=THUFIBER
  write(*,*) RunCommand
i= systemqq(RunCommand)
return
  end subroutine