#ifdef TIMING
time_vec=time_vec+MPI_Wtime()-time01
#endif
+C Introduction of shielding effect first for each peptide group
+C the shielding factor is set this factor is describing how each
+C peptide group is shielded by side-chains
+C the matrix - shield_fac(i) the i index describe the ith between i and i+1
+ if (shield_mode.gt.0) then
+ call set_shield_fac
+ endif
c print *,"Processor",myrank," left VEC_AND_DERIV"
if (ipot.lt.6) then
#ifdef SPLITELE