Fortran runtime error

General discussion about Elmer
Post Reply
jaime
Posts: 6
Joined: 02 Sep 2010, 00:35

Fortran runtime error

Post by jaime »

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
Juha
Site Admin
Posts: 357
Joined: 21 Aug 2009, 15:11

Re: Fortran runtime error

Post by Juha »

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
jaime
Posts: 6
Joined: 02 Sep 2010, 00:35

Re: Fortran runtime error

Post by jaime »

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
Juha
Site Admin
Posts: 357
Joined: 21 Aug 2009, 15:11

Re: Fortran runtime error

Post by Juha »

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
Juha
Site Admin
Posts: 357
Joined: 21 Aug 2009, 15:11

Re: Fortran runtime error

Post by Juha »

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
jaime
Posts: 6
Joined: 02 Sep 2010, 00:35

Re: Fortran runtime error

Post by jaime »

Hi,
Thank you very much, now it's working!!
Regards
Jaime
Post Reply