added some aliases for output to finally converge from 'shear_rate' to 'shearrate' and from 'accumulatedshear' to 'accumulated_shear'

code/config/Phase_Dislotwin_TWIP-Steel-FeMnC.config

@@ -6,18 +6,18 @@ plasticity          dislotwin
 
 #(output)           edge_density
 #(output)           dipole_density
-#(output)           shear_rate_slip
+#(output)           shearrate_slip
 #(output)           accumulated_shear_slip
 #(output)           mfp_slip
 #(output)           resolved_stress_slip
 #(output)           threshold_stress_slip
 #(output)           twin_fraction
-#(output)           shear_rate_twin
+#(output)           shearrate_twin
 #(output)           accumulated_shear_twin
 #(output)           mfp_twin
 #(output)           resolved_stress_twin
 #(output)           threshold_stress_twin
-#(output)           shear_rate_shearband	
+#(output)           shearrate_shearband	
 #(output)           resolved_stress_shearband	
 #(output)           sb_eigenvalues
 #(output)           sb_eigenvectors

code/config/Phase_Nonlocal_Aluminum.config

@@ -87,7 +87,7 @@ plasticity              nonlocal
 (output)                fluxDensity_screw_neg.z
 (output)                maximumDipoleHeight_edge
 (output)                maximumDipoleHeight_screw
-(output)                accumulatedshear
+(output)                accumulated_shear
 (output)                dislocationstress
 
 lattice_structure       fcc

code/config/Phase_Phenopowerlaw_Aluminum.config

@@ -6,12 +6,12 @@ plasticity              phenopowerlaw
 (output)                resistance_slip
 (output)                shearrate_slip
 (output)                resolvedstress_slip
-(output)                accumulatedshear_slip
+(output)                accumulated_shear_slip
 (output)                totalshear
 (output)                resistance_twin
 (output)                shearrate_twin
 (output)                resolvedstress_twin
-(output)                accumulatedshear_twin
+(output)                accumulated_shear_twin
 (output)                totalvolfrac_twin
 
 lattice_structure       fcc

code/plastic_disloKMC.f90

