====================================================================
 
      Description:
      This routine computes, with a jumping algorithm, the border of
      the stability (or acceptance loss) once a direction (from the
      origin of course) is fixed in the phase space.
 
      Arguments:
      r1   = radius of stability. r1 computed is still stable, but 1%
             close the border of stability. r1 needs to have an initial
             value which is out of the border of stability.
      n1   = are the angles varialbes (see gen_x subroutine)
      n2   = "
      n3   = "
      nmax = maximum value for n1,n2,n3
 
      NOTE:
           if r=0 is not stable, because COD is so deformed that r=0
           is not stable, then r1 is set equal to zero, and a message
           is prompt and written in the log file.
 
 ================================================
 
subroutine borderda_ex(r1,n1,n2,n3,nmax,lost_ext,gen_x_ext)
 
COMMONS: 
common-initi.f
common-2Dsc-simulation.f
 
SUBROUTINES USED:
  messages_n
  gen_x_ext
 
micromap-index
Mon Jul 17 08:57:31 CEST 2017