geom.f

Go to the documentation of this file.

      subroutine geom(a,b,c,r,g)                                        
c***********************************************************************
      implicit real*8 (a-h,o-z)
      g=c                                                               
      r=0.0d0                                                             
      if(g.le.1.0e-15)return
      dif1=b-a                                                          
      dif2=c-b                                                          
      dif2p=dabs(100.0d0*dif2/c)
      if(dif2p .lt. 1.0e-2) go to 10                                     
      r=dif2/dif1                                                       
      g=a+dif1/(1.0-r)                                                  
10    continue                                                          
      return                                                            
      end                                                               
c***********************************************************************