I’m trying to use magfld.f routine to set two region with different magnetic field, one Bx= 0.1T, another By=0.3T. But after compiling and creating executable, I get run error and I do not know why.
Here are my subroutine and flair file. Am I doing something wrong?