223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

somewhat better readable format without as much linebreaks

This commit is contained in:
Christoph Kords 2012-02-10 07:50:29 +00:00
parent 7a23833a39
commit ff53b9754b
1 changed files with 35 additions and 62 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

97
code/constitutive_nonlocal.f90
View File

@ -445,74 +445,47 @@ enddo
!*** sanity checks !*** sanity checks
if (myStructure < 1 .or. myStructure > 3)& if (myStructure < 1 .or. myStructure > 3) call IO_error(205)
call IO_error(205_pInt) if (sum(constitutive_nonlocal_Nslip(:,i)) <= 0_pInt) call IO_error(235,ext_msg='Nslip')
if (sum(constitutive_nonlocal_Nslip(:,i)) <= 0_pInt)&
call IO_error(235_pInt,ext_msg='Nslip')
do o = 1,maxval(phase_Noutput) do o = 1,maxval(phase_Noutput)
if(len(constitutive_nonlocal_output(o,i)) > 64)& if(len(constitutive_nonlocal_output(o,i)) > 64) call IO_error(666)
call IO_error(666_pInt)
enddo enddo
do f = 1,lattice_maxNslipFamily do f = 1,lattice_maxNslipFamily
if (constitutive_nonlocal_Nslip(f,i) > 0_pInt) then if (constitutive_nonlocal_Nslip(f,i) > 0_pInt) then
if (constitutive_nonlocal_rhoSglEdgePos0(f,i) < 0.0_pReal)& if (constitutive_nonlocal_rhoSglEdgePos0(f,i) < 0.0_pReal) call IO_error(235,ext_msg='rhoSglEdgePos0')
call IO_error(235_pInt,ext_msg='rhoSglEdgePos0') if (constitutive_nonlocal_rhoSglEdgeNeg0(f,i) < 0.0_pReal) call IO_error(235,ext_msg='rhoSglEdgeNeg0')
if (constitutive_nonlocal_rhoSglEdgeNeg0(f,i) < 0.0_pReal)& if (constitutive_nonlocal_rhoSglScrewPos0(f,i) < 0.0_pReal) call IO_error(235,ext_msg='rhoSglScrewPos0')
call IO_error(235_pInt,ext_msg='rhoSglEdgeNeg0') if (constitutive_nonlocal_rhoSglScrewNeg0(f,i) < 0.0_pReal) call IO_error(235,ext_msg='rhoSglScrewNeg0')
if (constitutive_nonlocal_rhoSglScrewPos0(f,i) < 0.0_pReal)& if (constitutive_nonlocal_rhoDipEdge0(f,i) < 0.0_pReal) call IO_error(235,ext_msg='rhoDipEdge0')
call IO_error(235_pInt,ext_msg='rhoSglScrewPos0') if (constitutive_nonlocal_rhoDipScrew0(f,i) < 0.0_pReal) call IO_error(235,ext_msg='rhoDipScrew0')
if (constitutive_nonlocal_rhoSglScrewNeg0(f,i) < 0.0_pReal)& if (constitutive_nonlocal_burgersPerSlipFamily(f,i) <= 0.0_pReal) call IO_error(235,ext_msg='burgers')
call IO_error(235_pInt,ext_msg='rhoSglScrewNeg0') if (constitutive_nonlocal_lambda0PerSlipFamily(f,i) <= 0.0_pReal) call IO_error(235,ext_msg='lambda0')
if (constitutive_nonlocal_rhoDipEdge0(f,i) < 0.0_pReal)& if (constitutive_nonlocal_minimumDipoleHeightPerSlipFamily(f,1,i) <= 0.0_pReal) &
call IO_error(235_pInt,ext_msg='rhoDipEdge0') call IO_error(235,ext_msg='minimumDipoleHeightEdge')
if (constitutive_nonlocal_rhoDipScrew0(f,i) < 0.0_pReal)& if (constitutive_nonlocal_minimumDipoleHeightPerSlipFamily(f,2,i) <= 0.0_pReal) &
call IO_error(235_pInt,ext_msg='rhoDipScrew0') call IO_error(235,ext_msg='minimumDipoleHeightScrew')
if (constitutive_nonlocal_burgersPerSlipFamily(f,i) <= 0.0_pReal)& if (constitutive_nonlocal_peierlsStressPerSlipFamily(f,1,i) <= 0.0_pReal) call IO_error(235,ext_msg='peierlsStressEdge')
call IO_error(235_pInt,ext_msg='burgers') if (constitutive_nonlocal_peierlsStressPerSlipFamily(f,2,i) <= 0.0_pReal) call IO_error(235,ext_msg='peierlsStressScrew')
if (constitutive_nonlocal_lambda0PerSlipFamily(f,i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='lambda0')
if (constitutive_nonlocal_minimumDipoleHeightPerSlipFamily(f,1,i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='minimumDipoleHeightEdge')
if (constitutive_nonlocal_minimumDipoleHeightPerSlipFamily(f,2,i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='minimumDipoleHeightScrew')
if (constitutive_nonlocal_peierlsStressPerSlipFamily(f,1,i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='peierlsStressEdge')
if (constitutive_nonlocal_peierlsStressPerSlipFamily(f,2,i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='peierlsStressScrew')
endif endif
enddo enddo
if (any(constitutive_nonlocal_interactionSlipSlip(1:maxval(lattice_interactionSlipSlip(:,:,myStructure)),i) < 0.0_pReal))& if (any(constitutive_nonlocal_interactionSlipSlip(1:maxval(lattice_interactionSlipSlip(:,:,myStructure)),i) < 0.0_pReal)) &
call IO_error(235_pInt,ext_msg='interaction_SlipSlip') call IO_error(235,ext_msg='interaction_SlipSlip')
if (constitutive_nonlocal_R(i) < 0.0_pReal)& if (constitutive_nonlocal_R(i) < 0.0_pReal) call IO_error(235,ext_msg='r')
call IO_error(235_pInt,ext_msg='r') if (constitutive_nonlocal_atomicVolume(i) <= 0.0_pReal) call IO_error(235,ext_msg='atomicVolume')
if (constitutive_nonlocal_atomicVolume(i) <= 0.0_pReal)& if (constitutive_nonlocal_Dsd0(i) <= 0.0_pReal) call IO_error(235,ext_msg='selfDiffusionPrefactor')
call IO_error(235_pInt,ext_msg='atomicVolume') if (constitutive_nonlocal_Qsd(i) <= 0.0_pReal) call IO_error(235,ext_msg='selfDiffusionEnergy')
if (constitutive_nonlocal_Dsd0(i) <= 0.0_pReal)& if (constitutive_nonlocal_aTolRho(i) <= 0.0_pReal) call IO_error(235,ext_msg='aTol_rho')
call IO_error(235_pInt,ext_msg='selfDiffusionPrefactor') if (constitutive_nonlocal_doublekinkwidth(i) <= 0.0_pReal) call IO_error(235,ext_msg='doublekinkwidth')
if (constitutive_nonlocal_Qsd(i) <= 0.0_pReal)& if (constitutive_nonlocal_solidSolutionEnergy(i) <= 0.0_pReal) call IO_error(235,ext_msg='solidSolutionEnergy')
call IO_error(235_pInt,ext_msg='selfDiffusionEnergy') if (constitutive_nonlocal_solidSolutionSize(i) <= 0.0_pReal) call IO_error(235,ext_msg='solidSolutionSize')
if (constitutive_nonlocal_aTolRho(i) <= 0.0_pReal)& if (constitutive_nonlocal_solidSolutionConcentration(i) <= 0.0_pReal) call IO_error(235,ext_msg='solidSolutionConcentration')
call IO_error(235_pInt,ext_msg='aTol_rho') if (constitutive_nonlocal_p(i) <= 0.0_pReal .or. constitutive_nonlocal_p(i) > 1.0_pReal) call IO_error(235,ext_msg='p')
if (constitutive_nonlocal_doublekinkwidth(i) <= 0.0_pReal)& if (constitutive_nonlocal_q(i) < 1.0_pReal .or. constitutive_nonlocal_q(i) > 2.0_pReal) call IO_error(235,ext_msg='q')
call IO_error(235_pInt,ext_msg='doublekinkwidth') if (constitutive_nonlocal_viscosity(i) <= 0.0_pReal) call IO_error(235,ext_msg='viscosity')
if (constitutive_nonlocal_solidSolutionEnergy(i) <= 0.0_pReal)& if (constitutive_nonlocal_fattack(i) <= 0.0_pReal) call IO_error(235,ext_msg='attackFrequency')
call IO_error(235_pInt,ext_msg='solidSolutionEnergy') if (constitutive_nonlocal_rhoSglScatter(i) < 0.0_pReal) call IO_error(235,ext_msg='rhoSglScatter')
if (constitutive_nonlocal_solidSolutionSize(i) <= 0.0_pReal)& if (constitutive_nonlocal_surfaceTransmissivity(i) < 0.0_pReal &
call IO_error(235_pInt,ext_msg='solidSolutionSize') .or. constitutive_nonlocal_surfaceTransmissivity(i) > 1.0_pReal) call IO_error(235,ext_msg='surfaceTransmissivity')
if (constitutive_nonlocal_solidSolutionConcentration(i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='solidSolutionConcentration')
if (constitutive_nonlocal_p(i) <= 0.0_pReal .or. constitutive_nonlocal_p(i) > 1.0_pReal)&
call IO_error(235_pInt,ext_msg='p')
if (constitutive_nonlocal_q(i) < 1.0_pReal .or. constitutive_nonlocal_q(i) > 2.0_pReal)&
call IO_error(235_pInt,ext_msg='q')
if (constitutive_nonlocal_viscosity(i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='viscosity')
if (constitutive_nonlocal_fattack(i) <= 0.0_pReal)&
call IO_error(235_pInt,ext_msg='attackFrequency')
if (constitutive_nonlocal_rhoSglScatter(i) < 0.0_pReal)&
call IO_error(235_pInt,ext_msg='rhoSglScatter')
if (constitutive_nonlocal_surfaceTransmissivity(i) < 0.0_pReal .or. constitutive_nonlocal_surfaceTransmissivity(i) > 1.0_pReal)&
call IO_error(235_pInt,ext_msg='surfaceTransmissivity')
!*** determine total number of active slip systems !*** determine total number of active slip systems