Bug #1155

pssolar: problem with filling of terminators

Added by Andreas over 2 years ago. Updated almost 2 years ago.

Status:ClosedStart date:2017-09-26
Priority:NormalDue 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

pssolar_fill.sh Magnifier (1.08 KB) Andreas, 2017-09-26 00:56

pssolar_fill_ok.pdf (240 KB) Andreas, 2017-09-26 00:56

pssolar_fill_problem.pdf (230 KB) Andreas, 2017-09-26 00:56

History

#1 Updated by Paul over 2 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

More on this:

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 2 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 over 2 years ago

Great! Revision? Still getting the same weird number as in the post above.

#5 Updated by Joaquim over 2 years ago

r19554

#6 Updated by Paul about 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.

#7 Updated by Paul almost 2 years ago

  • Status changed from Resolved to Closed

Closed as fixed and no further feedback.

Also available in: Atom PDF