preparing storage of output parameters

2019-02-13 08:06:22 +01:00 · 2019-02-13 08:06:22 +01:00 · a421525d15
parent f34c10a477
commit a421525d15
4 changed files with 121 additions and 30 deletions
--- a/src/source_damage_anisoBrittle.f90
+++ b/src/source_damage_anisoBrittle.f90
@ -56,6 +56,8 @@ module source_damage_anisoBrittle
     totalNcleavage
   integer(pInt), dimension(:), allocatable :: &
     Ncleavage
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID                                                                                       !< ID of each post result output
 end type tParameters

 type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@ -112,8 +114,6 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
   config_phase, &
   material_Nphase, &
   MATERIAL_partPhase
- use numerics,only: &
-   numerics_integrator
 use lattice, only: &
   lattice_maxNcleavageFamily, &
   lattice_NcleavageSystem
@ -123,17 +123,22 @@ subroutine source_damage_anisoBrittle_init(fileUnit)

 integer(pInt), allocatable, dimension(:) :: chunkPos
 integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
+ integer(pInt) :: NofMyPhase,p   ,i
 integer(pInt) :: Nchunks_CleavageFamilies = 0_pInt, j
- character(len=pStringLen) :: &
-   extmsg = ''
 character(len=65536) :: &
   tag  = '', &
   line = ''
 integer(pInt),          dimension(0), parameter :: emptyIntArray    = [integer(pInt)::]
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID

- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_anisoBrittle_LABEL//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ANISOBRITTLE_LABEL//' init  -+>>>'
 write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"

@ -185,6 +190,24 @@ subroutine source_damage_anisoBrittle_init(fileUnit)
   
   prm%Ncleavage = config%getInts('ncleavage',defaultVal=emptyIntArray)
   
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ANISOBRITTLE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('anisobrittle_drivingforce')
+
+     end select
+
+   enddo
+
   end associate
   
 enddo
--- a/src/source_damage_anisoDuctile.f90
+++ b/src/source_damage_anisoDuctile.f90
@ -57,6 +57,8 @@ module source_damage_anisoDuctile
     totalNslip
   integer(pInt), dimension(:), allocatable :: &
     Nslip
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID
 end type tParameters

 type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@ -113,8 +115,6 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
   config_phase, &
   material_Nphase, &
   MATERIAL_partPhase
- use numerics,only: &
-   numerics_integrator
 use lattice, only: &
   lattice_maxNslipFamily, &
   lattice_NslipSystem
@ -124,17 +124,22 @@ subroutine source_damage_anisoDuctile_init(fileUnit)

 integer(pInt), allocatable, dimension(:) :: chunkPos
 integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
+ integer(pInt) :: NofMyPhase,p ,i
 integer(pInt) :: Nchunks_SlipFamilies = 0_pInt, j
- character(len=pStringLen) :: &
-   extmsg = ''
 character(len=65536) :: &
   tag  = '', &
   line = ''
 integer(pInt),          dimension(0), parameter :: emptyIntArray    = [integer(pInt)::]
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID

- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_anisoDuctile_LABEL//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ANISODUCTILE_LABEL//' init  -+>>>'
 write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"

@ -184,6 +189,24 @@ subroutine source_damage_anisoDuctile_init(fileUnit)
   
   prm%Nslip  = config%getInts('nslip',defaultVal=emptyIntArray)
   
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ANISODUCTILE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('anisoductile_drivingforce')
+
+     end select
+
+   enddo
+
   end associate
   
 enddo
--- a/src/source_damage_isoBrittle.f90
+++ b/src/source_damage_isoBrittle.f90
@ -39,6 +39,8 @@ module source_damage_isoBrittle
     critStrainEnergy, &
     N, &
     aTol
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID
 end type tParameters

 type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@ -95,23 +97,26 @@ subroutine source_damage_isoBrittle_init(fileUnit)
   config_phase, &
   material_Nphase, &
   MATERIAL_partPhase
- use numerics,only: &
-   numerics_integrator

 implicit none
 integer(pInt), intent(in) :: fileUnit

 integer(pInt), allocatable, dimension(:) :: chunkPos
 integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
- character(len=pStringLen) :: &
-   extmsg = ''
+ integer(pInt) :: NofMyPhase,p,i   
 character(len=65536) :: &
   tag  = '', &
   line = ''
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID

- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_isoBrittle_label//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+ 
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ISOBRITTLE_LABEL//' init  -+>>>'
 write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"

@ -156,6 +161,24 @@ subroutine source_damage_isoBrittle_init(fileUnit)
   if (prm%N                  <= 0.0_pReal) extmsg = trim(extmsg)//' isobrittle_n'
   if (prm%critStrainEnergy   <= 0.0_pReal) extmsg = trim(extmsg)//' isobrittle_criticalstrainenergy'
   
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ISOBRITTLE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('isobrittle_drivingforce')
+       
+     end select
+
+   enddo
+
   end associate
   
 enddo
--- a/src/source_damage_isoDuctile.f90
+++ b/src/source_damage_isoDuctile.f90
@ -40,6 +40,8 @@ module source_damage_isoDuctile
     critPlasticStrain, &
     N, &
     aTol
+   integer(kind(undefined_ID)), allocatable, dimension(:) :: &
+     outputID
 end type tParameters

 type(tParameters), dimension(:), allocatable, private :: param                                     !< containers of constitutive parameters (len Ninstance)
@ -97,23 +99,25 @@ subroutine source_damage_isoDuctile_init(fileUnit)
   material_Nphase, &
   MATERIAL_partPhase

- use numerics,only: &
-   numerics_integrator
-
 implicit none
 integer(pInt), intent(in) :: fileUnit

 integer(pInt), allocatable, dimension(:) :: chunkPos
 integer(pInt) :: Ninstance,mySize=0_pInt,phase,instance,source,sourceOffset,o
- integer(pInt) :: sizeState, sizeDotState, sizeDeltaState
- integer(pInt) :: NofMyPhase,p   
- character(len=pStringLen) :: &
-   extmsg = ''
+ integer(pInt) :: NofMyPhase,p,i
 character(len=65536) :: &
   tag  = '', &
   line = ''
+ character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
+ integer(kind(undefined_ID)) :: &
+   outputID

- write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_damage_isoDuctile_label//' init  -+>>>'
+ character(len=pStringLen) :: &
+   extmsg = ''
+ character(len=65536), dimension(:), allocatable :: &
+   outputs
+
+ write(6,'(/,a)')   ' <<<+-  source_'//SOURCE_DAMAGE_ISODUCTILE_LABEL//' init  -+>>>'
 write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"

@ -158,6 +162,24 @@ subroutine source_damage_isoDuctile_init(fileUnit)
   if (prm%N                   <= 0.0_pReal) extmsg = trim(extmsg)//' isoductile_ratesensitivity'
   if (prm%critPlasticStrain   <= 0.0_pReal) extmsg = trim(extmsg)//' isoductile_criticalplasticstrain'
   
+!--------------------------------------------------------------------------------------------------
+!  exit if any parameter is out of range
+   if (extmsg /= '') &
+     call IO_error(211_pInt,ext_msg=trim(extmsg)//'('//SOURCE_DAMAGE_ISODUCTILE_LABEL//')')
+
+!--------------------------------------------------------------------------------------------------
+!  output pararameters
+   outputs = config%getStrings('(output)',defaultVal=emptyStringArray)
+   allocate(prm%outputID(0))
+   do i=1_pInt, size(outputs)
+     outputID = undefined_ID
+     select case(outputs(i))
+       case ('isoductile_drivingforce')
+       
+     end select
+
+   enddo
+
   end associate
   
 enddo