Changeset 39203

Timestamp:

Nov 30, 2015, 8:50:10 AM (11 years ago)

Author:

eugene

Message:

update galmotion demo

File:

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

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

@@ -17 +17 @@
 $W_sol  =   7.20; # km/sec
 
-$trans_phi = -62.8717488056*RAD_DEG;
-$trans_Xo  = 282.8594812080*RAD_DEG;
+$trans_phi = -62.8717488056;
+$trans_Xo  = 282.8594812080;
 $trans_xo  =  32.9319185700;
 
@@ -47 +47 @@
   set B = dasin(phi)
 
+  resize -n galactic 1600 800
   region 0 0 90 ait; cgrid -c red
-  # cplot L B
-  # break
+
+  resize -n celestial 1600 800
+  region 0 0 90 ait; cgrid -c red
 
   set R = L
@@ -64 +66 @@
   set uB = uB_gal + uB_sol
 
-  cneedles L B uL uB -scale $3
-  break
+  dev -n galactic
+  cneedles L B uL uB -scale $3 -c black -lw 0
+  cplot L B -pt 0 -sz 0.3 -c black
+  label +x Galactic -fn courier 24
 
   # C1, C2 are from http://arxiv.org/pdf/1306.2945v2.pdf
@@ -75 +79 @@
   set uR = cosBinv * (C1 * uL - C2 * uB);
   set uD = cosBinv * (C1 * uB + C2 * uL);
+
+  dev -n celestial
+  style -c red -lw 0 -pt 0 -sz 0.3; cggrid; style -c blue; galactic
+  cplot R D -pt 0 -sz 0.3 -c black
+  cneedles R D uR uD -scale $3 -c black -lw 0
+  label +x Celestial -fn courier 24
+
+  
 end
+
+macro cggrid
+
+ local delta off
+ $delta = 2.0
+
+ delete _lgrid _bgrid
+
+ create _tgrid 0 360 $delta
+ set _n = ramp(_tgrid) 
+ set _lgrid_t = _tgrid * (_n % 2) + (_tgrid - $delta) * (_n % 2 == 0)
+
+ for off -70 90 20
+   set _bgrid_t = zero(_lgrid_t) + $off
+   concat _lgrid_t _lgrid
+   concat _bgrid_t _bgrid
+ end
+
+ create _tgrid -90 90 $delta
+ set _n = ramp(_tgrid) 
+ set _bgrid_t = _tgrid * (_n % 2) + (_tgrid - $delta) * (_n % 2 == 0)
+
+ for off 0 360 30
+   set _lgrid_t = zero(_bgrid_t) + $off
+   concat _lgrid_t _lgrid
+   concat _bgrid_t _bgrid
+ end
+
+ csystem G C _lgrid _bgrid
+ cplot _lgrid _bgrid -pt 100
+end