Added HTRUDeep

Also minor change to the _powerlaw function

lib/python/populate.py

@@ -382,10 +382,9 @@ def _powerlaw( minval, maxval, power):
     c = -1.0 * logmax * power
     nmax = 10.0**(power*logmin + c)
 
-    # Dunc's code uses a goto statement
     # slightly worried about inf loops here...
     while True:
-        log = logmin + (logmax-logmin)*random.random()
+        log = random.uniform(logmin, logmax)
         n = 10.0**(power*log + c)
 
         if nmax*random.random() <= n:

lib/surveys/HTRUD

@@ -0,0 +1,23 @@
+###############################################################################
+# model of the hitrun-high multibeam survey - main parameters
+###############################################################################
+  1.0  ! survey degradation factor
+  0.6  ! antenna gain (K/Jy)
+   4320. ! integration time (s)
+  0.064 ! sampling time (ms)
+  23.  ! system temperature (K)
+  1390 ! centre frequency (MHz)
+  340  ! bandwidth (MHz)
+  0.039 ! channel bandwidth (MHz)
+  2    ! number of polarizations
+  14.  ! full-width half maximum (arcmin)
+  0.   ! minimum RA (deg)
+  360. ! maximum RA (deg)
+ -90.  ! minimum DEC (deg)
+  10.  ! maximum DEC (deg)
+ -80  ! minimum Galactic longitude (deg)
+  30.  ! maximum Galactic longitude (deg)
+  0.   ! minimum abs(Galactic latitude) (deg)
+  3.5   ! maximum abs(Galactic latitude) (deg)
+  .35  ! fractional survey coverage (0-1)
+  9.0  ! signal-to-noise ratio