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
#1 Updated by Paul almost 3 years ago
- 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.
#2 Updated by Andreas over 2 years ago
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
#6 Updated by Paul over 2 years ago
- 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.