anchor.diff:                call MC_updateState()
                 call MC_applyBC(BC)
                 call MC_printState()
  -                 
  +
                 !Write the state if a restart is also being written. 
                 if(mod(lstep,ntout).eq.0 ) call MC_writeState(lstep)
               endif
  
positive.code:             tcorecp2 = TMRC()
              write(6,*) 'T(core) cpu = ',tcorecp2-tcorecp1
           endif
          
        call wtime
  
        call destroyfncorp
  
   3000 continue !end of NTSEQ loop
  c
        print*, "DONE"
          call destruct_sum_vi_area
          call ifbc_mfree
  c     
  c.... ---------------------->  Post Processing  <----------------------
  c     
  c.... print out the last step
  c     
  !      if ( (irs .ge. 1) .and. ((mod(lstep, ntout) .ne. 0) .or.
  !     &    (nstp .eq. 0))) then
  !          if( (mod(lstep, ntoutv) .eq. 0) .and.
  
negative.code:               !restart is also being written. 
                if(mod(lstep, ntout) == 0) call BC_writePhase(lstep)
              endif
  
              !.... Yi Chen Duct geometry8
              if(isetBlowing_Duct.gt.0)then
                if(ifixBlowingVel_Duct.eq.0)then
  