Search Comments and Ratings

go

Comments and Ratings

   
Date File Comment by Comment Rating
15 Jul 2009 sun_position.m Sun position given observer time/location. Author: Vincent Roy E, K

Very helpful submission

23 Jan 2009 sun_position.m Sun position given observer time/location. Author: Vincent Roy 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. Author: Vincent Roy 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. Author: Vincent Roy RAMOS, EDGAR

Congrats guy!!

11 Nov 2008 sun_position.m Sun position given observer time/location. Author: Vincent Roy 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...
Would anyone want to give it a check?

24 Oct 2008 sun_position.m Sun position given observer time/location. Author: Vincent Roy 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. Author: Vincent Roy Xinan, Lin

Thank you for your perfect program!!!

24 Oct 2007 sun_position.m Sun position given observer time/location. Author: Vincent Roy Johnstone, P

This m-file is great! Very useful, well written.

03 Aug 2007 sun_position.m Sun position given observer time/location. Author: Vincent Roy 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. Author: Vincent Roy Nguyen, Son

Well report

22 Jan 2007 sun_position.m Sun position given observer time/location. Author: Vincent Roy 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. Author: Vincent Roy nazar, nadeem

15 Feb 2006 sun_position.m Sun position given observer time/location. Author: Vincent Roy TRAN THANH, DUNG

18 Nov 2005 sun_position.m Sun position given observer time/location. Author: Vincent Roy Wales, Steve

28 Aug 2005 sun_position.m Sun position given observer time/location. Author: Vincent Roy Steensen, frederik

Very high-quality. I have a translation to Delphi

27 Aug 2005 sun_position.m Sun position given observer time/location. Author: Vincent Roy Saenger, Mat

Cool!

10 Jul 2005 sun_position.m Sun position given observer time/location. Author: Vincent Roy Jin, Xin

Bingo! This is exactly what I am looking for.

07 Oct 2004 sun_position.m Sun position given observer time/location. Author: Vincent Roy 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. Author: Vincent Roy Bergami, Andrea

Very useful script!

 

MATLAB Central Terms of Use

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