anchor.diff:    if(Jok==1)then
       Jcur=0
     else
  -    write(*,*),'set'
       call setBC(x)
       call syncState(x)
       call findAdvJac(grid,rho,rhou,rhoE,p,JacP,JacC)
  
positive.code:       call writeState('rst_'//trim(adjustl(tmpStr))//'.vtk')
        tNext=tNext+tInt
      end if
      write(*,*),t
    end do
    call clear()
  end program
  
negative.code:   double precision::rhs(5)
    
    n=max(floor(-log10(pGamm)),0)
    ! adaptive number of Jacobi iteration
    if(n>=7)then
      z(1:5*grid%nC)=res(1:5*grid%nC)
      ier=0
      return
    else if(n>=6)then
      n=0
    else if(n>=5)then
      n=1
    else if(n>=4)then
      n=2
    else
      n=3
  