anchor.diff:                          y(:,7)  = yold(:,7)
                           ac(:,7) = zero 
                           if (ivconstraint .eq. 1) then
  -     &                       
   c ... applying the volume constraint
   c
                              call solvecon (y,    x,      iBC,  BC, 
  
positive.code:               call TD_bufferData()
                call TD_writeData(fvarts, .false.)
              endif
              
              !Update the Mach Control
              if(exts .and. exMC) then
                !Note: the function MC_updateState must be called after
  
negative.code: c              call write_field2(myrank,'a'//char(0),'ybar'//char(0),
  c     &                          4,ybar,'d'//char(0),nshg,ndof+8,
  c     &                         lstep,istep)
              endif  ! end write fields to restart files
  
   2000    continue  !end of NSTEP loop
   2001    continue
  
           ttim(1) = REAL(secs(0.0)) / 100. - ttim(1)
           ttim(2) = secs(0.0)              - ttim(2)
  