#macro radec2pos(h,m,dec)
  #local ra=h*15+m*0.25;
 <cos(radians(ra))*cos(radians(dec)), sin(radians(dec)),
sin(radians(ra))*cos(radians(dec))>
#end

/* VENUS
echov(radec2pos(18,11,67.16))
echov(radec2pos(18,11,67.16))
*/
// NorthPoles from RaDec Wiki

#declare Yp = array[11] {
  <0.122353441903628, 0.897797192115935, -0.423071905334163>
  <0.091274333318625, 0.87840046198445, -0.469128579882686>
  <0.018690819744392, 0.921592341274493, -0.387708924015789>
  <0.0, 1.0, 0.0>
  <0.446158464461141, 0.797442047751279, -0.406237375272049>
  <-0.014654507041372, 0.902510203469804, -0.430419072598169>
  <0.085438088017201, 0.993650637868801, 0.073203435567528>
  <0.212015998366143, 0.261768361853467, 0.941555384015917>
  <0.358536930454372, 0.681359937097356, -0.638122171389651>
  <-0.673685502106019, -0.15798530500782, 0.72193385268579>
  <0.0,1.0,0.0>
}

#declare Xp = array[11] {
  vnormalize(vcross(Yp[0],<0,1,0>))
  vnormalize(vcross(Yp[1],<0,1,0>))
  vnormalize(vcross(Yp[2],<0,1,0>))
  <1.0, 0.0, 0.0>
  vnormalize(vcross(Yp[4],<0,1,0>))
  vnormalize(vcross(Yp[5],<0,1,0>))
  vnormalize(vcross(Yp[6],<0,1,0>))
  vnormalize(vcross(Yp[7],<0,1,0>))
  vnormalize(vcross(Yp[8],<0,1,0>))
  vnormalize(vcross(Yp[9],<0,1,0>))
  <1.0,0.0,0.0>
}

#declare Zp = array[11] {
  vnormalize(vcross(Xp[0],Yp[0]))
  vnormalize(vcross(Xp[1],Yp[1]))
  vnormalize(vcross(Xp[2],Yp[2]))
  <0.0, 0.0, 1.0>
  vnormalize(vcross(Xp[4],Yp[4]))
  vnormalize(vcross(Xp[5],Yp[5]))
  vnormalize(vcross(Xp[6],Yp[6]))
  vnormalize(vcross(Xp[7],Yp[7]))
  vnormalize(vcross(Xp[8],Yp[8]))
  vnormalize(vcross(Xp[9],Yp[9]))
  <0.0, 0.0, 1.0>
}
// example
// rotate <0,rate[2]*days,0>
// Venusmatrix = matrix
<Xp[2].x,Xp[2].y,Xp[2].z,Yp[2].x,Yp[2].y,Yp[2].z,Zp[2].x,Zp[2].y,Zp[2].z,0,0,0>


/*
alias rates {
  ; (name,daily motion in host rads per day) return degrees per day rotation for
planets and moons
  if (!$1) return 0
  if ($1 == mercury) return 6.138505
  if ($1 == venus) return -1.481369
  if ($1 == earth) return 360.98564736
  if ($1 == earth-blued1024) return 360.985626
  if ($1 == mars) return 350.891899
  if ($1 == jupiter) return 870.536642
  if ($1 == saturn) return 801.1126565
  if ($1 == uranus) return -501.16242
  if ($1 == neptune) return 536.312849
  if ($1 == pluto) return -56.36246
  if ($1 == planet-x) return 187.0
  var %sync = $rad2deg($2)
  if ($1 == triton) return $calc(%sync * -1)
  return %sync
}


rates are *-1 because rotate <0,-n,0> lookin down = counterclockwise - (sun sets
in the west)

then we have deltaRot alignment: I have 2, earth and mars
%br.longs 0 0 5.0 40.0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0
thinking earth should be 0 and 5 is hardly detectable but seemed to look better
in my right-hand planetarium

I guess its up to me, since no one is going for it.
And anyone else has plenty of data to get started.


Peace,aQ
(unless you prefer violence)

Post a reply to this message