Just looked at this approach.
Funny, to calculate the rotation angle as the longitude difference, plus a correction.
Hadn't thought of that as a method, previously.
But if it works, it works.
A few questions (I like to be curious):
1. Did you discover, why they define os = rcentrum?
2. I see that you define sitelat and sitelon input in degrees and minutes, instead of decimals? Most people would use decimals, I daresay.
3. Any reason why you don't use radius=6371 instead of 6378? And use rcentrum=42164, directly (without horb)?
4. Did you discover, in what cases abs(landa*180/pi) > 90 ; so that the error value sets in?
Edit: just when the longitude difference is too large?
(In my calculator, I calculate till 'invisible/below horizon', from the calculated elevation angle. So that leads to my 'error message'.)
Greetz,
A33
Funny, to calculate the rotation angle as the longitude difference, plus a correction.
Hadn't thought of that as a method, previously.
But if it works, it works.
A few questions (I like to be curious):
1. Did you discover, why they define os = rcentrum?
2. I see that you define sitelat and sitelon input in degrees and minutes, instead of decimals? Most people would use decimals, I daresay.
3. Any reason why you don't use radius=6371 instead of 6378? And use rcentrum=42164, directly (without horb)?
4. Did you discover, in what cases abs(landa*180/pi) > 90 ; so that the error value sets in?
Edit: just when the longitude difference is too large?
(In my calculator, I calculate till 'invisible/below horizon', from the calculated elevation angle. So that leads to my 'error message'.)
Greetz,
A33