ELMER SOLVER (v 8.3) STARTED AT: 2017/07/04 13:42:29 ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.3 (Rev: 8068c86, Compiled: 2017-06-17) MAIN: HYPRE library linked in. MAIN: MUMPS library linked in. MAIN: ============================================================= MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: structure_mesh_mortar_V02.sif LoadInputFile: Scanning input file: structure_mesh_mortar_V02.sif LoadInputFile: Loading input file: structure_mesh_mortar_V02.sif Model Input: Unlisted keyword: [calculate magnetic field density] in section: [solver 3] Model Input: Unlisted keyword: [magnetic flux density {n}] in section: [boundary condition 1] Model Input: Unlisted keyword: [magnetic flux density {n}] in section: [boundary condition 2] Model Input: Unlisted keyword: [p {e}] in section: [boundary condition 3] Model Input: Unlisted keyword: [magnetic flux density {n}] in section: [boundary condition 4] Loading user function library: [StatCurrentSolve]...[StatCurrentSolver_Init0] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init0] Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init0] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0] Loading user function library: [SaveData]...[SaveScalars_Init0] Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy_Init0] LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./structured_extrusion_small_constant_element_size LoadMesh: Reading header info from file: ./structured_extrusion_small_constant_element_size/mesh.header LoadMesh: Number of nodes in mesh: 276777 LoadMesh: Number of bulk elements in mesh: 256000 LoadMesh: Number of boundary elements in mesh: 40960 LoadMesh: Initial number of max element nodes: 8 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./structured_extrusion_small_constant_element_size/mesh.nodes LoadMesh: Performing coordinate mapping LoadMesh: Dimension of model is: 3 LoadMesh: Dimension of mesh is: 3 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./structured_extrusion_small_constant_element_size/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./structured_extrusion_small_constant_element_size/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Remapping bodies LoadMesh: Minimum initial body index: 1 LoadMesh: Maximum initial body index: 4 LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 6 ElmerAsciiMesh: Performing step: 5 NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining faces in 3D mesh FindMeshFaces3D: Number of faces found: 788480 FindMeshEdges: Determining edges in 3D mesh FindMeshEdges3D: Number of edges found: 809256 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done LoadMesh: Elapsed REAL time: 2.8882 (s) WARNING:: CheckTimer: Requesting time from non-existing timer: LoadMesh MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.2125 (s) WARNING:: CheckTimer: Requesting time from non-existing timer: MeshStabParams MAIN: ------------------------------------- AddSolvers: Setting up 6 solvers AddSolvers: Setting up solver 1: static current conduction AddEquationBasics: Using procedure: StatCurrentSolve StatCurrentSolver AddEquationBasics: Setting up solver: static current conduction Loading user function library: [StatCurrentSolve]...[StatCurrentSolver_Init] Loading user function library: [StatCurrentSolve]...[StatCurrentSolver_bulk] Loading user function library: [StatCurrentSolve]...[StatCurrentSolver] AddEquationBasics: Creating standard variable: potential OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: static current conduction...done. OptimizeBandwidth: Half bandwidth without optimization: 25454 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 547 OptimizeBandwidth: --------------------------------------------------------- AddSolvers: Setting up solver 2: mgdynamics AddEquationBasics: Using procedure: MagnetoDynamics WhitneyAVSolver AddEquationBasics: Setting up solver: mgdynamics Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_bulk] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver] AddEquationBasics: Creating standard variable: p OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: mgdynamics...done. OptimizeBandwidth: Half bandwidth without optimization: 458402 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 12354 OptimizeBandwidth: --------------------------------------------------------- AddSolvers: Setting up solver 3: mgdynamicscalc AddEquationBasics: Using procedure: MagnetoDynamics MagnetoDynamicsCalcFields AddEquationBasics: Setting up solver: mgdynamicscalc Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init] Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_bulk] Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields] AddEquationBasics: Creating standard variable: hr_dummy OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: mgdynamicscalc...done. OptimizeBandwidth: Half bandwidth without optimization: 153818 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 2804 OptimizeBandwidth: --------------------------------------------------------- AddSolvers: Setting up solver 4: AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddSolvers: Setting up solver 5: AddEquationBasics: Using procedure: SaveData SaveScalars Loading user function library: [SaveData]...[SaveScalars_Init] Loading user function library: [SaveData]...[SaveScalars] AddSolvers: Setting up solver 6: never AddEquationBasics: Using procedure: MagnetoDynamics MagnetoDynamics_Dummy AddEquationBasics: Setting up solver: never Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy_Init] Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy_bulk] Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy] AddEquationBasics: Creating standard variable: cf_dummy AddMeshCoordinatesAndTime: Setting mesh coordinates and time SetInitialConditions: Setting up initial conditions (if any) MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: static current conduction StatCurrentSolve: ------------------------------------- StatCurrentSolve: STAT CURRENT SOLVER: StatCurrentSolve: ------------------------------------- DefaultStart: Starting solver: static current conduction StatElecSolve: Starting Assembly... DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions SetNodalLoads: Checking for nodal loads for variable: potential DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StatCurrentSolve: Assembly (s) : 0.76800000000000068 SolveSystem: Solving linear system IterSolver: Using iterative method: bicgstab CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 744997 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 744998 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.4726E-01 2 0.2885E-01 3 0.1380E-01 4 0.4563E-02 5 0.3705E-02 6 0.2785E-02 7 0.2226E-02 8 0.1810E-02 9 0.1537E-02 10 0.1360E-02 11 0.1206E-02 12 0.1053E-02 13 0.9444E-03 14 0.8742E-03 15 0.7989E-03 16 0.7358E-03 17 0.6782E-03 18 0.6331E-03 19 0.5924E-03 20 0.5507E-03 21 0.5041E-03 22 0.4592E-03 23 0.4172E-03 24 0.3609E-03 25 0.3001E-03 26 0.2388E-03 27 0.1770E-03 28 0.1201E-03 29 0.6692E-04 30 0.3042E-04 31 0.1096E-04 32 0.4279E-05 33 0.1455E-05 34 0.5046E-06 35 0.1954E-06 36 0.1070E-06 37 0.7724E-07 38 0.6275E-07 39 0.5139E-07 40 0.4268E-07 41 0.3604E-07 42 0.3046E-07 43 0.2576E-07 44 0.2068E-07 45 0.1656E-07 46 0.1204E-07 47 0.7921E-08 47 0.7921E-08 ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.34889838E-01 2.0000000 ) :: static current conduction StatCurrentSolve: Solve (s) : 0.33599999999999852 StatCurrentSolve: Total Heating Power : 81.713547615054750 StatCurrentSolve: Effective Resistance : 3.7019565155223780E-005 StatCurrentSolve: Total Volume : 3.9999999999891652E-007 StatCurrentSolve: Total Current : 1485.7008657282681 StatCurrentSolve: Control Scaling : 2.6923320112890021E-002 DefaultFinish: Finished solver: static current conduction ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.93935027E-03 2.0000000 ) :: static current conduction GenerateProjectors: Generating mortar projectors PeriodicProjector: ----------------------------------------------------- PeriodicProjector: Creating projector between BCs 5 and 6 CreateInterfaceMeshes: Making a list of elements at interface CreateInterfaceMeshes: Number of interface elements: 1600, 1600 CreateInterfaceMeshes: Number of interface nodes: 1717, 1717 PeriodicProjector: Enforcing > Level Projector = True < with dimensional reduction RadialInterfaceMeshes: Transformed extrema for this boundary (phi,r,z) RadialInterfaceMeshes: Coordinate 1: 0.000E+00 2.500E-02 RadialInterfaceMeshes: Coordinate 2: 1.800E+02 1.800E+02 RadialInterfaceMeshes: Coordinate 3: 0.000E+00 4.000E-03 RadialInterfaceMeshes: Discrepancy from constant angle (degs): 0.000E+00 RadialInterfaceMeshes: Transformed extrema for target boundary (phi,r,z) RadialInterfaceMeshes: Coordinate 1: 0.000E+00 2.500E-02 RadialInterfaceMeshes: Coordinate 2: 0.000E+00 0.000E+00 RadialInterfaceMeshes: Coordinate 3: 0.000E+00 4.000E-03 RadialInterfaceMeshes: Discrepancy from constant angle (degs): 0.000E+00 RadialInterfaceMeshes: Discrepancy in maximum radius: 0.000E+00 RadialInterfaceMeshes: Discrepancy in minimum radius: 0.000E+00 LevelProjector: Creating projector for a levelized mesh LevelProjector: Checking mesh skew LevelProjector: Maximum skew in this mesh: 0.000E+00 LevelProjector: Checking mesh skew LevelProjector: Maximum skew in target mesh: 0.000E+00 LevelProjector: Initial number of slave nodes 1717 out of 1717 LevelProjector: Number of potential nodes in projector: 1717 LevelProjector: Eliminating redundant nodes from projector: 17 LevelProjector: Number of active nodes in projector: 1700 LevelProjector: Number of potential edges in projector: 3316 LevelProjector: Eliminating redundant edges from projector: 16 LevelProjector: Number of active edges in projector: 3300 LevelProjector: Max number of rows in projector: 5000 LevelProjector: Creating strong stride projector for nodal dofs LevelProjector: Creating strong stride projector for edges assuming strides LevelProjector: Minimum absolute sum of edge weights: 1.00000E+00 LevelProjector: Maximum absolute sum of edge weights: 1.00000E+00 LevelProjector: Created strong constraints for edge dofs List_ToCRSMatrix: Number of entries in CRS matrix: 10000 ListToCRSMatrix: Matrix format changed from List to CRS LevelProjector: Projector created PeriodicProjector: Elapsed REAL time: 0.1149 (s) WARNING:: CheckTimer: Requesting time from non-existing timer: PeriodicProjector PeriodicProjector: Projector created, now exiting... SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: mgdynamics WhitneyAVSolver: ------------------------------------------- WhitneyAVSolver: Solving the AV equations with edge elements OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: mgdynamics...done. OptimizeBandwidth: Half bandwidth without optimization: 153818 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 2804 OptimizeBandwidth: --------------------------------------------------------- DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions SetNodalLoads: Checking for nodal loads for variable: jfix DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set GenerateConstraintMatrix: There are 5000 initial rows in constraint matrices GenerateConstraintMatrix: Adding projector for BC: 5 GenerateConstraintMatrix: Allocating 1700 rows and 3400 nonzeros SolveSystem: Solving linear system SolveSystem: Solving linear system with constraint matrix SolveWithLinearRestriction: SolveWithLinearRestriction: Adding ConstraintMatrix into CollectionMatrix SolveWithLinearRestriction: Adding Stiffness Matrix into CollectionMatrix SolveWithLinearRestriction: Reverting CollectionMatrix back to CRS matrix List_ToCRSMatrix: Number of entries in CRS matrix: 7105709 ListToCRSMatrix: Matrix format changed from List to CRS SolveWithLinearRestriction: CollectionMatrix done SolveWithLinearRestriction: CollectionVector done SolveWithLinearRestriction: Now going for the coupled linear system IterSolver: Using iterative method: bicgstab List_ToCRSMatrix: Number of entries in CRS matrix: 285277 ListToCRSMatrix: Matrix format changed from List to CRS CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 285277 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 285278 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01 1 0.8044E-01 2 0.3090E-01 3 0.1299E-01 4 0.7984E-02 5 0.5837E-02 6 0.4774E-02 7 0.3980E-02 8 0.3221E-02 9 0.2561E-02 10 0.1979E-02 11 0.1668E-02 12 0.1495E-02 13 0.1351E-02 14 0.1231E-02 15 0.1121E-02 16 0.1023E-02 17 0.9327E-03 18 0.8506E-03 19 0.7874E-03 20 0.7375E-03 21 0.6941E-03 22 0.6549E-03 23 0.6146E-03 24 0.5752E-03 25 0.5378E-03 26 0.5022E-03 27 0.4721E-03 28 0.4449E-03 29 0.4195E-03 30 0.3952E-03 31 0.3704E-03 32 0.3482E-03 33 0.3272E-03 34 0.3065E-03 35 0.2850E-03 36 0.2618E-03 37 0.2385E-03 38 0.2167E-03 39 0.1993E-03 40 0.1772E-03 41 0.1787E-03 42 0.1750E-03 43 0.1599E-03 44 0.1565E-03 45 0.1627E-03 46 0.1461E-03 47 0.1422E-03 48 0.1501E-03 49 0.1377E-03 50 0.1279E-03 51 0.1129E-03 52 0.9821E-04 53 0.9558E-04 54 0.2382E-03 55 0.1063E-03 56 0.8640E-04 57 0.8236E-04 58 0.7766E-04 59 0.8365E-04 60 0.7518E-04 61 0.7332E-04 62 0.7216E-04 63 0.6899E-04 64 0.6781E-04 65 0.6626E-04 66 0.6755E-04 67 0.5367E-04 68 0.6242E-04 69 0.6112E-04 70 0.5799E-04 71 0.5688E-04 72 0.4927E-04 73 0.4481E-04 74 0.4206E-04 75 0.4101E-04 76 0.4049E-04 77 0.3970E-04 78 0.3929E-04 79 0.3702E-04 80 0.3734E-04 81 0.3700E-04 82 0.3697E-04 83 0.3656E-04 84 0.3496E-04 85 0.3414E-04 86 0.3278E-04 87 0.3364E-04 88 0.3329E-04 89 0.2946E-04 90 0.3661E-04 91 0.2798E-04 92 0.2755E-04 93 0.2699E-04 94 0.2666E-04 95 0.2605E-04 96 0.2512E-04 97 0.2486E-04 98 0.2484E-04 99 0.2438E-04 100 0.2379E-04 101 0.2353E-04 102 0.2312E-04 103 0.2285E-04 104 0.2248E-04 105 0.2213E-04 106 0.2277E-04 107 0.2175E-04 108 0.2065E-04 109 0.2053E-04 110 0.1949E-04 111 0.1895E-04 112 0.1838E-04 113 0.1789E-04 114 0.1748E-04 115 0.1666E-04 116 0.1678E-04 117 0.1578E-04 118 0.1525E-04 119 0.1222E-04 120 0.2596E-04 121 0.1386E-04 122 0.1332E-04 123 0.1237E-04 124 0.1176E-04 125 0.1178E-04 126 0.1105E-04 127 0.1040E-04 128 0.1052E-04 129 0.9651E-05 130 0.9252E-05 131 0.8763E-05 132 0.7002E-05 133 0.1090E-04 134 0.7697E-05 135 0.7693E-05 136 0.7625E-05 137 0.7260E-05 138 0.7107E-05 139 0.7004E-05 140 0.6788E-05 141 0.6812E-05 142 0.6720E-05 143 0.6268E-05 144 0.6108E-05 145 0.6055E-05 146 0.5839E-05 147 0.5787E-05 148 0.5728E-05 149 0.5621E-05 150 0.5828E-05 151 0.5528E-05 152 0.5416E-05 153 0.5289E-05 154 0.5230E-05 155 0.5103E-05 156 0.5070E-05 157 0.4947E-05 158 0.4908E-05 159 0.4859E-05 160 0.4801E-05 161 0.4750E-05 162 0.4704E-05 163 0.4640E-05 164 0.4601E-05 165 0.4547E-05 166 0.4507E-05 167 0.4465E-05 168 0.4372E-05 169 0.4345E-05 170 0.4389E-05 171 0.4278E-05 172 0.4148E-05 173 0.4099E-05 174 0.3996E-05 175 0.3955E-05 176 0.3870E-05 177 0.3813E-05 178 0.3753E-05 179 0.3656E-05 180 0.3604E-05 181 0.3517E-05 182 0.3445E-05 183 0.3365E-05 184 0.3249E-05 185 0.1948E-05 186 0.1468E-05 187 0.1303E-05 188 0.1530E-05 189 0.1466E-05 190 0.1170E-05 191 0.1144E-05 192 0.1066E-05 193 0.1030E-05 194 0.9437E-06 195 0.1188E-05 196 0.9257E-06 197 0.9045E-06 198 0.8593E-06 199 0.8196E-06 200 0.8002E-06 201 0.7350E-06 202 0.7341E-06 203 0.7318E-06 204 0.6621E-06 205 0.6374E-06 206 0.5890E-06 207 0.5493E-06 208 0.4299E-06 209 0.5125E-06 210 0.4909E-06 211 0.4408E-06 212 0.4366E-06 213 0.4247E-06 214 0.3789E-06 215 0.3650E-06 216 0.3400E-06 217 0.3291E-06 218 0.3156E-06 219 0.2986E-06 220 0.2935E-06 221 0.2980E-06 222 0.2682E-06 223 0.2628E-06 224 0.2610E-06 225 0.2522E-06 226 0.2371E-06 227 0.2272E-06 228 0.2221E-06 229 0.1848E-06 230 0.2001E-06 231 0.1970E-06 232 0.3329E-06 233 0.1477E-06 234 0.1191E-06 235 0.1152E-06 236 0.1093E-06 237 0.1091E-06 238 0.1072E-06 239 0.1013E-06 240 0.9865E-07 241 0.9514E-07 242 0.9342E-07 243 0.8946E-07 244 0.8736E-07 245 0.7944E-07 246 0.8039E-07 247 0.7605E-07 248 0.7157E-07 249 0.1291E-06 250 0.6292E-07 251 0.6014E-07 252 0.5401E-07 253 0.5509E-07 254 0.4730E-07 255 0.4537E-07 256 0.3852E-07 257 0.3693E-07 258 0.3689E-07 259 0.3532E-07 260 0.2699E-07 261 0.2438E-07 262 0.1135E-07 263 0.7406E-08 264 0.1671E-07 265 0.1512E-07 266 0.8152E-08 267 0.1025E-07 268 0.8372E-08 269 0.6329E-08 270 0.6347E-08 271 0.5886E-08 272 0.4617E-08 273 0.9293E-08 274 0.5594E-08 275 0.3406E-08 276 0.2149E-08 277 0.2870E-08 278 0.2442E-08 279 0.1645E-08 280 0.1703E-08 281 0.1401E-08 282 0.1015E-08 283 0.1081E-08 284 0.1085E-08 285 0.7386E-09 286 0.6436E-09 287 0.5120E-09 288 0.4269E-09 289 0.3517E-09 290 0.2920E-09 291 0.2636E-09 292 0.2088E-09 293 0.1950E-09 294 0.2578E-09 295 0.2183E-09 296 0.1356E-09 297 0.1274E-09 298 0.1018E-09 299 0.9475E-10 300 0.6914E-10 301 0.6911E-10 302 0.8194E-10 303 0.4773E-10 304 0.4808E-10 305 0.4704E-10 306 0.3979E-10 307 0.3691E-10 308 0.3442E-10 309 0.3118E-10 310 0.2909E-10 311 0.2705E-10 312 0.2514E-10 313 0.2278E-10 314 0.2177E-10 315 0.1837E-10 316 0.1790E-10 317 0.1783E-10 318 0.1658E-10 319 0.1419E-10 320 0.1346E-10 321 0.1309E-10 322 0.1087E-10 323 0.1005E-10 324 0.7975E-11 324 0.7975E-11 SolveWithLinearRestriction: Picking solution from collection solution SolveWithLinearRestriction: All done DefaultStart: Starting solver: mgdynamics MGDynAssembly: Elapsed REAL time: 21.3402 (s) WARNING:: CheckTimer: Requesting time from non-existing timer: MGDynAssembly DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions SetNodalLoads: Checking for nodal loads for variable: p DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set Program received signal SIGSEGV: Segmentation fault - invalid memory reference. Backtrace for this error: #0 0x7F514D721E08 #1 0x7F514D720F90 #2 0x7F514D3724AF #3 0x7F513E4BAB74 #4 0x7F513E4BCD53 #5 0x7F513E4CFC90 #6 0x7F514DC18030 #7 0x7F514DC2A904 #8 0x7F514DC2BEA0 #9 0x7F514DC2D829 #10 0x7F514DDE1779 #11 0x401075 in MAIN__ at Solver.F90:69 Segmentation fault