anchor.diff:     RETURN
   END IF
   
  +WRITE(UnOut, '(A)')  'NNodes, Nelems, NProps, NCMass'
   WRITE(UnOut, '(4(1x, I5))' ) Init%NNode, Init%NElem, Init%NProp, Init%NCMass
  +WRITE(UnOut, '(A)')  'Node No.    X       Y      Z'
   WRITE(UnOut, '(F6.0, E15.6, E15.6, E15.6)') ((Init%Nodes(i, j), j = 1, Init%JointsCol), i = 1, Init%NNode)
  +WRITE(UnOut, '(A)')  'Elem No.,    Node_I,     Node_J,      Prop_I,      Prop_J'
   WRITE(UnOut, '(5(1x, I6))') ((p%Elems(i, j), j = 1, Init%MembersCol), i = 1, Init%NElem)
  +WRITE(UnOut, '(A)')  'Prop No.     YoungE       ShearG       MatDens     XsecD      XsecT'
   WRITE(UnOut, '(F6.0, E15.6, E15.6, E15.6, E15.6, E15.6 ) ') ((Init%Props(i, j), j = 1, 6), i = 1, Init%NProp)
  +WRITE(UnOut, '(A)')  'No. of Reaction DOFs'
   WRITE(UnOut, '(I6)')p%NReact*6
  +WRITE(UnOut, '(A)')  'Reaction DOF ID      LOCK'
   WRITE(UnOut, '(I6, I6)') ((Init%BCs(i, j), j = 1, 2), i = 1, p%NReact*6)
  +WRITE(UnOut, '(A)')  'No. of Interface DOFs'
   WRITE(UnOut, '(I6)')Init%NInterf*6
  +WRITE(UnOut, '(A)')  'Interface DOF ID      LOCK'
   WRITE(UnOut, '(I6, I6)') ((Init%IntFc(i, j), j = 1, 2), i = 1, Init%NInterf*6)
  +WRITE(UnOut, '(A)')  'Number of concentrated masses'
   WRITE(UnOut, '(I6)')Init%NCMass
  +WRITE(UnOut, '(A)')  'JointCMass     Mass         JXX             JYY             JZZ'
   WRITE(UnOut, '(F6.0, E15.6, E15.6, E15.6, E15.6)') ((Init%Cmass(i, j), j = 1, 5), i = 1, Init%NCMass)
  +WRITE(UnOut, '(A)')  'Number of members'
   WRITE(UnOut, '(I6)')p%NMembers
  -
  +WRITE(UnOut, '(A)')  'Number of nodes per member'
   WRITE(tempStr, '(I10)') Init%NDiv + 1 
  +WRITE(UnOut, '(A)')  'Member Nodes'
   OutFmt = ('('//trim(tempStr)//'(I6))')
   WRITE(UnOut, trim(OutFmt)) ((Init%MemberNodes(i, j), j = 1, Init%NDiv+1), i = 1, p%NMembers)
   
  
positive.code:    DirCos = 0
     
     IF ( EqualRealNos(xz, 0.0) ) THEN
        IF( y2 < y1) THEN
           DirCos(1, 1) = 1.0
           DirCos(2, 3) = -1.0
           DirCos(3, 2) = 1.0
  
negative.code:            CALL GetNewProp(kprop, TempProps(Prop1, 2), TempProps(Prop1, 3),&
                             TempProps(Prop1, 4), d1+dd, t1+dt, TempProps, TempNProp, Init%PropSetsCol)           
             kelem = kelem + 1
             CALL GetNewElem(kelem, Node1, knode, Prop1, kprop, p)  
             nprop = kprop              
        ELSE
             kelem = kelem + 1
             CALL GetNewElem(kelem, Node1, knode, Prop1, Prop1, p)                
             nprop = Prop1 
        ENDIF
        
  