2007-03-20 19:25:22 +05:30
|
|
|
|
|
|
|
!##############################################################
|
|
|
|
MODULE prec
|
|
|
|
!##############################################################
|
|
|
|
|
|
|
|
implicit none
|
2008-02-21 18:20:37 +05:30
|
|
|
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Precision of real and integer variables ***
|
2008-03-07 20:51:46 +05:30
|
|
|
integer, parameter :: pReal = selected_real_kind(15,300) ! 15 significant digits, up to 1e+-300
|
|
|
|
integer, parameter :: pInt = selected_int_kind(9) ! up to +- 1e9
|
2008-02-21 18:20:37 +05:30
|
|
|
|
|
|
|
! *** Strain increment considered significant ***
|
2008-02-18 23:41:24 +05:30
|
|
|
real(pReal), parameter :: relevantStrain = 1.0e-7_pReal
|
2008-02-15 15:34:49 +05:30
|
|
|
|
|
|
|
! *** Numerical parameters ***
|
|
|
|
integer(pInt), parameter :: ijaco = 1_pInt ! frequency of FEM Jacobi update
|
|
|
|
integer(pInt), parameter :: nCutback = 10_pInt ! cutbacks in time-step integration
|
|
|
|
integer(pInt), parameter :: nReg = 1_pInt ! regularization attempts for Jacobi inversion
|
2008-02-18 16:16:18 +05:30
|
|
|
real(pReal), parameter :: pert_Fg = 1.0e-5_pReal ! strain perturbation for FEM Jacobi
|
2008-02-15 15:34:49 +05:30
|
|
|
integer(pInt), parameter :: nOuter = 10_pInt ! outer loop limit
|
2008-02-21 18:20:37 +05:30
|
|
|
integer(pInt), parameter :: nInner = 1000_pInt ! inner loop limit
|
2008-02-15 15:34:49 +05:30
|
|
|
real(pReal), parameter :: reltol_Outer = 1.0e-4_pReal ! relative tolerance in outer loop (state)
|
|
|
|
real(pReal), parameter :: reltol_Inner = 1.0e-6_pReal ! relative tolerance in inner loop (Lp)
|
|
|
|
real(pReal), parameter :: abstol_Inner = 1.0e-8_pReal ! absolute tolerance in inner loop (Lp)
|
2007-03-20 19:25:22 +05:30
|
|
|
|
|
|
|
END MODULE prec
|