Hello,
I'm trying to solve a synthetic glacier but when I run the solver I got this error message:
At line 2158 of file Lists.f90
Fortran runtime error: Bad real number in item 1 of list input
Could anyone help me?
I enclose the sif file
Thanks
---- The sif file ----
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! !!
!! Cubo con deslizamiento !!
!! !!
!! !!
!! !!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
check keywords warn
echo on
$n = 3.0 !Exponente rel. const.
$rhoi = 9.1376e-19 !Densidad del hielo
$rhow = 10.0e-19 ! Densidad del agua
$g = 9.7696e15 !Gravedad
$AA = 1000.0
$b = 0.215443469
$HH = 0.5*AA*2.0*pow(2.0,0.5)/3.0
$fp = fopen("par_deslizamiento.dat","r")
$xl = fscanf(fp,"%lf")
$yl = fscanf(fp,"%lf")
$zl = fscanf(fp,"%lf")
$vvlx = fscanf(fp,"%lf")
$vvly = fscanf(fp,"%lf")
$vvlz = fscanf(fp,"%lf")
$function U(X) import xl,yl,zl,vvlx,vvly,vvlz,n,b_cons,rhoi,g,HH { _U = ((2.0/(n+1.0))*(pow((rhoi*g*(1.0/3.0)/b_cons),n))*pow(HH,n+1.0)-pow(HH-((X(2)-(zl-(vvlx*(X(0)-xl)+vvly*(X(1)-yl))/vvlz))*2.0*pow(2.0,0.5)/3.0),n+1.0))*(pow(2.
0,0.5)/2.0)*(2.0*pow(2.0,0.5)/3.0) }
$function V(X) import xl,yl,zl,vvlx,vvly,vvlz,n,b_cons,rhoi,g,HH { _V =
((2.0/(n+1.0))*(pow((rhoi*g*(1.0/3.0)/b_cons),n))*pow(HH,n+1.0)-pow(HH-((X(2)-(zl-(vvlx*(X(0)-xl)+vvly*(X(1)-yl))/vvlz))*2.0*pow(2.0,0.5)/3.0),n+1.0))*(pow(2.
0,0.5)/2.0)*(2.0*pow(2.0,0.5)/3.0) }
$function W(X) import xl,yl,zl,vvlx,vvly,vvlz,n,b_cons,rhoi,g,HH { _W = ((2.0/(n+1.0))*(pow((rhoi*g*(1.0/3.0)/b_cons),n))*pow(HH,n+1.0)-pow(HH-((X(2)-(zl-(vvlx*(X(0)-xl)+vvly*(X(1)-yl))/vvlz))*2.0*pow(2.0,0.5)/3.0),n+1.0))*(1.0/3.0) }
Header
Mesh DB "." "cubo"
End
Constants
! No constant needed
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Simulation
Coordinate System = Cartesian
Simulation Type = Steady State
Steady State Min Iterations = 1
Steady State Max Iterations = 1
Output File = "cubo_steady".result"
Post File = "cubo_steady".ep"
max output level = 100
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Body 1
Equation = 1
Body Force = 1
Material = 1
Initial Condition = 1
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Initial Condition 1
Pressure = Real 0.0
Velocity 1 = Real 0.0
Velocity 2 = Real 0.0
Velocity 3 = Real 0.0
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Body Force 1
Flow BodyForce 1 = Real 0.0
Flow BodyForce 2 = Real 0.0
Flow BodyForce 3 = Real -9.7696e15 !MPa - a - m Esto es el vector gravedad (0,0,-1)*g
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Material 1
Density = Real 9.1376e-19 ! MPa - a - m Densidad del hielo
Viscosity Model = String "power law"
Viscosity = Real 0.170997595 ! MPa - a - m Esto es (2A)^-1/n (o sea B*2^-1/n)(Tomo A=0.1 bar-3 a-1 del artículo de Vieli que es 0.1e3 MPa-3 a-1)
Viscosity Exponent = Real MATC 1/n
Critical Shear Rate = Real 1.0e-10
Cauchy = Logical False
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Solver 1
Equation = "Navier-Stokes"
Stabilization Method = String Bubbles
Linear System Solver = Direct
Linear System Direct Method = umfpack
Nonlinear System Max Iterations = 100
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Newton After Iterations = 5
Nonlinear System Newton After Tolerance = 1.0e-02
Nonlinear System Relaxation Factor = 1.00
Steady State Convergence Tolerance = Real 1.0e-3
End
Solver 2
Equation = String "StressSolver"
Procedure = File "./ComputeDevStressNS" "ComputeDevStress"
! this is just a dummy, hence no output is needed
!-----------------------------------------------------------------------
Variable = -nooutput "Sij"
Variable DOFs = 1
! the name of the variable containing the flow solution (U,V,W,Pressure)
!-----------------------------------------------------------------------
Flow Solver Name = String "Flow Solution"
Exported Variable 1 = "Stress" ! [Sxx, Syy, Szz, Sxy] in 2D
! [Sxx, Syy, Szz, Sxy, Syz, Szx] in 3D
Exported Variable 1 DOFs = 6 ! 4 in 2D, 6 in 3D
Linear System Solver = "Iterative"
Linear System Iterative Method = "BiCGStab"
Linear System Max Iterations = 300
Linear System Convergence Tolerance = 1.0E-09
Linear System Abort Not Converged = True
Linear System Preconditioning = "ILU0"
Linear System Residual Output = 1
End
Solver 3
Exec Solver = After All
Procedure = File "SaveData" "SaveLine"
Filename = "hansbreen_surface.dat"
File Append = Logical False
End
Solver 4
Exec Solver = After TimeStep
Procedure = "./MySaveData" "SaveScalars"
Filename = "hansbreen_scalars.dat"
File Append = Logical True
Variable 1 = String "flow solution"
Operator 1 = String "Volume"
Variable 2 = String "Velocity 1"
Operator 2 = String "max abs"
Variable 3 = String "flow solution"
Operator 3 = String "Convective flux"
Variable 4 = String "cpu time"
Variable 5 = String "cpu memory"
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Equation 1
Active Solvers(4) = 1 2 3 4
NS Convect = Logical False
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Bedrock
Boundary Condition 1
Target Boundaries = 1
Flow Force BC = Logical True
Normal-Tangential Velocity = Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
End
! Lateral
Boundary Condition 2
Target Boundaries = 2
Flow Force BC = Logical True
Normal-Tangential Velocity = Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True !Para salvar los datos de flujo por cada
superficie
End
! Front
Boundary Condition 3
Target Boundaries = 3
Flow Force BC = Logical True
Normal-Tangential Velocity = Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True !Para salvar los datos de flujo por cada
superficie
End
! Lateral
Boundary Condition 4
Target Boundaries = 4
Flow Force BC = Logical True
Normal-Tangential Velocity = Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True !Para salvar los datos de flujo por cada
superficie
End
! Divide
Boundary Condition 5
Target Boundaries = 5
Flow Force BC = Logical True
Normal-Tangential Velocity = Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True !Para salvar los datos de flujo por cada
superficie
End
! Upper Surface
Boundary Condition 6
Target Boundaries = 6
Save Line = Logical True !Para salvar los datos de la superficie
Flux Integrate = Logical True !Para salvar los datos de flujo por cada
superficie
End
Fortran runtime error
Re: Fortran runtime error
Hi Jaime,
hard to say anything without a testcase (including the input file par_deslizamiento.dat).
Some calculation might result in overflow or nan, for example. The sif-file syntax
seems OK.
Regards, Juha
hard to say anything without a testcase (including the input file par_deslizamiento.dat).
Some calculation might result in overflow or nan, for example. The sif-file syntax
seems OK.
Regards, Juha
Re: Fortran runtime error
Hi Juha,
Thanks for your fast answer.
This is the input file par_deslizamiento.dat:
1200.000 200.0000 900.0000
900000.0 900000.0 3600000.
and Thomas this is the complete output:
ELMER SOLVER (v 5.5.0) STARTED AT: 2010/09/02 11:57:00
MAIN:
MAIN: ==========================================
MAIN: E L M E R S O L V E R S T A R T I N G
MAIN: Library version: 5.5.0 (Rev: )
MAIN: HYPRE library linked in.
MAIN: ==========================================
MAIN:
MAIN:
MAIN: -----------------------
MAIN: Reading Model ...
Header
Mesh DB "." "cubo"
End
Constants
End
Simulation
Coordinate System = Cartesian
Cartesian
Simulation Type = Steady State
Steady State
Steady State Min Iterations = 1
1
Steady State Max Iterations = 1
1
Output File = "cubo_steady".result"
"cubo_steady".result"
Post File = "cubo_steady".ep"
"cubo_steady".ep"
max output level = 100
100
End
Body 1
Equation = 1
1
Body Force = 1
1
Material = 1
1
Initial Condition = 1
1
End
Initial Condition 1
Pressure = Real 0.0
Real 0.0
Velocity 1 = Real 0.0
Real 0.0
Velocity 2 = Real 0.0
Real 0.0
Velocity 3 = Real 0.0
Real 0.0
End
Body Force 1
Flow BodyForce 1 = Real 0.0
Real 0.0
Flow BodyForce 2 = Real 0.0
Real 0.0
Flow BodyForce 3 = Real -9.7696e15
Real -9.7696e15
End
Material 1
Density = Real 9.1376e-19
Real 9.1376e-19
Viscosity Model = String "power law"
String "power law"
Viscosity = Real 0.170997595
Real 0.170997595
Viscosity Exponent = Real MATC 1/n
Real MATC 1/n
Critical Shear Rate = Real 1.0e-10
Real 1.0e-10
Cauchy = Logical False
Logical False
End
Solver 1
Equation = "Navier-Stokes"
"Navier-Stokes"
Stabilization Method = String Bubbles
String Bubbles
Linear System Solver = Direct
Direct
Linear System Direct Method = umfpack
umfpack
Nonlinear System Max Iterations = 100
100
Nonlinear System Convergence Tolerance = 1.0e-5
1.0e-5
Nonlinear System Newton After Iterations = 5
5
Nonlinear System Newton After Tolerance = 1.0e-02
1.0e-02
Nonlinear System Relaxation Factor = 1.00
1.00
Steady State Convergence Tolerance = Real 1.0e-3
Real 1.0e-3
End
Solver 2
Equation = String "StressSolver"
String "StressSolver"
Procedure = File "./ComputeDevStressNS" "ComputeDevStress"
File "./ComputeDevStressNS" "ComputeDevStress"
Variable = -nooutput "Sij"
-nooutput "Sij"
Variable DOFs = 1
1
Flow Solver Name = String "Flow Solution"
String "Flow Solution"
Exported Variable 1 = "Stress"
"Stress"
Exported Variable 1 DOFs = 6
6
Linear System Solver = "Iterative"
"Iterative"
Linear System Iterative Method = "BiCGStab"
"BiCGStab"
Linear System Max Iterations = 300
300
Linear System Convergence Tolerance = 1.0E-09
1.0E-09
Linear System Abort Not Converged = True
True
Linear System Preconditioning = "ILU0"
"ILU0"
Linear System Residual Output = 1
1
End
Solver 3
Exec Solver = After All
After All
Procedure = File "SaveData" "SaveLine"
File "SaveData" "SaveLine"
Filename = "hansbreen_surface.dat"
"hansbreen_surface.dat"
File Append = Logical False
Logical False
End
Solver 4
Exec Solver = After TimeStep
After TimeStep
Procedure = "./MySaveData" "SaveScalars"
"./MySaveData" "SaveScalars"
Filename = "hansbreen_scalars.dat"
"hansbreen_scalars.dat"
File Append = Logical True
Logical True
Variable 1 = String "flow solution"
String "flow solution"
Operator 1 = String "Volume"
String "Volume"
Variable 2 = String "Velocity 1"
String "Velocity 1"
Operator 2 = String "max abs"
String "max abs"
Variable 3 = String "flow solution"
String "flow solution"
Operator 3 = String "Convective flux"
String "Convective flux"
Variable 4 = String "cpu time"
String "cpu time"
Variable 5 = String "cpu memory"
String "cpu memory"
End
Equation 1
Active Solvers(4) = 1 2 3 4
Size 4 = 1 2 3 4
1 2 3 4
NS Convect = Logical False
Logical False
End
Boundary Condition 1
Target Boundaries = 1
1
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
End
Boundary Condition 2
Target Boundaries = 2
2
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 3
Target Boundaries = 3
3
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 4
Target Boundaries = 4
4
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 5
Target Boundaries = 5
5
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 6
Target Boundaries = 6
6
Save Line = Logical True
Logical True
Flux Integrate = Logical True
Logical True
End
check keywords warn
echo on
Header
Mesh DB "." "cubo"
End
Constants
End
Simulation
Coordinate System = Cartesian
Cartesian
Simulation Type = Steady State
Steady State
Steady State Min Iterations = 1
1
Steady State Max Iterations = 1
1
Output File = "cubo_steady".result"
"cubo_steady".result"
Post File = "cubo_steady".ep"
"cubo_steady".ep"
max output level = 100
100
End
Body 1
Equation = 1
1
Body Force = 1
1
Material = 1
1
Initial Condition = 1
1
End
Initial Condition 1
Pressure = Real 0.0
Real 0.0
Velocity 1 = Real 0.0
Real 0.0
Velocity 2 = Real 0.0
Real 0.0
Velocity 3 = Real 0.0
Real 0.0
End
Body Force 1
Flow BodyForce 1 = Real 0.0
Real 0.0
Flow BodyForce 2 = Real 0.0
Real 0.0
Flow BodyForce 3 = Real -9.7696e15
Real -9.7696e15
End
Material 1
Density = Real 9.1376e-19
Real 9.1376e-19
Viscosity Model = String "power law"
String "power law"
Viscosity = Real 0.170997595
Real 0.170997595
Viscosity Exponent = Real MATC 1/n
Real MATC 1/n
Critical Shear Rate = Real 1.0e-10
Real 1.0e-10
Cauchy = Logical False
Logical False
Model Input: Unlisted keyword: [cauchy] in section: [material 1]
End
Solver 1
Equation = "Navier-Stokes"
"Navier-Stokes"
Stabilization Method = String Bubbles
String Bubbles
Linear System Solver = Direct
Direct
Linear System Direct Method = umfpack
umfpack
Nonlinear System Max Iterations = 100
100
Nonlinear System Convergence Tolerance = 1.0e-5
1.0e-5
Nonlinear System Newton After Iterations = 5
5
Nonlinear System Newton After Tolerance = 1.0e-02
1.0e-02
Nonlinear System Relaxation Factor = 1.00
1.00
Steady State Convergence Tolerance = Real 1.0e-3
Real 1.0e-3
End
Solver 2
Equation = String "StressSolver"
String "StressSolver"
Procedure = File "./ComputeDevStressNS" "ComputeDevStress"
File "./ComputeDevStressNS" "ComputeDevStress"
Variable = -nooutput "Sij"
-nooutput "Sij"
Variable DOFs = 1
1
Flow Solver Name = String "Flow Solution"
String "Flow Solution"
Model Input: Unlisted keyword: [flow solver name] in section: [solver 2]
Exported Variable 1 = "Stress"
"Stress"
Exported Variable 1 DOFs = 6
6
Linear System Solver = "Iterative"
"Iterative"
Linear System Iterative Method = "BiCGStab"
"BiCGStab"
Linear System Max Iterations = 300
300
Linear System Convergence Tolerance = 1.0E-09
1.0E-09
Linear System Abort Not Converged = True
True
Linear System Preconditioning = "ILU0"
"ILU0"
Linear System Residual Output = 1
1
End
Solver 3
Exec Solver = After All
After All
Procedure = File "SaveData" "SaveLine"
File "SaveData" "SaveLine"
Filename = "hansbreen_surface.dat"
"hansbreen_surface.dat"
File Append = Logical False
Logical False
End
Solver 4
Exec Solver = After TimeStep
After TimeStep
Procedure = "./MySaveData" "SaveScalars"
"./MySaveData" "SaveScalars"
Filename = "hansbreen_scalars.dat"
"hansbreen_scalars.dat"
File Append = Logical True
Logical True
Variable 1 = String "flow solution"
String "flow solution"
Operator 1 = String "Volume"
String "Volume"
Variable 2 = String "Velocity 1"
String "Velocity 1"
Operator 2 = String "max abs"
String "max abs"
Variable 3 = String "flow solution"
String "flow solution"
Operator 3 = String "Convective flux"
String "Convective flux"
Variable 4 = String "cpu time"
String "cpu time"
Variable 5 = String "cpu memory"
String "cpu memory"
End
Equation 1
Active Solvers(4) = 1 2 3 4
Size 4 = 1 2 3 4
1 2 3 4
NS Convect = Logical False
Logical False
End
Boundary Condition 1
Target Boundaries = 1
1
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
End
Boundary Condition 2
Target Boundaries = 2
2
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 3
Target Boundaries = 3
3
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 4
Target Boundaries = 4
4
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 5
Target Boundaries = 5
5
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 6
Target Boundaries = 6
6
Save Line = Logical True
Logical True
Flux Integrate = Logical True
Logical True
End
MAIN: Done
MAIN: -----------------------
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: navier-stokes...done.
OptimizeBandwidth: Half bandwidth without optimization: 3018
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 672
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [./ComputeDevStressNS]...[ComputeDevStress_Init]
Loading user function library: [./ComputeDevStressNS]...[ComputeDevStress]
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: stresssolver...done.
OptimizeBandwidth: Half bandwidth without optimization: 3018
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 672
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [SaveData]...[SaveLine_Init]
WARNING:: ListFind:
WARNING:: ListFind: Requested property: [Equation], not found
WARNING:: ListFind:
Loading user function library: [SaveData]...[SaveLine]
Loading user function library: [./MySaveData]...[SaveScalars_Init]
WARNING:: ListFind:
WARNING:: ListFind: Requested property: [Equation], not found
WARNING:: ListFind:
Loading user function library: [./MySaveData]...[SaveScalars]
At line 2158 of file Lists.f90
Fortran runtime error: Bad real number in item 1 of list input
Regards, Jaime
Thanks for your fast answer.
This is the input file par_deslizamiento.dat:
1200.000 200.0000 900.0000
900000.0 900000.0 3600000.
and Thomas this is the complete output:
ELMER SOLVER (v 5.5.0) STARTED AT: 2010/09/02 11:57:00
MAIN:
MAIN: ==========================================
MAIN: E L M E R S O L V E R S T A R T I N G
MAIN: Library version: 5.5.0 (Rev: )
MAIN: HYPRE library linked in.
MAIN: ==========================================
MAIN:
MAIN:
MAIN: -----------------------
MAIN: Reading Model ...
Header
Mesh DB "." "cubo"
End
Constants
End
Simulation
Coordinate System = Cartesian
Cartesian
Simulation Type = Steady State
Steady State
Steady State Min Iterations = 1
1
Steady State Max Iterations = 1
1
Output File = "cubo_steady".result"
"cubo_steady".result"
Post File = "cubo_steady".ep"
"cubo_steady".ep"
max output level = 100
100
End
Body 1
Equation = 1
1
Body Force = 1
1
Material = 1
1
Initial Condition = 1
1
End
Initial Condition 1
Pressure = Real 0.0
Real 0.0
Velocity 1 = Real 0.0
Real 0.0
Velocity 2 = Real 0.0
Real 0.0
Velocity 3 = Real 0.0
Real 0.0
End
Body Force 1
Flow BodyForce 1 = Real 0.0
Real 0.0
Flow BodyForce 2 = Real 0.0
Real 0.0
Flow BodyForce 3 = Real -9.7696e15
Real -9.7696e15
End
Material 1
Density = Real 9.1376e-19
Real 9.1376e-19
Viscosity Model = String "power law"
String "power law"
Viscosity = Real 0.170997595
Real 0.170997595
Viscosity Exponent = Real MATC 1/n
Real MATC 1/n
Critical Shear Rate = Real 1.0e-10
Real 1.0e-10
Cauchy = Logical False
Logical False
End
Solver 1
Equation = "Navier-Stokes"
"Navier-Stokes"
Stabilization Method = String Bubbles
String Bubbles
Linear System Solver = Direct
Direct
Linear System Direct Method = umfpack
umfpack
Nonlinear System Max Iterations = 100
100
Nonlinear System Convergence Tolerance = 1.0e-5
1.0e-5
Nonlinear System Newton After Iterations = 5
5
Nonlinear System Newton After Tolerance = 1.0e-02
1.0e-02
Nonlinear System Relaxation Factor = 1.00
1.00
Steady State Convergence Tolerance = Real 1.0e-3
Real 1.0e-3
End
Solver 2
Equation = String "StressSolver"
String "StressSolver"
Procedure = File "./ComputeDevStressNS" "ComputeDevStress"
File "./ComputeDevStressNS" "ComputeDevStress"
Variable = -nooutput "Sij"
-nooutput "Sij"
Variable DOFs = 1
1
Flow Solver Name = String "Flow Solution"
String "Flow Solution"
Exported Variable 1 = "Stress"
"Stress"
Exported Variable 1 DOFs = 6
6
Linear System Solver = "Iterative"
"Iterative"
Linear System Iterative Method = "BiCGStab"
"BiCGStab"
Linear System Max Iterations = 300
300
Linear System Convergence Tolerance = 1.0E-09
1.0E-09
Linear System Abort Not Converged = True
True
Linear System Preconditioning = "ILU0"
"ILU0"
Linear System Residual Output = 1
1
End
Solver 3
Exec Solver = After All
After All
Procedure = File "SaveData" "SaveLine"
File "SaveData" "SaveLine"
Filename = "hansbreen_surface.dat"
"hansbreen_surface.dat"
File Append = Logical False
Logical False
End
Solver 4
Exec Solver = After TimeStep
After TimeStep
Procedure = "./MySaveData" "SaveScalars"
"./MySaveData" "SaveScalars"
Filename = "hansbreen_scalars.dat"
"hansbreen_scalars.dat"
File Append = Logical True
Logical True
Variable 1 = String "flow solution"
String "flow solution"
Operator 1 = String "Volume"
String "Volume"
Variable 2 = String "Velocity 1"
String "Velocity 1"
Operator 2 = String "max abs"
String "max abs"
Variable 3 = String "flow solution"
String "flow solution"
Operator 3 = String "Convective flux"
String "Convective flux"
Variable 4 = String "cpu time"
String "cpu time"
Variable 5 = String "cpu memory"
String "cpu memory"
End
Equation 1
Active Solvers(4) = 1 2 3 4
Size 4 = 1 2 3 4
1 2 3 4
NS Convect = Logical False
Logical False
End
Boundary Condition 1
Target Boundaries = 1
1
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
End
Boundary Condition 2
Target Boundaries = 2
2
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 3
Target Boundaries = 3
3
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 4
Target Boundaries = 4
4
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 5
Target Boundaries = 5
5
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 6
Target Boundaries = 6
6
Save Line = Logical True
Logical True
Flux Integrate = Logical True
Logical True
End
check keywords warn
echo on
Header
Mesh DB "." "cubo"
End
Constants
End
Simulation
Coordinate System = Cartesian
Cartesian
Simulation Type = Steady State
Steady State
Steady State Min Iterations = 1
1
Steady State Max Iterations = 1
1
Output File = "cubo_steady".result"
"cubo_steady".result"
Post File = "cubo_steady".ep"
"cubo_steady".ep"
max output level = 100
100
End
Body 1
Equation = 1
1
Body Force = 1
1
Material = 1
1
Initial Condition = 1
1
End
Initial Condition 1
Pressure = Real 0.0
Real 0.0
Velocity 1 = Real 0.0
Real 0.0
Velocity 2 = Real 0.0
Real 0.0
Velocity 3 = Real 0.0
Real 0.0
End
Body Force 1
Flow BodyForce 1 = Real 0.0
Real 0.0
Flow BodyForce 2 = Real 0.0
Real 0.0
Flow BodyForce 3 = Real -9.7696e15
Real -9.7696e15
End
Material 1
Density = Real 9.1376e-19
Real 9.1376e-19
Viscosity Model = String "power law"
String "power law"
Viscosity = Real 0.170997595
Real 0.170997595
Viscosity Exponent = Real MATC 1/n
Real MATC 1/n
Critical Shear Rate = Real 1.0e-10
Real 1.0e-10
Cauchy = Logical False
Logical False
Model Input: Unlisted keyword: [cauchy] in section: [material 1]
End
Solver 1
Equation = "Navier-Stokes"
"Navier-Stokes"
Stabilization Method = String Bubbles
String Bubbles
Linear System Solver = Direct
Direct
Linear System Direct Method = umfpack
umfpack
Nonlinear System Max Iterations = 100
100
Nonlinear System Convergence Tolerance = 1.0e-5
1.0e-5
Nonlinear System Newton After Iterations = 5
5
Nonlinear System Newton After Tolerance = 1.0e-02
1.0e-02
Nonlinear System Relaxation Factor = 1.00
1.00
Steady State Convergence Tolerance = Real 1.0e-3
Real 1.0e-3
End
Solver 2
Equation = String "StressSolver"
String "StressSolver"
Procedure = File "./ComputeDevStressNS" "ComputeDevStress"
File "./ComputeDevStressNS" "ComputeDevStress"
Variable = -nooutput "Sij"
-nooutput "Sij"
Variable DOFs = 1
1
Flow Solver Name = String "Flow Solution"
String "Flow Solution"
Model Input: Unlisted keyword: [flow solver name] in section: [solver 2]
Exported Variable 1 = "Stress"
"Stress"
Exported Variable 1 DOFs = 6
6
Linear System Solver = "Iterative"
"Iterative"
Linear System Iterative Method = "BiCGStab"
"BiCGStab"
Linear System Max Iterations = 300
300
Linear System Convergence Tolerance = 1.0E-09
1.0E-09
Linear System Abort Not Converged = True
True
Linear System Preconditioning = "ILU0"
"ILU0"
Linear System Residual Output = 1
1
End
Solver 3
Exec Solver = After All
After All
Procedure = File "SaveData" "SaveLine"
File "SaveData" "SaveLine"
Filename = "hansbreen_surface.dat"
"hansbreen_surface.dat"
File Append = Logical False
Logical False
End
Solver 4
Exec Solver = After TimeStep
After TimeStep
Procedure = "./MySaveData" "SaveScalars"
"./MySaveData" "SaveScalars"
Filename = "hansbreen_scalars.dat"
"hansbreen_scalars.dat"
File Append = Logical True
Logical True
Variable 1 = String "flow solution"
String "flow solution"
Operator 1 = String "Volume"
String "Volume"
Variable 2 = String "Velocity 1"
String "Velocity 1"
Operator 2 = String "max abs"
String "max abs"
Variable 3 = String "flow solution"
String "flow solution"
Operator 3 = String "Convective flux"
String "Convective flux"
Variable 4 = String "cpu time"
String "cpu time"
Variable 5 = String "cpu memory"
String "cpu memory"
End
Equation 1
Active Solvers(4) = 1 2 3 4
Size 4 = 1 2 3 4
1 2 3 4
NS Convect = Logical False
Logical False
End
Boundary Condition 1
Target Boundaries = 1
1
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
End
Boundary Condition 2
Target Boundaries = 2
2
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 3
Target Boundaries = 3
3
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 4
Target Boundaries = 4
4
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 5
Target Boundaries = 5
5
Flow Force BC = Logical True
Logical True
Normal-Tangential Velocity = Logical False
Logical False
Velocity 1 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "U(tx)"
Velocity 2 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "V(tx)"
Velocity 3 = Variable Coordinate 1, Coordinate 2, Coordinate 3
Variable Coordinate 1, Coordinate 2, Coordinate 3
Real MATC "W(tx)"
Flux Integrate = Logical True
Logical True
End
Boundary Condition 6
Target Boundaries = 6
6
Save Line = Logical True
Logical True
Flux Integrate = Logical True
Logical True
End
MAIN: Done
MAIN: -----------------------
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: navier-stokes...done.
OptimizeBandwidth: Half bandwidth without optimization: 3018
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 672
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [./ComputeDevStressNS]...[ComputeDevStress_Init]
Loading user function library: [./ComputeDevStressNS]...[ComputeDevStress]
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: stresssolver...done.
OptimizeBandwidth: Half bandwidth without optimization: 3018
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 672
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [SaveData]...[SaveLine_Init]
WARNING:: ListFind:
WARNING:: ListFind: Requested property: [Equation], not found
WARNING:: ListFind:
Loading user function library: [SaveData]...[SaveLine]
Loading user function library: [./MySaveData]...[SaveScalars_Init]
WARNING:: ListFind:
WARNING:: ListFind: Requested property: [Equation], not found
WARNING:: ListFind:
Loading user function library: [./MySaveData]...[SaveScalars]
At line 2158 of file Lists.f90
Fortran runtime error: Bad real number in item 1 of list input
Regards, Jaime
Re: Fortran runtime error
Hi Jaime,
OK, i could finally reproduce the trouble and found out that the U, V & W MATC
functions try to import/use a undefined variable called 'b_cons'.
I'll try to do something about the diagnostics though.
Regards, Juha
OK, i could finally reproduce the trouble and found out that the U, V & W MATC
functions try to import/use a undefined variable called 'b_cons'.
I'll try to do something about the diagnostics though.
Regards, Juha
Re: Fortran runtime error
Hi again,
current svn trunk should now point out this kind of error directly.
OTOH it could be useful to test such, somewhat complicated, MATC constructions
using the standalone 'matc' executable in the Elmer distribution. (I could have spotted
the trouble immediately myself if i'd done that
Regards, Juha
current svn trunk should now point out this kind of error directly.
OTOH it could be useful to test such, somewhat complicated, MATC constructions
using the standalone 'matc' executable in the Elmer distribution. (I could have spotted
the trouble immediately myself if i'd done that
Regards, Juha
Re: Fortran runtime error
Hi,
Thank you very much, now it's working!!
Regards
Jaime
Thank you very much, now it's working!!
Regards
Jaime