2007-03-20 19:25:22 +05:30
|
|
|
|
|
|
|
!##############################################################
|
|
|
|
MODULE prec
|
|
|
|
!##############################################################
|
|
|
|
|
|
|
|
implicit none
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Precision of real and integer variables ***
|
2007-03-20 19:25:22 +05:30
|
|
|
integer, parameter :: pReal = 8
|
2007-03-20 21:28:53 +05:30
|
|
|
integer, parameter :: pInt = 4
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Numerical parameters ***
|
|
|
|
! *** How frequently the jacobian is recalculated ***
|
2007-12-14 19:06:04 +05:30
|
|
|
integer (pInt), parameter :: ijaco = 1_pInt
|
2007-04-10 16:52:06 +05:30
|
|
|
! *** Maximum number of internal cutbacks in time step ***
|
2007-04-11 15:35:08 +05:30
|
|
|
integer(pInt), parameter :: nCutback = 7_pInt
|
2007-04-10 16:52:06 +05:30
|
|
|
! *** Maximum number of regularization attempts for Jacobi inversion ***
|
2007-04-11 15:35:08 +05:30
|
|
|
integer(pInt), parameter :: nReg = 1_pInt
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Perturbation of strain array for numerical calculation of FEM Jacobi matrix ***
|
|
|
|
real(pReal), parameter :: pert_e=1.0e-5_pReal
|
2007-04-11 15:35:08 +05:30
|
|
|
! *** Maximum number of iterations in outer (state variables) loop ***
|
2008-01-11 00:23:57 +05:30
|
|
|
integer(pInt), parameter :: nState = 500_pInt
|
2007-04-11 15:35:08 +05:30
|
|
|
! *** Convergence criteria for outer (state variables) loop ***
|
2007-11-15 16:51:37 +05:30
|
|
|
real(pReal), parameter :: reltol_State = 1.0e-6_pReal
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Maximum number of iterations in inner (stress) loop ***
|
2008-01-11 00:23:57 +05:30
|
|
|
integer(pInt), parameter :: nStress = 1000_pInt
|
2007-11-15 16:51:37 +05:30
|
|
|
! *** Convergence criteria for inner (stress) loop ***
|
|
|
|
real(pReal), parameter :: reltol_Stress = 1.0e-6_pReal
|
|
|
|
! *** Convergence criteria for inner (stress) loop ***
|
|
|
|
real(pReal), parameter :: abstol_Stress = 1.0e3_pReal
|
|
|
|
! *** Convergence criteria for inner (stress) loop ***
|
|
|
|
real(pReal), parameter :: abstol_ResStress = 1.0_pReal
|
2007-03-20 19:25:22 +05:30
|
|
|
|
|
|
|
END MODULE prec
|