anchor.diff:  !DO i=1,p%elem_total
   !WRITE(*,*) p%Jacobian(:,i)
   !ENDDO
  +!STOP
      !CALL WrScr( "Finished reading input" )
      ! Allocate continuous states
      CALL AllocAry(x%q,p%dof_total,'x%q',ErrStat2,ErrMsg2)
  
positive.code:                            elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                hhx(i,igp)*Qe(m,n)*hhx(j,igp)*Jaco(igp)*gw(igp)
                             elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                hhx(i,igp)*Pe(m,n)*hpx(j,igp)*Jaco(igp)*gw(igp)
                             elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                hpx(i,igp)*Oe(m,n)*hhx(j,igp)*Jaco(igp)*gw(igp)
                             elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                hpx(i,igp)*Stif(m,n)*hpx(j,igp)*Jaco(igp)*gw(igp)
                             elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                hhx(i,igp)*Ki(m,n)*hhx(j,igp)*Jaco(igp)*gw(igp)
                             elm(temp_id1,temp_id2) = elm(temp_id1,temp_id2) + &
  
negative.code:                                    elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                        hhx(i,igp)*Qd(m,n)*hhx(j,igp)*Jaco(igp)*gw(igp)
                                     elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                        hhx(i,igp)*Pd(m,n)*hpx(j,igp)*gw(igp)
                                     elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                        hpx(i,igp)*Od(m,n)*hhx(j,igp)*gw(igp)
                                     elk(temp_id1,temp_id2) = elk(temp_id1,temp_id2) + &
                                        hpx(i,igp)*Sd(m,n)*hpx(j,igp)*gw(igp)/Jaco(igp)
                                     elg(temp_id1,temp_id2) = elg(temp_id1,temp_id2) + &
                                        hhx(i,igp)*Xd(m,n)*hhx(j,igp)*Jaco(igp)*gw(igp)
                                     elg(temp_id1,temp_id2) = elg(temp_id1,temp_id2) + &
                                        hhx(i,igp)*Yd(m,n)*hpx(j,igp)*gw(igp)
                                     elg(temp_id1,temp_id2) = elg(temp_id1,temp_id2) + &
                                        hpx(i,igp)*Gd(m,n)*hhx(j,igp)*gw(igp)
                                     elg(temp_id1,temp_id2) = elg(temp_id1,temp_id2) + &
                                        hpx(i,igp)*betaC(m,n)*hpx(j,igp)*gw(igp)/Jaco(igp)
                             ENDDO
                         ENDDO
                     ENDDO
  