How to use source.f to set up I-131?

Dear @fengfengxu,
You don’t need a source routine to model an isotope source, but you can (and should) use Fluka built-in functionalities.
Have a look at a similar question (albeit for different isotopes): How to define a isotope radioactive γ source?