1 | include(dom.inc)
   2 | 
   3 |       SUBROUTINE pdf(ph,cb,x,f)
   4 | 
   5 |       IMPLICIT NONE
   6 | 
   7 |       include 'dom_constants.h'
   8 | 
   9 | !     IN
  10 |       DOM_REAL   :: ph, cb, x
  11 | 
  12 | !     OUT
  13 |       DOM_REAL   :: f
  14 | 
  15 |       IF(x.ne.0.and.ph.ne.0.and.cb.ne.0) THEN
  16 |         f = sqrt((ph*cb)/(2*pi*x*x*x)) *                                &
  17 |      &      exp(-(ph*(x-cb)*(x-cb))/(2*x*cb))  
  18 |       ELSE
  19 |         f = 0.
  20 |       ENDIF     
  21 | 
  22 |       END SUBROUTINE pdf