1 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   2 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   3 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   4 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   5 | c
   6 | c.......................................................................
   7 | c
   8 | c  k-sampling for surface-surface exchanges
   9 | c
  10 | c  in :  * $\phi$                ---> phrass
  11 | c        * $\overline k$         ---> cbrass
  12 | c        * $l$                   ---> alrass
  13 | c
  14 | c  out : * $k$                   ---> rarass
  15 | c
  16 | c.......................................................................
  17 |       subroutine rass(phrass,cbrass,alrass,rarass)
  18 | c.......................................................................
  19 |       implicit double precision (a-h,o-z)
  20 | c.......................................................................
  21 |  
  22 |       rass1=dsqrt(1.d0+2.d0*cbrass*alrass/phrass)
  23 |       rass2=phrass*rass1
  24 |       rass3=cbrass/rass1
  25 |  
  26 |       call gira(rass2,rass3,rarass)
  27 |   
  28 | c.......................................................................
  29 |       return
  30 |       end