Rachit Soni UNiMORE

Thank you for replying. To elaborate further, I will say I am not using Fepest. I have written a script , in which i am producing a stage variation time series for each rainfall event from 1 hour to 48 hours for a return period of 100 years. I am using IDF curves and rating curves for this. Then I am adding these to the model in my script and running the simulation. I am performing feflow simulation for 3 days. It's taking a lot of time and causing errors, as explained previously. In fact, I am only changing the boundary condition in each simulation. But still, it's very slow to perform these many simulations. What can i do? 

I am getting 'unspecified error' message when running multiple simulations. I have to change boundary conditions in every simulations and then delete them after extracting results in each simulation . I have to perform 10000 such simulation. But after running some simulation it's getting stuck and then giving error message 'unspecified error'.  I am using 'SAMG' equation system solver. Number of threads for matrix assembly as single threaded. Number of threads for matrix solving as default. The another problem I am facing is computational time. It's taking a lot of computational time to solve a single simulation( approximately 5-10 minutes or sometimes more). How can I increase the computational speed? when I am using PETSc_ksp it's giving another error : Attempting to use an MPI routine after finalizing MPICH.

Thank you for your reply. But in my script: #There are 100 Monte carlo simulation #In each iteration of the Monte Carlo simulation, there are rainfall events lasting up to 24 hours with a step of 0.25   hours. #For each rainfall event, the stage time series is generated, which is added to the FEFLOW model for simulation and   analysis. #After each simulation, the maximum saturation at the observation node is to be extracted. Please help me with this, I will be thankful to you

I am using a Monte Carlo simulation. Due to this simulation, it is taking a lot of time (approximately 40 hours). Can I use parallel computing in my script? If yes, then how?

I want to use constant time stepping for simulation in feflow 8.0 . As per python API the function for constant time stepping accepts two arguments: initial increment and total number of time steps. Using python how can I set total number of time steps?

I want to import a .pow file using python, that contains time series with different IDs. I want to add this file in the .FEM model using python. I want to assign boundary conditions (hydraulic head = ID 1 from pow file. and max and min flow rates = rest of the IDs in alternate (like ID = 2 for min and ID = 3 for max flow rate etc.). This .pow file contain the data of time and stage for each ID(eg is shown below). How can I do that?