use of pssolar near equinox results in error
|Target version:||Candidate for next bugfix release|
|Affected version:||5.x-svn||Platform:||Mac OS X|
I am creating an animation of the whole Earth in orthographic projection showing the terminator using pssolar. When I try to render the frames for some times close to the (spring) equinox, I get the following cryptic error message:
% gmt pssolar -JG50/50/3c -Rg -Td+d2010-03-21T08:00:00.000 -Ggray > test.ps
pssolar: Error: Two points in input list are antipodal - great circle resampling is not unique!
pssolar: Error: Fix input data or use project -A to generate athe desired great circle by providing an azimuth.
If the date/time is set to slightly earlier or later, there is no problem. I have also tried some other projections, and they also provide no problem.
I would not know what to do with the suggestion in the second line of the error message.
I am using GMT 5.4.1, which is not yet listed in the "affected versions" dropdown.
- Status changed from New to In Progress
Our algorithm that converts the polygon to include polar caps gets fooled here. The first two points in the internal data are
so since it it touching the pole it is flagged as a polar cap. Later, this involves doing a detour to the pole but we end up going to the wrong pole... Will investigate.
You can try
gmt pssolar -Td+d2010-03-21T08:00:00.000 -M | gmt psxy -JH0/15c -Rg -Baf -Gred > t.ps
Does that give the right answer? I know it is not your projection of choice but curious if the polygon is correct.
- File pssolar_test.mp4 added
Thanks for your explanation. And yes, using other projections, such as your example, the polygon looks fine to me.
I have attached an example of what I'm trying to achieve. I'm using pssolar to create a semi-transparent overlay of the part of the Earth that is in shadow. But at solstice it seems the entire Earth is in sunlight, when using this projection. A couple of hours later (not included to keep the movie file small), the shadow reappears.