id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
809	IMPORTANT: Corrected bug in ana_srflux.h	arango		"Finally, I have time to browse the literature, Googling, and write a Matlab script to check a bug in '''ana_srflx.h''' as discussed in the [https://www.myroms.org/forum/viewtopic.php?f=19&t=3718&p=19840#p19840 Forum] Indeed, the '''15.0''' factor used to compute the '''zenith''' is incorrect. We need to have instead:
{{{
        zenith=cff1+cff2*COS(Hangle-lonr(i,j)*deg2rad)
}}}

In my Matlab script, I compared against a '''WRF-ROMS''' coupled solution for the US West Coast at a particular point in the grid:

[[Image(https://www.myroms.org/trac/swrad.png, center, 600)]]

Here, the '''X'''-axis is the time since initialization and it is in '''GMT''' time zone. We need to subtract 8 hours for local Pacific Standard Time ('''PST''') to analyze the results. As you can see, the '''WRF-ROMS''' coupled solution values agree well with the red curve that computes '''zenith''' with updated above equation. I am now satisfied that I proved to myself that it is a bug.

Many thanks to all that brought this bug to my attention."	bug	closed	major	Release ROMS/TOMS 3.7	Nonlinear	3.7	Fixed