anchor.diff:      call fkinmalloc(ioutFKIN,routFKIN,ier)
       call fkinspgmr(50,4,ier)
       !call fkindense(nEq,ier)
  +    call fkinsetrin('MAX_STEP',huge(1d0),ier)
       !call fkinsetrin('FNORM_TOL',1d-5,ier)
       !call fkinsetrin('SSTEP_TOL',1d-9,ier)
       call fkinsetiin('PRNT_LEVEL',1,ier)
  
positive.code:     do i=1,grid%nC
        j=(i-1)*5
        ! FIXME scale the variable
        x(j+1)=p(i)/xscale(j+1)
        x(j+2:j+4)=u(:,i)/xscale(j+2:j+4)
        x(j+5)=temp(i)/xscale(j+5)
      end do
      !$omp end parallel do
    end subroutine
  
negative.code:     dt=1d-4
      do i=1,grid%nC
        j=(i-1)*5
        xscale(j+1)=0.1d5
        xscale(j+2:j+4)=10d0
        xscale(j+5)=10d0
        rscale(j+1)=0.1d0
        rscale(j+2:j+4)=1d0
        rscale(j+5)=1d4
  