10
« on: March 26, 2024, 12:36:50 AM »
Hi Jean, Ed
can I find this calculation in PD as from Livio description? In Italy the call them direzioni placidiane
formula completa semiarchi placidus senza latitudine pianeta
complete formula semiarchs placidus without planet latitude
non traduco in inglese, le formule sono intuibili.
ob = obliquità 23.438130
dati forniti dal software
lob/lo[0][a] = longitudine pianeta promissore
lat = latitudine
dp = declinazione
dh = distanza oraria
ad = differenza ascensionale
a = significatore (0 = sole 1 = luna da tabella effemeridi)
aa = promissore
aspectangle[aaa] = aaa = tutti gli aspetti es 0,60,90,120 ecc
distance=command pd corrisponde all'angolo minore fra due punti
normalize= command pd riduce entro 360°
:psig0;`calcolo significatore xlat0 parte 1
u=1;
lob= lo[0][a];
los =lob;
ED = sin(los)*sin(ob);
ED = asin(ED);`distanza equinoziale Gansten
da = tan(lat)*tan(ED);da = asin(da);
Sad = normalize(90 + da);semiarco diurno
San = normalize(90 - da);semiarco notturno
dp = decp[a];
lp = 0;@cra; calcolo ar senza latitudine
dm_mc = normalize(ar-armc);if (dm_mc > 180) dm_mc -= 360;distanza meridiana mc
dm_ic = normalize(ar-aric);if (dm_ic > 180) dm_ic -= 360);distanza meridiana ic
return;
:ppro0;`promissore xlat0 parte 2
ora si calcola l'angolo del promissore alla congiunzione o agli aspetti col significatore
partendo dalla longitudine del promissore,aggiungendo o sottraendo l'angolo
lob = normalize(lo[0][aa]+/-aspectangle[aaa]);` nuova longitudine del promissore
'la direzione è sempre calcolata con latitudine 0 (lp = 0) proiettata sull' equatore
u=0;
lop = lob;
ED = sin(lop)*sin(ob);
ED = asin(ED);`gansten
da = tan(lat)*tan(ED);da = asin(da);
Sad = normalize(90 + da);
San = normalize(90 - da);
dp = decp[aa];
lp = 0;@cra;
dm_mc = normalize(ar-armc);if (dm_mc > 180) dm_mc -= 360;
dm_ic = normalize(ar-aric);if (dm_ic > 180) dm_ic -= 360);
if (normalize(-lo[0][14] + los) > 180)
{`se sopra orizzonte
dh = dm_mc[1] / Sad[1];
delta = normalize(dm_mc[0] - (dh * Sad[0]));
}
else
{`se sotto orizzonte
dh = dm_ic[1] / San[1];
delta = normalize(dm_ic[0] - (dh * San[0]));
}
delta = angolo direzione finale
si consiglia di usare 1 grado = 1 anno
return;