From 2458aa9c3ef65eb92d0c4cb3700c0ab8ed45910e Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Mon, 8 Dec 2014 16:27:23 +0000 Subject: [PATCH] BBar stabilisation now optional (default is off) --- code/numerics.f90 | 9 +++++++-- 1 file changed, 7 insertions(+), 2 deletions(-) diff --git a/code/numerics.f90 b/code/numerics.f90 index 93d243c1e..b981d3c6f 100644 --- a/code/numerics.f90 +++ b/code/numerics.f90 @@ -168,6 +168,8 @@ module numerics thermalOrder = 2_pInt, & damageOrder = 2_pInt, & vacancyDiffusionOrder = 2_pInt + logical, protected, public :: & + BBarStabilisation = .false. #endif public :: numerics_init @@ -432,10 +434,12 @@ subroutine numerics_init vacancyDiffusionOrder = IO_intValue(line,positions,2_pInt) case ('petsc_optionsfem') petsc_optionsFEM = trim(line(positions(4):)) + case ('bbarstabilisation') + BBarStabilisation = IO_intValue(line,positions,2_pInt) > 0_pInt #else case ('err_struct_tolabs','err_struct_tolrel','err_thermal_tol','err_damage_tol','err_vacancydiffusion_tol',& ! found FEM parameter for spectral/Abaqus/Marc build 'itmaxfem', 'itminfem','maxcutbackfem','maxstaggerediter','integrationorder',& - 'structorder','thermalorder', 'damageorder','petsc_optionsfem') + 'structorder','thermalorder', 'damageorder','petsc_optionsfem','bbarstabilisation') call IO_warning(40_pInt,ext_msg=tag) #endif case default ! found unknown keyword @@ -579,8 +583,9 @@ subroutine numerics_init write(6,'(a24,1x,es8.1)') ' err_struct_tolRel: ',err_struct_tolRel write(6,'(a24,1x,es8.1)') ' err_thermal_tol: ',err_thermal_tol write(6,'(a24,1x,es8.1)') ' err_damage_tol: ',err_damage_tol - write(6,'(a24,1x,es8.1)') ' err_vacancyDiffusion_tol:',err_vacancydiffusion_tol + write(6,'(a24,1x,es8.1)') ' err_vacancyDiff_tol: ',err_vacancyDiffusion_tol write(6,'(a24,1x,a)') ' PETSc_optionsFEM: ',trim(petsc_optionsFEM) + write(6,'(a24,1x,L8)') ' B-Bar stabilisation: ',BBarStabilisation #endif endif mainProcess3