Changeset 39215

Timestamp:

Dec 3, 2015, 3:52:08 PM (11 years ago)

Author:

eugene

Message:

add check.region command

File:

• branches/eam_branches/ipp-20151113/Ohana/src/uniphot/doc/galmotion.sh (modified) (3 diffs)

branches/eam_branches/ipp-20151113/Ohana/src/uniphot/doc/galmotion.sh

@@ -4 +4 @@
 
 ## FEAST-HIP
-if (0)
+if (1)
   $A_oort = +14.82; # km/sec/kpc
   $B_oort = -12.37; # km/sec/kpc
@@ -13 +13 @@
 
 ## ROESSER:
-if (1)
+if (0)
   $A_oort = +14.50; # km/sec/kpc
   $B_oort = -13.00; # km/sec/kpc
@@ -302 +302 @@
  cplot _lgrid _bgrid -pt 100
 end
+
+
+macro check.region
+  if ($0 != 6)
+    echo "USAGE: check.region (Rmin) (Rmax) (Dmin) (Dmax) (dist)"
+    break
+  end
+  
+  mkstars 10000 $5
+
+  set keep = (R > $1) && (R < $2) && (D > $3) && (D < $4)
+  subset  Rs =  R if keep
+  subset  Ds =  D if keep
+  subset uRs = uR if keep
+  subset uDs = uD if keep
+  lim -0.01 0.01 -0.01 0.01; clear; box; plot uRs uDs -pt 0 -sz 0.5
+  dot -0.003 -0.004 -pt 7 -c red -sz 4
+end