Bug #1155
pssolar: problem with filling of terminators
| Status: | Closed | Start date: | 2017-09-26 | |
|---|---|---|---|---|
| Priority: | Normal | Due date: | ||
| Assignee: |  Paul | % 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 1 year 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 1 year 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 over 1 year 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 over 1 year ago
Great! Revision? Still getting the same weird number as in the post above.
#5
Updated by Joaquim over 1 year ago
r19554
#6
Updated by Paul about 1 year 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.