Changeset 41733

Timestamp:

Jul 30, 2021, 10:31:04 AM (5 years ago)

Author:

eugene

Message:

working on galactic motion model

File:

• trunk/Ohana/src/uniphot/doc/galmotion.sh (modified) (2 diffs)

trunk/Ohana/src/uniphot/doc/galmotion.sh

@@ -138 +138 @@
 
   mkmotion_dmerr L B Dm $3 $4
+
+  # C1, C2 are from http://arxiv.org/pdf/1306.2945v2.pdf
+  set C1 =  dcos(D)*$trans_cos_phi + dsin(D)*(dcos(R)*$trans_sin_phi_sin_Xo - dsin(R)*$trans_sin_phi_cos_Xo);
+  set C2 =                                -1*(dcos(R)*$trans_sin_phi_cos_Xo + dsin(R)*$trans_sin_phi_sin_Xo);
+
+  set cosBinv = 1.0 / sqrt(C1*C1 + C2*C2);
+
+  set uR = cosBinv * (C1 * uL_o - C2 * uB_o);
+  set uD = cosBinv * (C1 * uB_o + C2 * uL_o);
 end
 
@@ -171 +180 @@
   set uL_o = uL_gal + uL_sol
   set uB_o = uB_gal + uB_sol
+
+  set uL = uL_o
+  set uB = uB_o
 
   set uL_sol_o = uL_sol