anchor.diff:      ! Subroutine parameter
       ! ---------------------------------------------------------------------------------------------- !
       integer :: l                                      ! Indice for x and y in loop
  -    integer :: max_i,max_j                            ! To prevent the SIGFPE
  +    integer :: max_i_l,max_j_l                        ! To prevent the SIGFPE
       real(kind=fPrec) :: xpow(0:this%n),ypow(0:this%m) ! Power of x and y
       real(kind=fPrec) :: r0                            ! Variable for the vector potential computation
       real(kind=fPrec) :: log_tmp
   
       ! Initialize vectors xpow and ypow
       ! ---------------------------------------------------------------------------------------------- !
  -    log_tmp=abs(log10_mb(abs(x)))
  -    if (log_tmp*(this%n)>230) then
  -      max_i=230/log_tmp
  -    else
  -      max_i=this%n
  -    endif
  -    log_tmp=abs(log10_mb(abs(y)))
  -    if (log_tmp*(this%m)>230) then
  -      max_j=230/log_tmp
  -    else
  -      max_j=this%m
  -    endif
  +    max_i_l=this%max_i;    max_j_l=this%max_j
   
       xpow(0)=one
  -    do l=1,max_i
  +    do l=1,max_i_l
         xpow(l)=xpow(l-1)*x 
       enddo
       ypow(0)=one
  -    do l=1,max_j
  +    do l=1,max_j_l
         ypow(l)=ypow(l-1)*y 
       enddo
   
  
positive.code:     do l=1,this%lx
        di=real(this%ij_TAx(1,l,z) + 1, fPrec)
        dj=real(this%ij_TAx(2,l,z)    , fPrec)
        if ((dj>0).and.(this%ij_TAx(1,l,z)+1<=max_i).and.(this%ij_TAx(2,l,z)-1<=max_j)) then
          r0 = r0 + (dj*(xpow(this%ij_TAx(1,l,z) + 1)*(ypow(this%ij_TAx(2,l,z) - 1)*this%TAx(l,z))))/di
        endif
      enddo
  
negative.code:     real(kind=fPrec) :: valA    ! Coefficient given by Horner subroutine
      real(kind=fPrec) :: dzover2 ! Dsigma/2
      real(kind=fPrec) :: g2d2inv ! 1/(delta+1)
      real(kind=fPrec) :: log_tmp
     
      ! Initialize the step size in z
      ! ---------------------------------------------------------------------------------------------- !
      dzover2=half*this%dz !loc
  
      ! Check size vectors xpow and ypow   (Prevent SIGFPE)
      ! ---------------------------------------------------------------------------------------------- !
      log_tmp=abs(log10_mb(x))
      if (log_tmp*(this%n)>230) then
        this%max_i=230/log_tmp-1
      else
        this%max_i=this%n
      endif
      log_tmp=abs(log10_mb(y))
      if (log_tmp*(this%m)>230) then
        this%max_j=230/log_tmp-1
      else
        this%max_j=this%m
      endif
  
      ! Initialize the step size in z
      ! ---------------------------------------------------------------------------------------------- !
  