anchor.diff:  program extended_examples
   
     ! modules
  +  use omp_lib
   
     use constants, only: ir
   
  
positive.code:   real (kind=ir), allocatable :: data_in(:,:,:), data_out(:,:,:)
    real (kind=ir), allocatable :: cflx(:,:,:), cfly(:,:,:)
    real (kind=ir), allocatable :: dcoeff, kcoeff
  
    integer :: ierr
    integer :: i, j, k, it
    integer :: istart, iend, jstart, jend, kstart, kend
    integer :: nincout
    integer :: iseed(80), nseed
  
  #ifdef _NETCDF
    character(len=255) :: fname = 'out.nc'
  #endif
  
    write(*,'(a)') 'Welcome to stencil3d!'
  
  
negative.code:     call swap_data()
  
    end do
  
    ! stop timer
    timespent = timespent + omp_get_wtime();
    PRINT *, nt ' time steps took total ', timespent, ' seconds, average ', timespent/nt , ' seconds'
  
    ! ****************** cleanup ******************
  
  