1 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   2 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   3 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   4 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   5 | c
   6 | c.......................................................................
   7 | c
   8 | c  k-sampling for gas-surface exchanges
   9 | c
  10 | c  in :  * $\phi$                ---> phrags
  11 | c        * $\overline k$         ---> cbrags
  12 | c        * $l$                   ---> alrags
  13 | c
  14 | c  out : * $k$                   ---> rarags
  15 | c
  16 | c.......................................................................
  17 |       subroutine rags(phrags,cbrags,alrags,rarags)
  18 | c.......................................................................
  19 |       implicit double precision (a-h,o-z)
  20 | c.......................................................................
  21 |       rags1=dsqrt(1.d0+2.d0/phrags*cbrags*alrags)
  22 |       rags2=phrags*rags1
  23 |       rags3=rags1/cbrags
  24 |       call gira(rags2,rags3,rags4)
  25 | 
  26 |       if (rags4.eq. 0.d+0) then
  27 |          rarags = cbrags
  28 |          ELSE
  29 |          rarags=1.d0/rags4
  30 |       endif
  31 | c.......................................................................
  32 |       return
  33 |       end