@@ -320,7 +320,7 @@ subroutine plastic_disloKMC_init(fileUnit)
              plastic_disloKMC_outputID(plastic_disloKMC_Noutput(instance),instance) = dipole_density_ID
              plastic_disloKMC_output(plastic_disloKMC_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_slip')
+           case ('shear_rate_slip','shearrate_slip')
              plastic_disloKMC_Noutput(instance) = plastic_disloKMC_Noutput(instance) + 1_pInt
              plastic_disloKMC_outputID(plastic_disloKMC_Noutput(instance),instance) = shear_rate_slip_ID
              plastic_disloKMC_output(plastic_disloKMC_Noutput(instance),instance) = &
@@ -360,7 +360,7 @@ subroutine plastic_disloKMC_init(fileUnit)
              plastic_disloKMC_outputID(plastic_disloKMC_Noutput(instance),instance) = twin_fraction_ID
              plastic_disloKMC_output(plastic_disloKMC_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_twin')
+           case ('shear_rate_twin','shearrate_twin')
              plastic_disloKMC_Noutput(instance) = plastic_disloKMC_Noutput(instance) + 1_pInt
              plastic_disloKMC_outputID(plastic_disloKMC_Noutput(instance),instance) = shear_rate_twin_ID
              plastic_disloKMC_output(plastic_disloKMC_Noutput(instance),instance) = &

code/plastic_disloUCLA.f90

@@ -328,7 +328,7 @@ subroutine plastic_disloUCLA_init(fileUnit)
              plastic_disloUCLA_outputID(plastic_disloUCLA_Noutput(instance),instance) = dipole_density_ID
              plastic_disloUCLA_output(plastic_disloUCLA_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_slip')
+           case ('shear_rate_slip','shearrate_slip')
              plastic_disloUCLA_Noutput(instance) = plastic_disloUCLA_Noutput(instance) + 1_pInt
              plastic_disloUCLA_outputID(plastic_disloUCLA_Noutput(instance),instance) = shear_rate_slip_ID
              plastic_disloUCLA_output(plastic_disloUCLA_Noutput(instance),instance) = &
@@ -368,7 +368,7 @@ subroutine plastic_disloUCLA_init(fileUnit)
              plastic_disloUCLA_outputID(plastic_disloUCLA_Noutput(instance),instance) = twin_fraction_ID
              plastic_disloUCLA_output(plastic_disloUCLA_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_twin')
+           case ('shear_rate_twin','shearrate_twin')
              plastic_disloUCLA_Noutput(instance) = plastic_disloUCLA_Noutput(instance) + 1_pInt
              plastic_disloUCLA_outputID(plastic_disloUCLA_Noutput(instance),instance) = shear_rate_twin_ID
              plastic_disloUCLA_output(plastic_disloUCLA_Noutput(instance),instance) = &

code/plastic_dislotwin.f90

@@ -376,7 +376,7 @@ subroutine plastic_dislotwin_init(fileUnit)
              plastic_dislotwin_outputID(plastic_dislotwin_Noutput(instance),instance) = dipole_density_ID
              plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_slip')
+           case ('shear_rate_slip','shearrate_slip')
              plastic_dislotwin_Noutput(instance) = plastic_dislotwin_Noutput(instance) + 1_pInt
              plastic_dislotwin_outputID(plastic_dislotwin_Noutput(instance),instance) = shear_rate_slip_ID
              plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &
@@ -416,7 +416,7 @@ subroutine plastic_dislotwin_init(fileUnit)
              plastic_dislotwin_outputID(plastic_dislotwin_Noutput(instance),instance) = twin_fraction_ID
              plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_twin')
+           case ('shear_rate_twin','shearrate_twin')
              plastic_dislotwin_Noutput(instance) = plastic_dislotwin_Noutput(instance) + 1_pInt
              plastic_dislotwin_outputID(plastic_dislotwin_Noutput(instance),instance) = shear_rate_twin_ID
              plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &
@@ -446,7 +446,7 @@ subroutine plastic_dislotwin_init(fileUnit)
              plastic_dislotwin_outputID(plastic_dislotwin_Noutput(instance),instance) = resolved_stress_shearband_ID
              plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &
                                                        IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('shear_rate_shearband')
+           case ('shear_rate_shearband','shearrate_shearband')
              plastic_dislotwin_Noutput(instance) = plastic_dislotwin_Noutput(instance) + 1_pInt
              plastic_dislotwin_outputID(plastic_dislotwin_Noutput(instance),instance) = shear_rate_shearband_ID
              plastic_dislotwin_output(plastic_dislotwin_Noutput(instance),instance) = &

code/plastic_nonlocal.f90

@@ -846,7 +846,7 @@ allocate(nonSchmidCoeff(lattice_maxNnonSchmid,maxNinstances),                  s
              plastic_nonlocal_outputID(plastic_nonlocal_Noutput(instance),instance) = maximumdipoleheight_screw_ID
              plastic_nonlocal_output(plastic_nonlocal_Noutput(instance),instance) = &
                IO_lc(IO_stringValue(line,positions,2_pInt))
-           case('accumulatedshear')
+           case ('accumulatedshear','accumulated_shear')
              plastic_nonlocal_Noutput(instance) = plastic_nonlocal_Noutput(instance) + 1_pInt
              plastic_nonlocal_outputID(plastic_nonlocal_Noutput(instance),instance) = accumulatedshear_ID
              plastic_nonlocal_output(plastic_nonlocal_Noutput(instance),instance) = &

code/plastic_phenopowerlaw.f90

@@ -269,7 +269,7 @@ subroutine plastic_phenopowerlaw_init(fileUnit)
              plastic_phenopowerlaw_outputID(plastic_phenopowerlaw_Noutput(instance),instance) = resistance_slip_ID
              plastic_phenopowerlaw_output(plastic_phenopowerlaw_Noutput(instance),instance) = &
                                                            IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('accumulatedshear_slip')
+           case ('accumulatedshear_slip','accumulated_shear_slip')
              plastic_phenopowerlaw_Noutput(instance) = plastic_phenopowerlaw_Noutput(instance) + 1_pInt
              plastic_phenopowerlaw_outputID(plastic_phenopowerlaw_Noutput(instance),instance) = accumulatedshear_slip_ID
              plastic_phenopowerlaw_output(plastic_phenopowerlaw_Noutput(instance),instance) = &
@@ -294,7 +294,7 @@ subroutine plastic_phenopowerlaw_init(fileUnit)
              plastic_phenopowerlaw_outputID(plastic_phenopowerlaw_Noutput(instance),instance) = resistance_twin_ID
              plastic_phenopowerlaw_output(plastic_phenopowerlaw_Noutput(instance),instance) = &
                                                            IO_lc(IO_stringValue(line,positions,2_pInt))
-           case ('accumulatedshear_twin')
+           case ('accumulatedshear_twin','accumulated_shear_twin')
              plastic_phenopowerlaw_Noutput(instance) = plastic_phenopowerlaw_Noutput(instance) + 1_pInt
              plastic_phenopowerlaw_outputID(plastic_phenopowerlaw_Noutput(instance),instance) = accumulatedshear_twin_ID
              plastic_phenopowerlaw_output(plastic_phenopowerlaw_Noutput(instance),instance) = &