Bug #1155
pssolar: problem with filling of terminators
Status: | Closed | Start date: | 2017-09-26 | |
---|---|---|---|---|
Priority: | Normal | Due date: | ||
Assignee: | % Done: | 100% | ||
Category: | - | |||
Target version: | Candidate for next minor release | |||
Affected version: | 6.x-svn | Platform: |
Description
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
History
#1
Updated by Paul over 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 about 3 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
#3
Updated by Joaquim about 3 years ago
OK, it now prints
Since you seam to be at "Harstad, Norway" wait another day before you can see the Sun
#4
Updated by Andreas about 3 years ago
Great! Revision? Still getting the same weird number as in the post above.
#5
Updated by Joaquim about 3 years ago
r19554
#6
Updated by Paul almost 3 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.
#7
Updated by Paul over 2 years ago
- Status changed from Resolved to Closed
Closed as fixed and no further feedback.