pssolar: problem with filling of terminators
|Target version:||Candidate for next minor release|
At certain dates/times the filling of terminators fails. Modified pssolar example, "(46) Day-light terminators and other sunlight parameters", included.
The time difference between the two plots is one hour (2017-09-26T11:40:58 vs. 2017-09-26T12:40:58).
Using Version 6.0.0_r18950
- Status changed from New to In Progress
- Assignee set to Paul
- Target version set to Candidate for next minor release
- Affected version changed from 5.x-svn to 6.x-svn
Fixing a few issues in pssolar but not fixed yet. Made a test script from this failure and added it to our test suite. r19268.
Tomorrow, on Jan 10th, the sun rises for the first time since early December in Harstad, Norway. See strange Sunrise/Sunset report for Jan 9th (i.e. Sun never comes up)
andreas:~$ gmt pssolar -I16.545833/68.8+d2018-01-09 Sun current position: long = 165.178329 lat = -22.135870 Azimuth = 18.7361 Elevation = -42.4641 Sunrise = -2147483648:00 Sunset = -2147483648:00 Noon = 11:01 Duration = -2147483648:00 andreas:~$ gmt pssolar -I16.545833/68.8+d2018-01-10 Sun current position: long = 165.281372 lat = -21.992271 Azimuth = 18.5835 Elevation = -42.3336 Sunrise = 10:46 Sunset = 11:16 Noon = 11:01 Duration = 00:30
- Status changed from In Progress to Resolved
- % Done changed from 0 to 100
Should be fixed. The terminator touched the north pole but did not get further south than -88 or so. Still, gmtlib_determine_pole determined it enclosed the S pole. Curiously, if you dump the path to file and plot with psxy it does the right thing. I added some more checking in the case a pole is touched and now it works. In r20187.