Merge branch 'per-system-initialization' into 'development'

enabling initialization per slip system

See merge request damask/DAMASK!937

Author: eisenlohr

src/phase_mechanical_plastic_dislotwin.f90

@@ -213,8 +213,15 @@ module function plastic_dislotwin_init() result(myPlasticity)
 
       f_edge       = math_expand(pl%get_as1dReal('f_edge',    requiredSize=size(N_sl), &
                                                  defaultVal=[(0.5_pREAL,i=1,size(N_sl))]),N_sl)
+
+#ifdef __GFORTRAN__
+      rho_mob_0    = pl%get_as1dReal('rho_mob_0', requiredChunks=N_sl)
+      rho_dip_0    = pl%get_as1dReal('rho_dip_0', requiredChunks=N_sl)
+#else
       rho_mob_0    = math_expand(pl%get_as1dReal('rho_mob_0', requiredSize=size(N_sl)),N_sl)
       rho_dip_0    = math_expand(pl%get_as1dReal('rho_dip_0', requiredSize=size(N_sl)),N_sl)
+#endif
+
       prm%v_0      = math_expand(pl%get_as1dReal('v_0',       requiredSize=size(N_sl)),N_sl)
       prm%b_sl     = math_expand(pl%get_as1dReal('b_sl',      requiredSize=size(N_sl)),N_sl)
       prm%Q_sl     = math_expand(pl%get_as1dReal('Q_sl',      requiredSize=size(N_sl)),N_sl)

src/phase_mechanical_plastic_phenopowerlaw.f90

@@ -144,8 +144,13 @@ module function plastic_phenopowerlaw_init() result(myPlasticity)
       prm%n_sl           = math_expand(pl%get_as1dReal('n_sl',          requiredSize=size(N_sl)), N_sl)
       prm%a_sl           = math_expand(pl%get_as1dReal('a_sl',          requiredSize=size(N_sl)), N_sl)
       prm%h_0_sl_sl      = math_expand(pl%get_as1dReal('h_0_sl-sl',     requiredSize=size(N_sl)), N_sl)
+#ifdef __GFORTRAN__
+      xi_0_sl            = pl%get_as1dReal('xi_0_sl',  requiredChunks=N_sl)
+      prm%xi_inf_sl      = pl%get_as1dReal('xi_inf_sl',requiredChunks=N_sl)
+#else
       xi_0_sl            = math_expand(pl%get_as1dReal('xi_0_sl',       requiredSize=size(N_sl)), N_sl)
       prm%xi_inf_sl      = math_expand(pl%get_as1dReal('xi_inf_sl',     requiredSize=size(N_sl)), N_sl)
+#endif
       prm%c_1            = math_expand(pl%get_as1dReal('c_1',           requiredSize=size(N_sl), &
                                                                         defaultVal=misc_zeros(size(N_sl))), N_sl)
       prm%c_2            = math_expand(pl%get_as1dReal('c_2',           requiredSize=size(N_sl), &