Marc中半椭球形热源、双椭球型热源子程

希望能得到一些在Marc中焊接热源的子程序，尤其是半椭球形热源和双椭球型热源，用Fortran编辑语言。