| Date | File | Comment by | Comment | Rating |
|---|---|---|---|---|
| 15 Jul 2009 | sun_position.m Sun position given observer time/location. | E, K | Very helpful submission |
|
| 23 Jan 2009 | sun_position.m Sun position given observer time/location. | Falofolio | Great tool! Good job! However, I wonder, if there will be an update on this script, after reading Jeff Gullett's remarks on the code? |
|
| 06 Jan 2009 | sun_position.m Sun position given observer time/location. | Gullett, Jeff | I have found a few other discrepancies in the algorithm. In the 2005 version of the document, the algorithm to compute the topocentric sun declination (eqn #39, line #825) contains the variable 'y' in the denominator, but the code presented later in the pdf document has the variable 'x' in its place. In the 2008 version of the algorithm document, this difference is resolved and the variable should be 'x' (as in the C-code, so it needs to be changed in the m-code presented on this page). Also, when computing the mean obliquity of the ecliptic (eqn #24, line #733), the documented steps have one of the terms as "-4680.93", but again, the C-code from the appendix used the term "-4680.96". Unfortunately, this discrepancy still exists in the 2008 document, so I don't know which term is correct. Using the wrong one will result in an error greater than the maximum error claimed by the algorithm designers. |
|
| 19 Nov 2008 | sun_position.m Sun position given observer time/location. | RAMOS, EDGAR | Congrats guy!! |
|
| 11 Nov 2008 | sun_position.m Sun position given observer time/location. | Salinas, Lu | Has anyone done this same script but in excel? I've been trying it and I still have some formula wrong... I guess I'm confuse with the t ut or the radians and degrees...
|
|
| 24 Oct 2008 | sun_position.m Sun position given observer time/location. | Gullett, Jeff | The algorithm has a flaw... The atmospheric refraction correction has a number of vertical asymptotes below -5 degrees true elevation to the sun. For elevation angles that hit these asymptotes, I was getting points that resulted in positive elevations, when the neighboring points would come out with the expected result of "sun well below the horizon". You need to add a conditional prior to applying the atmospheric correction to only add this in if the true elevation is greater than -5 degrees. As for the code itself, it is somewhat disorganized and difficult to read, but it does work. I had a specific application where I needed to improved the efficiency/readability of the code, so I ended up rewriting the thing from scratch based on the paper, except I store everything in radians instead of degrees. This is still a nice application, though. |
|
| 17 Oct 2008 | sun_position.m Sun position given observer time/location. | Xinan, Lin | Thank you for your perfect program!!! |
|
| 24 Oct 2007 | sun_position.m Sun position given observer time/location. | Johnstone, P | This m-file is great! Very useful, well written. |
|
| 03 Aug 2007 | sun_position.m Sun position given observer time/location. | Suwisto, Andora | Thanks for the program.....The reference from usno.navy is very well accepted and well known for scientist and engineer |
|
| 08 Feb 2007 | sun_position.m Sun position given observer time/location. | Nguyen, Son | Well report |
|
| 22 Jan 2007 | sun_position.m Sun position given observer time/location. | Young, George | The code works as advertise and is well documented. It even comes with a DOE report on the equations involved. The learning curve is extremely short. Thanks! |
|
| 29 May 2006 | sun_position.m Sun position given observer time/location. | nazar, nadeem | ||
| 15 Feb 2006 | sun_position.m Sun position given observer time/location. | TRAN THANH, DUNG | ||
| 18 Nov 2005 | sun_position.m Sun position given observer time/location. | Wales, Steve | ||
| 28 Aug 2005 | sun_position.m Sun position given observer time/location. | Steensen, frederik | Very high-quality. I have a translation to Delphi |
|
| 27 Aug 2005 | sun_position.m Sun position given observer time/location. | Saenger, Mat | Cool! |
|
| 10 Jul 2005 | sun_position.m Sun position given observer time/location. | Jin, Xin | Bingo! This is exactly what I am looking for. |
|
| 07 Oct 2004 | sun_position.m Sun position given observer time/location. | shengwang, lu | It is a grate work and very halpful for me.Thanks a lot. |
|
| 17 Mar 2004 | sun_position.m Sun position given observer time/location. | Bergami, Andrea | Very useful script! |
|
NOTICE: Any content you submit to MATLAB Central, including personal information, is not subject to the protections which may be afforded information collected under other sections of The MathWorks, Inc. Web site. You are entirely responsible for all content that you upload, post, e-mail, transmit or otherwise make available via MATLAB Central. The MathWorks does not control the content posted by visitors to MATLAB Central and, does not guarantee the accuracy, integrity, or quality of such content. Under no circumstances will The MathWorks be liable in any way for any content not authored by The MathWorks, or any loss or damage of any kind incurred as a result of the use of any content posted, e-mailed, transmitted or otherwise made available via MATLAB Central. Read the complete Terms prior to use.
Contact us at files@mathworks.com