1 | c subroutines ... for Malkmus !
   2 | c
   3 | c transmission functions,
   4 | c k-distributions,
   5 | c k-cumulatives,
   6 | c random k-generators.
   7 | c
   8 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   9 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  10 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  11 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  12 | c
  13 | c.......................................................................
  14 | c
  15 | c  transmission function for surface-surface exchanges
  16 | c
  17 | c  in :  * $\phi$                ---> phtrss
  18 | c        * $\overline k$         ---> cbtrss
  19 | c        * $l$                   ---> altrss
  20 | c
  21 | c  out : * ${\tilde \tau}^{ss}$  ---> trtrss
  22 | c
  23 | c.......................................................................
  24 |       subroutine trss(phtrss,cbtrss,altrss,trtrss)
  25 | c.......................................................................
  26 |       implicit double precision (a-h,o-z)
  27 | c.......................................................................
  28 |       trss1=phtrss*dsqrt(1.d0+2.d0*cbtrss*altrss/phtrss)
  29 |       trtrss=dexp(phtrss-trss1)
  30 | c.......................................................................
  31 |       return
  32 |       end