Marc有限元子程序二次开发

我有一个Marc有限元模型编写了一个子程序，但是在运行过程中子程序报错。子程序是一个蠕变本构方程，我的模型里包含三个不同区域，每个区域对应着不同的本构方程参数，请通过编写修改我的子程序实现